Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm

Epigenetic Modulation of GPER Expression in Gastric and Colonic Smooth Muscle of Male and Female Non-Obese Diabetic (NOD) Mice: Insights into H3K4me3 and H3K27ac Modifications

Type 1 diabetes (T1D) affects gastrointestinal (GI) motility, favoring gastroparesis, constipation, and fecal incontinence, which are more prevalent in women. The mechanisms are unknown. Given the G-protein-coupled estrogen receptor's (GPER) role in GI motility, we investigated sex-related diabetes-induced epigenetic changes in GPER. We assessed GPER mRNA and protein expression levels using qPCR and Western blot analyses, and quantified the changes in nuclear DNA methyltransferases and histone modifications (H3K4me3, H3Ac, and H3K27Ac) by ELISA kits. Targeted bisulfite and chromatin immunoprecipitation assays were used to evaluate DNA methylation and histone modifications around the GPER promoter by chromatin immunoprecipitation assays in gastric and colonic smooth muscle tissues of male and female control (CTR) and non-obese diabetic (NOD) mice. GPER expression was downregulated in NOD, with sex-dependent variations. In the gastric smooth muscle, not in colonic smooth muscle, downregulation coincided with differences in methylation ratios between regions 1 and 2 of the GPER promoter of NOD. DNA methylation was higher in NOD male colonic smooth muscle than in NOD females. H3K4me3 and H3ac enrichment decreased in NOD gastric smooth muscle. H3K4me3 levels diminished in the colonic smooth muscle of NOD. H3K27ac levels were unaffected, but enrichment decreased in NOD male gastric smooth muscle; however, it increased in the NOD male colonic smooth muscle and decreased in the female NOD colonic smooth muscle. Male NOD colonic smooth muscle exhibited decreased H3K27ac levels, not female, whereas female NOD colonic smooth muscle demonstrated diminished enrichment of H3ac at the GPER promoter, contrary to male NOD. Sex-specific epigenetic mechanisms contribute to T1D-mediated suppression of GPER expression in the GI tract. These insights advance our understanding of T1D complications and suggest promising avenues for targeted therapeutic interventions.

Understanding the role of DNA methylation in colorectal cancer: Mechanisms, detection, and clinical significance

Colorectal cancer (CRC) is one of the deadliest malignancies worldwide, ranking third in incidence and second in mortality. Remarkably, early stage localized CRC has a 5-year survival rate of over 90%; in stark contrast, the corresponding 5-year survival rate for metastatic CRC (mCRC) is only 14%. Compounding this problem is the staggering lack of effective therapeutic strategies. Beyond genetic mutations, which have been identified as critical instigators of CRC initiation and progression, the importance of epigenetic modifications, particularly DNA methylation (DNAm), cannot be underestimated, given that DNAm can be used for diagnosis, treatment monitoring and prognostic evaluation. This review addresses the intricate mechanisms governing aberrant DNAm in CRC and its profound impact on critical oncogenic pathways. In addition, a comprehensive review of the various techniques used to detect DNAm alterations in CRC is provided, along with an exploration of the clinical utility of cancer-specific DNAm alterations.

Associations in cell type-specific hydroxymethylation and transcriptional alterations of pediatric central nervous system tumors

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors, we utilize a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identify a preponderance differential Cytosine-phosphate-Guanine site hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like histone deacetylase 4 and insulin-like growth factor 1 receptor, are associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric central nervous system tumors.

Cellular and clinical impact of protein phosphatase enzyme epigenetic silencing in multiple cancer tissues

BACKGROUND

Protein Phosphatase Enzymes (PPE) and protein kinases simultaneously control phosphorylation mechanisms that tightly regulate intracellular signalling pathways and stimulate cellular responses. In human malignancies, PPE and protein kinases are frequently mutated resulting in uncontrolled kinase activity and PPE suppression, leading to cell proliferation, migration and resistance to anti-cancer therapies. Cancer associated DNA hypermethylation at PPE promoters gives rise to transcriptional silencing (epimutations) and is a hallmark of cancer. Despite recent advances in sequencing technologies, data availability and computational capabilities, only a fraction of PPE have been reported as transcriptionally inactive as a consequence of epimutations.

METHODS

In this study, we examined promoter-associated DNA methylation profiles in Protein Phosphatase Enzymes and their Interacting Proteins (PPEIP) in a cohort of 705 cancer patients in five tissues (Large intestine, Oesophagus, Lung, Pancreas and Stomach) in three cell models (primary tumours, cancer cell lines and 3D embedded cancer cell cultures). As a subset of PPEIP are known tumour suppressor genes, we analysed the impact of PPEIP promoter hypermethylation marks on gene expression, cellular networks and in a clinical setting.

RESULTS

Here, we report epimutations in PPEIP are a frequent occurrence in the cancer genome and manifest independent of transcriptional activity. We observed that different tumours have varying susceptibility to epimutations and identify specific cellular signalling networks that are primarily affected by epimutations. Additionally, RNA-seq analysis showed the negative impact of epimutations on most (not all) Protein Tyrosine Phosphatase transcription. Finally, we detected novel clinical biomarkers that inform on patient mortality and anti-cancer treatment sensitivity.

CONCLUSIONS

We propose that DNA hypermethylation marks at PPEIP frequently contribute to the pathogenesis of malignancies and within the precision medicine space, hold promise as biomarkers to inform on clinical features such as patient survival and therapeutic response.

DNA methylation dynamics during yak adipocyte differentiation

In mammals, epigenetic modifications involving DNA methylation are necessary for the completion of the cell differentiation process. However, the global DNA methylation landscape and its dynamics during yak adipocyte differentiation remain unexplored. Here, we performed whole-genome bisulfite sequencing (WGBS) to asses DNA methylation in yak preadipocytes and adipocytes, combining these results with those of our previous studies on changes in chromatin accessibility and gene expression during yak adipogenesis. The results showed that CG methylation levels were lower in promoter than in exon and intron, and gradually decreasing from the distal regions to transcription start site (TSS). There was a significant negative correlation between CG methylation levels located in promoter and gene expression levels. The 46 genes shared by CG differentially methylated regions (DMRs) and differential chromatin accessibility were significantly enriched in Hedgehog and PI3K-Akt signaling pathways. ATAC-seq peaks with high chromatin accessibility located in both promoter (≤ 2 kb from TSS) and distal (> 2 kb from TSS) regions corresponded to low methylation levels. Taken together, these findings demonstrated that DNA methylation and its interactions with chromatin accessibility regulate gene expression during yak adipocyte differentiation, contributing to the understanding of mechanisms of various epigenetic modifications and their interactions in adipogenesis.

Creation and Validation of the First Infinium DNA Methylation Array for the Human Imprintome

Background

Differentially methylated imprint control regions (ICRs) regulate the monoallelic expression of imprinted genes. Their epigenetic dysregulation by environmental exposures throughout life results in the formation of common chronic diseases. Unfortunately, existing Infinium methylation arrays lack the ability to profile these regions adequately. Whole genome bisulfite sequencing (WGBS) is the unique method able to profile these regions, but it is very expensive and it requires not only a high coverage but it is also computationally intensive to assess those regions.

Findings

To address this deficiency, we developed a custom methylation array containing 22,819 probes. Among them, 9,757 probes map to 1,088 out of the 1,488 candidate ICRs recently described. To assess the performance of the array, we created matched samples processed with the Human Imprintome array and WGBS, which is the current standard method for assessing the methylation of the Human Imprintome. We compared the methylation levels from the shared CpG sites and obtained a mean R 2 = 0.569. We also created matched samples processed with the Human Imprintome array and the Infinium Methylation EPIC v2 array and obtained a mean R 2 = 0.796. Furthermore, replication experiments demonstrated high reliability (ICC: 0.799-0.945).

Conclusions

Our custom array will be useful for replicable and accurate assessment, mechanistic insight, and targeted investigation of ICRs. This tool should accelerate the discovery of ICRs associated with a wide range of diseases and exposures, and advance our understanding of genomic imprinting and its relevance in development and disease formation throughout the life course.

