Saturation analysis for whole-genome bisulfite sequencing data

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

One-pot trimodal mapping of unmethylated, hydroxymethylated, and open chromatin sites unveils distinctive 5hmC roles at dynamic chromatin loci

We present a method, named Mx-TOP, for profiling of three epigenetic regulatory layers-chromatin accessibility, general DNA modification, and DNA hydroxymethylation-from a single library. The approach is based on chemo-enzymatic covalent tagging of unmodified CG sites and hydroxymethylated cytosine (5hmC) along with GC sites in chromatin, which are then mapped using tag-selective base-resolution TOP-seq sequencing. Our in-depth validation of the approach revealed its sensitivity and informativity in evaluating chromatin accessibility and DNA modification interactions that drive transcriptional regulation. We employed the technology in a study of chromatin and DNA demethylation dynamics during in vitro neuronal differentiation. The study highlighted the involvement of gene body 5hmC in modulating an extensive decoupling between promoter accessibility and transcription. The importance of 5hmC in chromatin remodeling was further demonstrated by the observed resistance of the developmentally acquired open loci to the global 5hmC erasure in neuronal progenitors.

Genetic impacts on DNA methylation: research findings and future perspectives

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

Hippocampal and cortical tissue-specific epigenetic clocks indicate an increased epigenetic age in a mouse model for Alzheimer's disease

Epigenetic clocks are based on age-associated changes in DNA methylation of CpG-sites, which can accurately measure chronological age in different species. Recently, several studies have indicated that the difference between chronological and epigenetic age, defined as the age acceleration, could reflect biological age indicating functional decline and age-associated diseases. In humans, an epigenetic clock associated Alzheimer's disease (AD) pathology with an acceleration of the epigenetic age. In this study, we developed and validated two mouse brain region-specific epigenetic clocks from the C57BL/6J hippocampus and cerebral cortex. Both clocks, which could successfully estimate chronological age, were further validated in a widely used mouse model for AD, the triple transgenic AD (3xTg-AD) mouse. We observed an epigenetic age acceleration indicating an increased biological age for the 3xTg-AD mice compared to non-pathological C57BL/6J mice, which was more pronounced in the cortex as compared to the hippocampus. Genomic region enrichment analysis revealed that age-dependent CpGs were enriched in genes related to developmental, aging-related, neuronal and neurodegenerative functions. Due to the limited access of human brain tissues, these epigenetic clocks specific for mouse cortex and hippocampus might be important in further unravelling the role of epigenetic mechanisms underlying AD pathology or brain aging in general.

DNA methylation study of Huntington's disease and motor progression in patients and in animal models

Although Huntington's disease (HD) is a well studied Mendelian genetic disorder, less is known about its associated epigenetic changes. Here, we characterize DNA methylation levels in six different tissues from 3 species: a mouse huntingtin (Htt) gene knock-in model, a transgenic HTT sheep model, and humans. Our epigenome-wide association study (EWAS) of human blood reveals that HD mutation status is significantly (p < 10^-7) associated with 33 CpG sites, including the HTT gene (p = 6.5 × 10^-26). These Htt/HTT associations were replicated in the Q175 Htt knock-in mouse model (p = 6.0 × 10^-8) and in the transgenic sheep model (p = 2.4 × 10^-88). We define a measure of HD motor score progression among manifest HD cases based on multiple clinical assessments. EWAS of motor progression in manifest HD cases exhibits significant (p < 10^-7) associations with methylation levels at three loci: near PEX14 (p = 9.3 × 10^-9), GRIK4 (p = 3.0 × 10^-8), and COX4I2 (p = 6.5 × 10^-8). We conclude that HD is accompanied by profound changes of DNA methylation levels in three mammalian species.

DNA methylation aging clocks: challenges and recommendations

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.

Systematic evaluation of library preparation methods and sequencing platforms for high-throughput whole genome bisulfite sequencing

