Search

Abstract Background DNA methylation showed notable potential to act as a diagnostic marker in many cancers. Many studies proposed DNA methylation biomarker in OSCC detection, while most of these studies are limited to specific cohorts or geographical location. However, the generalizability of DNA methylation as a diagnostic marker in oral cancer across different geographical locations is yet to be investigated. Methods We used genome-wide methylation data from 384 oral cavity cancer and normal tissues from TCGA HNSCC and eastern India. The common differentially methylated CpGs in these two cohorts were used to develop an Elastic-net model that can be used for the diagnosis of OSCC. The model was validated using 812 HNSCC and normal samples from different anatomical sites of oral cavity from seven countries. Droplet Digital PCR of methyl-sensitive restriction enzyme digested DNA (ddMSRE) was used for quantification of methylation and validation of the model with 22 OSCC and 22 contralateral normal samples. Additionally, pyrosequencing was used to validate the model using 46 OSCC and 25 adjacent normal and 21 contralateral normal tissue samples. Results With ddMSRE, our model showed 91% sensitivity, 100% specificity, and 95% accuracy in classifying OSCC from the contralateral normal tissues. Validation of the model with pyrosequencing also showed 96% sensitivity, 91% specificity, and 93% accuracy for classifying the OSCC from contralateral normal samples, while in case of adjacent normal samples we found similar sensitivity but with 20% specificity, suggesting the presence of early disease methylation signature at the adjacent normal samples. Methylation array data of HNSCC and normal tissues from different geographical locations and different anatomical sites showed comparable sensitivity, specificity, and accuracy in detecting oral cavity cancer with across. Similar results were also observed for different stages of oral cavity cancer. Conclusions Our model identified crucial genomic regions affected by DNA methylation in OSCC and showed similar accuracy in detecting oral cancer across different geographical locations. The high specificity of this model in classifying contralateral normal samples from the oral cancer compared to the adjacent normal samples suggested applicability of the model in early detection.

Abstract Background Restraining or slowing ageing hallmarks at the cellular level have been proposed as a route to increased organismal lifespan and healthspan. Consequently, there is great interest in anti-ageing drug discovery. However, this currently requires laborious and lengthy longevity analysis. Here, we present a novel screening readout for the expedited discovery of compounds that restrain ageing of cell populations in vitro and enable extension of in vivo lifespan. Methods Using Illumina methylation arrays, we monitored DNA methylation changes accompanying long-term passaging of adult primary human cells in culture. This enabled us to develop, test, and validate the CellPopAge Clock, an epigenetic clock with underlying algorithm, unique among existing epigenetic clocks for its design to detect anti-ageing compounds in vitro. Additionally, we measured markers of senescence and performed longevity experiments in vivo in Drosophila, to further validate our approach to discover novel anti-ageing compounds. Finally, we bench mark our epigenetic clock with other available epigenetic clocks to consolidate its usefulness and specialisation for primary cells in culture. Results We developed a novel epigenetic clock, the CellPopAge Clock, to accurately monitor the age of a population of adult human primary cells. We find that the CellPopAge Clock can detect decelerated passage-based ageing of human primary cells treated with rapamycin or trametinib, well-established longevity drugs. We then utilise the CellPopAge Clock as a screening tool for the identification of compounds which decelerate ageing of cell populations, uncovering novel anti-ageing drugs, torin2 and dactolisib (BEZ-235). We demonstrate that delayed epigenetic ageing in human primary cells treated with anti-ageing compounds is accompanied by a reduction in senescence and ageing biomarkers. Finally, we extend our screening platform in vivo by taking advantage of a specially formulated holidic medium for increased drug bioavailability in Drosophila. We show that the novel anti-ageing drugs, torin2 and dactolisib (BEZ-235), increase longevity in vivo. Conclusions Our method expands the scope of CpG methylation profiling to accurately and rapidly detecting anti-ageing potential of drugs using human cells in vitro, and in vivo, providing a novel accelerated discovery platform to test sought after anti-ageing compounds and geroprotectors.

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

Asteraceae is the world’s richest plant family and is found on all continents, in environments ranging from the coast to the highest mountains. The family shows all growth forms and, as in other angiosperm families, species richness is concentrated in tropical regions. South America has the highest diversity of Asteraceae in the world, yet taxonomic and distributional knowledge gaps remain. This study compiles an updated catalog of Asteraceae native to South America, based on national and regional checklists and ongoing large-scale flora projects. The resulting checklist includes a total of 6,940 species and 564 genera native to South America to date, which represent about a quarter of the family’s global diversity. Countries already considered to be megadiverse show the greatest diversity, such as Brazil with 2,095 species, followed by Peru (1,588), Argentina (1,377), and Colombia (1,244), with this diversity mainly focused on the Brazilian Highlands and the Andes. Species endemism also peaks in Brazil, but Sørensen distances reveal the Chilean flora to be eminently different from the rest of the continent. Tribes better represented in the continent are Eupatorieae, Senecioneae and Astereae, also with a remarkably presence of entirely South American subfamilies representing earliest diverging lineages of the Asteraceae, such as Barnadesioideae, Wunderlichioideae, Famatinanthoideae, and Stifftioideae. It is estimated that the discovery and description curves have not yet stabilized, and the number of species is likely to increase by 5 to 10% in the coming years, posing major challenges to continental-scale conservation.

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

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.

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

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.

