Search

The human placenta and its specialized cytotrophoblasts rapidly develop, have a compressed lifespan, govern pregnancy outcomes, and program the offspring's health. Understanding the molecular underpinnings of these behaviors informs development and disease. Profiling the extraembryonic epigenome and transcriptome during the 2nd and 3rd trimesters revealed H3K9 trimethylation overlapping deeply DNA hypomethylated domains with reduced gene expression and compartment-specific patterns that illuminated their functions. Cytotrophoblast DNA methylation increased, and several key histone modifications decreased across the genome as pregnancy advanced. Cytotrophoblasts from severe preeclampsia had substantially increased H3K27 acetylation globally and at genes that are normally downregulated at term but upregulated in this syndrome. In addition, some cases had an immature pattern of H3K27ac peaks, and others showed evidence of accelerated aging, suggesting subtype-specific alterations in severe preeclampsia. Thus, the cytotrophoblast epigenome dramatically reprograms during pregnancy, placental disease is associated with failures in this process, and H3K27 hyperacetylation is a feature of severe preeclampsia. • The epigenome of human placental cytotrophoblasts rapidly evolves over gestation • At genome scale, DNA methylation increases, and modified histones are lost • H3K9me3 occupancy of deeply hypomethylated DNA silences cytotrophoblast genes • Global increases in H3K27ac occupancy typify placental disease, for example, severe preeclampsia Zhang, Kim et al. show significant alterations in the cytotrophoblast epigenome over pregnancy, including increased DNA methylation and loss of modified histones. Throughout gestation, H3K9 trimethylation overlapping DNA hypomethylated domains reduces gene expression. In placental disease, for example, severe preeclampsia, cytotrophoblasts have global increases in H3K27 acetylation and at misexpressed genes. [ABSTRACT FROM AUTHOR]

Most genetic risk variants linked to ocular diseases are non-protein coding and presumably contribute to disease through dysregulation of gene expression, however, deeper understanding of their mechanisms of action has been impeded by an incomplete annotation of the transcriptional regulatory elements across different retinal cell types. To address this knowledge gap, we carried out single-cell multiomics assays to investigate gene expression, chromatin accessibility, DNA methylome and 3D chromatin architecture in human retina, macula, and retinal pigment epithelium (RPE)/choroid. We identified 420,824 unique candidate regulatory elements and characterized their chromatin states in 23 sub-classes of retinal cells. Comparative analysis of chromatin landscapes between human and mouse retina cells further revealed both evolutionarily conserved and divergent retinal gene-regulatory programs. Leveraging the rapid advancements in deep-learning techniques, we developed sequence-based predictors to interpret non-coding risk variants of retina diseases. Our study establishes retina-wide, single-cell transcriptome, epigenome, and 3D genome atlases, and provides a resource for studying the gene regulatory programs of the human retina and relevant diseases., Competing Interests: Competing interests: B.R. is a co-founder of Epigenome Technologies and has equity in Arima Genomics.

The dynamic three-dimensional (3D) organization of the human genome (the "4D Nucleome") is closely linked to genome function. Here, we integrate a wide variety of genomic data generated by the 4D Nucleome Project to provide a detailed view of human 3D genome organization in widely used embryonic stem cells (H1-hESCs) and immortalized fibroblasts (HFFc6). We provide extensive benchmarking of 3D genome mapping assays and integrate these diverse datasets to annotate spatial genomic features across scales. The data reveal a rich complexity of chromatin domains and their sub-nuclear positions, and over one hundred thousand structural loops and promoter-enhancer interactions. We developed 3D models of population-based and individual cell-to-cell variation in genome structure, establishing connections between chromosome folding, nuclear organization, chromatin looping, gene transcription, and DNA replication. We demonstrate the use of computational methods to predict genome folding from DNA sequence, uncovering potential effects of genetic variants on genome structure and function. Together, this comprehensive analysis contributes insights into human genome organization and enhances our understanding of connections between the regulation of genome function and 3D genome organization in general., Competing Interests: Conflicts of interest Job Dekker is a member of the scientific advisory board of Arima Genomics, San Diego, CA, USA and Omega Therapeutic, Cambridge, MA, USA. Sheng Zhong. is a founder and shareholder of Genemo, Inc., San Diego, CA, USA Bing Ren has equity in Arima Genomics Inc., San Diego, CA, USA and Epigenome Technologies, San Diego, CA, USA Clair Marchal is the director and founder of In silichrom ltd.