Whole genome bisulfite sequencing (WGBS), with its ability to interrogate methylation status at single CpG site resolution epigenome-wide, is a powerful technique for use in molecular experiments. Here, we aim to advance strategies for accurate and efficient WGBS for application in future large-scale epidemiological studies. We systematically compared the performance of three WGBS library preparation methods with low DNA input requirement (Swift Biosciences Accel-NGS, Illumina TruSeq and QIAGEN QIAseq) on two state-of-the-art sequencing platforms (Illumina NovaSeq and HiSeq X), and also assessed concordance between data generated by WGBS and methylation arrays. Swift achieved the highest proportion of CpG sites assayed and effective coverage at 26x (P < 0.001). TruSeq suffered from the highest proportion of PCR duplicates, while QIAseq failed to deliver across all quality metrics. There was little difference in performance between NovaSeq and HiSeq X, with the exception of higher read duplication rate on the NovaSeq (P < 0.05), likely attributable to the higher cluster densities on its flow cells. Systematic biases exist between WGBS and methylation arrays, with lower precision observed for WGBS across the range of depths investigated. To achieve a level of precision broadly comparable to the methylation array, a minimum coverage of 100x is recommended.

Integrating Signals from Sperm Methylome Analysis and Genome-Wide Association Study for a Better Understanding of Male Fertility in Cattle

Decreased male fertility is a big concern in both human society and the livestock industry. Sperm DNA methylation is commonly believed to be associated with male fertility. However, due to the lack of accurate male fertility records (i.e., limited mating times), few studies have investigated the comprehensive impacts of sperm DNA methylation on male fertility in mammals. In this study, we generated 10 sperm DNA methylomes and performed a preliminary correlation analysis between signals from sperm DNA methylation and signals from large-scale (n = 27,214) genome-wide association studies (GWAS) of 35 complex traits (including 12 male fertility-related traits). We detected genomic regions, which experienced DNA methylation alterations in sperm and were associated with aging and extreme fertility phenotypes (e.g., sire-conception rate or SCR). In dynamic hypomethylated regions (HMRs) and partially methylated domains (PMDs), we found genes (e.g., HOX gene clusters and microRNAs) that were involved in the embryonic development. We demonstrated that genomic regions, which gained rather than lost methylations during aging, and in animals with low SCR were significantly and selectively enriched for GWAS signals of male fertility traits. Our study discovered 16 genes as the potential candidate markers for male fertility, including SAMD5 and PDE5A. Collectively, this initial effort supported a hypothesis that sperm DNA methylation may contribute to male fertility in cattle and revealed the usefulness of functional annotations in enhancing biological interpretation and genomic prediction for complex traits and diseases.

Epigenetics and precision medicine in cardiovascular patients: from basic concepts to the clinical arena

Cardiovascular diseases (CVDs) remain the leading cause of mortality worldwide and also inflict major burdens on morbidity, quality of life, and societal costs. Considering that CVD preventive medications improve vascular outcomes in less than half of patients (often relative risk reductions range from 12% to 20% compared with placebo), precision medicine offers an attractive approach to refine the targeting of CVD medications to responsive individuals in a population and thus allocate resources more wisely and effectively. New tools furnished by advances in basic science and translational medicine could help achieve this goal. This approach could reach beyond the practitioners 'eyeball' assessment or venerable markers derived from the physical examination and standard laboratory evaluation. Advances in genetics have identified novel pathways and targets that operate in numerous diseases, paving the way for 'precision medicine'. Yet the inherited genome determines only part of an individual's risk profile. Indeed, standard genomic approaches do not take into account the world of regulation of gene expression by modifications of the 'epi'genome. Epigenetic modifications defined as 'heritable changes to the genome that do not involve changes in DNA sequence' have emerged as a new layer of biological regulation in CVD and could advance individualized risk assessment as well as devising and deploying tailored therapies. This review, therefore, aims to acquaint the cardiovascular community with the rapidly advancing and evolving field of epigenetics and its implications in cardiovascular precision medicine.

METHimpute: imputation-guided construction of complete methylomes from WGBS data

BACKGROUND

Whole-genome bisulfite sequencing (WGBS) has become the standard method for interrogating plant methylomes at base resolution. However, deep WGBS measurements remain cost prohibitive for large, complex genomes and for population-level studies. As a result, most published plant methylomes are sequenced far below saturation, with a large proportion of cytosines having either missing data or insufficient coverage.

RESULTS

Here we present METHimpute, a Hidden Markov Model (HMM) based imputation algorithm for the analysis of WGBS data. Unlike existing methods, METHimpute enables the construction of complete methylomes by inferring the methylation status and level of all cytosines in the genome regardless of coverage. Application of METHimpute to maize, rice and Arabidopsis shows that the algorithm infers cytosine-resolution methylomes with high accuracy from data as low as 6X, compared to data with 60X, thus making it a cost-effective solution for large-scale studies.

