An improved epigenetic counter to track mitotic age in normal and precancerous tissues.

The cumulative number of stem cell divisions in a tissue, known as mitotic age, is thought to be a major determinant of cancer-risk. Somatic mutational and DNA methylation (DNAm) clocks are promising tools to molecularly track mitotic age, yet their relationship is underexplored and their potential for cancer risk prediction in normal tissues remains to be demonstrated. Here we build and validate an improved pan-tissue DNAm counter of total mitotic age called stemTOC. We demonstrate that stemTOC's mitotic age proxy increases with the tumor cell-of-origin fraction in each of 15 cancer-types, in precancerous lesions, and in normal tissues exposed to major cancer risk factors. Extensive benchmarking against 6 other mitotic counters shows that stemTOC compares favorably, specially in the preinvasive and normal-tissue contexts. By cross-correlating stemTOC to two clock-like somatic mutational signatures, we confirm the mitotic-like nature of only one of these. Our data points towards DNAm as a promising molecular substrate for detecting mitotic-age increases in normal tissues and precancerous lesions, and hence for developing cancer-risk prediction strategies.

Donor whole blood DNA methylation is not a strong predictor of acute graft versus host disease in unrelated donor allogeneic haematopoietic cell transplantation.

Allogeneic hematopoietic cell transplantation (HCT) is used to treat many blood-based disorders and malignancies, however it can also result in serious adverse events, such as the development of acute graft-versus-host disease (aGVHD). This study aimed to develop a donor-specific epigenetic classifier to reduce incidence of aGVHD by improving donor selection. Genome-wide DNA methylation was assessed in a discovery cohort of 288 HCT donors selected based on recipient aGVHD outcome; this cohort consisted of 144 cases with aGVHD grades III-IV and 144 controls with no aGVHD. We applied a machine learning algorithm to identify CpG sites predictive of aGVHD. Receiver operating characteristic (ROC) curve analysis of these sites resulted in a classifier with an encouraging area under the ROC curve (AUC) of 0.91. To test this classifier, we used an independent validation cohort (n = 288) selected using the same criteria as the discovery cohort. Attempts to validate the classifier failed with the AUC falling to 0.51. These results indicate that donor DNA methylation may not be a suitable predictor of aGVHD in an HCT setting involving unrelated donors, despite the initial promising results in the discovery cohort. Our work highlights the importance of independent validation of machine learning classifiers, particularly when developing classifiers intended for clinical use.

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.

Open access-enabled evaluation of epigenetic age acceleration in colorectal cancer and development of a classifier with diagnostic potential.

Aberrant DNA methylation (DNAm) is known to be associated with the aetiology of cancer, including colorectal cancer (CRC). In the past, the availability of open access data has been the main driver of innovative method development and research training. However, this is increasingly being eroded by the move to controlled access, particularly of medical data, including cancer DNAm data. To rejuvenate this valuable tradition, we leveraged DNAm data from 1,845 samples (535 CRC tumours, 522 normal colon tissues adjacent to tumours, 72 colorectal adenomas, and 716 normal colon tissues from healthy individuals) from 14 open access studies deposited in NCBI GEO and ArrayExpress. We calculated each sample's epigenetic age (EA) using eleven epigenetic clock models and derived the corresponding epigenetic age acceleration (EAA). For EA, we observed that most first- and second-generation epigenetic clocks reflect the chronological age in normal tissues adjacent to tumours and healthy individuals [e.g., Horvath (r = 0.77 and 0.79), Zhang elastic net (EN) (r = 0.70 and 0.73)] unlike the epigenetic mitotic clocks (EpiTOC, HypoClock, MiAge) (r < 0.3). For EAA, we used PhenoAge, Wu, and the above mitotic clocks and found them to have distinct distributions in different tissue types, particularly between normal colon tissues adjacent to tumours and cancerous tumours, as well as between normal colon tissues adjacent to tumours and normal colon tissue from healthy individuals. Finally, we harnessed these associations to develop a classifier using elastic net regression (with lasso and ridge regularisations) that predicts CRC diagnosis based on a patient's sex and EAAs calculated from histologically normal controls (i.e., normal colon tissues adjacent to tumours and normal colon tissue from healthy individuals). The classifier demonstrated good diagnostic potential with ROC-AUC = 0.886, which suggests that an EAA-based classifier trained on relevant data could become a tool to support diagnostic/prognostic decisions in CRC for clinical professionals. Our study also reemphasises the importance of open access clinical data for method development and training of young scientists. Obtaining the required approvals for controlled access data would not have been possible in the timeframe of this study.