Age-related cognitive decline is associated with altered physiology of the hippocampus. While changes in gene expression have been observed in aging brain, the regulatory mechanisms underlying these changes remain underexplored. We generated single-nucleus gene expression, chromatin accessibility, DNA methylation, and 3D genome data from 40 human hippocampal tissues spanning adult lifespan. We observed a striking loss of astrocytes, OPC, and endothelial cells during aging, including astrocytes that play a role in regulating synapses. Microglia undergo a dramatic switch from a homeostatic state to a primed inflammatory state through DNA methylome and 3D genome reprogramming. Aged cells experience erosion of their 3D genome architecture. Our study identifies age-associated changes in cell types/states and gene regulatory features that provide insight into cognitive decline during human aging., Competing Interests: Competing interests: BR is a co-founder and consultant of Arima Genomics and co-founder of Epigenome Technologies.

Inhibiting epigenetic modulators can transcriptionally reactivate transposable elements (TEs). These TE transcripts often generate unique peptides that can serve as immunogenic antigens for immunotherapy. Here, we ask whether TEs activated by epigenetic therapy could appreciably increase the antigen repertoire in glioblastoma, an aggressive brain cancer with low mutation and neoantigen burden. We treated patient-derived primary glioblastoma stem cell lines, an astrocyte cell line and primary fibroblast cell lines with epigenetic drugs, and identified treatment-induced, TE-derived transcripts that are preferentially expressed in cancer cells. We verified that these transcripts could produce human leukocyte antigen class I-presented antigens using liquid chromatography with tandem mass spectrometry pulldown experiments. Importantly, many TEs were also transcribed, even in proliferating nontumor cell lines, after epigenetic therapy, which suggests that targeted strategies like CRISPR-mediated activation could minimize potential side effects of activating unwanted genomic regions. The results highlight both the need for caution and the promise of future translational efforts in harnessing treatment-induced TE-derived antigens for targeted immunotherapy., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Motivation: With single-cell DNA methylation studies yielding vast datasets, existing data formats struggle with the unique challenges of storage and efficient operations, highlighting a need for improved solutions., Results: BAllC (Binary All Cytosines) emerges as a tailored format for methylation data, addressing these challenges. BAllCools, its complementary software toolkit, enhances parsing, indexing, and querying capabilities, promising superior operational speeds and reduced storage needs., Availability and Implementation: https://github.com/jksr/ballcools., (© The Author(s) 2024. Published by Oxford University Press.)

Motivation: Unraveling the transcriptional programs that control how cells divide, differentiate, and respond to their environments requires a precise understanding of transcription factors' (TFs) DNA-binding activities. Calling cards (CC) technology uses transposons to capture transient TF binding events at one instant in time and then read them out at a later time. This methodology can also be used to simultaneously measure TF binding and mRNA expression from single-cell CC and to record and integrate TF binding events across time in any cell type of interest without the need for purification. Despite these advantages, there has been a lack of dedicated bioinformatics tools for the detailed analysis of CC data., Results: We introduce Pycallingcards, a comprehensive Python module specifically designed for the analysis of single-cell and bulk CC data across multiple species. Pycallingcards introduces two innovative peak callers, CCcaller and MACCs, enhancing the accuracy and speed of pinpointing TF binding sites from CC data. Pycallingcards offers a fully integrated environment for data visualization, motif finding, and comparative analysis with RNA-seq and ChIP-seq datasets. To illustrate its practical application, we have reanalyzed previously published mouse cortex and glioblastoma datasets. This analysis revealed novel cell-type-specific binding sites and potential sex-linked TF regulators, furthering our understanding of TF binding and gene expression relationships. Thus, Pycallingcards, with its user-friendly design and seamless interface with the Python data science ecosystem, stands as a critical tool for advancing the analysis of TF functions via CC data., Availability and Implementation: Pycallingcards can be accessed on the GitHub repository: https://github.com/The-Mitra-Lab/pycallingcards., (© The Author(s) 2024. Published by Oxford University Press.)