CONCLUSIONS

METHimpute provides methylation status calls and levels for all cytosines in the genome regardless of coverage, thus yielding complete methylomes even with low-coverage WGBS datasets. The method has been extensively tested in plants, but should also be applicable to other species. An implementation is available on Bioconductor.

Epigenetic modification of nucleic acids: from basic studies to medical applications

The epigenetic modification of nucleic acids represents one of the most significant areas of study in the field of nucleic acids because it makes gene regulation more complex and heredity more complicated, thus indicating its profound impact on aspects of heredity, growth, and diseases. The recent characterization of epigenetic modifications of DNA and RNA using chemical labelling strategies has promoted the discovery of these modifications, and the newly developed single-base or single-cell resolution mapping strategies have enabled large-scale epigenetic studies in eukaryotes. Due to these technological breakthroughs, several new epigenetic marks have been discovered that have greatly extended the scope and impact of epigenetic modifications in nucleic acids over the past few years. Because epigenetics is reversible and susceptible to environmental factors, it could potentially be a promising direction for clinical medicine research. In this review, we have comprehensively discussed how these epigenetic marks are involved in disease, including the pathogenesis, prevention, diagnosis and treatment of disease. These findings have revealed that the epigenetic modification of nucleic acids has considerable significance in various areas from methodology to clinical medicine and even in biomedical applications.

Statistical and integrative system-level analysis of DNA methylation data

Epigenetics plays a key role in cellular development and function. Alterations to the epigenome are thought to capture and mediate the effects of genetic and environmental risk factors on complex disease. Currently, DNA methylation is the only epigenetic mark that can be measured reliably and genome-wide in large numbers of samples. This Review discusses some of the key statistical challenges and algorithms associated with drawing inferences from DNA methylation data, including cell-type heterogeneity, feature selection, reverse causation and system-level analyses that require integration with other data types such as gene expression, genotype, transcription factor binding and other epigenetic information.

Enabling personalized medicine in Europe by the European Commission's funding activities

Personalized medicine (PM) is an emerging approach to prevention, diagnosis, treatment and care. It helps to address the challenge of the aging of the population, an increase in chronic disease and increasing healthcare costs. The EU is developing policies to move toward PM. This is underpinned by a sustained and significant investment starting in 2010. So far, a total of €3.2 billion has been invested in PM research across the medical innovation cycle ‘from bench to bedside’. This investment has come from the research framework programs FP7 and Horizon 2020. About a third of the total investment has been made in the context of the Innovative Medicines Initiative, the largest public–private partnership in life sciences globally.

The International Human Epigenome Consortium: A Blueprint for Scientific Collaboration and Discovery

