Search

Pancreatic ductal adenocarcinoma (PDAC) is believed to arise from the accumulation of a series of somatic mutations and is also frequently associated with pancreatic intraepithelial neoplasia (PanIN) lesions. However, there is still debate as to whether the cell type-of-origin of PanINs and PDACs in humans is acinar or ductal. As cell type identity is maintained epigenetically, DNA methylation changes during pancreatic neoplasia can provide a compelling perspective to examine this question. Here, we performed laser-capture microdissection on surgically resected specimens from 18 patients to isolate, with high purity, DNA for whole-genome bisulfite sequencing from four relevant cell types: acini, nonneoplastic ducts, PanIN lesions, and PDAC lesions. Differentially methylated regions (DMR) were identified using two complementary analytical approaches: bsseq, which identifies any DMRs but is particularly useful for large block-like DMRs, and informME, which profiles the potential energy landscape across the genome and is particularly useful for identifying differential methylation entropy. Both global methylation profiles and block DMRs clearly implicated an acinar origin for PanINs. At the gene level, PanIN lesions exhibited an intermediate acinar-ductal phenotype resembling acinar-to-ductal metaplasia. In 97.6% of PanIN-specific DMRs, PanIN lesions had an intermediate methylation level between normal and PDAC, which suggests from an information theory perspective that PanIN lesions are epigenetically primed to progress to PDAC. Thus, epigenomic analysis complements histopathology to define molecular progression toward PDAC. The shared epigenetic lineage between PanIN and PDAC lesions could provide an opportunity for prevention by targeting aberrantly methylated progression-related genes. Significance: Analysis of DNA methylation landscapes provides insights into the cell-of-origin of PanIN lesions, clarifies the role of PanIN lesions as metaplastic precursors to human PDAC, and suggests potential targets for chemoprevention.

Epigenetic information defines tissue identity and is largely inherited in development through DNA methylation. While studied mostly for mean differences, methylation also encodes stochastic change, defined as entropy in information theory. Analyzing allele-specific methylation in 49 human tissue sample datasets, we find that methylation entropy is associated with specific DNA binding motifs, regulatory DNA, and CpG density. Then applying information theory to 42 mouse embryo methylation datasets, we find that the contribution of methylation entropy to time- and tissue-specific patterns of development is comparable to the contribution of methylation mean, and methylation entropy is associated with sequence and chromatin features conserved with human. Moreover, methylation entropy is directly related to gene expression variability in development, suggesting a role for epigenetic entropy in developmental plasticity.

Blood lead (Pb) level (BLL) is a commonly used biomarker to evaluate associations with health effects. However, interventions to reduce the adverse effects of Pb require relating BLL to external exposure. Moreover, risk mitigation actions need to ensure protection of more susceptible individuals with a greater tendency to accumulate Pb. Because little data is available to quantify inter-individual variability in biokinetics of Pb, we investigated the influence of genetics and diet on BLL in the genetically diverse Collaborative Cross (CC) mouse population. Adult female mice from 49 CC strains received either a standard mouse chow or a chow mimicking the American diet while being provided water ad libitum with 1000 ppm Pb for 4 weeks. In both arms of the study, inter-strain variability was observed; however, in American diet-fed animals, the BLL was greater and more variable. Importantly, the degree of variation in BLL among strains on the American diet was greater (2.3) than the default variability estimate (1.6) used in setting the regulatory standards. Genetic analysis identified suggestive diet-associated haplotypes that were associated with variation in BLL, largely contributed by the PWK/PhJ strain. This study quantified the variation in BLL that is due to genetic background, diet, and their interactions, and observed that it may be greater than that assumed for current regulatory standards for Pb in drinking water. Moreover, this work highlights the need of characterizing inter-individual variation in BLL to ensure adequate public health interventions aimed at reducing human health risks from Pb.

Pairwise comparisons of promoter jsDMRs

Pairwise mean-based differentially methylated regions

Pairwise jsDMRs

Motif analysis