Individualized tumor-informed circulating tumor DNA analysis for postoperative monitoring of non-small cell lung cancer.

We report a personalized tumor-informed technology, Patient-specific pROgnostic and Potential tHErapeutic marker Tracking (PROPHET) using deep sequencing of 50 patient-specific variants to detect molecular residual disease (MRD) with a limit of detection of 0.004%. PROPHET and state-of-the-art fixed-panel assays were applied to 760 plasma samples from 181 prospectively enrolled early stage non-small cell lung cancer patients. PROPHET shows higher sensitivity of 45% at baseline with circulating tumor DNA (ctDNA). It outperforms fixed-panel assays in prognostic analysis and demonstrates a median lead-time of 299 days to radiologically confirmed recurrence. Personalized non-canonical variants account for 98.2% with prognostic effects similar to canonical variants. The proposed tumor-node-metastasis-blood (TNMB) classification surpasses TNM staging for prognostic prediction at the decision point of adjuvant treatment. PROPHET shows potential to evaluate the effect of adjuvant therapy and serve as an arbiter of the equivocal radiological diagnosis. These findings highlight the potential advantages of personalized cancer techniques in MRD detection.

Individualized dynamic methylation-based analysis of cell-free DNA in postoperative monitoring of lung cancer.

BACKGROUND: The feasibility of DNA methylation-based assays in detecting minimal residual disease (MRD) and postoperative monitoring remains unestablished. We aim to investigate the dynamic characteristics of cancer-related methylation signals and the feasibility of methylation-based MRD detection in surgical lung cancer patients. METHODS: Matched tumor, tumor-adjacent tissues, and longitudinal blood samples from a cohort (MEDAL) were analyzed by ultra-deep targeted sequencing and bisulfite sequencing. A tumor-informed methylation-based MRD (timMRD) was employed to evaluate the methylation status of each blood sample. Survival analysis was performed in the MEDAL cohort (n = 195) and validated in an independent cohort (DYNAMIC, n = 36). RESULTS: Tumor-informed methylation status enabled an accurate recurrence risk assessment better than the tumor-naïve methylation approach. Baseline timMRD-scores were positively correlated with tumor burden, invasiveness, and the existence and abundance of somatic mutations. Patients with higher timMRD-scores at postoperative time-points demonstrated significantly shorter disease-free survival in the MEDAL cohort (HR: 3.08, 95% CI: 1.48-6.42; P = 0.002) and the independent DYNAMIC cohort (HR: 2.80, 95% CI: 0.96-8.20; P = 0.041). Multivariable regression analysis identified postoperative timMRD-score as an independent prognostic factor for lung cancer. Compared to tumor-informed somatic mutation status, timMRD-scores yielded better performance in identifying the relapsed patients during postoperative follow-up, including subgroups with lower tumor burden like stage I, and was more accurate among relapsed patients with baseline ctDNA-negative status. Comparing to the average lead time of ctDNA mutation, timMRD-score yielded a negative predictive value of 97.2% at 120 days prior to relapse. CONCLUSIONS: The dynamic methylation-based analysis of peripheral blood provides a promising strategy for postoperative cancer surveillance. TRIAL REGISTRATION: This study (MEDAL, MEthylation based Dynamic Analysis for Lung cancer) was registered on ClinicalTrials.gov on 08/05/2018 (NCT03634826). https://clinicaltrials.gov/ct2/show/NCT03634826 .