The International Human Epigenome Consortium (IHEC) coordinates the generation of a catalog of high-resolution reference epigenomes of major primary human cell types. The studies now presented (see the Cell Press IHEC web portal at http://www.cell.com/consortium/IHEC) highlight the coordinated achievements of IHEC teams to gather and interpret comprehensive epigenomic datasets to gain insights in the epigenetic control of cell states relevant for human health and disease. PAPERCLIP.

Stress and the Emerging Roles of Chromatin Remodeling in Signal Integration and Stable Transmission of Reversible Phenotypes

The influence of early life experience and degree of parental-infant attachment on emotional development in children and adolescents has been comprehensively studied. Structural and mechanistic insight into the biological foundation and maintenance of mammalian defensive systems (metabolic, immune, nervous and behavioral) is slowly advancing through the emerging field of developmental molecular (epi)genetics. Initial evidence revealed that differential nurture early in life generates stable differences in offspring hypothalamic-pituitary-adrenal (HPA) regulation, in part, through chromatin remodeling and changes in DNA methylation of specific genes expressed in the brain, revealing physical, biochemical and molecular paths for the epidemiological concept of gene-environment interactions. Herein, a primary molecular mechanism underpinning the early developmental programming and lifelong maintenance of defensive (emotional) responses in the offspring is the alteration of chromatin domains of specific genomic regions from a condensed state (heterochromatin) to a transcriptionally accessible state (euchromatin). Conversely, DNA methylation promotes the formation of heterochromatin, which is essential for gene silencing, genomic integrity and chromosome segregation. Therefore, inter-individual differences in chromatin modifications and DNA methylation marks hold great potential for assessing the impact of both early life experience and effectiveness of intervention programs—from guided psychosocial strategies focused on changing behavior to pharmacological treatments that target chromatin remodeling and DNA methylation enzymes to dietary approaches that alter cellular pools of metabolic intermediates and methyl donors to affect nutrient bioavailability and metabolism. In this review article, we discuss the potential molecular mechanism(s) of gene regulation associated with chromatin modeling and programming of endocrine (e.g., HPA and metabolic or cardiovascular) and behavioral (e.g., fearfulness, vigilance) responses to stress, including alterations in DNA methylation and the role of DNA repair machinery. From parental history (e.g., drugs, housing, illness, nutrition, socialization) to maternal-offspring exchanges of nutrition, microbiota, antibodies and stimulation, the nature of nurture provides not only mechanistic insight into how experiences propagate from external to internal variables, but also identifies a composite therapeutic target, chromatin modeling, for gestational/prenatal stress, adolescent anxiety/depression and adult-onset neuropsychiatric disease.

Cell-type deconvolution in epigenome-wide association studies: a review and recommendations

A major challenge faced by epigenome-wide association studies (EWAS) is cell-type heterogeneity. As many EWAS have already demonstrated, adjusting for changes in cell-type composition can be critical when analyzing and interpreting findings from such studies. Because of their importance, a great number of different statistical algorithms, which adjust for cell-type composition, have been proposed. Some of the methods are 'reference based' in that they require a priori defined reference DNA methylation profiles of cell types that are present in the tissue of interest, while other algorithms are 'reference free.' At present, however, it is unclear how best to adjust for cell-type heterogeneity, as this may also largely depend on the type of tissue and phenotype being considered. Here, we provide a critical review of the major existing algorithms for correcting cell-type composition in the context of Illumina Infinium Methylation Beadarrays, with the aim of providing useful recommendations to the EWAS community.

Clinical implications of genome-wide DNA methylation studies in acute myeloid leukemia

Acute myeloid leukemia (AML) is the most common type of acute leukemia in adults. AML is a heterogeneous malignancy characterized by distinct genetic and epigenetic abnormalities. Recent genome-wide DNA methylation studies have highlighted an important role of dysregulated methylation signature in AML from biological and clinical standpoint. In this review, we will outline the recent advances in the methylome study of AML and overview the impacts of DNA methylation on AML diagnosis, treatment, and prognosis.

Tethered Oligonucleotide-Primed Sequencing, TOP-Seq: A High-Resolution Economical Approach for DNA Epigenome Profiling

Modification of CG dinucleotides in DNA is part of epigenetic regulation of gene function in vertebrates and is associated with complex human disease. Bisulfite sequencing permits high-resolution analysis of cytosine modification in mammalian genomes; however, its utility is often limited due to substantial cost. Here, we describe an alternative epigenome profiling approach, named TOP-seq, which is based on covalent tagging of individual unmodified CG sites followed by non-homologous priming of the DNA polymerase action at these sites to directly produce adjoining regions for their sequencing and precise genomic mapping. Pilot TOP-seq analyses of bacterial and human genomes showed a better agreement of TOP-seq with published bisulfite sequencing maps as compared to widely used MBD-seq and MRE-seq and permitted identification of long-range and gene-level differential methylation among human tissues and neuroblastoma cell types. Altogether, we propose an affordable single CG-resolution technique well suited for large-scale epigenome studies.

Information recovery from low coverage whole-genome bisulfite sequencing

The cost of whole-genome bisulfite sequencing (WGBS) remains a bottleneck for many studies and it is therefore imperative to extract as much information as possible from a given dataset. This is particularly important because even at the recommend 30X coverage for reference methylomes, up to 50% of high-resolution features such as differentially methylated positions (DMPs) cannot be called with current methods as determined by saturation analysis. To address this limitation, we have developed a tool that dynamically segments WGBS methylomes into blocks of comethylation (COMETs) from which lost information can be recovered in the form of differentially methylated COMETs (DMCs). Using this tool, we demonstrate recovery of ∼30% of the lost DMP information content as DMCs even at very low (5X) coverage. This constitutes twice the amount that can be recovered using an existing method based on differentially methylated regions (DMRs). In addition, we explored the relationship between COMETs and haplotypes in lymphoblastoid cell lines of African and European origin. Using best fit analysis, we show COMETs to be correlated in a population-specific manner, suggesting that this type of dynamic segmentation may be useful for integrated (epi)genome-wide association studies in the future.

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.