Pancreatic ductal adenocarcinoma (PDAC) is believed to arise from the accumulation of a series of somatic mutations and is also frequently associated with pancreatic intraepithelial neoplasia (PanIN) lesions. However, there is still debate as to whether the cell type-of-origin of PanINs and PDACs in humans is acinar or ductal. As cell type identity is maintained epigenetically, DNA methylation changes during pancreatic neoplasia can provide a compelling perspective to examine this question. Here, we performed laser-capture microdissection on surgically resected specimens from 18 patients to isolate, with high purity, DNA for whole-genome bisulfite sequencing from four relevant cell types: acini, nonneoplastic ducts, PanIN lesions, and PDAC lesions. Differentially methylated regions (DMR) were identified using two complementary analytical approaches: bsseq, which identifies any DMRs but is particularly useful for large block-like DMRs, and informME, which profiles the potential energy landscape across the genome and is particularly useful for identifying differential methylation entropy. Both global methylation profiles and block DMRs clearly implicated an acinar origin for PanINs. At the gene level, PanIN lesions exhibited an intermediate acinar-ductal phenotype resembling acinar-to-ductal metaplasia. In 97.6% of PanIN-specific DMRs, PanIN lesions had an intermediate methylation level between normal and PDAC, which suggests from an information theory perspective that PanIN lesions are epigenetically primed to progress to PDAC. Thus, epigenomic analysis complements histopathology to define molecular progression toward PDAC. The shared epigenetic lineage between PanIN and PDAC lesions could provide an opportunity for prevention by targeting aberrantly methylated progression-related genes.Significance:Analysis of DNA methylation landscapes provides insights into the cell-of-origin of PanIN lesions, clarifies the role of PanIN lesions as metaplastic precursors to human PDAC, and suggests potential targets for chemoprevention.

Supplementary figures and captions for supplementary tables

Variant and mutation calling

Progression DMRs

Diet-related metabolic syndrome is the largest contributor to adverse health in the United States. However, the study of gene-environment interactions and their epigenomic and transcriptomic integration is complicated by the lack of environmental and genetic control in humans that is possible in mouse models. Here we exposed three mouse strains, C57BL/6J (BL6), A/J, and NOD/ShiLtJ (NOD), to a high-fat high-carbohydrate diet, leading to varying degrees of metabolic syndrome. We then performed transcriptomic and genomic DNA methylation analyses and found overlapping but also highly divergent changes in gene expression and methylation upstream of the discordant metabolic phenotypes. Strain-specific pathway analysis of dietary effects reveals a dysregulation of cholesterol biosynthesis common to all three strains but distinct regulatory networks driving this dysregulation. This suggests a strategy for strain-specific targeted pharmacologic intervention of these upstream regulators informed by transcriptional regulation. As a pilot study, we administered the drug GW4064 to target one of these genotype-dependent networks, the Farnesoid X receptor pathway, and found that GW4064 exerts genotype-specific protection against dietary effects in BL6, as predicted by our transcriptomic analysis, as well as increased inflammatory-related gene expression changes in NOD. This pilot study demonstrates the potential efficacy of precision therapeutics for genotype-informed dietary metabolic intervention, and a mouse platform for guiding this approach.

Supplementary Figures 1-4, Methods and Materials from DNA Methyltransferase 1 and 3B Activate BAG-1 Expression via Recruitment of CTCFL/BORIS and Modulation of Promoter Histone Methylation

Supplementary Figures 1-6, Methods and Materials from CTCFL/BORIS Is a Methylation-Independent DNA-Binding Protein That Preferentially Binds to the Paternal H19 Differentially Methylated Region

In a previous genomic analysis, using somatic methyltransferase (DNMT) knockout cells, we showed that hypomethylation decreased the expression of as many genes as were observed to increase, suggesting a previously unknown mechanism for epigenetic regulation. To address this idea, the expression of the BAG family genes was used as a model. These genes were used because their expression was decreased in DNMT1−/−, DNMT3B−/−, and double knockout cells and increased in DNMT1-overexpressing and DNMT3B-overexpressing cells. Chromatin immunoprecipitation analysis of the BAG-1 promoter in DNMT1-overexpressing or DNMT3B-overexpressing cells showed a permissive dimethyl-H3-K4/dimethyl-H3-K9 chromatin status associated with DNA-binding of CTCFL/BORIS, as well as increased BAG-1 expression. In contrast, a nonpermissive dimethyl-H3-K4/dimethyl-H3-K9 chromatin status was associated with CTCF DNA-binding and decreased BAG-1 expression in the single and double DNMT knockout cells. BORIS short hairpin RNA knockdown decreased both promoter DNA-binding, as well as BAG-1 expression, and changed the dimethyl-H3-K4/dimethyl-H3-K9 ratio to that characteristic of a nonpermissive chromatin state. These results suggest that DNMT1 and DNMT3B regulate BAG-1 expression via insulator protein DNA-binding and chromatin dynamics by regulating histone dimethylation. [Cancer Res 2008;68(8):2726–35]