Recent advances in long-read sequencing technologies have not only dramatically increased sequencing read length but also have improved the accuracy of detecting chemical modifications to the canonical nucleotide bases, thus opening exciting venues to investigate the epigenome. Currently, the ability to visualize modified bases from long-read sequencing data in genome browsers is still limited, preventing users from easily and fully exploring these type of data. To address this limitation, the WashU Epigenome Browser introduces the modbed track type, which provides visualization of modification details in each single read as well as aggregated modifications of individual or multiple molecules across a dynamic range of resolutions. The modbed file can be uploaded for visualization as a local track or viewed with an accessible URL freely on the WashU Epigenome Browser at https://epigenomegateway.wustl.edu/., Competing Interests: The authors declare no competing interests., (© 2023 The Author(s).)

Divergence of cis-regulatory elements drives species-specific traits 1 , but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits., (© 2023. The Author(s).)

Maintaining healthy adipose tissue is crucial for metabolic health, requiring a deeper understanding of adipocyte development and response to high-calorie diets. This study highlights the importance of TET3 during white adipose tissue (WAT) development and expansion. Selective depletion of Tet3 in adipose precursor cells (APCs) reduces adipogenesis, protects against diet-induced adipose expansion, and enhances whole-body metabolism. Transcriptomic analysis of wild-type and Tet3 knockout (KO) APCs unveiled TET3 target genes, including Pparg and several genes linked to the extracellular matrix, pivotal for adipogenesis and remodeling. DNA methylation profiling and functional studies underscore the importance of DNA demethylation in gene regulation. Remarkably, targeted DNA demethylation at the Pparg promoter restored its transcription. In conclusion, TET3 significantly governs adipogenesis and diet-induced adipose expansion by regulating key target genes in APCs., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Environmental factors may alter the fetal genome to cause metabolic diseases. It is unknown whether embryonic immune cell programming impacts the risk of type 2 diabetes in later life. We demonstrate that transplantation of fetal hematopoietic stem cells (HSCs) made vitamin D deficient in utero induce diabetes in vitamin D-sufficient mice. Vitamin D deficiency epigenetically suppresses Jarid2 expression and activates the Mef2/PGC1a pathway in HSCs, which persists in recipient bone marrow, resulting in adipose macrophage infiltration. These macrophages secrete miR106-5p, which promotes adipose insulin resistance by repressing PIK3 catalytic and regulatory subunits and down-regulating AKT signaling. Vitamin D-deficient monocytes from human cord blood have comparable Jarid2/Mef2/PGC1a expression changes and secrete miR-106b-5p, causing adipocyte insulin resistance. These findings suggest that vitamin D deficiency during development has epigenetic consequences impacting the systemic metabolic milieu., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Genome browsers have become an intuitive and critical tool to visualize and analyze genomic features and data. Conventional genome browsers display data/annotations on a single reference genome/assembly; there are also genomic alignment viewer/browsers that help users visualize alignment, mismatch, and rearrangement between syntenic regions. However, there is a growing need for a comparative epigenome browser that can display genomic and epigenomic data sets across different species and enable users to compare them between syntenic regions. Here, we present the WashU Comparative Epigenome Browser. It allows users to load functional genomic data sets/annotations mapped to different genomes and display them over syntenic regions simultaneously. The browser also displays genetic differences between the genomes from single-nucleotide variants (SNVs) to structural variants (SVs) to visualize the association between epigenomic differences and genetic differences. Instead of anchoring all data sets to the reference genome coordinates, it creates independent coordinates of different genome assemblies to faithfully present features and data mapped to different genomes. It uses a simple, intuitive genome-align track to illustrate the syntenic relationship between different species. It extends the widely used WashU Epigenome Browser infrastructure and can be expanded to support multiple species. This new browser function will greatly facilitate comparative genomic/epigenomic research, as well as support the recent growing needs to directly compare and benchmark the T2T CHM13 assembly and other human genome assemblies., (© 2023 Zhuo et al.; Published by Cold Spring Harbor Laboratory Press.)