Subclassification of epithelioid sarcoma with potential therapeutic impact.

Epithelioid sarcoma is a rare and aggressive mesenchymal tumour, the genetic hallmark of which is the loss of expression of SMARCB1, a key member of the SWItch/Sucrose Non-Fermentable (SWI/SNF) chromatin remodelling complex. Hampered by its rarity, epithelioid sarcoma has received little research attention and therapeutic options for this disease remain limited. SMARCB1-deficient tumours also include malignant rhabdoid tumour, atypical teratoid and rhabdoid tumour, epithelioid malignant peripheral nerve sheath tumour, and poorly differentiated chordoma. Histologically, it can be challenging to distinguish epithelioid sarcoma from malignant rhabdoid tumour and other SMARCB1-deficient tumours, whereas methylation profiling shows that they represent distinct entities and facilitates their classification. Methylation studies on SMARCB1-deficient tumours, although not including epithelioid sarcomas, reported methylation subgroups which resulted in new clinical stratification and therapeutic approaches. In addition, emerging evidence indicates that immunotherapy, including immune checkpoint inhibitors, represents a promising therapeutic strategy for SMARCB1-deficient tumours. Here, we show that some epithelioid sarcomas share methylation patterns of malignant rhabdoid tumours indicating that this could help to distinguish these entities and guide treatment. Using gene expression data, we also showed that the immune environment of epithelioid sarcoma is characterised by a predominance of CD8+ lymphocytes and M2 macrophages. These findings have potential implications for the management of patients with epithelioid sarcoma. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Genomic-transcriptomic evolution in lung cancer and metastasis.

Intratumour heterogeneity (ITH) fuels lung cancer evolution, which leads to immune evasion and resistance to therapy1. Here, using paired whole-exome and RNA sequencing data, we investigate intratumour transcriptomic diversity in 354 non-small cell lung cancer tumours from 347 out of the first 421 patients prospectively recruited into the TRACERx study2,3. Analyses of 947 tumour regions, representing both primary and metastatic disease, alongside 96 tumour-adjacent normal tissue samples implicate the transcriptome as a major source of phenotypic variation. Gene expression levels and ITH relate to patterns of positive and negative selection during tumour evolution. We observe frequent copy number-independent allele-specific expression that is linked to epigenomic dysfunction. Allele-specific expression can also result in genomic-transcriptomic parallel evolution, which converges on cancer gene disruption. We extract signatures of RNA single-base substitutions and link their aetiology to the activity of the RNA-editing enzymes ADAR and APOBEC3A, thereby revealing otherwise undetected ongoing APOBEC activity in tumours. Characterizing the transcriptomes of primary-metastatic tumour pairs, we combine multiple machine-learning approaches that leverage genomic and transcriptomic variables to link metastasis-seeding potential to the evolutionary context of mutations and increased proliferation within primary tumour regions. These results highlight the interplay between the genome and transcriptome in influencing ITH, lung cancer evolution and metastasis.

Osteosarcoma: Novel prognostic biomarkers using circulating and cell-free tumour DNA.