Transcription factors regulating vasculogenesis and angiogenesis

Transcription factors (TFs) play a crucial role in regulating the dynamic and precise patterns of gene expression required for the initial specification of endothelial cells (ECs), and during endothelial growth and differentiation. While sharing many core features, ECs can be highly heterogeneous. Differential gene expression between ECs is essential to pattern the hierarchical vascular network into arteries, veins and capillaries, to drive angiogenic growth of new vessels, and to direct specialization in response to local signals. Unlike many other cell types, ECs have no single master regulator, instead relying on differing combinations of a necessarily limited repertoire of TFs to achieve tight spatial and temporal activation and repression of gene expression. Here, we will discuss the cohort of TFs known to be involved in directing gene expression during different stages of mammalian vasculogenesis and angiogenesis, with a primary focus on development.

Epigenetic changes to gene pathways linked to male fertility in ex situ black-footed ferrets

Environmental variation can influence the reproductive success of species managed under human care and in the wild, yet the mechanisms underlying this phenomenon remain largely mysterious. Molecular mechanisms such as epigenetic modifiers are important in mediating the timing and progression of reproduction in humans and model organisms, but few studies have linked epigenetic variation to reproductive fitness in wildlife. Here, we investigated epigenetic variation in black-footed ferrets (Mustela nigripes), an endangered North American mammal reliant on ex situ management for survival and persistence in the wild. Despite similar levels of genetic diversity in human-managed and wild-born populations, individuals in ex situ facilities exhibit reproductive problems, such as poor sperm quality. Differences across these settings suggest that an environmentally driven decline in reproductive capacity may be occurring in this species. We examined the role of DNA methylation, one well-studied epigenetic modifier, in this emergent condition. We leveraged blood, testes, and semen samples from male black-footed ferrets bred in ex situ facilities and found tissue-type specificity in DNA methylation across the genome, although 1360 Gene Ontology terms associated with male average litter size shared functions across tissues. We then constructed gene networks of differentially methylated genomic sites associated with three different reproductive phenotypes to explore the putative biological impact of variation in DNA methylation. Sperm gene networks associated with average litter size and sperm count were functionally enriched for candidate genes involved in reproduction, development, and its regulation through transcriptional repression. We propose that DNA methylation plays an important role in regulating these reproductive phenotypes, thereby impacting the fertility of male ex situ individuals. Our results provide information into how DNA methylation may function in the alteration of reproductive pathways and phenotypes in artificial environments. These findings provide early insights to conservation hurdles faced in the protection of this rare species.

Mouse methylation profiles for leukocyte cell types, and estimation of leukocyte fractions in inflamed gastrointestinal DNA samples

Precise analysis of tissue DNA and RNA samples is often hampered by contaminating non-target cells whose amounts are highly variable. DNA methylation profiles are specific to cell types, and can be utilized for assessment of the fraction of such contaminating non-target cells. Here, we aimed 1) to identify methylation profiles specific to multiple types of mouse leukocytes, and 2) to estimate the fraction of leukocytes infiltrating inflamed tissues using DNA samples. First, genome-wide DNA methylation analysis was conducted for three myeloid-lineage cells and four lymphoid-lineage cells isolated by fluorescence-activated cell sorting after magnetic-activated cell sorting from leukocytes in the spleen. Clustering analysis using CpG sites within enhancers separated the three myeloid-lineage cells and four lymphoid-lineage cells while that using promoter CpG islands (TSS200CGIs) did not. Among the 266,108 CpG sites analyzed, one CpG site was specifically hypermethylated (β value ≥ 0.7) in B cells, and four, seven, 183, and 34 CpG sites were specifically hypomethylated (β value < 0.2) in CD4+ T cells, CD8+ T cells, B cells, and NK cells, respectively. Importantly, cell type-specific hypomethylated CpG sites were located at genes involved in cell type-specific biological functions. Then, marker CpG sites to estimate the leukocyte fraction in a tissue with leukocyte infiltration were selected, and an estimation algorithm was established. The fractions of infiltrating leukocytes were estimated to be 1.6-12.4% in the stomach (n = 10) with Helicobacter pylori-induced inflammation and 1.5-4.3% in the colon with dextran sulfate sodium-induced colitis (n = 4), and the fractions were highly correlated with those estimated histologically using Cd45-stained tissue sections [R = 0.811 (p = 0.004)]. These results showed that mouse methylation profiles at CpG sites within enhancers reflected leukocyte cell lineages, and the use of marker CpG sites successfully estimated the leukocyte fraction in inflamed gastric and colon tissues.

Non-small Cell Lung Cancer Epigenomes Exhibit Altered DNA Methylation in Smokers and Never-smokers

Epigenetic alterations are widespread in cancer and can complement genetic alterations to influence cancer progression and treatment outcome. To determine the potential contribution of DNAmethylation alterations to tumor phenotype in non-small cell lung cancer (NSCLC) in both smoker and never-smoker patients, we performed genome-wide profiling of DNA methylation in 17 primary NSCLC tumors and 10 matched normal lung samples using the complementary assays, methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylation sensitive restriction enzyme sequencing (MRE-seq). We reported recurrent methylation changes in the promoters of several genes, many previously implicated in cancer, including FAM83A and SEPT9 (hypomethylation), as well as PCDH7, NKX2-1, and SOX17 (hypermethylation). Although many methylation changes between tumors and their paired normal samples were shared across patients, several were specific to a particular smoking status. For example, never-smokers displayed a greater proportion of hypomethylated differentially methylated regions (hypoDMRs) and a greater number of recurrently hypomethylated promoters, including those of ASPSCR1, TOP2A, DPP9, and USP39, all previously linked to cancer. Changes outside of promoters were also widespread and often recurrent, particularly methylation loss over repetitive elements, highly enriched for ERV1 subfamilies. Recurrent hypoDMRs were enriched for several transcription factor binding motifs, often for genes involved in signaling and cell proliferation. For example, 71% of recurrent promoter hypoDMRs contained a motif for NKX2-1. Finally, the majority of DMRs were located within an active chromatin state in tissues profiled using the Roadmap Epigenomics data, suggesting that methylation changes may contribute to altered regulatory programs through the adaptation of cell type-specific expression programs.

Genome-Wide DNA Methylation Profiles in Whole-Blood and Buccal Samples-Cross-Sectional, Longitudinal, and across Platforms

The field of DNA methylation research is rapidly evolving, focusing on disease and phenotype changes over time using methylation measurements from diverse tissue sources and multiple array platforms. Consequently, identifying the extent of longitudinal, inter-tissue, and inter-platform variation in DNA methylation is crucial for future advancement. DNA methylation was measured in 375 individuals, with 197 of those having 2 blood sample measurements ~10 years apart. Whole-blood samples were measured on Illumina Infinium 450K and EPIC methylation arrays, and buccal samples from a subset of 58 participants were measured on EPIC array. The data were analyzed with the aims to examine the correlation between methylation levels in longitudinal blood samples in 197 individuals, examine the correlation between methylation levels in the blood and buccal samples in 58 individuals, and examine the correlation between blood methylation profiles assessed on the EPIC and 450K arrays in 83 individuals. We identified 136,833, 7674, and 96,891 CpGs significantly and strongly correlated (>0.50) longitudinally, across blood and buccal samples as well as array platforms, respectively. A total of 3674 of these CpGs were shared across all three sets. Analysis of these shared CpGs identified previously found associations with aging, ancestry, and 7016 mQTLs as well.

Accurate prediction of pan-cancer types using machine learning with minimal number of DNA methylation sites

DNA methylation analysis has been applied to determine the primary site of cancer; however, robust and accurate prediction of cancer types with a minimum number of sites is still a significant scientific challenge. To build an accurate and robust cancer type prediction tool with a minimum number of DNA methylation sites, we internally benchmarked different DNA methylation site selection and ranking procedures, as well as different classification models. We used The Cancer Genome Atlas dataset (26 cancer types with 8296 samples) to train and test models and used an independent dataset (17 cancer types with 2738 samples) for model validation. A deep neural network model using a combined feature selection procedure (named MethyDeep) can predict 26 cancer types using 30 methylation sites with superior performance compared with the known methods for both primary and metastatic cancers in independent validation datasets. In conclusion, MethyDeep is an accurate and robust cancer type predictor with the minimum number of DNA methylation sites; it could help the cost-effective clarification of cancer of unknown primary patients and the liquid biopsy-based early screening of cancers.