Cryptic promoters within transposable elements (TEs) can be transcriptionally reactivated in tumors to create new TE-chimeric transcripts, which can produce immunogenic antigens. We performed a comprehensive screen for these TE exaptation events in 33 TCGA tumor types, 30 GTEx adult tissues and 675 cancer cell lines, and identified 1,068 TE-exapted candidates with the potential to generate shared tumor-specific TE-chimeric antigens (TS-TEAs). Whole-lysate and HLA-pulldown mass spectrometry data confirmed that TS-TEAs are presented on the surface of cancer cells. In addition, we highlight tumor-specific membrane proteins transcribed from TE promoters that constitute aberrant epitopes on the extracellular surface of cancer cells. Altogether, we showcase the high pan-cancer prevalence of TS-TEAs and atypical membrane proteins that could potentially be therapeutically exploited and targeted., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Understanding the impact of DNA methylation within different disease contexts often requires accurate assessment of these modifications in a genome-wide fashion. Frequently, patient-derived tissues stored in long-term hospital tissue banks have been preserved using formalin-fixation paraffin-embedding (FFPE). While these samples can comprise valuable resources for studying disease, the fixation process ultimately compromises the DNA's integrity and leads to degradation. Degraded DNA can complicate CpG methylome profiling using traditional techniques, particularly when performing methylation-sensitive restriction enzyme sequencing (MRE-seq), yielding high backgrounds and resulting in lowered library complexity. Here, we describe Capture MRE-seq, a new MRE-seq protocol tailored to preserving unmethylated CpG information when using samples with highly degraded DNA. The results using Capture MRE-seq correlate well (0.92) with traditional MRE-seq calls when profiling non-degraded samples, and can recover unmethylated regions in highly degraded samples when traditional MRE-seq fails, which we validate using bisulfite sequencing-based data (WGBS) as well as methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq)., (© 2023. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society 1,2 . However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals 3,4 . Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome 5 . To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity 6 . Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements., (© 2022. The Author(s).)

The human reference genome is the most widely used resource in human genetics and is due for a major update. Its current structure is a linear composite of merged haplotypes from more than 20 people, with a single individual comprising most of the sequence. It contains biases and errors within a framework that does not represent global human genomic variation. A high-quality reference with global representation of common variants, including single-nucleotide variants, structural variants and functional elements, is needed. The Human Pangenome Reference Consortium aims to create a more sophisticated and complete human reference genome with a graph-based, telomere-to-telomere representation of global genomic diversity. Here we leverage innovations in technology, study design and global partnerships with the goal of constructing the highest-possible quality human pangenome reference. Our goal is to improve data representation and streamline analyses to enable routine assembly of complete diploid genomes. With attention to ethical frameworks, the human pangenome reference will contain a more accurate and diverse representation of global genomic variation, improve gene-disease association studies across populations, expand the scope of genomics research to the most repetitive and polymorphic regions of the genome, and serve as the ultimate genetic resource for future biomedical research and precision medicine., (© 2022. Springer Nature Limited.)

Innate lymphocytes encompass a diverse array of phenotypic identities with specialized functions. DNA methylation and hydroxymethylation are essential for epigenetic fidelity and fate commitment. The landscapes of these modifications are unknown in innate lymphocytes. Here, we characterized the whole-genome distribution of methyl-CpG and 5-hydroxymethylcytosine (5hmC) in mouse innate lymphoid cell 3 (ILC3), ILC2 and natural killer (NK) cells. We identified differentially methylated regions (DMRs) and differentially hydroxymethylated regions (DHMRs) between ILC and NK cell subsets and correlated them with transcriptional signatures. We associated lineage-determining transcription factors (LDTFs) with demethylation and demonstrated unique patterns of DNA methylation/hydroxymethylation in relationship to open chromatin regions (OCRs), histone modifications and TF-binding sites. We further identified an association between hydroxymethylation and NK cell superenhancers (SEs). Using mice lacking the DNA hydroxymethylase TET2, we showed the requirement for TET2 in optimal production of hallmark cytokines by ILC3s and interleukin-17A (IL-17A) by inflammatory ILC2s. These findings provide a powerful resource for studying innate lymphocyte epigenetic regulation and decode the regulatory logic governing their identity., (© 2022. Springer Nature America, Inc.)