The concept of an epigenetic landscape describing potential cellular fates arising from pluripotent cells, first advanced by Conrad Waddington, has evolved in light of experiments showing nondeterministic outcomes of regulatory processes and mathematical methods for quantifying stochasticity. In this Review, we discuss modern approaches to epigenetic and gene regulation landscapes and the associated ideas of entropy and attractor states, illustrating how their definitions are both more precise and relevant to understanding cancer etiology and the plasticity of cancerous states. We address the interplay between different types of regulatory landscapes and how their changes underlie cancer progression. We also consider the roles of cellular aging and intrinsic and extrinsic stimuli in modulating cellular states and how landscape alterations can be quantitatively mapped onto phenotypic outcomes and thereby used in therapy development.

High-throughput third-generation nanopore sequencing devices have enormous potential for simultaneously observing epigenetic modifications in human cells over large regions of the genome. However, signals generated by these devices are subject to considerable noise that can lead to unsatisfactory detection performance and hamper downstream analysis. Here we develop a statistical method, CpelNano, for the quantification and analysis of 5mC methylation landscapes using nanopore data. CpelNano takes into account nanopore noise by means of a hidden Markov model (HMM) in which the true but unknown (“hidden”) methylation state is modeled through an Ising probability distribution that is consistent with methylation means and pairwise correlations, whereas nanopore current signals constitute the observed state. It then estimates the associated methylation potential energy function by employing the expectation-maximization (EM) algorithm and performs differential methylation analysis via permutation-based hypothesis testing. Using simulations and analysis of published data obtained from three human cell lines (GM12878, MCF-10A, and MDA-MB-231), we show that CpelNano can faithfully estimate DNA methylation potential energy landscapes, substantially improving current methods and leading to a powerful tool for the modeling and analysis of epigenetic landscapes using nanopore sequencing data.

In the search for alternatives to 6-aminonicotinamide (6AN), a series of 6-aminonicotinic acid esters were designed and synthesized as precursors of 6-amino-NADP

In cancer, linking epigenetic alterations to drivers of transformation has been difficult, in part because DNA methylation analyses must capture epigenetic variability, which is central to tumour heterogeneity and tumour plasticity. Here, by conducting a comprehensive analysis, based on information theory, of differences in methylation stochasticity in samples from patients with paediatric acute lymphoblastic leukaemia (ALL), we show that ALL epigenomes are stochastic and marked by increased methylation entropy at specific regulatory regions and genes. By integrating DNA methylation and single-cell gene-expression data, we arrived at a relationship between methylation entropy and gene-expression variability, and found that epigenetic changes in ALL converge on a shared set of genes that overlap with genetic drivers involved in chromosomal translocations across the disease spectrum. Our findings suggest that an epigenetically driven gene-regulation network, with UHRF1 (ubiquitin-like with PHD and RING finger domains 1) as a central node, links genetic drivers and epigenetic mediators in ALL. An information-theoretic analysis of DNA methylation in samples from patients with paediatric acute lymphoblastic leukaemia reveals that a regulatory set of driver genes harbour the greatest differences in methylation stochasticity.

In heterozygous genomes, allele-specific measurements can reveal biologically significant differences in DNA methylation between homologous alleles associated with local changes in genetic sequence. Current approaches for detecting such events from whole-genome bisulfite sequencing (WGBS) data perform statistically independent marginal analysis at individual cytosine-phosphate-guanine (CpG) sites, thus ignoring correlations in the methylation state, or carry-out a joint statistical analysis of methylation patterns at four CpG sites producing unreliable statistical evidence. Here, we employ the one-dimensional Ising model of statistical physics and develop a method for detecting allele-specific methylation (ASM) events within segments of DNA containing clusters of linked single-nucleotide polymorphisms (SNPs), called haplotypes. Comparisons with existing approaches using simulated and real WGBS data show that our method provides an improved fit to data, especially when considering large haplotypes. Importantly, the method employs robust hypothesis testing for detecting statistically significant imbalances in mean methylation level and methylation entropy, as well as for identifying haplotypes for which the genetic variant carries significant information about the methylation state. As such, our ASM analysis approach can potentially lead to biological discoveries with important implications for the genetics of complex human diseases., Allele-specific measurements can reveal differences in DNA methylation between homologous alleles associated with changes in genetic sequence. Here, the authors develop a method for detecting allele specific methylation events within haplotypes of linked SNPs, compare it with existing methods, and show it identifies haplotypes for which the genetic variant carries significant information about the methylation state of the allele of origin.