5-Aza-2′-Deoxycytidine (5-Aza-dC, Decitabine) Inhibits Collagen Type I and III Expression in TGF-β1-Treated Equine Endometrial Fibroblasts

Endometrosis negatively affects endometrial function and fertility in mares, due to excessive deposition of type I (COL1) and type III (COL3) collagens. The pro-fibrotic transforming growth factor (TGF-β1) induces myofibroblast differentiation, characterized by α-smooth muscle actin (α-SMA) expression, and collagen synthesis. In humans, fibrosis has been linked to epigenetic mechanisms. To the best of our knowledge, this has not been described in mare endometrium. Therefore, this study aimed to investigate the in vitro epigenetic regulation in TGF-β1-treated mare endometrial fibroblasts and the use of 5-aza-2′-deoxycytidine (5-aza-dC), an epigenetic modifier, as a putative treatment option for endometrial fibrosis. Methods and Results: The in vitro effects of TGF-β1 and of 5-aza-dC on DNA methyltransferases (DNMT1, DNMT3A, and DNMT3B), COL1A1, COL3A1, and α-SMA transcripts were analyzed in endometrial fibroblasts, and COL1 and COL3 secretion in a co-culture medium. TGF-β1 upregulated DNMT3A transcripts and collagen secretion. In TGF-β1-treated endometrial fibroblasts, DNA methylation inhibitor 5-aza-dC decreased collagen transcripts and secretion, but not α-SMA transcripts. Conclusion: These findings suggest a possible role of epigenetic mechanisms during equine endometrial fibrogenesis. The in vitro effect of 5-aza-dC on collagen reduction in TGF-β1-treated fibroblasts highlights this epigenetic involvement. This may pave the way to different therapeutic approaches for endometrosis.

Hydroxymethylation alterations in progenitor-like cell types of pediatric central nervous system tumors are associated with cell type-specific transcriptional changes

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors we utilized a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identified a preponderance differential CpG hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like HDAC4 and IGF1R, were associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric CNS tumors.

Environmental DNA methylation of Lymnaea stagnalis varies with age and is hypermethylated compared to tissue DNA

Environmental DNA (eDNA) approaches contributing to species identifications are quickly becoming the new norm in biomonitoring and ecosystem assessments. Yet, information such as age and health state of the population, which is vital to species biomonitoring, has not been accessible from eDNA. DNA methylation has the potential to provide such information on the state of a population. Here, we measured the methylation of eDNA along with tissue DNA (tDNA) of Lymnaea stagnalis at four life stages. We demonstrate that eDNA methylation varies with age and allows distinguishing among age classes. Moreover, eDNA was globally hypermethylated in comparison to tDNA. This difference was age-specific and connected to a limited number of eDNA sites. This differential methylation pattern suggests that eDNA release with age is partially regulated through DNA methylation. Our findings help to understand mechanisms involved in eDNA release and shows the potential of eDNA methylation analysis to assess age classes. Such age class assessments will encourage future eDNA studies to assess fundamental processes of population dynamics and functioning in ecology, biodiversity conservation and impact assessments.

Stress induced aging in mouse eye

Aging, a universal process that affects all cells in an organism, is a major risk factor for a group of neuropathies called glaucoma, where elevated intraocular pressure is one of the known stresses affecting the tissue. Our understanding of molecular impact of aging on response to stress in retina is very limited; therefore, we developed a new mouse model to approach this question experimentally. Here we show that susceptibility to response to stress increases with age and is primed on chromatin level. We demonstrate that ocular hypertension activates a stress response that is similar to natural aging and involves activation of inflammation and senescence. We show that multiple instances of pressure elevation cause aging of young retina as measured on transcriptional and DNA methylation level and are accompanied by local histone modification changes. Our data show that repeated stress accelerates appearance of aging features in tissues and suggest chromatin modifications as the key molecular components of aging. Lastly, our work further emphasizes the importance of early diagnosis and prevention as well as age-specific management of age-related diseases, including glaucoma.

Association between genome-wide epigenetic and genetic alterations in breast cancer tissue and response to HER2-targeted therapies in HER2-positive breast cancer patients: new findings and a systematic review

Recent evidence suggests that genetic and epigenetic mechanisms might be associated with acquired resistance to cancer therapies. The aim of this study was to assess the association of genome-wide genetic and epigenetic alterations with the response to anti-HER2 agents in HER2-positive breast cancer patients. PubMed was screened for articles published until March 2021 on observational studies investigating the association of genome-wide genetic and epigenetic alterations, measured in breast cancer tissues or blood, with the response to targeted treatment in HER2-positive breast cancer patients. Sixteen studies were included in the review along with ours, in which we compared the genome-wide DNA methylation pattern in breast tumor tissues of patients who acquired resistance to treatment (case group, n = 6) to that of patients who did not develop resistance (control group, n = 6). Among genes identified as differentially methylated between the breast cancer tissue of cases and controls, one of them, PRKACA, was also reported as differentially expressed in two studies included in the review. Although included studies were heterogeneous in terms of methodology and study population, our review suggests that genes of the PI3K pathway may play an important role in developing resistance to anti-HER2 agents in breast cancer patients. Genome-wide genetic and epigenetic alterations measured in breast cancer tissue or blood might be promising markers of resistance to anti-HER2 agents in HER2-positive breast cancer patients. Further studies are needed to confirm these data.

The Mutagenic Consequences of DNA Methylation within and across Generations

DNA methylation is an epigenetic modification with wide-ranging consequences across the life of an organism. This modification can be stable, persisting through development despite changing environmental conditions. However, in other contexts, DNA methylation can also be flexible, underlying organismal phenotypic plasticity. One underappreciated aspect of DNA methylation is that it is a potent mutagen; methylated cytosines mutate at a much faster rate than other genetic motifs. This mutagenic property of DNA methylation has been largely ignored in eco-evolutionary literature, despite its prevalence. Here, we explore how DNA methylation induced by environmental and other factors could promote mutation and lead to evolutionary change at a more rapid rate and in a more directed manner than through stochastic genetic mutations alone. We argue for future research on the evolutionary implications of DNA methylation driven mutations both within the lifetime of organisms, as well as across timescales.

Characterization of Transposon-Derived Accessible Chromatin Regions in Rice (Oryza Sativa)

Growing evidence indicates that transposons or transposable elements (TEs)-derived accessible chromatin regions (ACRs) play essential roles in multiple biological processes by interacting with trans-acting factors. However, the function of TE-derived ACRs in the regulation of gene expression in the rice genome has not been well characterized. In this study, we examined the chromatin dynamics in six types of rice tissues and found that ~8% of ACRs were derived from TEs and exhibited distinct levels of accessibility and conservation as compared to those without TEs. TEs exhibited a TE subtype-dependent impact on ACR formation, which can be mediated by changes in the underlying DNA methylation levels. Moreover, we found that tissue-specific TE-derived ACRs might function in the tissue development through the modulation of nearby gene expression. Interestingly, many genes in domestication sweeps were found to overlap with TE-derived ACRs, suggesting their potential functions in the rice domestication. In addition, we found that the expression divergence of 1070 duplicate gene pairs were associated with TE-derived ACRs and had distinct distributions of TEs and ACRs around the transcription start sites (TSSs), which may experience different selection pressures. Thus, our study provides some insights into the biological implications of TE-derived ACRs in the rice genome. Our results imply that these ACRs are likely involved in the regulation of tissue development, rice domestication and functional divergence of duplicated genes.

First Steps towards the Development of Epigenetic Biomarkers in Female Cheetahs (Acinonyx jubatus)