Purpose: Rare disease diagnosis is challenging in medical image-based artificial intelligence due to a natural class imbalance in datasets, leading to biased prediction models. Inherited retinal diseases (IRDs) are a research domain that particularly faces this issue. This study investigates the applicability of synthetic data in improving artificial intelligence-enabled diagnosis of IRDs using generative adversarial networks (GANs). Design: Diagnostic study of gene-labeled fundus autofluorescence (FAF) IRD images using deep learning. Participants: Moorfields Eye Hospital (MEH) dataset of 15 692 FAF images obtained from 1800 patients with confirmed genetic diagnosis of 1 of 36 IRD genes. Methods: A StyleGAN2 model is trained on the IRD dataset to generate 512 × 512 resolution images. Convolutional neural networks are trained for classification using different synthetically augmented datasets, including real IRD images plus 1800 and 3600 synthetic images, and a fully rebalanced dataset. We also perform an experiment with only synthetic data. All models are compared against a baseline convolutional neural network trained only on real data. Main Outcome Measures: We evaluated synthetic data quality using a Visual Turing Test conducted with 4 ophthalmologists from MEH. Synthetic and real images were compared using feature space visualization, similarity analysis to detect memorized images, and Blind/Referenceless Image Spatial Quality Evaluator (BRISQUE) score for no-reference-based quality evaluation. Convolutional neural network diagnostic performance was determined on a held-out test set using the area under the receiver operating characteristic curve (AUROC) and Cohen’s Kappa (κ). Results: An average true recognition rate of 63% and fake recognition rate of 47% was obtained from the Visual Turing Test. Thus, a considerable proportion of the synthetic images were classified as real by clinical experts. Similarity analysis showed that the synthetic images were not copies of the real images, indicating that copied real images, meaning the GAN was able to generalize. However, BRISQUE score analysis indicated that synthetic images were of significantly lower quality overall than real images (P < 0.05). Comparing the rebalanced model (RB) with the baseline (R), no significant change in the average AUROC and κ was found (R-AUROC = 0.86[0.85-88], RB-AUROC = 0.88[0.86-0.89], R-κ = 0.51[0.49-0.53], and RB-κ = 0.52[0.50-0.54]). The synthetic data trained model (S) achieved similar performance as the baseline (S-AUROC = 0.86[0.85-87], S-κ = 0.48[0.46-0.50]). Conclusions: Synthetic generation of realistic IRD FAF images is feasible. Synthetic data augmentation does not deliver improvements in classification performance. However, synthetic data alone deliver a similar performance as real data, and hence may be useful as a proxy to real data.Financial Disclosure(s): Proprietary or commercial disclosure may be found after the references.

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

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

Chronic kidney disease (CKD) related cardiovascular disease (CVD) is characterized by vascular remodelling with well-established structural and functional changes in the vascular wall such as arterial stiffness, matrix deposition, and calcification. These phenotypic changes resemble pathology seen in ageing, and are likely to be mediated by sustained alterations in gene expression, which may be caused by epigenetic changes such as tissue-specific DNA methylation. We aimed to investigate tissue specific changes in DNA methylation that occur in CKD-related CVD. Genome-wide DNA methylation changes were examined in bisulphite converted genomic DNA isolated from the vascular media of CKD and healthy arteries. Methylation-specific PCR was used to validate the array data, and the association between DNA methylation and gene and protein expression was examined. The DNA methylation age was compared to the chronological age in both cases and controls. Three hundred and nineteen differentially methylated regions (DMR) were identified spread across the genome. Pathway analysis revealed that DMRs associated with genes were involved in embryonic and vascular development, and signalling pathways such as TGFβ and FGF. Expression of top differentially methylated gene HOXA5 showed a significant negative correlation with DNA methylation. Interestingly, DNA methylation age and chronological age were highly correlated, but there was no evidence of accelerated age-related DNA methylation in the arteries of CKD patients. In conclusion, we demonstrated that differential DNA methylation in the arterial tissue of CKD patients represents a potential mediator of arterial pathology and may be used to uncover novel pathways in the genesis of CKD-associated complications.

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

Reliably documenting plant diversity is necessary to protect and sustainably benefit from it. At the heart of this documentation lie species concepts and the practical methods used to delimit taxa. Here, we apply a total-evidence, iterative methodology to delimit and document species in the South American genus Victoria (Nymphaeaceae). The systematics of Victoria has thus far been poorly characterized due to difficulty in attributing species identities to biological collections. This research gap stems from an absence of type material and biological collections, also the confused diagnosis of V. cruziana. With the goal of improving systematic knowledge of the genus, we compiled information from historical records, horticulture and geography and assembled a morphological dataset using citizen science and specimens from herbaria and living collections. Finally, we generated genomic data from a subset of these specimens. Morphological and geographical observations suggest four putative species, three of which are supported by nuclear population genomic and plastid phylogenomic inferences. We propose these three confirmed entities as robust species, where two correspond to the currently recognized V. amazonica and V. cruziana, the third being new to science, which we describe, diagnose and name here as V. boliviana Magdalena and L. T. Sm. Importantly, we identify new morphological and molecular characters which serve to distinguish the species and underpin their delimitations. Our study demonstrates how combining different types of character data into a heuristic, total-evidence approach can enhance the reliability with which biological diversity of morphologically challenging groups can be identified, documented and further studied.

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

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