Translocations of the KMT2A (MLL) gene define a biologically distinct and clinically aggressive subtype of acute myeloid leukaemia (AML), marked by a characteristic gene expression profile and few cooperating mutations. Although dysregulation of the epigenetic landscape in this leukaemia is particularly interesting given the low mutation frequency, its comprehensive analysis using whole genome bisulphite sequencing (WGBS) has not been previously performed. Here we investigated epigenetic dysregulation in nine MLL-rearranged (MLL-r) AML samples by comparing them to six normal myeloid controls, using a computational method that encapsulates mean DNA methylation measurements along with analyses of methylation stochasticity. We discovered a dramatically altered epigenetic profile in MLL-r AML, associated with genome-wide hypomethylation and a markedly increased DNA methylation entropy reflecting an increasingly disordered epigenome. Methylation discordance mapped to key genes and regulatory elements that included bivalent promoters and active enhancers. Genes associated with significant changes in methylation stochasticity recapitulated known MLL-r AML expression signatures, suggesting a role for the altered epigenetic landscape in the transcriptional programme initiated by MLL translocations. Accordingly, we established statistically significant associations between discordances in methylation stochasticity and gene expression in MLL-r AML, thus providing a link between the altered epigenetic landscape and the phenotype.

Data de publicació electrònica: 02-11-2022 Cells respond to physical stimuli, such as stiffness1, fluid shear stress2 and hydraulic pressure3,4. Extracellular fluid viscosity is a key physical cue that varies under physiological and pathological conditions, such as cancer5. However, its influence on cancer biology and the mechanism by which cells sense and respond to changes in viscosity are unknown. Here we demonstrate that elevated viscosity counterintuitively increases the motility of various cell types on two-dimensional surfaces and in confinement, and increases cell dissemination from three-dimensional tumour spheroids. Increased mechanical loading imposed by elevated viscosity induces an actin-related protein 2/3 (ARP2/3)-complex-dependent dense actin network, which enhances Na+/H+ exchanger 1 (NHE1) polarization through its actin-binding partner ezrin. NHE1 promotes cell swelling and increased membrane tension, which, in turn, activates transient receptor potential cation vanilloid 4 (TRPV4) and mediates calcium influx, leading to increased RHOA-dependent cell contractility. The coordinated action of actin remodelling/dynamics, NHE1-mediated swelling and RHOA-based contractility facilitates enhanced motility at elevated viscosities. Breast cancer cells pre-exposed to elevated viscosity acquire TRPV4-dependent mechanical memory through transcriptional control of the Hippo pathway, leading to increased migration in zebrafish, extravasation in chick embryos and lung colonization in mice. Cumulatively, extracellular viscosity is a physical cue that regulates both short- and long-term cellular processes with pathophysiological relevance to cancer biology. This work was supported in part by R01 CA257647 (to K.K. and D.M.G.), R01 GM134542 (to S.X.S. and K.K.), NSF 2045715 (to Y.L.), R01 AR071976 (to C.-M.F. and J.T.), R01 AR072644 (to C.-M.F. and J.T.) and R01 CA054358 (to A.P.F.), the Spanish Ministry of Science, Education and Universities through grants RTI2018 099718-B-100 (to M.A.V.) and an institutional “Maria de Maeztu” Programme for Units of Excellence in R&D and FEDER funds (to M.A.V.), and postdoctoral fellowships from the Fonds de recherche du Quebec—Nature et technologies and the Natural Sciences and Engineering Research Council of Canada (to A.K.). The opinions, findings and conclusions, or recommendations expressed are those of the authors and do not necessarily reflect the views of any of the funding agencies.