Free-ranging cheetahs (Acinonyx jubatus) are generally healthy, whereas cheetahs under human care, such as those in zoological gardens, suffer from ill-defined infectious and degenerative pathologies. These differences are only partially explained by husbandry management programs because both groups share low genetic diversity. However, mounting evidence suggests that physiological differences between populations in different environments can be tracked down to differences in epigenetic signatures. Here, we identified differentially methylated regions (DMRs) between free-ranging cheetahs and conspecifics in zoological gardens and prospect putative links to pathways relevant to immunity, energy balance and homeostasis. Comparing epigenomic DNA methylation profiles obtained from peripheral blood mononuclear cells (PBMCs) from eight free-ranging female cheetahs from Namibia and seven female cheetahs living in zoological gardens within Europe, we identified DMRs of which 22 were hypermethylated and 23 hypomethylated. Hypermethylated regions in cheetahs under human care were located in the promoter region of a gene involved in host-pathogen interactions (KLC1) and in an intron of a transcription factor relevant for the development of pancreatic β-cells, liver, and kidney (GLIS3). The most canonical mechanism of DNA methylation in promoter regions is assumed to repress gene transcription. Taken together, this could indicate that hypermethylation at the promoter region of KLC1 is involved in the reduced immunity in cheetahs under human care. This approach can be generalized to characterize DNA methylation profiles in larger cheetah populations under human care with a more granular longitudinal data collection, which, in the future, could be used to monitor the early onset of pathologies, and ultimately translate into the development of biomarkers with prophylactic and/or therapeutic potential.

Methylation quantitative trait loci are largely consistent across disease states in Crohn’s disease

Recently, we identified 1,189 CpG sites whose DNA methylation level in blood associated with Crohn’s disease. Here, we examined associations between DNA methylation and genetic variants to identify methylation quantitative trait loci across disease states in (1) 402 blood samples from 164 newly diagnosed pediatric Crohn’s disease patients taken at 2 time points (diagnosis and follow-up), and 74 non-inflammatory bowel disease controls, (2) 780 blood samples from a non-Crohn’s disease adult population, and (3) 40 ileal biopsies (17 Crohn’s disease cases and 23 non-inflammatory bowel disease controls) from group (1). Genome-wide DNAm profiling and genotyping were performed using the Illumina MethylationEPIC and Illumina Multi-Ethnic arrays. SNP-CpG associations were identified via linear models adjusted for age, sex, disease status, disease subtype, estimated cell proportions, and genotype-based principal components. In total, we observed 535,448 SNP-CpG associations between 287,881 SNPs and 12,843 CpG sites (P &lt; 8.21 × 10−14). Associations were highly consistent across different ages, races, disease states, and tissue types, suggesting that the majority of these methylation quantitative trait loci participate in common gene regulation. However, genes near CpGs associated with inflammatory bowel disease SNPs were enriched for 18 KEGG pathways relevant to inflammatory bowel disease-linked immune function and inflammatory responses. We observed suggestive evidence for a small number of tissue-specific associations and disease-specific associations in ileum, though larger studies will be needed to confirm these results. Our study concludes that the vast majority of blood-derived methylation quantitative trait loci are common across individuals, though a subset may be involved in processes related to Crohn’s disease. Independent cohort studies will be required to validate these findings.

Physiological tissue-specific and age-related reduction of mouse TDP-43 levels is regulated by epigenetic modifications

The cellular level of TDP-43 (also known as TARDBP) is tightly regulated; increases or decreases in TDP-43 have deleterious effects in cells. The predominant mechanism responsible for the regulation of the level of TDP-43 is an autoregulatory negative feedback loop. In this study, we identified an in vivo cause-effect relationship between Tardbp gene promoter methylation and specific histone modification and the TDP-43 level in tissues of mice at two different ages. Furthermore, epigenetic control was observed in mouse and human cultured cell lines. In amyotrophic lateral sclerosis, the formation of TDP-43-containing brain inclusions removes functional protein from the system. This phenomenon is continuous but compensated by newly synthesized protein. The balance between sequestration and new synthesis might become critical with ageing, if accompanied by an epigenetic modification-regulated decrease in newly synthesized TDP-43. Sequestration by aggregates would then decrease the amount of functional TDP-43 to a level lower than those needed by the cell and thereby trigger the onset of symptoms.

Summary: Identification of a cause-effect relationship between epigenetic modifications that occur on the promoter and histones of mouse TARDBP and the level of TDP-43 both in tissues and in cell culture.

Whole-genome profiling of DNA methylation and hydroxymethylation identifies distinct regulatory programs among innate lymphocytes

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 in mouse ILC3, ILC2, and NK cells. We identified differentially methylated and hydroxymethylated DNA regions between ILC-NK subsets and correlated them with transcriptional signatures. We associated lineage-determining transcription factors with demethylation and demonstrated unique patterns of DNA methylation/hydroxymethylation in relationship to open chromatin regions, histone modifications, and transcription factor binding sites. We further discovered a novel association between hydroxymethylation and NK cell super-enhancers. Using mice lacking DNA hydroxymethylase TET2, we showed its requirement for optimal production of hallmark cytokines by ILC3 and IL-17A by inflammatory ILC2. These findings provide a powerful resource for studying innate lymphocyte epigenetic regulation and decode the regulatory logic governing their identity.

Global DNA Methylation in Dental Implant Failure Due to Peri-Implantitis: An Exploratory Clinical Pilot Study

Background: Peri-implantitis (PIT) is highly prevalent in patients with dental implants and is a challenging condition to treat due to the limited outcomes reported for non-surgical and surgical therapies. Therefore, epigenetic therapeutics might be of key importance to treat PIT. However, developing epigenetic therapeutics is based on understanding the relationship between epigenetics and disease. To date, there is still scarce knowledge about the relationship between epigenetic modifications and PIT, which warrants further investigations. Aim: The purpose of this study was to evaluate the level of global DNA methylation associated with implant failure (IF) due to PIT compared to periodontally healthy (PH) patients. Material and Methods: A total of 20 participants were initially enrolled in this pilot, exploratory, single-blinded, cross-sectional clinical human study in two groups: 10 in the PH group and 10 in the IF group. In the participants who have completed the study, gingival tissue and bone samples were harvested from each participant and were used to perform global DNA methylation analysis. The percentage of global DNA methylation (5-mC%) was compared (1) between groups (PH and IF); (2) between the subgroups of gingival tissue and bone separately; (3) in the whole sample, comparing gingival tissue and bone; (4) within groups, comparing gingival tissue and bone. Demographic, periodontal, and peri-implant measurements as well as periodontal staging, were also recorded. All statistical comparisons were made at the 0.05 significance level. Results: Out of the initially enrolled 20 patients, only 19 completed the study and, thus, were included in the final analysis; 10 patients in the PH group and 9 patients in the IF group, contributing to a total of 38 samples. One patient from the IF group was excluded from the study due to systemic disease. The mean implant survival time was 10.8 years (2.17–15.25 years). Intergroup comparison, stratified by group, indicated a similar 5-mC% between the PH and IF groups in both gingival tissue and bone (p = 0.599), only in bone (p = 0.414), and only in gingival tissue (p = 0.744). Intragroup comparison, stratified by the type of sample, indicated a significantly higher 5-mC% in gingival tissue samples compared to bone in both the PH and IF groups (p = 0.001), in the PH group (p = 0.019), and in the IF group (p = 0.009). Conclusions: Within the limitations of this study, higher global DNA methylation levels were found in gingival tissue samples compared to bone, regardless of the study groups. However, similar global DNA methylation levels were observed overall between the IF and PH groups. Yet, differences in the global DNA methylation levels between gingival tissues and bone, regardless of the study group, could reflect a different epigenetic response between various tissues within the same microenvironment. Further studies are necessary to elucidate the present findings and to evaluate the role of epigenetic modifications in IF due to PIT.

The role of DNA methylation in genome-wide gene regulation during development

Although it is well known that DNA methylation serves to repress gene expression, precisely how it functions during the process of development remains unclear. Here, we propose that the overall pattern of DNA methylation established in the early embryo serves as a sophisticated mechanism for maintaining a genome-wide network of gene regulatory elements in an inaccessible chromatin structure throughout the body. As development progresses, programmed demethylation in each cell type then provides the specificity for maintaining select elements in an open structure. This allows these regulatory elements to interact with a large range of transcription factors and thereby regulate the gene expression profiles that define cell identity.