Background: Zebrafish pigment cell differentiation provides an attractive model for studying cell fate progression as a neural crest progenitor engenders diverse cell types, including two morphologically distinct pigment cells: black melanophores and reflective iridophores. Nontrivial classical genetic and transcriptomic approaches have revealed essential molecular mechanisms and gene regulatory circuits that drive neural crest-derived cell fate decisions. However, how the epigenetic landscape contributes to pigment cell differentiation, especially in the context of iridophore cell fate, is poorly understood., Results: We chart the global changes in the epigenetic landscape, including DNA methylation and chromatin accessibility, during neural crest differentiation into melanophores and iridophores to identify epigenetic determinants shaping cell type-specific gene expression. Motif enrichment in the epigenetically dynamic regions reveals putative transcription factors that might be responsible for driving pigment cell identity. Through this effort, in the relatively uncharacterized iridophores, we validate alx4a as a necessary and sufficient transcription factor for iridophore differentiation and present evidence on alx4a's potential regulatory role in guanine synthesis pathway., Conclusions: Pigment cell fate is marked by substantial DNA demethylation events coupled with dynamic chromatin accessibility to potentiate gene regulation through cis-regulatory control. Here, we provide a multi-omic resource for neural crest differentiation into melanophores and iridophores. This work led to the discovery and validation of iridophore-specific alx4a transcription factor., (© 2021. The Author(s).)

Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) is a technique widely used to investigate genome-wide chromatin accessibility. The recently published Omni-ATAC-seq protocol substantially improves the signal/noise ratio and reduces the input cell number. High-quality data are critical to ensure accurate analysis. Several tools have been developed for assessing sequencing quality and insertion size distribution for ATAC-seq data; however, key quality control (QC) metrics have not yet been established to accurately determine the quality of ATAC-seq data. Here, we optimized the analysis strategy for ATAC-seq and defined a series of QC metrics for ATAC-seq data, including reads under peak ratio (RUPr), background (BG), promoter enrichment (ProEn), subsampling enrichment (SubEn), and other measurements. We incorporated these QC tests into our recently developed ATAC-seq Integrative Analysis Package (AIAP) to provide a complete ATAC-seq analysis system, including quality assurance, improved peak calling, and downstream differential analysis. We demonstrated a significant improvement of sensitivity (20%-60%) in both peak calling and differential analysis by processing paired-end ATAC-seq datasets using AIAP. AIAP is compiled into Docker/Singularity, and it can be executed by one command line to generate a comprehensive QC report. We used ENCODE ATAC-seq data to benchmark and generate QC recommendations, and developed qATACViewer for the user-friendly interaction with the QC report. The software, source code, and documentation of AIAP are freely available at https://github.com/Zhang-lab/ATAC-seq_QC_analysis., (Copyright © 2021 The Authors. Published by Elsevier B.V. All rights reserved.)

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species 1 . However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

Transposable elements (TEs) are an abundant and rich genetic resource of regulatory sequences 1-3 . Cryptic regulatory elements within TEs can be epigenetically reactivated in cancer to influence oncogenesis in a process termed onco-exaptation 4 . However, the prevalence and impact of TE onco-exaptation events across cancer types are poorly characterized. Here, we analyzed 7,769 tumors and 625 normal datasets from 15 cancer types, identifying 129 TE cryptic promoter-activation events involving 106 oncogenes across 3,864 tumors. Furthermore, we interrogated the AluJb-LIN28B candidate: the genetic deletion of the TE eliminated oncogene expression, while dynamic DNA methylation modulated promoter activity, illustrating the necessity and sufficiency of a TE for oncogene activation. Collectively, our results characterize the global profile of TE onco-exaptation and highlight this prevalent phenomenon as an important mechanism for promiscuous oncogene activation and ultimately tumorigenesis.