AIM: Osteosarcoma (OS) is the most common primary bone tumour in children and adolescents. Circulating free (cfDNA) and circulating tumour DNA (ctDNA) are promising biomarkers for disease surveillance and prognostication in several cancer types; however, few such studies are reported for OS. The purpose of this study was to discover and validate methylation-based biomarkers to detect plasma ctDNA in patients with OS and explore their utility as prognostic markers. METHODS: Candidate CpG markers were selected through analysis of methylation array data for OS, non-OS tumours and germline samples. Candidates were validated in two independent OS datasets (n = 162, n = 107) and the four top-performing markers were selected. Methylation-specific digital droplet PCR (ddPCR) assays were designed and experimentally validated in OS tumour samples (n = 20) and control plasma samples. Finally, ddPCR assays were applied to pre-operative plasma and where available post-operative plasma from 72 patients with OS, and findings correlated with outcome. RESULTS: Custom ddPCR assays detected ctDNA in 69% and 40% of pre-operative plasma samples (n = 72), based on thresholds of one or two positive markers respectively. ctDNA was detected in 5/17 (29%) post-operative plasma samples from patients, which in four cases were associated with or preceded disease relapse. Both pre-operative cfDNA levels and ctDNA detection independently correlated with overall survival (p = 0.0015 and p = 0.0096, respectively). CONCLUSION: Our findings illustrate the potential of mutation-independent methylation-based ctDNA assays for OS. This study lays the foundation for multi-institutional collaborative studies to explore the utility of plasma-derived biomarkers in the management of OS.

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.

Novel epigenetic network biomarkers for early detection of esophageal cancer.

BACKGROUND: Early detection of esophageal cancer is critical to improve survival. Whilst studies have identified biomarkers, their interpretation and validity is often confounded by cell-type heterogeneity. RESULTS: Here we applied systems-epigenomic and cell-type deconvolution algorithms to a discovery set encompassing RNA-Seq and DNA methylation data from esophageal adenocarcinoma (EAC) patients and matched normal-adjacent tissue, in order to identify robust biomarkers, free from the confounding effect posed by cell-type heterogeneity. We identify 12 gene-modules that are epigenetically deregulated in EAC, and are able to validate all 12 modules in 4 independent EAC cohorts. We demonstrate that the epigenetic deregulation is present in the epithelial compartment of EAC-tissue. Using single-cell RNA-Seq data we show that one of these modules, a proto-cadherin module centered around CTNND2, is inactivated in Barrett's Esophagus, a precursor lesion to EAC. By measuring DNA methylation in saliva from EAC cases and controls, we identify a chemokine module centered around CCL20, whose methylation patterns in saliva correlate with EAC status. CONCLUSIONS: Given our observations that a CCL20 chemokine network is overactivated in EAC tissue and saliva from EAC patients, and that in independent studies CCL20 has been found to be overactivated in EAC tissue infected with the bacterium F. nucleatum, a bacterium that normally inhabits the oral cavity, our results highlight the possibility of using DNAm measurements in saliva as a proxy for changes occurring in the esophageal epithelium. Both the CTNND2/CCL20 modules represent novel promising network biomarkers for EAC that merit further investigation.

Evaluation of Epigenetic Age Acceleration Scores and Their Associations with CVD-Related Phenotypes in a Population Cohort.

We evaluated associations between nine epigenetic age acceleration (EAA) scores and 18 cardiometabolic phenotypes using an Eastern European ageing population cohort richly annotated for a diverse set of phenotypes (subsample, n = 306; aged 45-69 years). This was implemented by splitting the data into groups with positive and negative EAAs. We observed strong association between all EAA scores and sex, suggesting that any analysis of EAAs should be adjusted by sex. We found that some sex-adjusted EAA scores were significantly associated with several phenotypes such as blood levels of gamma-glutamyl transferase and low-density lipoprotein, smoking status, annual alcohol consumption, multiple carotid plaques, and incident coronary heart disease status (not necessarily the same phenotypes for different EAAs). We demonstrated that even after adjusting EAAs for sex, EAA-phenotype associations remain sex-specific, which should be taken into account in any downstream analysis involving EAAs. The obtained results suggest that in some EAA-phenotype associations, negative EAA scores (i.e., epigenetic age below chronological age) indicated more harmful phenotype values, which is counterintuitive. Among all considered epigenetic clocks, GrimAge was significantly associated with more phenotypes than any other EA scores in this Russian sample.

Chemotherapy induces canalization of cell state in childhood B-cell precursor acute lymphoblastic leukemia.