An Epigenetic Aging Clock for Cattle Using Portable Sequencing Technology

Extensively grazed cattle are often mustered only once a year. Therefore, birthdates are typically unknown or inaccurate. Birthdates would be useful for deriving important traits (growth rate; calving interval), breed registrations, and making management decisions. Epigenetic clocks use methylation of DNA to predict an individual's age. An epigenetic clock for cattle could provide a solution to the challenges of industry birthdate recording. Here we derived the first epigenetic clock for tropically adapted cattle using portable sequencing devices from tail hair, a tissue which is widely used in industry for genotyping. Cattle (n = 66) with ages ranging from 0.35 to 15.7 years were sequenced using Oxford Nanopore Technologies MinION and methylation was called at CpG sites across the genome. Sites were then filtered and used to calculate a covariance relationship matrix based on methylation state. Best linear unbiased prediction was used with 10-fold cross validation to predict age. A second methylation relationship matrix was also calculated that contained sites associated with genes used in the dog and human epigenetic clocks. The correlation between predicted age and actual age was 0.71 for all sites and 0.60 for dog and human gene epigenetic clock sites. The mean absolute deviation was 1.4 years for animals aged less than 3 years of age, and 1.5 years for animals aged 3-10 years. This is the first reported epigenetic clock using industry relevant samples in cattle.

Epigenetics of alcohol use disorder-A review of recent advances in DNA methylation profiling

Alcohol use disorder (AUD) is a major contributor to morbidity and mortality worldwide. Although there is a heritable component, the etiology of AUD is complex and can involve environmental exposures like trauma and can be associated with many different patterns of alcohol consumption. Epigenetic modifications, which can mediate the influence of genetic variants and environmental variables on gene expression, have emerged as an important area of AUD research. Over the past decade, the number of studies investigating AUD and DNA methylation, a form of epigenetic modification, has grown rapidly. Yet we are still far from understanding how DNA methylation contributes to or reflects aspects of AUD. In this paper, we reviewed studies of DNA methylation and AUD and discussed how the field has evolved. We found that global DNA and candidate DNA methylation studies did not produce replicable results. To assess whether findings of epigenome-wide association studies (EWAS) were replicated, we aggregated significant findings across studies and identified 184 genes and 15 gene ontological pathways that were differentially methylated in at least two studies and four genes and three gene ontological pathways that were differentially methylated in three studies. These genes and pathways repeatedly found enrichment of immune processes, which is in line with recent developments suggesting that the immune system may be altered in AUD. Finally, we assess the current limitations of studies of DNA methylation and AUD and make recommendations on how to design future studies to resolve outstanding questions.

Integrated multiomics analysis identifies molecular landscape perturbations during hyperammonemia in skeletal muscle and myotubes

Ammonia is a cytotoxic molecule generated during normal cellular functions. Dysregulated ammonia metabolism, which is evident in many chronic diseases such as liver cirrhosis, heart failure, and chronic obstructive pulmonary disease, initiates a hyperammonemic stress response in tissues including skeletal muscle and in myotubes. Perturbations in levels of specific regulatory molecules have been reported, but the global responses to hyperammonemia are unclear. In this study, we used a multiomics approach to vertically integrate unbiased data generated using an assay for transposase-accessible chromatin with high-throughput sequencing, RNA-Seq, and proteomics. We then horizontally integrated these data across different models of hyperammonemia, including myotubes and mouse and human muscle tissues. Changes in chromatin accessibility and/or expression of genes resulted in distinct clusters of temporal molecular changes including transient, persistent, and delayed responses during hyperammonemia in myotubes. Known responses to hyperammonemia, including mitochondrial and oxidative dysfunction, protein homeostasis disruption, and oxidative stress pathway activation, were enriched in our datasets. During hyperammonemia, pathways that impact skeletal muscle structure and function that were consistently enriched were those that contribute to mitochondrial dysfunction, oxidative stress, and senescence. We made several novel observations, including an enrichment in antiapoptotic B-cell leukemia/lymphoma 2 family protein expression, increased calcium flux, and increased protein glycosylation in myotubes and muscle tissue upon hyperammonemia. Critical molecules in these pathways were validated experimentally. Human skeletal muscle from patients with cirrhosis displayed similar responses, establishing translational relevance. These data demonstrate complex molecular interactions during adaptive and maladaptive responses during the cellular stress response to hyperammonemia.

The methylome of the model arbuscular mycorrhizal fungus, Rhizophagus irregularis, shares characteristics with early diverging fungi and Dikarya

Early-diverging fungi (EDF) are distinct from Dikarya and other eukaryotes, exhibiting high N6-methyldeoxyadenine (6mA) contents, rather than 5-methylcytosine (5mC). As plants transitioned to land the EDF sub-phylum, arbuscular mycorrhizal fungi (AMF; Glomeromycotina) evolved a symbiotic lifestyle with 80% of plant species worldwide. Here we show that these fungi exhibit 5mC and 6mA methylation characteristics that jointly set them apart from other fungi. The model AMF, R. irregularis, evolved very high levels of 5mC and greatly reduced levels of 6mA. However, unlike the Dikarya, 6mA in AMF occurs at symmetrical ApT motifs in genes and is associated with their transcription. 6mA is heterogeneously distributed among nuclei in these coenocytic fungi suggesting functional differences among nuclei. While far fewer genes are regulated by 6mA in the AMF genome than in EDF, most strikingly, 6mA methylation has been specifically retained in genes implicated in components of phosphate regulation; the quintessential hallmark defining this globally important symbiosis.

Anurag Chaturvedi et al. use long-read PacBio sequencing and LC-MS to profile 5mC and 6mA DNA methylation in the model arbuscular mycorrhizal fungus, Rhizophagus irregularis. Their results suggest that R. irregularis shows methylation profiles distinct from other early-diverging fungi, and Dikarya and provide further insight into how these fungi may have adapted to form symbiotic relationships with important plant species.

Nuclear Organization during Hepatogenesis in Zebrafish Requires Uhrf1

Acquisition of cellular fate during development is initiated and maintained by well-coordinated patterns of gene expression that are dictated by the epigenetic landscape and genome organization in the nucleus. While the epigenetic marks that mediate developmental gene expression patterns during organogenesis have been well studied, less is known about how epigenetic marks influence nuclear organization during development. This study examines the relationship between nuclear structure, chromatin accessibility, DNA methylation, and gene expression during hepatic outgrowth in zebrafish larvae. We investigate the relationship between these features using mutants that lack DNA methylation. Hepatocyte nuclear morphology was established coincident with hepatocyte differentiation at 80 h post-fertilization (hpf), and nuclear shape and size continued to change until the conclusion of outgrowth and morphogenesis at 120 hpf. Integrating ATAC-Seq analysis with DNA methylation profiling of zebrafish livers at 120 hpf showed that closed and highly methylated chromatin occupies most transposable elements and that open chromatin correlated with gene expression. DNA hypomethylation, due to mutation of genes encoding ubiquitin-like, containing PHD and RING Finger Domains 1 (uhrf1) and DNA methyltransferase (dnmt1), did not block hepatocyte differentiation, but had dramatic effects on nuclear organization. Hepatocytes in uhrf1 mutants have large, deformed nuclei with multiple nucleoli, downregulation of nucleolar genes, and a complete lack of the nuclear lamina. Loss of lamin B2 staining was phenocopied by dnmt1 mutation. Together, these data show that hepatocyte nuclear morphogenesis coincides with organ morphogenesis and outgrowth, and that DNA methylation directs chromatin organization, and, in turn, hepatocyte nuclear shape and size during liver development.

AIAP: A Quality Control and Integrative Analysis Package to Improve ATAC-seq Data Analysis

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.

Nasal DNA methylation differentiates severe from non-severe asthma in African-American children

BACKGROUND