Abstract Epigenetic clocks comprise a set of CpG sites whose DNA methylation levels measure subject age. These clocks are acknowledged as a highly accurate molecular correlate of chronological age in humans and other vertebrates. Also, extensive research is aimed at their potential to quantify biological aging rates and test longevity or rejuvenating interventions. Here, we discuss key challenges to understand clock mechanisms and biomarker utility. This requires dissecting the drivers and regulators of age-related changes in single-cell, tissue- and disease-specific models, as well as exploring other epigenomic marks, longitudinal and diverse population studies, and non-human models. We also highlight important ethical issues in forensic age determination and predicting the trajectory of biological aging in an individual.

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.

Abstract As epigenetic studies become more common and lead to new insights into health and disease, the return of individual epigenetic results to research participants, in particular in large-scale epigenomic studies, will be of growing importance. Members of the International Human Epigenome Consortium (IHEC) Bioethics Workgroup considered the potential ethical, legal, and social issues (ELSI) involved in returning epigenetic research results and incidental findings in order to produce a set of ‘Points-to-consider’ (P-t-C) for the epigenetics research community. These P-t-C draw on existing guidance on the return of genetic research results, while also integrating the IHEC Bioethics Workgroup’s ELSI research on and discussion of the issues associated with epigenetic data as well as the experience of a return of results pilot study by the Personal Genome Project UK (PGP-UK). Major challenges include how to determine the clinical validity and actionability of epigenetic results, and considerations related to environmental exposures and epigenetic marks, including circumstances warranting the sharing of results with family members and third parties. Interdisciplinary collaboration and good public communication regarding epigenetic risk will be important to advance the return of results framework for epigenetic science.