Comparison of intratumor genetic heterogeneity in cancer at diagnosis and relapse suggests that chemotherapy induces bottleneck selection of subclonal genotypes. However, evolutionary events subsequent to chemotherapy could also explain changes in clonal dominance seen at relapse. We, therefore, investigated the mechanisms of selection in childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL) during induction chemotherapy where maximal cytoreduction occurs. To distinguish stochastic versus deterministic events, individual leukemias were transplanted into multiple xenografts and chemotherapy administered. Analyses of the immediate post-treatment leukemic residuum at single-cell resolution revealed that chemotherapy has little impact on genetic heterogeneity. Rather, it acts on extensive, previously unappreciated, transcriptional and epigenetic heterogeneity in BCP-ALL, dramatically reducing the spectrum of cell states represented, leaving a genetically polyclonal but phenotypically uniform population with hallmark signatures relating to developmental stage, cell cycle and metabolism. Hence, canalization of cell state accounts for a significant component of bottleneck selection during induction chemotherapy.

DNA Methylome Alterations Are Associated with Airway Macrophage Differentiation and Phenotype during Lung Fibrosis.

Rationale: Airway macrophages (AMs) are key regulators of the lung environment and are implicated in the pathogenesis of idiopathic pulmonary fibrosis (IPF), a fatal respiratory disease with no cure. However, knowledge about the epigenetics of AMs in IPF is limited. Objectives: To assess the role of epigenetic regulation of AMs during lung fibrosis. Methods: We undertook DNA methylation (DNAm) profiling by using Illumina EPIC (850k) arrays in sorted AMs from healthy donors (n = 14) and donors with IPF (n = 30). Cell-type deconvolution was performed by using reference myeloid-cell DNA methylomes. Measurements and Main Results: Our analysis revealed that epigenetic heterogeneity was a key characteristic of IPF AMs. DNAm "clock" analysis indicated that epigenetic alterations in IPF AMs were not associated with accelerated aging. In differential DNAm analysis, we identified numerous differentially methylated positions (n = 11) and differentially methylated regions (n = 49) between healthy and IPF AMs, respectively. Differentially methylated positions and differentially methylated regions encompassed genes involved in lipid (LPCAT1 [lysophosphatidylcholine acyltransferase 1]) and glucose (PFKFB3 [6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3]) metabolism, and importantly, the DNAm status was associated with disease severity in IPF. Conclusions: Collectively, our data identify that changes in the epigenome are associated with the development and function of AMs in the IPF lung.

DNA methylation-based profiling of bone and soft tissue tumours: a validation study of the 'DKFZ Sarcoma Classifier'.

Diagnosing bone and soft tissue neoplasms remains challenging because of the large number of subtypes, many of which lack diagnostic biomarkers. DNA methylation profiles have proven to be a reliable basis for the classification of brain tumours and, following this success, a DNA methylation-based sarcoma classification tool from the Deutsches Krebsforschungszentrum (DKFZ) in Heidelberg has been developed. In this study, we assessed the performance of their classifier on DNA methylation profiles of an independent data set of 986 bone and soft tissue tumours and controls. We found that the 'DKFZ Sarcoma Classifier' was able to produce a diagnostic prediction for 55% of the 986 samples, with 83% of these predictions concordant with the histological diagnosis. On limiting the validation to the 820 cases with histological diagnoses for which the DKFZ Classifier was trained, 61% of cases received a prediction, and the histological diagnosis was concordant with the predicted methylation class in 88% of these cases, findings comparable to those reported in the DKFZ Classifier paper. The classifier performed best when diagnosing mesenchymal chondrosarcomas (CHSs, 88% sensitivity), chordomas (85% sensitivity), and fibrous dysplasia (83% sensitivity). Amongst the subtypes least often classified correctly were clear cell CHSs (14% sensitivity), malignant peripheral nerve sheath tumours (27% sensitivity), and pleomorphic liposarcomas (29% sensitivity). The classifier predictions resulted in revision of the histological diagnosis in six of our cases. We observed that, although a higher tumour purity resulted in a greater likelihood of a prediction being made, it did not correlate with classifier accuracy. Our results show that the DKFZ Classifier represents a powerful research tool for exploring the pathogenesis of sarcoma; with refinement, it has the potential to be a valuable diagnostic tool.

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.