Understanding the role of DNA methylation often requires accurate assessment and comparison of these modifications in a genome-wide fashion. Sequencing-based DNA methylation profiling provides an unprecedented opportunity to map and compare complete DNA CpG methylomes. These include whole genome bisulfite sequencing (WGBS), Reduced-Representation Bisulfite-Sequencing (RRBS), and enrichment-based methods such as MeDIP-seq, MBD-seq, and MRE-seq. An investigator needs a method that is flexible with the quantity of input DNA, provides the appropriate balance among genomic CpG coverage, resolution, quantitative accuracy, and cost, and comes with robust bioinformatics software for analyzing the data. In this chapter, we describe four protocols that combine state-of-the-art experimental strategies with state-of-the-art computational algorithms to achieve this goal. We first introduce two experimental methods that are complementary to each other. MeDIP-seq, or methylation-dependent immunoprecipitation followed by sequencing, uses an anti-methylcytidine antibody to enrich for methylated DNA fragments, and uses massively parallel sequencing to reveal identity of enriched DNA. MRE-seq, or methylation-sensitive restriction enzyme digestion followed by sequencing, relies on a collection of restriction enzymes that recognize CpG containing sequence motifs, but only cut when the CpG is unmethylated. Digested DNA fragments enrich for unmethylated CpGs at their ends, and these CpGs are revealed by massively parallel sequencing. The two computational methods both implement advanced statistical algorithms that integrate MeDIP-seq and MRE-seq data. M&M is a statistical framework to detect differentially methylated regions between two samples. methylCRF is a machine learning framework that predicts CpG methylation levels at single CpG resolution, thus raising the resolution and coverage of MeDIP-seq and MRE-seq to a comparable level of WGBS, but only incurring a cost of less than 5% of WGBS. Together these methods form an effective, robust, and affordable platform for the investigation of genome-wide DNA methylation.

Uterine carcinosarcoma (UCS) is a form of endometrial cancer simultaneously exhibiting carcinomatous and sarcomatous elements, but the underlying molecular and epigenetic basis of this disease is poorly understood. We generated complete DNA methylomes for both the carcinomatous and the sarcomatous components of three UCS samples separated by laser capture microdissection and compared DNA methylomes of UCS with those of normal endometrium as well as methylomes derived from endometrioid carcinoma, serous endometrial carcinoma, and endometrial stromal sarcoma. We identified epigenetic lesions specific to carcinosarcoma and specific to its two components. Hallmarks of DNA methylation abnormalities in UCS included global hypomethylation, especially in repetitive elements, and hypermethylation of tumor suppressor gene promoters. Among these, aberrant DNA methylation of MIR200 genes is a key feature of UCS. The carcinoma component of UCS was characterized by hypermethylation of promoters of EMILIN1, NEFM, and CLEC14A, genes that are associated with tumor vascularization. In contrast, DNA methylation changes of PKP3, FAM83F, and TCP11 were more characteristic of the sarcoma components. Our findings highlight the epigenetic signatures that distinguish the two components of UCS, providing a valuable resource for investigation of this disease., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

DNA methylation is an important epigenetic modification involved in many biological processes and diseases. Many studies have mapped DNA methylation changes associated with embryogenesis, cell differentiation, and cancer at a genome-wide scale. Our understanding of genome-wide DNA methylation changes in a developmental or disease-related context has been steadily growing. However, the investigation of which CpGs are variably methylated in different normal cell or tissue types is still limited. Here, we present an in-depth analysis of 54 single-CpG-resolution DNA methylomes of normal human cell types by integrating high-throughput sequencing-based methylation data. We found that the ratio of methylated to unmethylated CpGs is relatively constant regardless of cell type. However, which CpGs made up the unmethylated complement was cell-type specific. We categorized the 26,000,000 human autosomal CpGs based on their methylation levels across multiple cell types to identify variably methylated CpGs and found that 22.6% exhibited variable DNA methylation. These variably methylated CpGs formed 660,000 variably methylated regions (VMRs), encompassing 11% of the genome. By integrating a multitude of genomic data, we found that VMRs enrich for histone modifications indicative of enhancers, suggesting their role as regulatory elements marking cell type specificity. VMRs enriched for transcription factor binding sites in a tissue-dependent manner. Importantly, they enriched for GWAS variants, suggesting that VMRs could potentially be implicated in disease and complex traits. Taken together, our results highlight the link between CpG methylation variation, genetic variation, and disease risk for many human cell types., (Copyright © 2016 Gu et al.)