Asthma is highly heterogeneous, and severity evaluation is key to asthma management. DNA methylation (DNAm) contributes to asthma pathogenesis. This study aimed to identify nasal epithelial DNAm differences between severe and nonsevere asthmatic children and evaluate the impact of environmental exposures.

METHODS

Thirty-three nonsevere and 22 severe asthmatic African American children were included in an epigenome-wide association study. Genome-wide nasal epithelial DNAm and gene expression were measured. CpG sites associated with asthma severity and environmental exposures and predictive of severe asthma were identified. DNAm was correlated with gene expression. Enrichment for transcription factor (TF) binding sites or histone modifications surrounding DNAm differences were determined.

RESULTS

We identified 816 differentially methylated CpG positions (DMPs) and 10 differentially methylated regions (DMRs) associated with asthma severity. Three DMPs exhibited discriminatory ability for severe asthma. Intriguingly, six DMPs were simultaneously associated with asthma, allergic asthma, total IgE, environmental IgE, and FeNO in an independent cohort of children. Twenty-seven DMPs were associated with traffic-related air pollution or secondhand smoke. DNAm at 22 DMPs was altered by diesel particles or allergen in human bronchial epithelial cells. DNAm levels at 39 DMPs were correlated with mRNA expression. Proximal to 816 DMPs, three histone marks and several TFs involved in asthma pathogenesis were enriched.

CONCLUSIONS

Significant differences in nasal epithelial DNAm were observed between nonsevere and severe asthma in African American children, a subset of which may be useful to predict disease severity. These CpG sites are subjected to the influences of environmental exposures and may regulate gene expression.

Common DNA methylation dynamics in endometriod adenocarcinoma and glioblastoma suggest universal epigenomic alterations in tumorigenesis

Trends in altered DNA methylation have been defined across human cancers, revealing global loss of methylation (hypomethylation) and focal gain of methylation (hypermethylation) as frequent cancer hallmarks. Although many cancers share these trends, little is known about the specific differences in DNA methylation changes across cancer types, particularly outside of promoters. Here, we present a comprehensive comparison of DNA methylation changes between two distinct cancers, endometrioid adenocarcinoma (EAC) and glioblastoma multiforme (GBM), to elucidate common rules of methylation dysregulation and changes unique to cancers derived from specific cells. Both cancers exhibit significant changes in methylation over regulatory elements. Notably, hypermethylated enhancers within EAC samples contain several transcription factor binding site clusters with enriched disease ontology terms highlighting uterine function, while hypermethylated enhancers in GBM are found to overlap active enhancer marks in adult brain. These findings suggest that loss of original cellular identity may be a shared step in tumorigenesis.

Characterizing Genetic Regulatory Elements in Ovine Tissues

The Ovine Functional Annotation of Animal Genomes (FAANG) project, part of the broader livestock species FAANG initiative, aims to identify and characterize gene regulatory elements in domestic sheep. Regulatory element annotation is essential for identifying genetic variants that affect health and production traits in this important agricultural species, as greater than 90% of variants underlying genetic effects are estimated to lie outside of transcribed regions. Histone modifications that distinguish active or repressed chromatin states, CTCF binding, and DNA methylation were used to characterize regulatory elements in liver, spleen, and cerebellum tissues from four yearling sheep. Chromatin immunoprecipitation with sequencing (ChIP-seq) was performed for H3K4me3, H3K27ac, H3K4me1, H3K27me3, and CTCF. Nine chromatin states including active promoters, active enhancers, poised enhancers, repressed enhancers, and insulators were characterized in each tissue using ChromHMM. Whole-genome bisulfite sequencing (WGBS) was performed to determine the complement of whole-genome DNA methylation with the ChIP-seq data. Hypermethylated and hypomethylated regions were identified across tissues, and these locations were compared with chromatin states to better distinguish and validate regulatory elements in these tissues. Interestingly, chromatin states with the poised enhancer mark H3K4me1 in the spleen and cerebellum and CTCF in the liver displayed the greatest number of hypermethylated sites. Not surprisingly, active enhancers in the liver and spleen, and promoters in the cerebellum, displayed the greatest number of hypomethylated sites. Overall, chromatin states defined by histone marks and CTCF occupied approximately 22% of the genome in all three tissues. Furthermore, the liver and spleen displayed in common the greatest percent of active promoter (65%) and active enhancer (81%) states, and the liver and cerebellum displayed in common the greatest percent of poised enhancer (53%), repressed enhancer (68%), hypermethylated sites (75%), and hypomethylated sites (73%). In addition, both known and de novo CTCF-binding motifs were identified in all three tissues, with the highest number of unique motifs identified in the cerebellum. In summary, this study has identified the regulatory regions of genes in three tissues that play key roles in defining health and economically important traits and has set the precedent for the characterization of regulatory elements in ovine tissues using the Rambouillet reference genome.

Exploring epigenetic biomarkers of universal specificities and commonalities among pan-cancer cohorts in The Cancer Genome Atlas

Aim: To explore the mechanism of cancer by employing a comprehensive analysis of DNA methylation patterns and variations among pan-cancer cohorts. Materials & methods: This research focused on the discovery of universally specific or common biomarkers by mathematical statistics and machine learning methods in The Cancer Genome Atlas. Results: We found 138 differently methylated CpGs (DMCs) with a common methylation trend and eight common differently methylated regions in different cancer cohorts. Additionally, we found 99 DMCs to distinguish 32 different cancer cohorts in random forest analysis because of the specificity mechanism, but each DMC still had high instability. Conclusion: Our results could facilitate the development of biomarkers that are universally specific and common features across pan-cancer cohorts.

A workflow for simultaneous DNA copy number and methylome analysis of inner cell mass and trophectoderm cells from human blastocysts

OBJECTIVE

To establish a workflow for isolating single trophectoderm (TE) and inner cell mass (ICM) cells and to simultaneously evaluate these cells for copy number variation (CNV) as well as methylome development.

DESIGN

Experimental.

SETTING

Academic medical center.

PATIENT(S)

Donated genetically abnormal blastocysts.

INTERVENTION(S)

Single cells were isolated, followed by bisulfite conversion and sequencing to identify CNV and methylome profiles.

MAIN OUTCOME MEASURE(S)

CNV and methylation profiling.

RESULT(S)

Two embryos were dissociated, isolating 46 single cells, with 17 ICM and 12 TE cells selected for further downstream analysis. Chromosome ploidies and embryo sex were concordant with the results from conventional aneuploidy testing. In 3 of the 29 cells, additional aneuploidies were discovered, indicating possible mosaicism undetected by routine preimplantation genetic testing for aneuploidy. CpG methylation frequency was higher in ICM cells compared with TE cells (44.3% vs. 32.4%), respectively, while non-CpG methylation frequency was similar among both cell types. CpG methylation levels accurately distinguished ICM from TE cells epigenetically.

CONCLUSION(S)

We describe an effective workflow for isolating and sequencing single ICM and TE cells from human blastocysts. The use of methylation profiling can help distinguish these two cell populations better then morphologic identification alone. TE cells had significantly lower levels of DNA methylation, which may be explained in part by the fact that these cells have begun the process of differentiation and are transcriptionally more active than ICM. This approach may be used to explore the genetic complexities within human embryos, specifically among the two primary cell types seen at this stage of development.

Temporal changes in DNA methylation and RNA expression in a small song bird: within- and between-tissue comparisons

Background

DNA methylation is likely a key mechanism regulating changes in gene transcription in traits that show temporal fluctuations in response to environmental conditions. To understand the transcriptional role of DNA methylation we need simultaneous within-individual assessment of methylation changes and gene expression changes over time. Within-individual repeated sampling of tissues, which are essential for trait expression is, however, unfeasible (e.g. specific brain regions, liver and ovary for reproductive timing). Here, we explore to what extend between-individual changes in DNA methylation in a tissue accessible for repeated sampling (red blood cells (RBCs)) reflect such patterns in a tissue unavailable for repeated sampling (liver) and how these DNA methylation patterns are associated with gene expression in such inaccessible tissues (hypothalamus, ovary and liver). For this, 18 great tit (Parus major) females were sacrificed at three time points (n = 6 per time point) throughout the pre-laying and egg-laying period and their blood, hypothalamus, ovary and liver were sampled.