Glioblastomas acquire myeloid-affiliated transcriptional programs via epigenetic immunoediting to elicit immune evasion.

Glioblastoma multiforme (GBM) is an aggressive brain tumor for which current immunotherapy approaches have been unsuccessful. Here, we explore the mechanisms underlying immune evasion in GBM. By serially transplanting GBM stem cells (GSCs) into immunocompetent hosts, we uncover an acquired capability of GSCs to escape immune clearance by establishing an enhanced immunosuppressive tumor microenvironment. Mechanistically, this is not elicited via genetic selection of tumor subclones, but through an epigenetic immunoediting process wherein stable transcriptional and epigenetic changes in GSCs are enforced following immune attack. These changes launch a myeloid-affiliated transcriptional program, which leads to increased recruitment of tumor-associated macrophages. Furthermore, we identify similar epigenetic and transcriptional signatures in human mesenchymal subtype GSCs. We conclude that epigenetic immunoediting may drive an acquired immune evasion program in the most aggressive mesenchymal GBM subtype by reshaping the tumor immune microenvironment.

EPISCORE: cell type deconvolution of bulk tissue DNA methylomes from single-cell RNA-Seq data.

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.

Epigenome-Wide Association Study of Incident Type 2 Diabetes in a British Population: EPIC-Norfolk Study.

Epigenetic changes may contribute substantially to risks of diseases of aging. Previous studies reported seven methylation variable positions (MVPs) robustly associated with incident type 2 diabetes mellitus (T2DM). However, their causal roles in T2DM are unclear. In an incident T2DM case-cohort study nested within the population-based European Prospective Investigation into Cancer and Nutrition (EPIC)-Norfolk cohort, we used whole blood DNA collected at baseline, up to 11 years before T2DM onset, to investigate the role of methylation in the etiology of T2DM. We identified 15 novel MVPs with robust associations with incident T2DM and robustly confirmed three MVPs identified previously (near to TXNIP, ABCG1, and SREBF1). All 18 MVPs showed directionally consistent associations with incident and prevalent T2DM in independent studies. Further conditional analyses suggested that the identified epigenetic signals appear related to T2DM via glucose and obesity-related pathways acting before the collection of baseline samples. We integrated genome-wide genetic data to identify methylation-associated quantitative trait loci robustly associated with 16 of the 18 MVPs and found one MVP, cg00574958 at CPT1A, with a possible direct causal role in T2DM. None of the implicated genes were previously highlighted by genetic association studies, suggesting that DNA methylation studies may reveal novel biological mechanisms involved in tissue responses to glycemia.

Neoantigen-directed immune escape in lung cancer evolution.

The interplay between an evolving cancer and a dynamic immune microenvironment remains unclear. Here we analyse 258 regions from 88 early-stage, untreated non-small-cell lung cancers using RNA sequencing and histopathology-assessed tumour-infiltrating lymphocyte estimates. Immune infiltration varied both between and within tumours, with different mechanisms of neoantigen presentation dysfunction enriched in distinct immune microenvironments. Sparsely infiltrated tumours exhibited a waning of neoantigen editing during tumour evolution, indicative of historical immune editing, or copy-number loss of previously clonal neoantigens. Immune-infiltrated tumour regions exhibited ongoing immunoediting, with either loss of heterozygosity in human leukocyte antigens or depletion of expressed neoantigens. We identified promoter hypermethylation of genes that contain neoantigenic mutations as an epigenetic mechanism of immunoediting. Our results suggest that the immune microenvironment exerts a strong selection pressure in early-stage, untreated non-small-cell lung cancers that produces multiple routes to immune evasion, which are clinically relevant and forecast poor disease-free survival.