Thyroid hormone binds to nuclear receptors and regulates gene transcription. Here we report that in mice, at around the first day of life, there is a transient surge in hepatocyte type 2 deiodinase (D2) that activates the prohormone thyroxine to the active hormone triiodothyronine, modifying the expression of ∼165 genes involved in broad aspects of hepatocyte function, including lipid metabolism. Hepatocyte-specific D2 inactivation (ALB-D2KO) is followed by a delay in neonatal expression of key lipid-related genes and a persistent reduction in peroxisome proliferator-activated receptor-γ expression. Notably, the absence of a neonatal D2 peak significantly modifies the baseline and long-term hepatic transcriptional response to a high-fat diet (HFD). Overall, changes in the expression of approximately 400 genes represent the HFD response in control animals toward the synthesis of fatty acids and triglycerides, whereas in ALB-D2KO animals, the response is limited to a very different set of only approximately 200 genes associated with reverse cholesterol transport and lipase activity. A whole genome methylation profile coupled to multiple analytical platforms indicate that 10-20% of these differences can be related to the presence of differentially methylated local regions mapped to sites of active/suppressed chromatin, thus qualifying as epigenetic modifications occurring as a result of neonatal D2 inactivation. The resulting phenotype of the adult ALB-D2KO mouse is dramatic, with greatly reduced susceptibility to diet-induced steatosis, hypertriglyceridemia, and obesity.

Summary: We present methylC track, an efficient mechanism for visualizing single-base resolution DNA methylation data on a genome browser. The methylC track dynamically integrates the level of methylation, the position and context of the methylated cytosine (i.e. CG, CHG and CHH), strand and confidence level (e.g. read coverage depth in the case of whole-genome bisulfite sequencing data). Investigators can access and integrate these information visually at specific locus or at the genome-wide level on the WashU EpiGenome Browser in the context of other rich epigenomic datasets., Availability and Implementation: The methylC track is part of the WashU EpiGenome Browser, which is open source and freely available at http://epigenomegateway.wustl.edu/browser/. The most up-to-date instructions and tools for preparing methylC track are available at http://epigenomegateway.wustl.edu/+/cmtk., Contact: twang@genetics.wustl.edu, Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author 2014. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Aberrant DNA hypomethylation may play an important role in the growth rate of glioblastoma (GBM), but the functional impact on transcription remains poorly understood. We assayed the GBM methylome with MeDIP-seq and MRE-seq, adjusting for copy number differences, in a small set of non-glioma CpG island methylator phenotype (non-G-CIMP) primary tumors. Recurrent hypomethylated loci were enriched within a region of chromosome 5p15 that is specified as a cancer amplicon and also encompasses TERT, encoding telomerase reverse transcriptase, which plays a critical role in tumorigenesis. Overall, 76 gene body promoters were recurrently hypomethylated, including TERT and the oncogenes GLI3 and TP73. Recurring hypomethylation also affected previously unannotated alternative promoters, and luciferase reporter assays for three of four of these promoters confirmed strong promoter activity in GBM cells. Histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq on tissue from the GBMs uncovered peaks that coincide precisely with tumor-specific decrease of DNA methylation at 200 loci, 133 of which are in gene bodies. Detailed investigation of TP73 and TERT gene body hypomethylation demonstrated increased expression of corresponding alternate transcripts, which in TP73 encodes a truncated p73 protein with oncogenic function and in TERT encodes a putative reverse transcriptase-null protein. Our findings suggest that recurring gene body promoter hypomethylation events, along with histone H3K4 trimethylation, alter the transcriptional landscape of GBM through the activation of a limited number of normally silenced promoters within gene bodies, in at least one case leading to expression of an oncogenic protein.