Epigenetic modifications confer stable transcriptional patterns in the brain, and both normal and abnormal brain function involve specialized brain regions. We examined DNA methylation by whole-genome bisulfite sequencing in neuronal and non-neuronal populations from four brain regions (anterior cingulate gyrus, hippocampus, prefrontal cortex, and nucleus accumbens) as well as chromatin accessibility in the latter two. We find pronounced differences in both CpG and non-CpG methylation (CG-DMRs and CH-DMRs) only in neuronal cells across brain regions. Neuronal CH-DMRs were highly associated with differential gene expression, whereas CG-DMRs were consistent with chromatin accessibility and enriched for regulatory regions. These CG-DMRs comprise ~12 Mb of the genome that is highly enriched for genomic regions associated with heritability of neuropsychiatric traits including addictive behavior, schizophrenia, and neuroticism, thus suggesting a mechanistic link between pathology and differential neuron-specific epigenetic regulation in distinct brain regions. An extensive profile of DNA methylation in neuronal and non-neuronal cells across four brain regions is reported, showing that differential epigenetic marks are enriched for DNA variants associated with neuropsychiatric traits.

Single-cell omics is transforming our understanding of cell biology and disease, yet the systems-level analysis and interpretation of single-cell data faces many challenges. In this perspective, we describe the impact that fundamental concepts from statistical mechanics, notably entropy, stochastic processes and critical phenomena, are having on single-cell data analysis. We further advocate the need for more bottom-up modelling of single-cell data, and to embrace a statistical mechanics analysis paradigm to help attain a deeper understanding of single-cell systems biology.

Fetal sex-specific differences in placental morphology and physiology have been associated with sexually dimorphic health outcomes. However, the molecular mechanisms underlying these sex differences are not well understood. We performed whole genome bisulfite sequencing in 133 placenta samples and discovered a significant difference in DNA methylation (DNAm) at the ZNF300 gene locus between male and female offspring and replicated this result in 6 independent datasets. Additionally, the sex-specific pattern appears to be placenta-specific, is robust to a wide range of gestational ages and adverse health outcomes and is present in sorted placenta villous cytotrophoblast cells. Integration of DNAm, genetic, and placental morphology data from the same individuals revealed ZNF300 methylation is also associated with placenta area, perimeter, and max diameter, genetic variants on chromosomes 5 and X, and may mediate the effects of genetic variation on placental area.

High-throughput third-generation sequencing devices, such as the Oxford Nanopore Technologies (ONT) MinION sequencer, can generate long reads that span thousands of bases. This new technology opens the possibility of considering a wide range of epigenetic modifications and provides the capability of interrogating previously inaccessible regions of the genome, such as highly repetitive regions, as well as of performing comprehensive allele-specific methylation analysis, among other applications. It is well-known, however, that detection of DNA methylation from nanopore data results in a substantially reduced per-read accuracy when comparing to WGBS, due to noise introduced by the sequencer and its underlying chemistry. It is therefore imperative that methods are developed for the reliable modeling and analysis of the DNA methylation landscape using nanopore data. Here we introduce such method that takes into account the presence of noise introduced by the ONT sequencer and, by using simulations, we provide evidence of its potential. The proposed approach establishes a solid foundation for the development of a comprehensive framework for the statistical analysis of DNA methylation, and possibly of other epigenetic marks, using third-generation sequencing.