Results

We simultaneously assessed DNA methylation changes (via reduced representation bisulfite sequencing) and changes in gene expression (via RNA-seq and qPCR) over time. In general, we found a positive correlation between changes in CpG site methylation in RBCs and liver across timepoints. For CpG sites in close proximity to the transcription start site, an increase in RBC methylation over time was associated with a decrease in the expression of the associated gene in the ovary. In contrast, no such association with gene expression was found for CpG site methylation within the gene body or the 10 kb up- and downstream regions adjacent to the gene body.

Conclusion

Temporal changes in DNA methylation are largely tissue-general, indicating that changes in RBC methylation can reflect changes in DNA methylation in other, often less accessible, tissues such as the liver in our case. However, associations between temporal changes in DNA methylation with changes in gene expression are mostly tissue- and genomic location-dependent. The observation that temporal changes in DNA methylation within RBCs can relate to changes in gene expression in less accessible tissues is important for a better understanding of how environmental conditions shape traits that temporally change in expression in wild populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07329-9.

Machine learning-based integrative analysis of methylome and transcriptome identifies novel prognostic DNA methylation signature in uveal melanoma

Uveal melanoma (UVM) is the most common primary intraocular human malignancy with a high mortality rate. Aberrant DNA methylation has rapidly emerged as a diagnostic and prognostic signature in many cancers. However, such DNA methylation signature available in UVM remains limited. In this study, we performed a genome-wide integrative analysis of methylome and transcriptome and identified 40 methylation-driven prognostic genes (MDPGs) associated with the tumorigenesis and progression of UVM. Then, we proposed a machine-learning-based discovery and validation strategy to identify a DNA methylation-driven signature (10MeSig) composing of 10 MDPGs (AZGP1, BAI1, CCDC74A, FUT3, PLCD1, S100A4, SCN8A, SEMA3B, SLC25A38 and SLC44A3), which stratified 80 patients of the discovery cohort into two risk subtypes with significantly different overall survival (HR = 29, 95% CI: 6.7-126, P < 0.001). The 10MeSig was validated subsequently in an independent cohort with 57 patients and yielded a similar prognostic value (HR = 2.1, 95% CI: 1.2-3.7, P = 0.006). Multivariable Cox regression analysis showed that the 10MeSig is an independent predictive factor for the survival of patients with UVM. With a prospective validation study, this 10MeSig will improve clinical decisions and provide new insights into the pathogenesis of UVM.

Methylation of the glucocorticoid receptor promoter in children: Links with parents as teachers, early life stress, and behavior problems

The present study examined the effect of early life stress (ELS) on the glucocorticoid receptor gene (NR3C1) methylation, the associations between NR3C1 methylation and behavior problems, and the effect of the program Parents as Teachers (PAT) on NR3C1 methylation. Participants included 132 children, 72 assigned to the PAT intervention group and 60 to the PAT control group. Children were aged 3 years, and were living in psychosocially at-risk families. We assessed NR3C1 methylation of the NGFI-A binding regions of exon 1F via sodium bisulfite sequencing from saliva DNA. Results indicated that (a) children living in families receiving PAT had decreased methylation at one single cytosine–guanine dinucleotides (CpG) site; (b) current maternal depressive symptoms and parental disagreement were predictive of increased methylation of mean NGFI-A and three single CpG sites; and (c) increased methylation of mean NGFI-A and one single CpG site was significantly associated with increased internalizing and externalizing symptoms. In addition, mean NGFI-A was a mediator of the association between parental disagreement and a child's affective problems. These results suggest that PAT may contribute to preventing NR3C1 methylation in preschool children living in psychosocially at-risk situations, and confirm previous findings on the associations between ELS, NR3C1 methylation, and behavior problems.

Biological Aging Measures Based on Blood DNA Methylation and Risk of Cancer: A Prospective Study

Background

We previously investigated the association between 5 "first-generation" measures of epigenetic aging and cancer risk in the Melbourne Collaborative Cohort Study. This study assessed cancer risk associations for 3 recently developed methylation-based biomarkers of aging: PhenoAge, GrimAge, and predicted telomere length.

Methods

We estimated rate ratios (RRs) for the association between these 3 age-adjusted measures and risk of colorectal (N = 813), gastric (N = 165), kidney (N = 139), lung (N = 327), mature B-cell (N = 423), prostate (N = 846), and urothelial (N = 404) cancer using conditional logistic regression models. We also assessed associations by time since blood draw and by cancer subtype, and we investigated potential nonlinearity.

Results

We observed relatively strong associations of age-adjusted PhenoAge with risk of colorectal, kidney, lung, mature B-cell, and urothelial cancers (RR per SD was approximately 1.2-1.3). Similar findings were obtained for age-adjusted GrimAge, but the association with lung cancer risk was much larger (RR per SD = 1.82, 95% confidence interval [CI] = 1.44 to 2.30), after adjustment for smoking status, pack-years, starting age, time since quitting, and other cancer risk factors. Most associations appeared linear, larger than for the first-generation measures, and were virtually unchanged after adjustment for a large set of sociodemographic, lifestyle, and anthropometric variables. For cancer overall, the comprehensively adjusted rate ratio per SD was 1.13 (95% CI = 1.07 to 1.19) for PhenoAge and 1.12 (95% CI = 1.05 to 1.20) for GrimAge and appeared larger within 5 years of blood draw (RR = 1.29, 95% CI = 1.15 to 1.44 and 1.19, 95% CI = 1.06 to 1.33, respectively).

Conclusions

The methylation-based measures PhenoAge and GrimAge may provide insights into the relationship between biological aging and cancer and be useful to predict cancer risk, particularly for lung cancer.

Tissue-specific usage of transposable element-derived promoters in mouse development

BACKGROUND

Transposable elements (TEs) are a significant component of eukaryotic genomes and play essential roles in genome evolution. Mounting evidence indicates that TEs are highly transcribed in early embryo development and contribute to distinct biological functions and tissue morphology.

RESULTS

We examine the epigenetic dynamics of mouse TEs during the development of five tissues: intestine, liver, lung, stomach, and kidney. We found that TEs are associated with over 20% of open chromatin regions during development. Close to half of these accessible TEs are only activated in a single tissue and a specific developmental stage. Most accessible TEs are rodent-specific. Across these five tissues, 453 accessible TEs are found to create the transcription start sites of downstream genes in mouse, including 117 protein-coding genes and 144 lincRNA genes, 93.7% of which are mouse-specific. Species-specific TE-derived transcription start sites are found to drive the expression of tissue-specific genes and change their tissue-specific expression patterns during evolution.

CONCLUSION

Our results suggest that TE insertions increase the regulatory potential of the genome, and some TEs have been domesticated to become a crucial component of gene and regulate tissue-specific expression during mouse tissue development.

TCDD-induced multi- and transgenerational changes in the methylome of male zebrafish gonads

The legacy endocrine disrupting chemical and aryl hydrocarbon receptor agonist, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), is produced as a byproduct of industrial processes and causes adverse health effects ranging from skin irritation to cancer. TCDD endpoints are also observed in subsequent, unexposed generations; however, the mechanisms of these multi- and transgenerational effects are unknown. We hypothesized an epigenetic mechanism, specifically DNA methylation for the transgenerational, male-mediated reproductive effects of developmental TCDD exposure. Using whole genome bisulfite sequencing, we evaluated DNA methylation changes in three generations of zebrafish, the first of which was exposed to TCDD during sexual development at 50 ppt for 1 h at both 3- and 7-week post-fertilization. We discovered that TCDD induces multi- and transgenerational methylomic changes in testicular tissue from zebrafish with decreased reproductive capacity, but most significantly in the indirectly exposed F1 generation. In comparing differentially methylated genes to concurrent transcriptomic changes, we identified several genes and pathways through which transgenerational effects of low level TCDD exposure are likely inherited. These include significant differential methylation of genes involved in reproduction, endocrine function, xenobiotic metabolism, and epigenetic processing. Notably, a number of histone modification genes were both differentially methylated and expressed in all generations, and many differentially methylated genes overlapped between multiple generations. Collectively, our results suggest that DNA methylation is a promising mechanism to explain male-mediated transgenerational reproductive effects of TCDD exposure in zebrafish, and these effects are likely inherited through integration of multiple epigenetic pathways.