Abstract Background Due to limited access to brain tissue, the precise mechanisms underlying neuro-axonal dysfunction in neurological disorders such as multiple sclerosis (MS) are largely unknown. In that context, profiling DNA methylation, which is a stable and cell type-specific regulatory epigenetic mark of genome activity, offers a unique opportunity to characterize the molecular mechanisms underpinning brain pathology in situ. We examined DNA methylation patterns of neuronal nuclei isolated from post-mortem brain tissue to infer processes that occur in neurons of MS patients. Results We isolated subcortical neuronal nuclei from post-mortem white matter tissue of MS patients and non-neurological controls using flow cytometry. We examined bulk DNA methylation changes (total n = 29) and further disentangled true DNA methylation (5mC) from neuron-specific DNA hydroxymethylation (5hmC) (n = 17), using Illumina Infinium 450K arrays. We performed neuronal sub-type deconvolution using glutamate and GABA methylation profiles to further reduce neuronal sample heterogeneity. In total, we identified 2811 and 1534 significant (genome-wide adjusted P value

In recent years, there has been a significant increase in whole genome sequencing data of individual genomes produced by research projects as well as direct to consumer service providers. While many of these sources provide their users with an interpretation of the data, there is a lack of free, open tools for generating reports exploring the data in an easy to understand manner. GenomeChronicler was developed as part of the Personal Genome Project UK (PGP-UK) to address this need. PGP-UK provides genomic, transcriptomic, epigenomic and self-reported phenotypic data under an open-access model with full ethical approval. As a result, the reports generated by GenomeChronicler are intended for research purposes only and include information relating to potentially beneficial and potentially harmful variants, but without clinical curation. GenomeChronicler can be used with data from whole genome or whole exome sequencing, producing a genome report containing information on variant statistics, ancestry and known associated phenotypic traits. Example reports are available from the PGP-UK data page (personalgenomes.org.uk/data). The objective of this method is to leverage existing resources to find known phenotypes associated with the genotypes detected in each sample. The provided trait data is based primarily upon information available in SNPedia, but also collates data from ClinVar, GETevidence, and gnomAD to provide additional details on potential health implications, presence of genotype in other PGP participants and population frequency of each genotype. The analysis can be run in a self-contained environment without requiring internet access, making it a good choice for cases where privacy is essential or desired: any third party project can embed GenomeChronicler within their off-line safe-haven environments. GenomeChronicler can be run for one sample at a time, or in parallel making use of the Nextflow workflow manager. The source code is available from GitHub (https://github.com/PGP-UK/GenomeChronicler), container recipes are available for Docker and Singularity, as well as a pre-built container from SingularityHub (https://singularity-hub.org/collections/3664) enabling easy deployment in a variety of settings. Users without access to computational resources to run GenomeChronicler can access the software from the Lifebit CloudOS platform (https://lifebit.ai/cloudos) enabling the production of reports and variant calls from raw sequencing data in a scalable fashion.

ABSTRACT: Background: Alopecia areata (AA) is considered a highly heritable, T-cell-mediated autoimmune disease of the hair follicle. However, no convincing susceptibility gene has yet been pinpointed in the major histocompatibility complex (MHC), a genome region known to be associated with AA as compared to other regions. Methods: We engineered mice carrying AA risk allele identified by haplotype sequencing for the MHC region using allele-specific genome editing with the CRISPR/Cas9 system. Finally, we performed functional evaluations in the mice and AA patients with and without the risk allele. Findings: We identified a variant (rs142986308, p.Arg587Trp) in the coiled-coil alpha-helical rod protein 1 (CCHCR1) gene as the only non-synonymous variant in the AA risk haplotype. Furthermore, mice engineered to carry the risk allele displayed a hair loss phenotype. Transcriptomics further identified CCHCR1 as a novel component interacting with hair cortex keratin in hair shafts. Both, these alopecic mice and AA patients with the risk allele displayed morphologically impaired hair and comparable differential expression of hair-related genes, including hair keratin and keratin-associated proteins (KRTAPs). Interpretation: Our results implicate CCHCR1 with the risk allele in a previously unidentified subtype of AA based on aberrant keratinization in addition to autoimmune events. Funding: This work was supported by JSPS KAKENHI (JP16K10177) and the NIHR UCLH Biomedical Research center (BRC84/CN/SB/5984).

Summary: Although thousands of solid tumors have been sequenced to date, a fundamental under-sampling bias is inherent in current methodologies. This is caused by a tissue sample input of fixed dimensions (e.g., 6 mm biopsy), which becomes grossly under-powered as tumor volume scales. Here, we demonstrate representative sequencing (Rep-Seq) as a new method to achieve unbiased tumor tissue sampling. Rep-Seq uses fixed residual tumor material, which is homogenized and subjected to next-generation sequencing. Analysis of intratumor tumor mutation burden (TMB) variability shows a high level of misclassification using current single-biopsy methods, with 20% of lung and 52% of bladder tumors having at least one biopsy with high TMB but low clonal TMB overall. Misclassification rates by contrast are reduced to 2% (lung) and 4% (bladder) when a more representative sampling methodology is used. Rep-Seq offers an improved sampling protocol for tumor profiling, with significant potential for improved clinical utility and more accurate deconvolution of clonal structure.

Abstract Background Radiation therapy is an indispensable part of various treatment modalities for breast cancer. Specifically, for non-inflammatory locally advanced breast cancer (LABC) patients, preoperative radiotherapy (pRT) is currently indicated as a second line therapy in the event of lack of response to neoadjuvant chemotherapy. Still approximately one third of patients fails to respond favourably to pRT. The aim of this study was to explore molecular mechanisms underlying differential response to radiotherapy (RT) to identify predictive biomarkers and potential targets for increasing radiosensitivity. Methods The study was based on a cohort of 134 LABC patients, treated at the Institute of Oncology and Radiology of Serbia (IORS) with pRT, without previous or concomitant systemic therapy. Baseline transcriptional profiles were established using Agilent 60 K microarray platform in a subset of 23 formalin-fixed paraffin-embedded (FFPE) LABC tumour samples of which 11 radiotherapy naïve and 3 post-radiotherapy samples passed quality control and were used for downstream analysis. Biological networks and signalling pathways underlying differential response to RT were identified using Ingenuity Pathways Analysis software. Predictive value of candidate genes in the preoperative setting was further validated by qRT-PCR in an independent subset of 60 LABC samples of which 42 had sufficient quality for data analysis, and in postoperative setting using microarray data from 344 node-negative breast cancer patients (Erasmus cohort, GSE2034 and GSE5327) treated either with surgery only (20%) or surgery with RT (80%). Results We identified 192 significantly differentially expressed genes (FDR

Age-related changes in DNA methylation were observed in cross-sectional studies, but longitudinal evidence is still limited. Here, we aimed to characterize longitudinal age-related methylation patterns using 1011 blood samples collected from 385 Swedish twins (age at entry: mean 69 and standard deviation 9.7, 73 monozygotic and 96 dizygotic pairs) up to five times (mean 2.6) over 20 years (mean 8.7). We identified 1316 age-associated methylation sites (P

Abstract Background Rheumatoid arthritis is a common autoimmune disorder influenced by both genetic and environmental factors. Epigenome-wide association studies can identify environmentally mediated epigenetic changes such as altered DNA methylation, which may also be influenced by genetic factors. To investigate possible contributions of DNA methylation to the aetiology of rheumatoid arthritis with minimum confounding genetic heterogeneity, we investigated genome-wide DNA methylation in disease-discordant monozygotic twin pairs. Methods Genome-wide DNA methylation was assessed in 79 monozygotic twin pairs discordant for rheumatoid arthritis using the HumanMethylation450 BeadChip array (Illumina). Discordant twins were tested for both differential DNA methylation and methylation variability between rheumatoid arthritis and healthy twins. The methylation variability signature was then compared with methylation variants from studies of other autoimmune diseases and with an independent healthy population. Results We have identified a differentially variable DNA methylation signature that suggests multiple stress response pathways may be involved in the aetiology of the disease. This methylation variability signature also highlighted potential epigenetic disruption of multiple RUNX3 transcription factor binding sites as being associated with disease development. Comparison with previously performed epigenome-wide association studies of rheumatoid arthritis and type 1 diabetes identified shared pathways for autoimmune disorders, suggesting that epigenetics plays a role in autoimmunity and offering the possibility of identifying new targets for intervention. Conclusions Through genome-wide analysis of DNA methylation in disease-discordant monozygotic twins, we have identified a differentially variable DNA methylation signature, in the absence of differential methylation in rheumatoid arthritis. This finding supports the importance of epigenetic variability as an emerging component in autoimmune disorders.

Abstract DNA methylation is the most studied epigenetic modification due to its role in regulating gene expression, and its involvement in the pathogenesis of cancer and several diseases upon aberrations in methylation. The method of choice to evaluate genome-wide methylation has been the Illumina HumanMethylation450 BeadChip (450K), but it was recently replaced with the MethylationEPIC BeadChip (EPIC). We therefore sought to validate the EPIC array in comparison to the 450K array for both fresh-frozen (FF) and formalin-fixed paraffin-embedded (FFPE) tumours. We also performed analysis on the EPIC array with paired FF and FFPE samples to adapt to a clinical setting where FFPE is routinely used. Further, we compared two restoration methods, REPLI-g and Infinium, for FFPE-derived DNA on the EPIC array. The Pearson correlation of β values for common probes on the 450K and EPIC array was high for both FF (mean: 0.992) and FFPE (mean: 0.984) samples. The β values generated from the EPIC array for FFPE samples correlated well with the paired FF tumours, but varied between 0.901 and 0.987. We did note that sample pairs with lower correlation had less bimodal density distributions of β values and displayed higher noise in the copy number alteration plots (generated from the methylation array data) in the FFPE sample. Both REPLI-g and the Infinium restoration for FFPE samples performed well on the EPIC array and generated equivalent correlation scores to the paired FF sample.

The heterogeneity of colorectal cancer has important clinical and therapeutic implications. Here the authors analysed the responses of a large biobank of organoids and xenografts derived from colorectal patients to a panel of clinically relevant therapeutic agents to identify genes signatures associated with drug response.

The incidence of type 1 diabetes is increasing, potentially implicating non-genetic factors. Here the authors conduct an epigenome-wide association study in disease-discordant twins and find increased DNA methylation variability at genes associated with immune cell metabolism and the cell cycle.

Abstract Stephan Beck discusses recent developments in sharing personal genomes as part of the Personal Genome Project in the UK and globally, and how these efforts are advancing research.

Here, Libertini and colleagues devise a computation tool that can analyze whole-genome bisulfite sequencing (WGBS) data to recover of ∼30% of the lost differential methylation position information. They use COMETgazer and COMETvintage to analyze 13 diffferent methylome data to demonstrate their performance.

Women with germline variants in BRCA genes are predisposed to ovarian cancer. In this study, the authors demonstrate that fimbrial tissue from the ovary, the site of ovarian cancer, in BRCAmutant carriers contains marked DNA methylation changes compared with the proximal region of the ovary.

Differential impact of genetic and environmental influences on DNA methylation may result in sex- and age-related physiological variation and disease susceptibility. By analysing DNA methylome of 2,603 individuals from twin families, here, the authors establish a catalogue of between-individual variation in DNA methylation.

The methylation of non-CpG residues is a poorly understood marker of pluripotent cells, gradually lost as cells differentiate. Here the authors show non-CG methylation can be used as a marker of differentiation potential.

Glioblastoma (GBM) is an aggressive brain tumor whose growth is driven by stem cell-like cells. BMP signaling triggers cell-cycle exit and differentiation of GBM stem cells (GSCs) and, therefore, might have therapeutic value. However, the epigenetic mechanisms that accompany differentiation remain poorly defined. It is also unclear whether cell-cycle arrest is terminal. Here we find only a subset of GSC cultures exhibit astrocyte differentiation in response to BMP. Although overtly differentiated non-cycling astrocytes are generated, they remain vulnerable to cell-cycle re-entry and fail to appropriately reconfigure DNA methylation patterns. Chromatin accessibility mapping identified loci that failed to alter in response to BMP and these were enriched in SOX transcription factor-binding motifs. SOX transcription factors, therefore, may limit differentiation commitment. A similar propensity for cell-cycle re-entry and de-differentiation was observed in GSC-derived oligodendrocyte-like cells. These findings highlight significant obstacles to BMP-induced differentiation as therapy for GBM.

Human genetic studies have emphasised the dominant contribution of pancreatic islet dysfunction to development of Type 2 Diabetes (T2D). However, limited annotation of the islet epigenome has constrained efforts to define the molecular mechanisms mediating the, largely regulatory, signals revealed by Genome-Wide Association Studies (GWAS). We characterised patterns of chromatin accessibility (ATAC-seq, n = 17) and DNA methylation (whole-genome bisulphite sequencing, n = 10) in human islets, generating high-resolution chromatin state maps through integration with established ChIP-seq marks. We found enrichment of GWAS signals for T2D and fasting glucose was concentrated in subsets of islet enhancers characterised by open chromatin and hypomethylation, with the former annotation predominant. At several loci (including CDC123, ADCY5, KLHDC5) the combination of fine-mapping genetic data and chromatin state enrichment maps, supplemented by allelic imbalance in chromatin accessibility pinpointed likely causal variants. The combination of increasingly-precise genetic and islet epigenomic information accelerates definition of causal mechanisms implicated in T2D pathogenesis.

Following publication of the original article [1], the authors reported an error in Additional file 1.

Recent studies have demonstrated that the DNA methylome changes with age. This epigenetic drift may have deep implications for cellular differentiation and disease development. However, it remains unclear how much of this drift is functional or caused by underlying changes in cell subtype composition. Moreover, no study has yet comprehensively explored epigenetic drift at different genomic length scales and in relation to regulatory elements. Here we conduct an in-depth analysis of epigenetic drift in blood tissue. We demonstrate that most of the age-associated drift is independent of the increase in the granulocyte to lymphocyte ratio that accompanies aging and that enrichment of age-hypermethylated CpG islands increases upon adjustment for cellular composition. We further find that drift has only a minimal impact on in-cis gene expression, acting primarily to stabilize pre-existing baseline expression levels. By studying epigenetic drift at different genomic length scales, we demonstrate the existence of mega-base scale age-associated hypomethylated blocks, covering approximately 14% of the human genome, and which exhibit preferential hypomethylation in age-matched cancer tissue. Importantly, we demonstrate the feasibility of integrating Illumina 450k DNA methylation with ENCODE data to identify transcription factors with key roles in cellular development and aging. Specifically, we identify REST and regulatory factors of the histone methyltransferase MLL complex, whose function may be disrupted in aging. In summary, most of the epigenetic drift seen in blood is independent of changes in blood cell type composition, and exhibits patterns at different genomic length scales reminiscent of those seen in cancer. Integration of Illumina 450k with appropriate ENCODE data may represent a fruitful approach to identify transcription factors with key roles in aging and disease.

The Infinium 450K Methylation array is an established tool for measuring methylation. However, the bisulfite (BS) reaction commonly used with the 450K array cannot distinguish between 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC). The oxidative-bisulfite assay disambiguates 5mC and 5hmC. We describe the use of oxBS in conjunction with the 450K array (oxBS-array) to analyse 5hmC/5mC in cerebellum DNA. The "methylation" level derived by the BS reaction is the combined level of 5mC and 5hmC at a given base, while the oxBS reaction gives the level of 5mC alone. The level of 5hmC is derived by subtracting the oxBS level from the BS level. Here we present an analysis method that distinguishes genuine positive levels of 5hmC at levels as low as 3%. We performed four replicates of the same sample of cerebellum and found a high level of reproducibility (average r for BS = 98.3, and average r for oxBS = 96.8). In total, 114,734 probes showed a significant positive measurement for 5hmC. The range at which we were able to distinguish 5hmC occupancy was between 3% and 42%. In order to investigate the effects of multiple replicates on 5hmC detection we also simulated fewer replicates and found that decreasing the number of replicates to two reduced the number of positive probes identified by > 50%. We validated our results using qPCR in conjunction with glucosylation of 5hmC sites followed by MspI digestion and we found good concordance with the array estimates (r = 0.94). This experiment provides a map of 5hmC in the cerebellum and a robust dataset for use as a standard in future 5hmC analyses. We also provide a novel method for validating the presence of 5hmC at low levels, and highlight some of the pitfalls associated with measuring 5hmC and 5mC.

Recently it has been observed that cancer tissue is characterised by an increased variability in DNA methylation patterns. However, how the correlative patterns in genome-wide DNA methylation change during the carcinogenic progress has not yet been explored. Here we study genome-wide inter-CpG correlations in DNA methylation, in addition to single site variability, during cervical carcinogenesis. We demonstrate how the study of changes in DNA methylation covariation patterns across normal, intra-epithelial neoplasia and invasive cancer allows the identification of CpG sites that indicate the risk of neoplastic transformation in stages prior to neoplasia. Importantly, we show that the covariation in DNA methylation at these risk CpG loci is maximal immediately prior to the onset of cancer, supporting the view that high epigenetic diversity in normal cells increases the risk of cancer. Consistent with this, we observe that invasive cancers exhibit increased covariation in DNA methylation at the risk CpG sites relative to normal tissue, but lower levels relative to pre-cancerous lesions. We further show that the identified risk CpG sites undergo preferential DNA methylation changes in relation to human papilloma virus infection and age. Results are validated in independent data including prospectively collected samples prior to neoplastic transformation. Our data are consistent with a phase transition model of carcinogenesis, in which epigenetic diversity is maximal prior to the onset of cancer. The model and algorithm proposed here may allow, in future, network biomarkers predicting the risk of neoplastic transformation to be identified.

Many genes important in immunity are found as multigene families. The butyrophilin genes are members of the B7 family, playing diverse roles in co-regulation and perhaps in antigen presentation. In humans, a fixed number of butyrophilin genes are found in and around the major histocompatibility complex (MHC), and show striking association with particular autoimmune diseases. In chickens, BG genes encode homologues with somewhat different domain organisation. Only a few BG genes have been characterised, one involved in actin-myosin interaction in the intestinal brush border, and another implicated in resistance to viral diseases. We characterise all BG genes in B12 chickens, finding a multigene family organised as tandem repeats in the BG region outside the MHC, a single gene in the MHC (the BF-BL region), and another single gene on a different chromosome. There is a precise cell and tissue expression for each gene, but overall there are two kinds, those expressed by haemopoietic cells and those expressed in tissues (presumably non-haemopoietic cells), correlating with two different kinds of promoters and 5' untranslated regions (5'UTR). However, the multigene family in the BG region contains many hybrid genes, suggesting recombination and/or deletion as major evolutionary forces. We identify BG genes in the chicken whole genome shotgun sequence, as well as by comparison to other haplotypes by fibre fluorescence in situ hybridisation, confirming dynamic expansion and contraction within the BG region. Thus, the BG genes in chickens are undergoing much more rapid evolution compared to their homologues in mammals, for reasons yet to be understood.

Schickendantziella trichosepala (Speg.) Speg., a new record from Bolivian fora. The taxonomic history, description, iconography and geographical distribution of this monotypic South American genus are presented.

Monozygotic (MZ) twin pair discordance for childhood-onset Type 1 Diabetes (T1D) is ∼50%, implicating roles for genetic and non-genetic factors in the aetiology of this complex autoimmune disease. Although significant progress has been made in elucidating the genetics of T1D in recent years, the non-genetic component has remained poorly defined. We hypothesized that epigenetic variation could underlie some of the non-genetic component of T1D aetiology and, thus, performed an epigenome-wide association study (EWAS) for this disease. We generated genome-wide DNA methylation profiles of purified CD14+ monocytes (an immune effector cell type relevant to T1D pathogenesis) from 15 T1D-discordant MZ twin pairs. This identified 132 different CpG sites at which the direction of the intra-MZ pair DNA methylation difference significantly correlated with the diabetic state, i.e. T1D-associated methylation variable positions (T1D-MVPs). We confirmed these T1D-MVPs display statistically significant intra-MZ pair DNA methylation differences in the expected direction in an independent set of T1D-discordant MZ pairs (P = 0.035). Then, to establish the temporal origins of the T1D-MVPs, we generated two further genome-wide datasets and established that, when compared with controls, T1D-MVPs are enriched in singletons both before (P = 0.001) and at (P = 0.015) disease diagnosis, and also in singletons positive for diabetes-associated autoantibodies but disease-free even after 12 years follow-up (P = 0.0023). Combined, these results suggest that T1D-MVPs arise very early in the etiological process that leads to overt T1D. Our EWAS of T1D represents an important contribution toward understanding the etiological role of epigenetic variation in type 1 diabetes, and it is also the first systematic analysis of the temporal origins of disease-associated epigenetic variation for any human complex disease.

Cellular differentiation involves widespread epigenetic reprogramming, including modulation of DNA methylation patterns. Using Differential Methylation Hybridization (DMH) in combination with a custom DMH array containing 51,243 features covering more than 16,000 murine genes, we carried out a genome-wide screen for cell- and tissue-specific differentially methylated regions (tDMRs) in undifferentiated embryonic stem cells (ESCs), in in-vitro induced neural stem cells (NSCs) and 8 differentiated embryonic and adult tissues. Unsupervised clustering of the generated data showed distinct cell- and tissue-specific DNA methylation profiles, revealing 202 significant tDMRs (p1.96) enrichment for genes involved in neural differentiation, including, for example, Jag1 and Tcf4. Our results provide robust evidence for the relevance of DNA methylation in early neural development and identify novel marker candidates for neural cell differentiation.

DNA methylation plays an important role in biological processes in human health and disease. Recent technological advances allow unbiased whole-genome DNA methylation (methylome) analysis to be carried out on human cells. Using whole-genome bisulfite sequencing at 24.7-fold coverage (12.3-fold per strand), we report a comprehensive (92.62%) methylome and analysis of the unique sequences in human peripheral blood mononuclear cells (PBMC) from the same Asian individual whose genome was deciphered in the YH project. PBMC constitute an important source for clinical blood tests world-wide. We found that 68.4% of CpG sites and 80% displayed allele-specific expression (ASE). These data demonstrate that ASM is a recurrent phenomenon and is highly correlated with ASE in human PBMCs. Together with recently reported similar studies, our study provides a comprehensive resource for future epigenomic research and confirms new sequencing technology as a paradigm for large-scale epigenomics studies.

Recent multi-dimensional approaches to the study of complex disease have revealed powerful insights into how genetic and epigenetic factors may underlie their aetiopathogenesis. We examined genotype-epigenotype interactions in the context of Type 2 Diabetes (T2D), focussing on known regions of genomic susceptibility. We assayed DNA methylation in 60 females, stratified according to disease susceptibility haplotype using previously identified association loci. CpG methylation was assessed using methylated DNA immunoprecipitation on a targeted array (MeDIP-chip) and absolute methylation values were estimated using a Bayesian algorithm (BATMAN). Absolute methylation levels were quantified across LD blocks, and we identified increased DNA methylation on the FTO obesity susceptibility haplotype, tagged by the rs8050136 risk allele A (p = 9.40×10(-4), permutation p = 1.0×10(-3)). Further analysis across the 46 kb LD block using sliding windows localised the most significant difference to be within a 7.7 kb region (p = 1.13×10(-7)). Sequence level analysis, followed by pyrosequencing validation, revealed that the methylation difference was driven by the co-ordinated phase of CpG-creating SNPs across the risk haplotype. This 7.7 kb region of haplotype-specific methylation (HSM), encapsulates a Highly Conserved Non-Coding Element (HCNE) that has previously been validated as a long-range enhancer, supported by the histone H3K4me1 enhancer signature. This study demonstrates that integration of Genome-Wide Association (GWA) SNP and epigenomic DNA methylation data can identify potential novel genotype-epigenotype interactions within disease-associated loci, thus providing a novel route to aid unravelling common complex diseases.

BACKGROUND:Recent studies have shown that DNA methylation (DNAm) markers in peripheral blood may hold promise as diagnostic or early detection/risk markers for epithelial cancers. However, to date no study has evaluated the diagnostic and predictive potential of such markers in a large case control cohort and on a genome-wide basis. PRINCIPAL FINDINGS:By performing genome-wide DNAm profiling of a large ovarian cancer case control cohort, we here demonstrate that active ovarian cancer has a significant impact on the DNAm pattern in peripheral blood. Specifically, by measuring the methylation levels of over 27,000 CpGs in blood cells from 148 healthy individuals and 113 age-matched pre-treatment ovarian cancer cases, we derive a DNAm signature that can predict the presence of active ovarian cancer in blind test sets with an AUC of 0.8 (95% CI (0.74-0.87)). We further validate our findings in another independent set of 122 post-treatment cases (AUC = 0.76 (0.72-0.81)). In addition, we provide evidence for a significant number of candidate risk or early detection markers for ovarian cancer. Furthermore, by comparing the pattern of methylation with gene expression data from major blood cell types, we here demonstrate that age and cancer elicit common changes in the composition of peripheral blood, with a myeloid skewing that increases with age and which is further aggravated in the presence of ovarian cancer. Finally, we show that most cancer and age associated methylation variability is found at CpGs located outside of CpG islands. SIGNIFICANCE:Our results underscore the potential of DNAm profiling in peripheral blood as a tool for detection or risk-prediction of epithelial cancers, and warrants further in-depth and higher CpG coverage studies to further elucidate this role.

There are two main classes of natural killer (NK) cell receptors in mammals, the killer cell immunoglobulin-like receptors (KIR) and the structurally unrelated killer cell lectin-like receptors (KLR). While KIR represent the most diverse group of NK receptors in all primates studied to date, including humans, apes, and Old and New World monkeys, KLR represent the functional equivalent in rodents. Here, we report a first digression from this rule in lemurs, where the KLR (CD94/NKG2) rather than KIR constitute the most diverse group of NK cell receptors. We demonstrate that natural selection contributed to such diversification in lemurs and particularly targeted KLR residues interacting with the peptide presented by MHC class I ligands. We further show that lemurs lack a strict ortholog or functional equivalent of MHC-E, the ligands of non-polymorphic KLR in "higher" primates. Our data support the existence of a hitherto unknown system of polymorphic and diverse NK cell receptors in primates and of combinatorial diversity as a novel mechanism to increase NK cell receptor repertoire.

Numerous genetic association studies have implicated the KIAA0319 gene on human chromosome 6p22 in dyslexia susceptibility. The causative variant(s) remains unknown but may modulate gene expression, given that (1) a dyslexia-associated haplotype has been implicated in the reduced expression of KIAA0319, and (2) the strongest association has been found for the region spanning exon 1 of KIAA0319. Here, we test the hypothesis that variant(s) responsible for reduced KIAA0319 expression resides on the risk haplotype close to the gene's transcription start site. We identified seven single-nucleotide polymorphisms on the risk haplotype immediately upstream of KIAA0319 and determined that three of these are strongly associated with multiple reading-related traits. Using luciferase-expressing constructs containing the KIAA0319 upstream region, we characterized the minimal promoter and additional putative transcriptional regulator regions. This revealed that the minor allele of rs9461045, which shows the strongest association with dyslexia in our sample (max p-value = 0.0001), confers reduced luciferase expression in both neuronal and non-neuronal cell lines. Additionally, we found that the presence of this rs9461045 dyslexia-associated allele creates a nuclear protein-binding site, likely for the transcriptional silencer OCT-1. Knocking down OCT-1 expression in the neuronal cell line SHSY5Y using an siRNA restores KIAA0319 expression from the risk haplotype to nearly that seen from the non-risk haplotype. Our study thus pinpoints a common variant as altering the function of a dyslexia candidate gene and provides an illustrative example of the strategic approach needed to dissect the molecular basis of complex genetic traits.