The Genotype-Tissue Expression (GTEx) project dissects how genetic variation affects gene expression and splicing. Some human genetic variants affect the amount of RNA produced and the splicing of gene transcripts, crucial steps in development and maintaining a healthy individual. However, some of these changes only occur in a small number of tissues within the body. The Genotype-Tissue Expression (GTEx) project has been expanded over time, and, looking at the final data in version 8, Aguet et al. present a deep characterization of genetic associations and gene expression and splicing in 838 individuals over 49 tissues (see the Perspective by Wilson). This large study was able to characterize the details underlying many aspects of gene expression and provides a resource with which to better understand the fundamental molecular mechanisms of how genetic variants affect gene regulation and complex traits in humans. Science, this issue p. 1318; see also p. 1298 The Genotype-Tissue Expression (GTEx) project was established to characterize genetic effects on the transcriptome across human tissues and to link these regulatory mechanisms to trait and disease associations. Here, we present analyses of the version 8 data, examining 15,201 RNA-sequencing samples from 49 tissues of 838 postmortem donors. We comprehensively characterize genetic associations for gene expression and splicing in cis and trans, showing that regulatory associations are found for almost all genes, and describe the underlying molecular mechanisms and their contribution to allelic heterogeneity and pleiotropy of complex traits. Leveraging the large diversity of tissues, we provide insights into the tissue specificity of genetic effects and show that cell type composition is a key factor in understanding gene regulatory mechanisms in human tissues. We thank the donors and their families for their generous gifts of organ donation for transplantation and tissue donations for the GTEx research project; the Genomics Platform at the Broad Institute for data generation; J. Struewing for support and leadership of the GTEx project; M. Khan and C. Stolte for the illustrations in Fig. 1; and R. Do, D. Jordan, and M. Verbanck for providing GWAS pleiotropy scores. Funding: This work was supported by the Common Fund of the Office of the Director, U.S. National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, NIA, NIAID, and NINDS through NIH contracts HHSN261200800001E (Leidos Prime contract with NCI: A.M.S., D.E.T., N.V.R., J.A.M., L.S., M.E.B., L.Q., T.K., D.B., K.R., and A.U.), 10XS170 (NDRI: W.F.L., J.A.T., G.K., A.M., S.S., R.H., G.Wa., M.J., M.Wa., L.E.B., C.J., J.W., B.R., M.Hu., K.M., L.A.S., H.M.G., M.Mo., and L.K.B.), 10XS171 (Roswell Park Cancer Institute: B.A.F., M.T.M., E.K., B.M.G., K.D.R., and J.B.), 10X172 (Science Care Inc.), 12ST1039 (IDOX), 10ST1035 (Van Andel Institute: S.D.J., D.C.R., and D.R.V.), HHSN268201000029C (Broad Institute: F.A., G.G., K.G.A., A.V.S., X.Li., E.T., S.G., A.G., S.A., K.H.H., D.T.N., K.H., S.R.M., and J.L.N.), 5U41HG009494 (F.A., G.G., and K.G.A.), and through NIH grants R01 DA006227-17 (University of Miami Brain Bank: D.C.M. and D.A.D.), Supplement to University of Miami grant DA006227 (D.C.M. and D.A.D.), R01 MH090941 (University of Geneva), R01 MH090951 and R01 MH090937 (University of Chicago), R01 MH090936 (University of North Carolina–Chapel Hill), R01MH101814 (M.M.-A., V.W., S.B.M., R.G., E.T.D., D.G.-M., and A.V.), U01HG007593 (S.B.M.), R01MH101822 (C.D.B.), U01HG007598 (M.O. and B.E.S.), U01MH104393 (A.P.F.), extension H002371 to 5U41HG002371 (W.J.K.), as well as other funding sources: R01MH106842 (T.L., P.M., E.F., and P.J.H.), R01HL142028 (T.L., Si.Ka., and P.J.H.), R01GM122924 (T.L. and S.E.C.), R01MH107666 (H.K.I.), P30DK020595 (H.K.I.), UM1HG008901 (T.L.), R01GM124486 (T.L.), R01HG010067 (Y.Pa.), R01HG002585 (G.Wa. and M.St.), Gordon and Betty Moore Foundation GBMF 4559 (G.Wa. and M.St.), 1K99HG009916-01 (S.E.C.), R01HG006855 (Se.Ka. and R.E.H.), BIO2015-70777-P, Ministerio de Economia y Competitividad and FEDER funds (M.M.-A., V.W., R.G., and D.G.-M.), la Caixa Foundation ID 100010434 under agreement LCF/BQ/SO15/52260001 (D.G.-M.), NIH CTSA grant UL1TR002550-01 (P.M.), Marie-Skłodowska Curie fellowship H2020 Grant 706636 (S.K.-H.), R35HG010718 (E.R.G.), FPU15/03635, Ministerio de Educación, Cultura y Deporte (M.M.-A.),R01MH109905, 1R01HG010480 (A.Ba.), Searle Scholar Program (A.Ba.), R01HG008150 (S.B.M.), 5T32HG000044-22, NHGRI Institutional Training Grant in Genome Science (N.R.G.), EU IMI program (UE7-DIRECT-115317-1) (E.T.D. and A.V.), FNS funded project RNA1 (31003A_149984) (E.T.D. and A.V.), DK110919 (F.H.), F32HG009987 (F.H.), Massachusetts Lions Eye Research Fund Grant (A.R.H.), Wellcome grant WT108749/Z/15/Z (P.F.), and European Molecular Biology Laboratory (P.F. and D.Z.). Peer Reviewed "Article signat per 1 autors/es del BSC membres del THE GTEX CONSORTIUM: Marta Mele Messeguer"