Characterizing the Causal Pathway for Genetic Variants Associated with Neurological Phenotypes Using Human Brain-Derived Proteome Data

Leveraging high-dimensional molecular datasets can help us develop mechanistic insight into associations between genetic variants and complex traits. In this study, we integrated human proteome data derived from brain tissue to evaluate whether targeted proteins putatively mediate the effects of genetic variants on seven neurological phenotypes (Alzheimer disease, amyotrophic lateral sclerosis, depression, insomnia, intelligence, neuroticism, and schizophrenia). Applying the principles of Mendelian randomization (MR) systematically across the genome highlighted 43 effects between genetically predicted proteins derived from the dorsolateral prefrontal cortex and these outcomes. Furthermore, genetic colocalization provided evidence that the same causal variant at 12 of these loci was responsible for variation in both protein and neurological phenotype. This included genes such as DCC, which encodes the netrin-1 receptor and has an important role in the development of the nervous system (p = 4.29 × 10-11 with neuroticism), as well as SARM1, which has been previously implicated in axonal degeneration (p = 1.76 × 10-08 with amyotrophic lateral sclerosis). We additionally conducted a phenome-wide MR study for each of these 12 genes to assess potential pleiotropic effects on 700 complex traits and diseases. Our findings suggest that genes such as SNX32, which was initially associated with increased risk of Alzheimer disease, may potentially influence other complex traits in the opposite direction. In contrast, genes such as CTSH (which was also associated with Alzheimer disease) and SARM1 may make worthwhile therapeutic targets because they did not have genetically predicted effects on any of the other phenotypes after correcting for multiple testing.

Glucose-induced, duration-dependent genome-wide DNA methylation changes in human endothelial cells

DNA methylation, a critical epigenetic mechanism, plays an important role in governing gene expressions during biological processes such as aging, which is well known to be accelerated in hyperglycemia (diabetes). In the present study, we investigated the effects of glucose on whole genome DNA methylation in small [human retinal microvascular endothelial cells (HRECs)] and large [human umbilical vein endothelial cells (HUVECs)] vessel endothelial cell (EC) lines exposed to basal or high glucose-containing media for variable lengths of time. Using the Infinium EPIC array, we obtained 773,133 CpG sites (probes) for analysis. Unsupervised clustering of the top 5% probes identified four distinct clusters within EC groups, with significant methylation differences attributed to EC types and the duration of cell culture rather than glucose stimuli alone. When comparing the ECs incubated for 2 days versus 7 days, hierarchical clustering analyses [methylation change >10% and false discovery rate (FDR) <0.05] identified 17,354 and 128 differentially methylated CpGs for HUVECs and HRECs, respectively. Predominant DNA hypermethylation was associated with the length of culture and was enriched for gene enhancer elements and regions surrounding CpG shores and shelves. We identified 88 differentially methylated regions (DMRs) for HUVECs and 8 DMRs for HRECs (all FDR <0.05). Pathway enrichment analyses of DMRs highlighted involvement of regulators of embryonic development (i.e., HOX genes) and cellular differentiation [transforming growth factor-β (TGF-β) family members]. Collectively, our findings suggest that DNA methylation is a complex process that involves tightly coordinated, cell-specific mechanisms. Such changes in methylation overlap genes critical for cellular differentiation and embryonic development.

Intraspecific and interspecific investigations of skeletal DNA methylation and femur morphology in primates

OBJECTIVES

Epigenetic mechanisms influence the development and maintenance of complex phenotypes and may also contribute to the evolution of species-specific phenotypes. With respect to skeletal traits, little is known about the gene regulation underlying these hard tissues or how tissue-specific patterns are associated with bone morphology or vary among species. To begin exploring these topics, this study evaluates one epigenetic mechanism, DNA methylation, in skeletal tissues from five nonhuman primate species which display anatomical and locomotor differences representative of their phylogenetic groups.

MATERIALS AND METHODS

First, we test whether intraspecific variation in skeletal DNA methylation is associated with intraspecific variation in femur morphology. Second, we identify interspecific differences in DNA methylation and assess whether these lineage-specific patterns may have contributed to species-specific morphologies. Specifically, we use the Illumina Infinium MethylationEPIC BeadChip to identify DNA methylation patterns in femur trabecular bone from baboons (n = 28), macaques (n = 10), vervets (n = 10), chimpanzees (n = 4), and marmosets (n = 6).

RESULTS

Significant differentially methylated positions (DMPs) were associated with a subset of morphological variants, but these likely have small biological effects and may be confounded by other variables associated with morphological variation. Conversely, several species-specific DMPs were identified, and these are found in genes enriched for functions associated with complex skeletal traits.

DISCUSSION

Overall, these findings reveal that while intraspecific epigenetic variation is not readily associated with skeletal morphology differences, some interspecific epigenetic differences in skeletal tissues exist and may contribute to evolutionarily distinct phenotypes. This work forms a foundation for future explorations of gene regulation and skeletal trait evolution in primates.

Regenerating zebrafish fin epigenome is characterized by stable lineage-specific DNA methylation and dynamic chromatin accessibility

BACKGROUND

Zebrafish can faithfully regenerate injured fins through the formation of a blastema, a mass of proliferative cells that can grow and develop into the lost body part. After amputation, various cell types contribute to blastema formation, where each cell type retains fate restriction and exclusively contributes to regeneration of its own lineage. Epigenetic changes that are associated with lineage restriction during regeneration remain underexplored.

RESULTS

We produce epigenome maps, including DNA methylation and chromatin accessibility, as well as transcriptomes, of osteoblasts and other cells in uninjured and regenerating fins. This effort reveals regeneration as a process of highly dynamic and orchestrated transcriptomic and chromatin accessibility changes, coupled with stably maintained lineage-specific DNA methylation. The epigenetic signatures also reveal many novel regeneration-specific enhancers, which are experimentally validated. Regulatory networks important for regeneration are constructed through integrative analysis of the epigenome map, and a knockout of a predicted upstream regulator disrupts normal regeneration, validating our prediction.

CONCLUSION

Our study shows that lineage-specific DNA methylation signatures are stably maintained during regeneration, and regeneration enhancers are preset as hypomethylated before injury. In contrast, chromatin accessibility is dynamically changed during regeneration. Many enhancers driving regeneration gene expression as well as upstream regulators of regeneration are identified and validated through integrative epigenome analysis.

Population-level inferences from environmental DNA-Current status and future perspectives

Environmental DNA (eDNA) extracted from water samples has recently shown potential as a valuable source of population genetic information for aquatic macroorganisms. This approach offers several potential advantages compared with conventional tissue-based methods, including the fact that eDNA sampling is noninvasive and generally more cost-efficient. Currently, eDNA approaches have been limited to single-marker studies of mitochondrial DNA (mtDNA), and the relationship between eDNA haplotype composition and true haplotype composition still needs to be thoroughly verified. This will require testing of bioinformatic and statistical software to correct for erroneous sequences, as well as biases and random variation in relative sequence abundances. However, eDNA-based population genetic methods have far-reaching potential for both basic and applied research. In this paper, we present a brief overview of the achievements of eDNA-based population genetics to date, and outline the prospects for future developments in the field, including the estimation of nuclear DNA (nuDNA) variation and epigenetic information. We discuss the challenges associated with eDNA samples as opposed to those of individual tissue samples and assess whether eDNA might offer additional types of information unobtainable with tissue samples. Lastly, we provide recommendations for determining whether an eDNA approach would be a useful and suitable choice in different research settings. We limit our discussion largely to contemporary aquatic systems, but the advantages, challenges, and perspectives can to a large degree be generalized to eDNA studies with a different spatial and temporal focus.