DNA methylation plays key roles in diverse biological processes such as X chromosome inactivation, transposable element repression, genomic imprinting, and tissue-specific gene expression. Sequencing-based DNA methylation profiling provides an unprecedented opportunity to map and compare complete DNA methylomes. This includes one of the most widely applied technologies for measuring DNA methylation: methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq), coupled with a complementary method, methylation-sensitive restriction enzyme sequencing (MRE-seq). A computational approach that integrates data from these two different but complementary assays and predicts methylation differences between samples has been unavailable. Here, we present a novel integrative statistical framework M&M (for integration of MeDIP-seq and MRE-seq) that dynamically scales, normalizes, and combines MeDIP-seq and MRE-seq data to detect differentially methylated regions. Using sample-matched whole-genome bisulfite sequencing (WGBS) as a gold standard, we demonstrate superior accuracy and reproducibility of M&M compared to existing analytical methods for MeDIP-seq data alone. M&M leverages the complementary nature of MeDIP-seq and MRE-seq data to allow rapid comparative analysis between whole methylomes at a fraction of the cost of WGBS. Comprehensive analysis of nineteen human DNA methylomes with M&M reveals distinct DNA methylation patterns among different tissue types, cell types, and individuals, potentially underscoring divergent epigenetic regulation at different scales of phenotypic diversity. We find that differential DNA methylation at enhancer elements, with concurrent changes in histone modifications and transcription factor binding, is common at the cell, tissue, and individual levels, whereas promoter methylation is more prominent in reinforcing fundamental tissue identities.

Recent advancements in sequencing-based DNA methylation profiling methods provide an unprecedented opportunity to map complete DNA methylomes. These include whole-genome bisulfite sequencing (WGBS, MethylC-seq, or BS-seq), reduced-representation bisulfite sequencing (RRBS), and enrichment-based methods such as MeDIP-seq, MBD-seq, and MRE-seq. These methods yield largely comparable results but differ significantly in extent of genomic CpG coverage, resolution, quantitative accuracy, and cost, at least while using current algorithms to interrogate the data. None of these existing methods provides single-CpG resolution, comprehensive genome-wide coverage, and cost feasibility for a typical laboratory. We introduce methylCRF, a novel conditional random fields-based algorithm that integrates methylated DNA immunoprecipitation (MeDIP-seq) and methylation-sensitive restriction enzyme (MRE-seq) sequencing data to predict DNA methylation levels at single-CpG resolution. Our method is a combined computational and experimental strategy to produce DNA methylomes of all 28 million CpGs in the human genome for a fraction (<10%) of the cost of whole-genome bisulfite sequencing methods. methylCRF was benchmarked for accuracy against Infinium arrays, RRBS, WGBS sequencing, and locus-specific bisulfite sequencing performed on the same human embryonic stem cell line. methylCRF transformation of MeDIP-seq/MRE-seq was equivalent to a biological replicate of WGBS in quantification, coverage, and resolution. We used conventional bisulfite conversion, PCR, cloning, and sequencing to validate loci where our predictions do not agree with whole-genome bisulfite data, and in 11 out of 12 cases, methylCRF predictions of methylation level agree better with validated results than does whole-genome bisulfite sequencing. Therefore, methylCRF transformation of MeDIP-seq/MRE-seq data provides an accurate, inexpensive, and widely accessible strategy to create full DNA methylomes.

The density of epilithic algae and concentration of chlorophyll a, Shannon-Wiener index and primary productivity were chosen to study the river continuum characteristics of Xiangxi River. The results showed that in general, the algal density, chlorophyll a concentration and Shannon-Wiener index tended to be increased with increasing stream order, but this trend was not very evident in the first three orders and even reversed. The primary productivity was also higher in the higher than that in the lower stream orders, and reached the maximum in the fourth order, then decreased slightly. The ratio of gross primary productivity to community respiration was the highest in the first stream orders, which means autochthonous productivity was dominant in the river. The results did not absolutely tally with the RCC forecast, therefore, we can conclude that Xiangxi River had been undergone anthropogenic influence.

Cluster analysis and principal component analysis were used to analyze the physical and chemical characteristics of water samples collected from 19 sampling sites along Xiangxi River. The result of cluster analysis showed that Xiangxi River could be divided into three reaches, which belonged to different sub-basins. The principle component analysis on each reach showed that the main information of water quality was different from each other. In the upper reach, the main information of river water quality was hardness and total alkalinity; in the middle reach, the main information was dissolved phosphorus, total phosphorus and Cl-, in the lower reach, the main information was pH, NO2(-)-N, total nitrogen and COD. Furthermore, considering the geographical and environmental backgrounds of the Xiangxi River basin, the reason that caused the different distribution pattern of water quality in Xiangxi River was discussed. Hence, a simple approach to analyze the spatial distribution of water quality was provided.