BackgroundDNA methylation dynamics in the brain are associated with normal development and neuropsychiatric disease and differ across functionally distinct brain regions. Previous studies of genome-wide methylation differences among human brain regions focus on limited numbers of individuals and one to two brain regions.ResultsUsing GTEx samples, we generate a resource of DNA methylation in purified neuronal nuclei from 8 brain regions as well as lung and thyroid tissues from 12 to 23 donors. We identify differentially methylated regions between brain regions among neuronal nuclei in both CpG (181,146) and non-CpG (264,868) contexts, few of which were unique to a single pairwise comparison. This significantly expands the knowledge of differential methylation across the brain by 10-fold. In addition, we present the first differential methylation analysis among neuronal nuclei from basal ganglia tissues and identify unique CpG differentially methylated regions, many associated with ion transport. We also identify 81,130 regions of variably CpG methylated regions, i.e., variable methylation among individuals in the same brain region, which are enriched in regulatory regions and in CpG differentially methylated regions. Many variably methylated regions are unique to a specific brain region, with only 202 common across all brain regions, as well as lung and thyroid. Variably methylated regions identified in the amygdala, anterior cingulate cortex, and hippocampus are enriched for heritability of schizophrenia.ConclusionsThese data suggest that epigenetic variation in these particular human brain regions could be associated with the risk for this neuropsychiatric disorder.

Allele expression (AE) analysis robustly measures cis-regulatory effects. Here, we present and demonstrate the utility of a vast AE resource generated from the GTEx v8 release, containing 15,253 samples spanning 54 human tissues for a total of 431 million measurements of AE at the SNP level and 153 million measurements at the haplotype level. In addition, we develop an extension of our tool phASER that allows effect sizes of cis-regulatory variants to be estimated using haplotype-level AE data. This AE resource is the largest to date, and we are able to make haplotype-level data publicly available. We anticipate that the availability of this resource will enable future studies of regulatory variation across human tissues.

Epigenetics in Disease Prevention and Mitigation Epigenetics is the regulation of gene expression through alterations in DNA or associated factors (other than the DNA sequence). These factors control the diverse manifestations of diseases. Insights into epigenetic modification may lead to new therapies for common diseases.

Background DNA methylation is a stable form of epigenetic memory used by cells to control gene expression. Whole genome bisulfite sequencing (WGBS) has emerged as a gold-standard experimental technique for studying DNA methylation by producing high resolution genome-wide methylation profiles. Statistical modeling and analysis is employed to computationally extract and quantify information from these profiles in an effort to identify regions of the genome that demonstrate crucial or aberrant epigenetic behavior. However, the performance of most currently available methods for methylation analysis is hampered by their inability to directly account for statistical dependencies between neighboring methylation sites, thus ignoring significant information available in WGBS reads. Results We present a powerful information-theoretic approach for genome-wide modeling and analysis of WGBS data based on the 1D Ising model of statistical physics. This approach takes into account correlations in methylation by utilizing a joint probability model that encapsulates all information available in WGBS methylation reads and produces accurate results even when applied on single WGBS samples with low coverage. Using the Shannon entropy, our approach provides a rigorous quantification of methylation stochasticity in individual WGBS samples genome-wide. Furthermore, it utilizes the Jensen-Shannon distance to evaluate differences in methylation distributions between a test and a reference sample. Differential performance assessment using simulated and real human lung normal/cancer data demonstrate a clear superiority of our approach over DSS, a recently proposed method for WGBS data analysis. Critically, these results demonstrate that marginal methods become statistically invalid when correlations are present in the data. Conclusions This contribution demonstrates clear benefits and the necessity of modeling joint probability distributions of methylation using the 1D Ising model of statistical physics and of quantifying methylation stochasticity using concepts from information theory. By employing this methodology, substantial improvement of DNA methylation analysis can be achieved by effectively taking into account the massive amount of statistical information available in WGBS data, which is largely ignored by existing methods. Electronic supplementary material The online version of this article (10.1186/s12859-018-2086-5) contains supplementary material, which is available to authorized users.