Intragenic DNA methylation alters chromatin structure and elongation efficiency in mammalian cells

p53 and TDG are dominant in regulating the activity of the human de novo DNA methyltransferase DNMT3A on nucleosomes

DNA methylation and histone tail modifications are interrelated mechanisms involved in a wide range of biological processes, and disruption of this crosstalk is linked to diseases such as acute myeloid leukemia. In addition, DNA methyltransferase 3A (DNMT3A) activity is modulated by several regulatory proteins, including p53 and thymine DNA glycosylase (TDG). However, the relative role of histone tails and regulatory proteins in the simultaneous coordination of DNMT3A activity remains obscure. We observed that DNMT3A binds H3 tails and p53 or TDG at distinct allosteric sites to form DNMT3A–H3 tail-p53 or –TDG multiprotein complexes. Functional characterization of DNMT3A–H3 tail-p53 or –TDG complexes on human-derived synthetic histone H3 tails, mononucleosomes, or polynucleosomes shows p53 and TDG play dominant roles in the modulation of DNMT3A activity. Intriguingly, this dominance occurs even when DNMT3A is actively methylating nucleosome substrates. The activity of histone modifiers is influenced by their ability to sense modifications on histone tails within the same nucleosome or histone tails on neighboring nucleosomes. In contrast, we show here that DNMT3A acts on DNA within a single nucleosome, on nucleosomal DNA within adjacent nucleosomes, and DNA not associated with the DNMT3A–nucleosome complex. Our findings have direct bearing on how the histone code drives changes in DNA methylation and highlight the complex interplay between histone tails, epigenetic enzymes, and modulators of enzymatic activity.

The myonuclear DNA methylome in response to an acute hypertrophic stimulus

In addition to multi-nucleated muscle fibres, numerous resident and infiltrating mononuclear cells populate the muscle compartment. As most epigenetic assays in skeletal muscle are conducted on whole tissue homogenates, essentially nothing is known about regulatory processes exclusively within muscle fibres in vivo. Utilizing a novel genetically modified mouse model developed by our laboratory, we (1) outline a simple and rapid workflow for isolating pure myonuclei from small tissue samples via fluorescent activated cell sorting and extracting high-quality large-fragment DNA for downstream analyses, and (2) provide information on myonuclear and interstitial cell nuclear CpG DNA methylation via reduced representation bisulphite sequencing (RRBS) using mice that were subjected to an acute mechanical overload of the plantaris muscle. In 3-month-old mice, myonuclei are ~50% of total nuclei in sham and ~30% in 3-d overloaded muscle, the difference being attributable to mononuclear cell infiltration and proliferation with overload. In purified myonuclei, pathway analysis of hypomethylated promoter regions following overload was distinct from interstitial nuclei and revealed marked regulation of factors that converge on the master regulator of muscle growth mTOR, and on autophagy. Specifically, acute hypomethylation of Rheb, Rictor, Hdac1, and Hdac2, in addition to a major driver of ribosome biogenesis Myc, reveals the epigenetic regulation of hypertrophic signalling within muscle fibres that may underpin the long-term growth response to loading. This study provides foundational information on global myonuclear epigenetics in vivo using RRBS, and demonstrates the importance of isolating specific nuclear populations to study the epigenetic regulation of skeletal muscle fibre adaptation.

Genome-Wide Assessment of DNA Methylation in Chicken Cardiac Tissue Exposed to Different Incubation Temperatures and CO2 Levels

Temperature and CO₂ concentration during incubation have profound effects on broiler chick development, and numerous studies have identified significant effects on hatch heart weight (HW) as a result of differences in these parameters. Early life environment has also been shown to affect broiler performance later in life; it has thus been suggested that epigenetic mechanisms may mediate long-term physiological changes induced by environmental stimuli. DNA methylation is an epigenetic modification that can confer heritable changes in gene expression. Using reduced-representation bisulfite sequencing (RRBS), we assessed DNA methylation patterns in cardiac tissue of 84 broiler hatchlings incubated at two egg shell temperatures (EST; 37.8°C and 38.9°C) and three CO₂ concentrations (0.1%, 0.4%, and 0.8%) from day 8 of incubation onward. We assessed differential methylation between EST treatments and identified 2,175 differentially methylated (DM) CpGs (1,121 hypermethylated, 1,054 hypomethylated at 38.9° vs. 37.8°) in 269 gene promoters and 949 intragenic regions. DM genes (DMGs) were associated with heart developmental processes, including cardiomyocyte proliferation and differentiation. We identified enriched binding motifs among DM loci, including those for transcription factors associated with cell proliferation and heart development among hypomethylated CpGs that suggest increased binding ability at higher EST. We identified 9,823 DM CpGs between at least two CO₂ treatments, with the greatest difference observed between 0.8 and 0.1% CO₂ that disproportionately impacted genes involved in cardiac muscle development and response to low oxygen levels. Using HW measurements from the same chicks, we performed an epigenome-wide association study (EWAS) for HW, and identified 23 significantly associated CpGs, nine of which were also DM between ESTs. We found corresponding differences in transcript abundance between ESTs in three DMGs (ABLIM2, PITX2, and THRSP). Hypomethylation of an exonic CpG in PITX2 at 38.9°C was associated with increased expression, and suggests increased cell proliferation in broiler hatchlings incubated at higher temperatures. Overall, these results identified numerous epigenetic associations between chick incubation factors and heart development that may manifest in long-term differences in animal performance.

Histone modifications, DNA methylation, and the epigenetic code of alcohol use disorder

Alcohol use disorder (AUD) is a leading cause of morbidity and mortality. Despite AUD's substantial contributions to lost economic productivity and quality of life, there are only a limited number of approved drugs for treatment of AUD in the United States. This chapter will update progress made on the epigenetic basis of AUD, with particular focus on histone post-translational modifications and DNA methylation and how these two epigenetic mechanisms interact to contribute to neuroadaptive processes leading to initiation, maintenance and progression of AUD pathophysiology. We will also evaluate epigenetic therapeutic strategies that have arisen from preclinical models of AUD and epigenetic biomarkers that have been discovered in human populations with AUD.

Comparative analysis of the transcriptome, methylome, and metabolome during pollen abortion of a seedless citrus mutant

Pollen abortion could be mainly attributed to abnormal meiosis in the mutant. Multiomics analysis uncovered significant epigenetic variations between the mutant and its wild type during the pollen abortion process. Male sterility caused by aborted pollen can result in seedless fruit. A seedless Ponkan mandarin mutant (bud sport) was used to compare the transcriptome, methylome, and metabolome with its progenitor to understand the mechanism of citrus pollen abortion. Cytological observations showed that the anther of the mutant could form microspore mother cells, although the microspores failed to develop fertile pollen at the anther dehiscence stage. Based on pollen phenotypic analysis, pollen abortion could be mainly attributed to abnormal meiosis in the mutant. A transcriptome analysis uncovered the molecular mechanisms underlying pollen abortion between the mutant and its wild type. A total of 5421 differentially expressed genes were identified, and some of these genes were involved in the meiosis, hormone biosynthesis and signaling, carbohydrate, and flavonoid pathways. A total of 50,845 differentially methylated regions corresponding to 15,426 differentially methylated genes in the genic region were found between the mutant and its wild type by the methylome analysis. The expression level of these genes was negatively correlated with their methylation level, especially in the promoter regions. In addition, 197 differential metabolites were identified between the mutant and its wild type based on the metabolome analysis. The transcription and metabolome analysis further indicated that the expression of genes in the flavonoid, carbohydrate, and hormone metabolic pathways was significantly modulated in the pollen of the mutant. These results indicated that demethylation may alleviate the silencing of carbohydrate genes in the mutant, resulting in excessive starch and sugar hydrolysis and thereby causing pollen abortion in the mutant.

Ctcf haploinsufficiency mediates intron retention in a tissue-specific manner

CTCF is a master regulator of gene transcription and chromatin organisation with occupancy at thousands of DNA target sites genome-wide. While CTCF is essential for cell survival, CTCF haploinsufficiency is associated with tumour development and hypermethylation. Increasing evidence demonstrates CTCF as a key player in several mechanisms regulating alternative splicing (AS), however, the genome-wide impact of Ctcf dosage on AS has not been investigated.

We examined the effect of Ctcf haploinsufficiency on gene expression and AS in five tissues from Ctcf hemizygous (Ctcf^+/-) mice. Reduced Ctcf levels caused distinct tissue-specific differences in gene expression and AS in all tissues. An increase in intron retention (IR) was observed in Ctcf^+/- liver and kidney. In liver, this specifically impacted genes associated with cytoskeletal organisation, splicing and metabolism. Strikingly, most differentially retained introns were short, with a high GC content and enriched in Ctcf binding sites in their proximal upstream genomic region. This study provides new insights into the effects of CTCF haploinsufficiency on organ transcriptomes and the role of CTCF in AS regulation.

Distinct Expression and Methylation Patterns for Genes with Different Fates following a Single Whole-Genome Duplication in Flowering Plants

For most sequenced flowering plants, multiple whole-genome duplications (WGDs) are found. Duplicated genes following WGD often have different fates that can quickly disappear again, be retained for long(er) periods, or subsequently undergo small-scale duplications. However, how different expression, epigenetic regulation, and functional constraints are associated with these different gene fates following a WGD still requires further investigation due to successive WGDs in angiosperms complicating the gene trajectories. In this study, we investigate lotus (Nelumbo nucifera), an angiosperm with a single WGD during the K-pg boundary. Based on improved intraspecific-synteny identification by a chromosome-level assembly, transcriptome, and bisulfite sequencing, we explore not only the fundamental distinctions in genomic features, expression, and methylation patterns of genes with different fates after a WGD but also the factors that shape post-WGD expression divergence and expression bias between duplicates. We found that after a WGD genes that returned to single copies show the highest levels and breadth of expression, gene body methylation, and intron numbers, whereas the long-retained duplicates exhibit the highest degrees of protein-protein interactions and protein lengths and the lowest methylation in gene flanking regions. For those long-retained duplicate pairs, the degree of expression divergence correlates with their sequence divergence, degree in protein-protein interactions, and expression level, whereas their biases in expression level reflecting subgenome dominance are associated with the bias of subgenome fractionation. Overall, our study on the paleopolyploid nature of lotus highlights the impact of different functional constraints on gene fate and duplicate divergence following a single WGD in plant.

Experimental Parasite Infection Causes Genome-Wide Changes in DNA Methylation

Parasites are arguably among the strongest drivers of natural selection, constraining hosts to evolve resistance and tolerance mechanisms. Although, the genetic basis of adaptation to parasite infection has been widely studied, little is known about how epigenetic changes contribute to parasite resistance and eventually, adaptation. Here, we investigated the role of host DNA methylation modifications to respond to parasite infections. In a controlled infection experiment, we used the three-spined stickleback fish, a model species for host-parasite studies, and their nematode parasite Camallanus lacustris. We showed that the levels of DNA methylation are higher in infected fish. Results furthermore suggest correlations between DNA methylation and shifts in key fitness and immune traits between infected and control fish, including respiratory burst and functional trans-generational traits such as the concentration of motile sperm. We revealed that genes associated with metabolic, developmental, and regulatory processes (cell death and apoptosis) were differentially methylated between infected and control fish. Interestingly, genes such as the neuropeptide FF receptor 2 and the integrin alpha 1 as well as molecular pathways including the Th1 and Th2 cell differentiation were hypermethylated in infected fish, suggesting parasite-mediated repression mechanisms of immune responses. Altogether, we demonstrate that parasite infection contributes to genome-wide DNA methylation modifications. Our study brings novel insights into the evolution of vertebrate immunity and suggests that epigenetic mechanisms are complementary to genetic responses against parasite-mediated selection.

Genome-Wide DNA Methylation and Gene Expression Profiles in Cows Subjected to Different Stress Level as Assessed by Cortisol in Milk

Dairy cattle health, wellbeing and productivity are deeply affected by stress. Its influence on metabolism and immune response is well known, but the underlying epigenetic mechanisms require further investigation. In this study, we compared DNA methylation and gene expression signatures between two dairy cattle populations falling in the high- and low-variant tails of the distribution of milk cortisol concentration (MC), a neuroendocrine marker of stress in dairy cows. Reduced Representation Bisulfite Sequencing was used to obtain a methylation map from blood samples of these animals. The high and low groups exhibited similar amounts of methylated CpGs, while we found differences among non-CpG sites. Significant methylation changes were detected in 248 genes. We also identified significant fold differences in the expression of 324 genes. KEGG and Gene Ontology (GO) analysis showed that genes of both groups act together in several pathways, such as nervous system activity, immune regulatory functions and glucocorticoid metabolism. These preliminary results suggest that, in livestock, cortisol secretion could act as a trigger for epigenetic regulation and that peripheral changes in methylation can provide an insight into central nervous system functions.

Comparative whole genome DNA methylation profiling across cattle tissues reveals global and tissue-specific methylation patterns

BACKGROUND

Efforts to improve animal health, and understand genetic bases for production, may benefit from a comprehensive analysis of animal genomes and epigenomes. Although DNA methylation has been well studied in humans and other model species, its distribution patterns and regulatory impacts in cattle are still largely unknown. Here, we present the largest collection of cattle DNA methylation epigenomic data to date.

RESULTS

Using Holstein cattle, we generated 29 whole genome bisulfite sequencing (WGBS) datasets for 16 tissues, 47 corresponding RNA-seq datasets, and 2 whole genome sequencing datasets. We did read mapping and DNA methylation calling based on two different cattle assemblies, demonstrating the high quality of the long-read-based assembly markedly improved DNA methylation results. We observed large differences across cattle tissues in the methylation patterns of global CpG sites, partially methylated domains (PMDs), hypomethylated regions (HMRs), CG islands (CGIs), and common repeats. We detected that each tissue had a distinct set of PMDs, which showed tissue-specific patterns. Similar to human PMD, cattle PMDs were often linked to a general decrease of gene expression and a decrease in active histone marks and related to long-range chromatin organizations, like topologically associated domains (TADs). We tested a classification of the HMRs based on their distributions relative to transcription start sites (TSSs) and detected tissue-specific TSS-HMRs and genes that showed strong tissue effects. When performing cross-species comparisons of paired genes (two opposite strand genes with their TSS located in the same HMR), we found out they were more consistently co-expressed among human, mouse, sheep, goat, yak, pig, and chicken, but showed lower consistent ratios in more divergent species. We further used these WGBS data to detect 50,023 experimentally supported CGIs across bovine tissues and found that they might function as a guard against C-to-T mutations for TSS-HMRs. Although common repeats were often heavily methylated, some young Bov-A2 repeats were hypomethylated in sperm and could affect the promoter structures by exposing potential transcription factor binding sites.

CONCLUSIONS

This study provides a comprehensive resource for bovine epigenomic research and enables new discoveries about DNA methylation and its role in complex traits.

Genome-wide DNA methylation profiles may reveal new possible epigenetic pathogenesis of sporadic congenital cataract

Aim: To investigate the possible epigenetic pathogenesis of sporadic congenital cataract. Materials & methods: We conducted whole genome bisulfite sequencing on peripheral blood from sporadic binocular or monocular congenital cataract patients and cataract-free participants. Results: We found massive differentially methylated regions within the whole genomes between any two groups. Meanwhile, we identified five genes (ACTN4, ACTG1, TUBA1A, TUBA1C, TUBB4B) for the binocular and control groups and TUBA1A for the monocular and control groups as the core differentially methylated region-related genes. The proteins encoded by these core genes are involved in building cytoskeleton and intercellular junctions. Conclusion: Changes in the methylation levels of core genes may disturb the function of cytoskeleton and intercellular junctions, eventually leading to sporadic congenital cataract.

The Landscape of DNA Methylation Associated With the Transcriptomic Network of Intramuscular Adipocytes Generates Insight Into Intramuscular Fat Deposition in Chicken

Intramuscular fat (IMF), which regulated by genetics, nutrition and environment is an important factor that influencing meat quality. Up to now, the epigenetic regulation mechanism underlying poultry IMF deposition remains poorly understood. Here, we focused on the DNA methylation, which usually regulate genes in transcription level. To look into the essential role of DNA methylation on the IMF deposition, chicken intramuscular preadipocytes were isolated and cultured in vitro, and a model of intramuscular adipocyte differentiation was constructed. Combined the whole genome bisulfite sequencing (WGBS) and RNA-Seq technologies, we identified several methylated genes, which mainly affecting fatty acid metabolism and muscle development. Furthermore, we reported that DNA methylation regulate intramuscular adipogenesis by regulating the genes, such as collagen, type VI, alpha 1 (COL6A1) thus affecting IMF deposition. Overexpression of COL6A1 increases the lipid droplet and inhibits cell proliferation by regulating CHAD and CAMK2 in intramuscular adipocytes, while knockdown of COL6A1 shows the opposite effect. Taken together, our results reveal that DNA methylation plays an important role in poultry IMF deposition.

Transcriptional regulation of genes bearing intronic heterochromatin in the rice genome

Intronic regions of eukaryotic genomes accumulate many Transposable Elements (TEs). Intronic TEs often trigger the formation of transcriptionally repressive heterochromatin, even within transcription-permissive chromatin environments. Although TE-bearing introns are widely observed in eukaryotic genomes, their epigenetic states, impacts on gene regulation and function, and their contributions to genetic diversity and evolution, remain poorly understood. In this study, we investigated the genome-wide distribution of intronic TEs and their epigenetic states in the Oryza sativa genome, where TEs comprise 35% of the genome. We found that over 10% of rice genes contain intronic heterochromatin, most of which are associated with TEs and repetitive sequences. These heterochromatic introns are longer and highly enriched in promoter-proximal positions. On the other hand, introns also accumulate hypomethylated short TEs. Genes with heterochromatic introns are implicated in various biological functions. Transcription of genes bearing intronic heterochromatin is regulated by an epigenetic mechanism involving the conserved factor OsIBM2, mutation of which results in severe developmental and reproductive defects. Furthermore, we found that heterochromatic introns evolve rapidly compared to non-heterochromatic introns. Our study demonstrates that heterochromatin is a common epigenetic feature associated with actively transcribed genes in the rice genome.

Discovery, Validation, and Application of Novel Methylated DNA Markers for Detection of Esophageal Cancer in Plasma

PURPOSE

The burden of esophageal cancer continues to rise, and noninvasive screening tools are needed. Methylated DNA markers (MDM) assayed from plasma show promise in detection of other cancers. For esophageal cancer detection, we aimed to discover and validate MDMs in tissue, and determine their feasibility when assayed from plasma.

EXPERIMENTAL DESIGN

Whole-methylome sequencing was performed on DNA extracted from 37 tissues (28 EC; 9 normal esophagus) and 8 buffy coat samples. Top MDMs were validated by methylation specific PCR on tissue from 76 EC (41 adeno, 35 squamous cell) and 17 normal esophagus. Quantitative allele-specific real-time target and signal amplification was used to assay MDMs in plasma from 183 patients (85 EC, 98 controls). Recursive partitioning (rPART) identified MDM combinations predictive of esophageal cancer. Validation was performed in silico by bootstrapping.

RESULTS

From discovery, 23 candidate MDMs were selected for independent tissue validation; median area under the receiver operating curve (AUC) for individual MDMs was 0.93. Among 12 MDMs advanced to plasma testing, rPART modeling selected a 5 MDM panel (FER1L4, ZNF671, ST8SIA1, TBX15, ARHGEF4) which achieved an AUC of 0.93 (95% CI, 0.89-0.96) on best-fit and 0.81 (95% CI, 0.75-0.88) on cross-validation. At 91% specificity, the panel detected 74% of esophageal cancer overall, and 43%, 64%, 77%, and 92% of stages I, II, III, and IV, respectively. Discrimination was not affected by age, sex, smoking, or body mass index.

CONCLUSIONS

Novel MDMs assayed from plasma detect esophageal cancer with moderate accuracy. Further optimization and clinical testing are warranted.

DNA Methylation and Susceptibility to Autism Spectrum Disorder

The prevalence of autism spectrum disorder (ASD) has been increasing steadily over the last 20 years; however, the molecular basis for the majority of ASD cases remains unknown. Recent advances in next-generation sequencing and detection of DNA modifications have made methylation-dependent regulation of transcription an attractive hypothesis for being a causative factor in ASD etiology. Evidence for abnormal DNA methylation in ASD can be seen on multiple levels, from genetic mutations in epigenetic machinery to loci-specific and genome-wide changes in DNA methylation. Epimutations in DNA methylation can be acquired throughout life, as global DNA methylation reprogramming is dynamic during embryonic development and the early postnatal period that corresponds to the peak time of synaptogenesis. However, technical advances and causative evidence still need to be established before abnormal DNA methylation and ASD can be confidently associated.

Differentially methylated gene patterns between age-matched sarcopenic and non-sarcopenic women

BACKGROUND

Sarcopenia is characterized by progressive decreases in muscle mass, muscle strength, and muscle function with ageing. Although many studies have investigated the mechanisms of sarcopenia, its connection with epigenetic factors, such as DNA methylation, still remains poorly understood. The aim of this study was to explore sarcopenia-related DNA methylation differences in blood samples between age-matched sarcopenic and non-sarcopenic older women.

METHODS

A sarcopenic group (n = 24) was identified and selected from a set of 247 older Caucasian women (aged 65-80 years) based on cut-off points of skeletal muscle index at 6.75 kg/m2 and grip strength at 26 kg (the lower quintile of grip strength in the set). A non-sarcopenic group (n = 24) was created with a similar age distribution as that of the sarcopenic group. DNA methylation patterns of whole blood samples from both groups were analysed using Infinium MethylationEPIC BeadChip arrays. Differentially methylated cytosin-phosphate-guanine sites (dmCpGs) were identified at a P value threshold of 0.01 by comparing methylation levels between the sarcopenic and non-sarcopenic groups at each CpG site. dmCpG-related genes were annotated based on Homo sapiens hg19 genome build. The functions of these genes were further examined by GO and KEGG pathway enrichment analysis.

RESULTS

The global methylation level of all analysed CpG sites (n = 788 074) showed no significant difference between the sarcopenic and non-sarcopenic groups (0.812), while the average methylation level of dmCpGs (n = 6258) was significantly lower in the sarcopenic group (0.004). The sarcopenic group had significantly higher methylation levels in TSS200 (the region from transcription start site to 200 nucleotides upstream of the site) and lower methylation levels in gene body and 3'UTR regions. In respect of CpG regions, CpG islands in promoters and some intragenic regions showed greater levels of methylation in the sarcopenic group. dmCpG-related KEGG pathways were mainly associated with muscle function, actin cytoskeleton regulation, and energy metabolism. Seven genes (HSPB1, PBX4, CNKSR3, ORMDL3, MIR10A, ZNF619, and CRADD) were found with the same methylation direction as previous studies of blood sample methylation during ageing. Fifty-four genes were shared with previous studies of resistance training.

CONCLUSIONS

Our results improve understanding of epigenetic mechanisms of sarcopenia by identifying sarcopenia-related DNA methylation differences in blood samples of older women. These methylation differences suggest underlying alterations of gene expression and pathway function, which can partially explain sarcopenia-related muscular changes.

Genetic and epigenetic effects on centromere establishment

Endogenous chromosomes contain centromeres to direct equal chromosomal segregation in mitosis and meiosis. The location and function of existing centromeres is usually maintained through cell cycles and generations. Recent studies have investigated how the centromere-specific histone H3 variant CENP-A is assembled and replenished after DNA replication to epigenetically propagate the centromere identity. However, existing centromeres occasionally become inactivated, with or without change in underlying DNA sequences, or lost after chromosomal rearrangements, resulting in acentric chromosomes. New centromeres, known as neocentromeres, may form on ectopic, non-centromeric chromosomal regions to rescue acentric chromosomes from being lost, or form dicentric chromosomes if the original centromere is still active. In addition, de novo centromeres can form after chromatinization of purified DNA that is exogenously introduced into cells. Here, we review the phenomena of naturally occurring and experimentally induced new centromeres and summarize the genetic (DNA sequence) and epigenetic features of these new centromeres. We compare the characteristics of new and native centromeres to understand whether there are different requirements for centromere establishment and propagation. Based on our understanding of the mechanisms of new centromere formation, we discuss the perspectives of developing more stably segregating human artificial chromosomes to facilitate gene delivery in therapeutics and research.

Systematic screening identifies a 2-gene signature as a high-potential prognostic marker of undifferentiated pleomorphic sarcoma/myxofibrosarcoma

The Cancer Genome Atlas (TCGA) Research Network confirmed that undifferentiated pleomorphic sarcoma (UPS) and myxofibrosarcoma (MFS) share a high level of genomic similarities and fall into a single spectrum of tumour. However, no molecular prognostic biomarkers have been identified in UPS/MFS. In this study, by extracting data from TCGA-Sarcoma (SARC), we explored relapse-related genes, their prognostic value and possible mechanisms of the dysregulations. After systematic screening, ITGA10 and PPP2R2B were included to construct a 2-gene signature. The 2-gene signature had an AUC value of 0.83 and had an independent prognostic value in relapse-free survival (RFS) (HR: 2.966, 95%CI: 1.995-4.410 P < .001), and disease-specific survival (DSS) (HR: 2.283, 95%CI: 1.358-3.835, P = .002), as a continuous variable. Gene-level copy number alterations (CNAs) were irrelevant to their dysregulation. Two CpG sites (cg15585341 and cg04126335) around the promoter of ITGA10 showed strong negative correlations with ITGA10 expression (Pearson's r < -0.6). Transcript preference was observed in PPP2R2B expression. The methylation of some CpG sites in two gene body regions showed at least moderate positive correlations (Pearson's r > .4) with PPP2R2B expression. Besides, the 2-gene signature showed a moderate negative correlation with CD4 + T cell infiltration. High-level CD4 + T cell infiltration and neutrophil infiltration were associated with significantly better RFS. Based on these findings, we infer that the 2-gene signature might be a potential prognostic marker in patients with UPS/MFS. Considering the potential benefits of immunotherapy for UPS/MFS patients, it is imperative to explore the predictive value of this signature in immunotherapeutic responses in the future.

DNA methylation, microRNA expression profiles and their relationships with transcriptome in grass-fed and grain-fed Angus cattle rumen tissue

Rumen is an organ for supplying nutrients for the growth and production of bovine, which might function differently under grass-fed and grain-fed regimens considering the association of gene expression, DNA methylation, and microRNA expression. The objective of this study was to explore the potential mechanism influencing rumen function of grass-fed and grain-fed animals. Methylated DNA binding domain sequencing (MBD-Seq) and microRNA-Seq were respectively utilized to detect the DNA methylation and microRNA expression in rumen tissue of grass-fed and grain-fed Angus cattle. Combined analysis revealed that the expression of the differentially expressed genes ADAMTS3 and ENPP3 was correlated with the methylation abundance of the corresponding differentially methylated regions (DMRs) inside these two genes, and these two genes were reported to be respectively involved in biosynthesis and regulation of glycosyltransferase activity; the differentially expressed microRNA bta-mir-122 was predicted to possibly target the differentially expressed genes OCLN and RBM47, potentially affecting the rumen function; the microRNA bta-mir-655 was exclusively detected in grain-fed group; its targets were significantly enriched in insulin and TGF-beta signaling pathways, which might worked together to regulate the function of rumen, resulting in different characteristics between grass-fed and grain-fed cattle. Collectively, our results provided insights into understanding the mechanisms determining rumen function and unraveled the biological basis underlying the economic traits to improve the productivity of animals.

Epimutations are associated with CHROMOMETHYLASE 3-induced de novo DNA methylation

In many plant species, a subset of transcribed genes are characterized by strictly CG-context DNA methylation, referred to as gene body methylation (gbM). The mechanisms that establish gbM are unclear, yet flowering plant species naturally without gbM lack the DNA methyltransferase, CMT3, which maintains CHG (H = A, C, or T) and not CG methylation at constitutive heterochromatin. Here, we identify the mechanistic basis for gbM establishment by expressing CMT3 in a species naturally lacking CMT3. CMT3 expression reconstituted gbM through a progression of de novo CHG methylation on expressed genes, followed by the accumulation of CG methylation that could be inherited even following loss of the CMT3 transgene. Thus, gbM likely originates from the simultaneous targeting of loci by pathways that promote euchromatin and heterochromatin, which primes genes for the formation of stably inherited epimutations in the form of CG DNA methylation.

Epimutations are associated with CHROMOMETHYLASE 3-induced de novo DNA methylation

In many plant species, a subset of transcribed genes are characterized by strictly CG-context DNA methylation, referred to as gene body methylation (gbM). The mechanisms that establish gbM are unclear, yet flowering plant species naturally without gbM lack the DNA methyltransferase, CMT3, which maintains CHG (H = A, C, or T) and not CG methylation at constitutive heterochromatin. Here, we identify the mechanistic basis for gbM establishment by expressing CMT3 in a species naturally lacking CMT3. CMT3 expression reconstituted gbM through a progression of de novo CHG methylation on expressed genes, followed by the accumulation of CG methylation that could be inherited even following loss of the CMT3 transgene. Thus, gbM likely originates from the simultaneous targeting of loci by pathways that promote euchromatin and heterochromatin, which primes genes for the formation of stably inherited epimutations in the form of CG DNA methylation.

Aberrant DNA methylation defines isoform usage in cancer, with functional implications

Alternative transcript isoforms are common in tumors and act as potential drivers of cancer. Mechanisms determining altered isoform expression include somatic mutations in splice regulatory sites or altered splicing factors. However, since DNA methylation is known to regulate transcriptional isoform activity in normal cells, we predicted the highly dysregulated patterns of DNA methylation present in cancer also affect isoform activity. We analyzed DNA methylation and RNA-seq isoform data from 18 human cancer types and found frequent correlations specifically within 11 cancer types. Examining the top 25% of variable methylation sites revealed that the location of the methylated CpG site in a gene determined which isoform was used. In addition, the correlated methylation-isoform patterns classified tumors into known subtypes and predicted distinct protein functions between tumor subtypes. Finally, methylation-correlated isoforms were enriched for oncogenes, tumor suppressors, and cancer-related pathways. These findings provide new insights into the functional impact of dysregulated DNA methylation in cancer and highlight the relationship between the epigenome and transcriptome.

In eukaryotes, one gene can be transcribed into multiple RNA sequences (or isoforms) that are subsequently translated into proteins with different functions in response to specific cellular needs. Recent studies showed that cancer cells can obtain abnormal functions via expressing different isoforms. In normal cells, isoform expression can be regulated by DNA methylation–a molecular signature with attached methyl groups on DNA sequences. Given that dysregulation of DNA methylation is a cancer hallmark, we suspect the same regulation holds in cancer and contributes to cancer progression. In this study, we analyzed data from 18 human cancer types and found frequent correlations in 11 cancer types between specific isoform usage and DNA methylation depending on the location of the methylated site in a gene. These correlation patterns can classify heterogeneous tumors in a cancer type into homogeneous subtypes and are predicted to change protein functions via isoform switching between subtypes. Finally, we found cancer-related genes often harbored more DNA methylation-isoform correlations than genes not implicated in cancer. This finding could help us to better understand the functional impact of DNA methylation alterations via regulation of isoform expression in tumorigenesis and to further improve the cancer treatment.

The landscape of DNA methylation associated with the transcriptomic network in layers and broilers generates insight into embryonic muscle development in chicken

Scope: As DNA methylation is one of the key epigenetic mechanisms involved in embryonic muscle development, elucidating its relationship with non-coding RNAs and genes is essential for understanding early muscle development. The methylome profiles of pre-hatching chicken across multiple developmental stages remain incomplete although several related studies have been reported. Methods: In this study, we performed single-base-resolution bisulfite sequencing together with RNA-seq of broilers and layers in different embryonic development points (E10, E13, E16 and E19) to explore the genetic basis of embryonic muscle development in chicken. The differential methylated regions and novel lncRNAs were identified for association analyses. Through genomic position and correlation analysis between DMRs and lncRNAs, the target lncRNAs were detected to participate in the embryonic muscle formation and the results were then verified in vitro experiments. Results: Comparison of methylome profiles between two chicken lines revealed that lower methylation in broilers might contribute to muscle development in embryonic period. Differential methylated region analysis showed that the majority of differential methylated regions were hypo-DMRs for broilers. Differential methylated genes were significantly enriched in muscle development-related terms at E13 and E19. Furthermore, we identified a long non-coding RNA MyH1-AS that potentially regulated embryonic muscle development, proved by the regulatory network construction and further in vitro experiments. Conclusion: Our study revealed an integrative landscape of middle- to late-stage of embryonic myogenesis in chicken, gave rise to a comprehensive understanding of epigenetic and transcriptional regulation in muscle development. Moreover, we provided a reliable data resource for further embryonic muscle development studies.

Pan-cancer genomic analysis links 3'UTR DNA methylation with increased gene expression in T cells

Background

Investigations into the function of non-promoter DNA methylation have yielded new insights into the epigenetic regulation of gene expression. However, integrated genome-wide non-promoter DNA methylation and gene expression analyses across a wide number of tumour types and corresponding normal tissues have not been performed.

Methods

To investigate the impact of non-promoter DNA methylation on cancer pathogenesis, we performed a large-scale analysis of gene expression and DNA methylation profiles, finding enrichment in the 3’UTR DNA methylation positively correlated with gene expression. Filtering for genes in which 3’UTR DNA methylation strongly correlated with gene expression yielded a list of genes enriched for functions involving T cell activation.

Findings

The important immune checkpoint gene Havcr2 showed a substantial increase in 3’UTR DNA methylation upon T cell activation and subsequent upregulation of gene expression in mice. Furthermore, this increase in Havcr2 gene expression was abrogated by treatment with decitabine.

Interpretation

These findings indicate that the 3’UTR is a functionally relevant DNA methylation site. Additionally, we show a potential novel mechanism of HAVCR2 regulation in T cells, providing new insights for modulating immune checkpoint blockade.

Linking endotypes to omics profiles in difficult-to-control asthma using the diagnostic Chinese medicine syndrome differentiation algorithm

Objective: Patients with difficult-to-control asthma have difficulty breathing almost all of the time, even leading to life-threatening asthma attacks. However, only few diagnostic markers for this disease have been identified. We aimed to take advantage of unique Chinese medicine theories for phenotypic classification and to explore molecular signatures in difficult-to-control asthma. Methods: The Chinese medicine syndrome differentiation algorithm (CMSDA) is a syndrome-scoring classification method based on the Chinese medicine overall observation theory. Patients with difficult-to-control asthma were classified into Cold- and Hot-pattern groups according to the CMSDA. DNA methylation and metabolomic profiles were obtained using Infinium Human Methylation 450 BeadChip and gas chromatography-mass spectrometer. Subsequently, an integrated bioinformatics analysis was performed to compare those two patterns and identify Cold/Hot-associated candidates, followed by functional validation studies. Results: A total of 20 patients with difficult-to-control asthma were enrolled in the study. Ten were grouped as Cold and 10 as Hot according to the CMSDA. We identified distinct whole-genome DNA methylation and metabolomic profiles between Cold- and Hot-pattern groups. ALDH3A1 gene exhibited variations in the DNA methylation probe cg10791966, while two metabolic pathways were associated with those two patterns. Conclusions: Our study introduced a novel diagnostic classification approach, the CMSDA, for difficult-to-control asthma. This is an alternative way to categorize diverse syndromes and link endotypes with omics profiles of this disease. ALDH3A1 might be a potential biomarker for precision diagnosis of difficult-to-control asthma.

Hypomethylation in the promoter region of ZPBP as a potential litter size indicator in Berkshire pigs

In pigs, litter size is typically defined as the total number of piglets born (TNB) or the number of piglets born alive (NBA). Increasing pig litter size is of great economic interest as a means to increase productivity. The capacity of the uterus is a critical component of litter size and may play a central role in prolificacy. In this study, we investigated litter-size-related epigenetic markers in uterine tissue from Berkshire pigs with smaller litter size groups (SLGs) and larger litter size groups (LLGs) using genome-wide bisulfite sequencing (GWBS). A total of 3269 differentially methylated regions (DMRs) were identified: 1566 were hypermethylated and 1703 hypomethylated in LLG compared to SLG. The zona pellucida binding protein (ZPBP) gene was significantly hypomethylated in the LLG promoter region, and its expression was significantly upregulated in uterine tissue. Thus, the methylation status of ZPBP gene was identified as a potential indicator of litter size. Furthermore, we verified its negative correlation with litter size traits (TNB and NBA) in whole blood samples from 172 Berkshire sows as a blood-based biomarker by a porcine methylation-specific restriction enzyme polymerase chain reaction (PMP) assay. The results suggest that the methylation status of the ZPBP gene can serve as a valuable epigenetic biomarker for hyperprolific sows.

New guidelines for DNA methylome studies regarding 5-hydroxymethylcytosine for understanding transcriptional regulation

Many DNA methylome profiling methods cannot distinguish between 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC). Because 5mC typically acts as a repressive mark whereas 5hmC is an intermediate form during active demethylation, the inability to separate their signals could lead to incorrect interpretation of the data. Is the extra information contained in 5hmC signals worth the additional experimental and computational costs? Here we combine whole-genome bisulfite sequencing (WGBS) and oxidative WGBS (oxWGBS) data in various human tissues to investigate the quantitative relationships between gene expression and the two forms of DNA methylation at promoters, transcript bodies, and immediate downstream regions. We find that 5mC and 5hmC signals correlate with gene expression in the same direction in most samples. Considering both types of signals increases the accuracy of expression levels inferred from methylation data by a median of 18.2% as compared to having only WGBS data, showing that the two forms of methylation provide complementary information about gene expression. Differential analysis between matched tumor and normal pairs is particularly affected by the superposition of 5mC and 5hmC signals in WGBS data, with at least 25%–40% of the differentially methylated regions (DMRs) identified from 5mC signals not detected from WGBS data. Our results also confirm a previous finding that methylation signals at transcript bodies are more indicative of gene expression levels than promoter methylation signals. Overall, our study provides data for evaluating the cost-effectiveness of some experimental and analysis options in the study of DNA methylation in normal and cancer samples.

The long non-coding RNA Kcnq1ot1 controls maternal p57 expression in muscle cells by promoting H3K27me3 accumulation to an intragenic MyoD-binding region

BACKGROUND

The cell-cycle inhibitor p57^kip2 plays a critical role in mammalian development by coordinating cell proliferation and differentiation in many cell types. p57^kip2 expression is finely regulated by several epigenetic mechanisms, including paternal imprinting. Kcnq1ot1, a long non-coding RNA (LncRNA), whose gene maps to the p57^Kip2 imprinting domain, is expressed exclusively from the paternal allele and participates in the cis-silencing of the neighboring imprinted genes through chromatin-level regulation. In light of our previous evidence of a functional interaction between myogenic factors and imprinting control elements in the regulation of the maternal p57^Kip2 allele during muscle differentiation, we examined the possibility that also Kcnq1ot1 could play an imprinting-independent role in the control of p57^Kip2 expression in muscle cells.

RESULTS

We found that Kcnq1ot1 depletion by siRNA causes the upregulation of the maternal and functional p57^Kip2 allele during differentiation, suggesting a previously undisclosed role for this LncRNA. Consistently, Chromatin Oligo-affinity Precipitation assays showed that Kcnq1ot1 physically interacts not only with the paternal imprinting control region of the locus, as already known, but also with both maternal and paternal alleles of a novel p57^Kip2 regulatory region, located intragenically and containing two binding sites for the muscle-specific factor MyoD. Moreover, chromatin immunoprecipitation assays after Kcnq1ot1 depletion demonstrated that the LncRNA is required for the accumulation of H3K27me3, a chromatin modification catalyzed by the histone-methyl-transferase EZH2, at the maternal p57^kip2 intragenic region. Finally, upon differentiation, the binding of MyoD to this region and its physical interaction with Kcnq1ot1, analyzed by ChIP and RNA immunoprecipitation assays, correlate with the loss of EZH2 and H3K27me3 from chromatin and with p57^Kip2 de-repression.

CONCLUSIONS

These findings highlight the existence of an imprinting-independent role of Kcnq1ot1, adding new insights into the biology of a still mysterious LncRNA. Moreover, they expand our knowledge about the molecular mechanisms underlying the tight and fine regulation of p57^Kip2 during differentiation and, possibly, its aberrant silencing observed in several pathologic conditions.

Can epigenetics translate environmental cues into phenotypes?

Living organisms are constantly exposed to wide ranges of environmental cues. They react to these cues by undergoing a battery of phenotypic responses, such as by altering their physiological and behavioral traits, in order to adapt and survive in the changed environments. The adaptive response of a species induced by environmental cues is typically thought to be associated with its genetic diversity such that higher genetic diversity provides increased adaptive potential. This originates from the general consensus that phenotypic traits have a genetic basis and are subject to Darwinian natural selection and Mendelian inheritance. There is no doubt about the validity of these principles, supported by the successful introgression of specific traits during (selective) breeding. However, a range of recent studies provided fascinating evidences suggesting that environmental effects experienced by an organism during its lifetime can have marked influences on its phenotype, and additionally the organism can pass on the acquired phenotypes to its subsequent generations through non-genetic mechanisms (also termed as epigenetic mechanism) - a notion that dates back to Lamarck and has been controversial ever since. In this review, we describe how the epigenetics has reshaped our long perception about the inheritance/development of phenotypes within organisms, contrasting with the classical gene-based view of inheritance. We particularly highlighted recent developments in our understanding of inheritance of parental environmental induced phenotypic traits in multicellular organisms under different environmental conditions, and discuss how modifications of the epigenome contribute to the determination of the adult phenotype of future generations.

Expression of Sclerostin in Osteoporotic Fracture Patients Is Associated with DNA Methylation in the CpG Island of the SOST Gene

PURPOSE

SOST gene is one of the key factors in regulating bone absorption. Although there are reports showing diverse transcription factors, epigenetic modification could be responsible for regulating SOST gene expression. There is still little exploration on promoter methylation status of SOST gene in osteoporotic bone tissues. The aim of this study is to investigate the involvement of CpG methylation in regulation of SOST expression in patients with primary osteoporosis.

METHODS

The diagnosis of osteoporosis was established on the basis of dual energy X-ray absorptiometry to measure BMD. All femoral bone tissues were separated in surgeries. After extracting total RNA and protein, we checked the relative expression levels of SOST by quantitative real-time PCR and western blot. Also, immunohistochemical staining was performed to observe the expression of SOST protein in the bone samples. The genomic DNA of non-OPF (non-osteoporotic fracture bone tissues) and OPF (osteoporotic fracture bone tissues) were treated by bisulfite modification, and methylation status of CpG sites in the CpG island of SOST gene promoter was determined by DNA sequencing.

RESULTS

SOST gene expression in the non-OPF group was lower than that in OPF group. Bisulfite sequencing result showed that SOST gene promoter was slightly demethylated in the OPF group, as compared with non-OPF group.

CONCLUSION

Our study demonstrated that DNA methylation influenced the transcriptional expression of SOST gene, which probably may play an important role in the pathogenesis of primary osteoporosis.

Comparative whole genome DNA methylation profiling of cattle sperm and somatic tissues reveals striking hypomethylated patterns in sperm

Background

Although sperm DNA methylation has been studied in humans and other species, its status in cattle is largely unknown.

Results

Using whole-genome bisulfite sequencing (WGBS), we profiled the DNA methylome of cattle sperm through comparison with three somatic tissues (mammary gland, brain, and blood). Large differences between cattle sperm and somatic cells were observed in the methylation patterns of global CpGs, pericentromeric satellites, partially methylated domains (PMDs), hypomethylated regions (HMRs), and common repeats. As expected, we observed low methylation in the promoter regions and high methylation in the bodies of active genes. We detected selective hypomethylation of megabase domains of centromeric satellite clusters, which may be related to chromosome segregation during meiosis and their rapid transcriptional activation upon fertilization. We found more PMDs in sperm cells than in somatic cells and identified meiosis-related genes such asKIF2B and REPIN1, which are hypomethylated in sperm but hypermethylated in somatic cells. In addition to the common HMRs around gene promoters, which showed substantial differences between sperm and somatic cells, the sperm-specific HMRs also targeted to distinct spermatogenesis-related genes, including BOLL, MAEL, ASZ1, SYCP3, CTCFL, MND1, SPATA22, PLD6, DDX4, RBBP8, FKBP6, and SYCE1. Although common repeats were heavily methylated in both sperm and somatic cells, some young Bov-A2 repeats, which belong to the SINE family, were hypomethylated in sperm and could affect the promoter structures by introducing new regulatory elements.

Conclusions

Our study provides a comprehensive resource for bovine sperm epigenomic research and enables new discoveries about DNA methylation and its role in male fertility.

The methylome of the marbled crayfish links gene body methylation to stable expression of poorly accessible genes

Background

The parthenogenetic marbled crayfish (Procambarus virginalis) is a novel species that has rapidly invaded and colonized various different habitats. Adaptation to different environments appears to be independent of the selection of genetic variants, but epigenetic programming of the marbled crayfish genome remains to be understood.

Results

Here, we provide a comprehensive analysis of DNA methylation in marbled crayfish. Whole-genome bisulfite sequencing of multiple replicates and different tissues revealed a methylation pattern that is characterized by gene body methylation of housekeeping genes. Interestingly, this pattern was largely tissue invariant, suggesting a function that is unrelated to cell fate specification. Indeed, integrative analysis of DNA methylation, chromatin accessibility and mRNA expression patterns revealed that gene body methylation correlated with limited chromatin accessibility and stable gene expression, while low-methylated genes often resided in chromatin with higher accessibility and showed increased expression variation. Interestingly, marbled crayfish also showed reduced gene body methylation and higher gene expression variability when compared with their noninvasive mother species, Procambarus fallax.

Conclusions

Our results provide novel insights into invertebrate gene body methylation and its potential role in adaptive gene regulation.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0229-6) contains supplementary material, which is available to authorized users.

Glia-specific APOE epigenetic changes in the Alzheimer's disease brain

The apolipoprotein E gene (APOE) is the strongest genetic risk factor for developing Alzheimer's disease (AD). Our recent identification of altered APOE DNA methylation in AD postmortem brain (PMB) prompted this follow-up study. Our goals were to (i) validate the AD-differential methylation of APOE in an independent PMB study cohort and (ii) determine the cellular populations (i.e., neuronal vs. non-neuronal) of AD PMB that contribute to this differential methylation. Here, we obtained an independent cohort of 57 PMB (42 AD and 15 controls) and quantified their APOE methylation levels from frontal lobe and cerebellar tissue. We also applied fluorescence-activated nuclei sorting (FANS) to separate neuronal nuclei from non-neuronal nuclei within the tissue of 15 AD and 14 control subjects. Bisulfite pyrosequencing was used to generate DNA methylation profiles of APOE from both bulk PMB and FANS nuclei. Our results provide independent validation that the APOE CGI holds lower DNA methylation levels in AD compared to control in frontal lobe but not cerebellar tissue. Our data also indicate that the non-neuronal cells of the AD brain, which are mainly composed of glia, are the main contributors to the lower APOE DNA methylation observed in AD PMB. Given that astrocytes are the primary producers of ApoE in the brain our results suggest that alteration of epigenetically regulated APOE expression in glia could be an important part of APOE's strong effect on AD risk.

Investigation of epigenetic regulatory networks associated with autism spectrum disorder (ASD) by integrated global LINE-1 methylation and gene expression profiling analyses

BACKGROUND

The exact cause and mechanisms underlying the pathobiology of autism spectrum disorder (ASD) remain unclear. Dysregulation of long interspersed element-1 (LINE-1) has been reported in the brains of ASD-like mutant mice and ASD brain tissues. However, the role and methylation of LINE-1 in individuals with ASD remain unclear. In this study, we aimed to investigate whether LINE-1 insertion is associated with differentially expressed genes (DEGs) and to assess LINE-1 methylation in ASD.

METHODS

To identify DEGs associated with LINE-1 in ASD, we reanalyzed previously published transcriptome profiles and overlapped them with the list of LINE-1-containing genes from the TranspoGene database. An Ingenuity Pathway Analysis (IPA) of DEGs associated with LINE-1 insertion was conducted. DNA methylation of LINE-1 was assessed via combined bisulfite restriction analysis (COBRA) of lymphoblastoid cell lines from ASD individuals and unaffected individuals, and the methylation levels were correlated with the expression levels of LINE-1 and two LINE-1-inserted DEGs, C1orf27 and ARMC8.

RESULTS

We found that LINE-1 insertion was significantly associated with DEGs in ASD. The IPA showed that LINE-1-inserted DEGs were associated with ASD-related mechanisms, including sex hormone receptor signaling and axon guidance signaling. Moreover, we observed that the LINE-1 methylation level was significantly reduced in lymphoblastoid cell lines from ASD individuals with severe language impairment and was inversely correlated with the transcript level. The methylation level of LINE-1 was also correlated with the expression of the LINE-1-inserted DEG C1orf27 but not ARMC8.

CONCLUSIONS

In ASD individuals with severe language impairment, LINE-1 methylation was reduced and correlated with the expression levels of LINE-1 and the LINE-1-inserted DEG C1orf27. Our findings highlight the association of LINE-1 with DEGs in ASD blood samples and warrant further investigation. The molecular mechanisms of LINE-1 and the effects of its methylation in ASD pathobiology deserve further study.

Integrated analysis of methylome, transcriptome and miRNAome of three pig breeds

AIM

Integrated analysis of methylome and transcriptome may help understand the molecular basis of the different breeds with different traits of commercial interest.

MATERIALS & METHODS

We obtained the first genome-wide methylome with single-base resolution, miRNAome and transcriptome from three swine breeds.

RESULTS

We displayed the landscape of the three omics in the whole-genome level. Integrated outcomes of methylome with genetic selection, miRNAome and transcriptome are also provided. Finally, we identified 11 candidate differentially methylated genes associated with phenotype variance in pigs.

CONCLUSION

DNA methylation not only suppresses transcriptome but also miRNAome. The different -omics data have complicated interaction in directly or indirectly and exhibited close relations with the distinct phenotypic traits of growth, disease resistance and energy metabolism.

DNA methylation profiles correlated to striped bass sperm fertility

Background

Striped bass (Morone saxatilis) spermatozoa are used to fertilize in vitro the eggs of white bass (M. chrysops) to produce the preferred hybrid for the striped bass aquaculture industry. Currently, only one source of domestic striped bass juveniles is available to growers that is not obtained from wild-caught parents and is thus devoid of any genetic improvement in phenotypic traits of importance to aquaculture. Sperm epigenetic modification has been predicted to be associated with fertility, which could switch genes on and off without changing the DNA sequence itself. DNA methylation is one of the most common epigenetic modification types and changes in sperm epigenetics can be correlated to sub-fertility or infertility in male striped bass. The objective of this study was to find the differentially methylated regions (DMRs) between high-fertility and sub-fertility male striped bass, which could potentially regulate the fertility performance.

Results

In our present study, we performed DNA methylation analysis of high-fertility and sub-fertility striped bass spermatozoa through MBD-Seq methods. A total of 171 DMRs were discovered in striped bass sperm correlated to fertility. Based on the annotation of these DMRs, we conducted a functional classification analysis and two important groups of genes including the WDR3/UTP12 and GPCR families, were discovered to be related to fertility performance of striped bass. Proteins from the WDR3/UTP12 family are involved in forming the sperm flagella apparatus in vertebrates and GPCRs are involved in hormonal signaling and regulation of tissue development, proliferation and differentiation.

Conclusions

Our results contribute insights into understanding the mechanism of fertility in striped bass, which will provide powerful tools to maximize reproductive efficiencies and to identify those males with superior gametes for this important aquaculture species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4548-6) contains supplementary material, which is available to authorized users.

Non-homologous end joining induced alterations in DNA methylation: A source of permanent epigenetic change

In addition to genetic mutations, epigenetic revision plays a major role in the development and progression of cancer; specifically, inappropriate DNA methylation or demethylation of CpG residues may alter the expression of genes that promote tumorigenesis. We hypothesize that DNA repair, specifically the repair of DNA double strand breaks (DSB) by Non-Homologous End Joining (NHEJ) may play a role in this process. Using a GFP reporter system inserted into the genome of HeLa cells, we are able to induce targeted DNA damage that enables the cells, after successfully undergoing NHEJ repair, to express WT GFP. These GFP+ cells were segregated into two expression classes, one with robust expression (Bright) and the other with reduced expression (Dim). Using a DNA hypomethylating drug (AzadC) we demonstrated that the different GFP expression levels was due to differential methylation statuses of CpGs in regions on either side of the break site. Deep sequencing analysis of this area in sorted Bright and Dim populations revealed a collection of different epi-alleles that display patterns of DNA methylation following repair by NHEJ. These patterns differ between Bright and Dim cells which are hypo- and hypermethylated, respectively, and between the post-repair populations and the original, uncut cells. These data suggest that NHEJ repair facilitates a rewrite of the methylation landscape in repaired genes, elucidating a potential source for the altered methylation patterns seen in cancer cells, and understanding the mechanism by which this occurs could provide new therapeutic targets for preventing this process from contributing to tumorigenesis.

Whole genome bisulfite sequencing reveals unique adaptations to high-altitude environments in Tibetan chickens

BACKGROUND

Tibetan chickens living at high altitudes show specific adaptations to high-altitude conditions, but the epigenetic modifications associated with these adaptations have not been characterized.

RESULTS

We investigated the genome-wide DNA methylation patterns in Tibetan chicken blood by using whole genome bisulfite sequencing. Generally, Tibetan chickens exhibited analogous methylation patterns to that of lowland chickens. A total of 3.92% of genomic cytosines were methylcytosines and 51.22% of cytosines in CG contexts were methylated, which was less than those in lowland chicken (55.69%). Moreover, the base adjacent to the methylcytosines of mCHGs in Tibetan chickens had a preference for T, which was different from that in lowland chickens. In Tibetan chickens, the methylation levels in the promoter were relatively low, while the gene body was also maintained in a hypomethylated state. DNA methylation levels in regions upstream of the transcription start site of genes were negatively correlated with the level of gene expression, and DNA methylation of gene body regions was also negatively related to gene expression.

CONCLUSIONS

We generated the genome-wide DNA methylation patterns in Tibetan chickens and our results will be helpful for future epigenetic studies related to adaptations to high-altitude conditions.

The long non-coding RNA HOTAIR is transcriptionally activated by HOXA9 and is an independent prognostic marker in patients with malignant glioma

The lncRNA HOTAIR has been implicated in several human cancers. Here, we evaluated the molecular alterations and upstream regulatory mechanisms of HOTAIR in glioma, the most common primary brain tumors, and its clinical relevance. HOTAIR gene expression, methylation, copy-number and prognostic value were investigated in human gliomas integrating data from online datasets and our cohorts. High levels of HOTAIR were associated with higher grades of glioma, particularly IDH wild-type cases. Mechanistically, HOTAIR was overexpressed in a gene dosage-independent manner, while DNA methylation levels of particular CpGs in HOTAIR locus were associated with HOTAIR expression levels in GBM clinical specimens and cell lines. Concordantly, the demethylating agent 5-Aza-2'-deoxycytidine affected HOTAIR transcriptional levels in a cell line-dependent manner. Importantly, HOTAIR was frequently co-expressed with HOXA9 in high-grade gliomas from TCGA, Oncomine, and our Portuguese and French datasets. Integrated in silico analyses, chromatin immunoprecipitation, and qPCR data showed that HOXA9 binds directly to the promoter of HOTAIR. Clinically, GBM patients with high HOTAIR expression had a significantly reduced overall survival, independently of other prognostic variables. In summary, this work reveals HOXA9 as a novel direct regulator of HOTAIR, and establishes HOTAIR as an independent prognostic marker, providing new therapeutic opportunities to treat this highly aggressive cancer.

Epigenetic modulation of Fgf21 in the perinatal mouse liver ameliorates diet-induced obesity in adulthood

The nutritional environment to which animals are exposed in early life can lead to epigenetic changes in the genome that influence the risk of obesity in later life. Here, we demonstrate that the fibroblast growth factor-21 gene (Fgf21) is subject to peroxisome proliferator-activated receptor (PPAR) α-dependent DNA demethylation in the liver during the postnatal period. Reductions in Fgf21 methylation can be enhanced via pharmacologic activation of PPARα during the suckling period. We also reveal that the DNA methylation status of Fgf21, once established in early life, is relatively stable and persists into adulthood. Reduced DNA methylation is associated with enhanced induction of hepatic FGF21 expression after PPARα activation, which may partly explain the attenuation of diet-induced obesity in adulthood. We propose that Fgf21 methylation represents a form of epigenetic memory that persists into adulthood, and it may have a role in the developmental programming of obesity.

Genome-Wide Analysis on the Landscape of Transcriptomes and Their Relationship With DNA Methylomes in the Hypothalamus Reveals Genes Related to Sexual Precocity in Jining Gray Goats

The Jining Gray goat is famous for its sexual precocity; however, the exact regulatory mechanism is still unknown. The hypothalamus is the key centrum in the process of animal reproduction, especially in signal transduction, and the initiation of puberty. The identification of potential genes and pathways in the hypothalamus of Jining Gray goat is critical to understanding the regulatory mechanism of sexual precocity in these goats. In this study, mRNA transcriptome analysis of the hypothalamus of juvenile and pubertal goats revealed eight genes (NTS, ADORA1, CRH, UCN3, E2F2, PDGFRB, GNRH1, and CACNA1C) and three pathways [neuroactive ligand-receptor interaction; gonadotropin-releasing hormone (GnRH) signal; melanoma] that are involved in this regulation. Subsequent methylation analysis on differentially methylated region (DMR) genes revealed the potential regulation network that influences pubertal onset. Correlation analysis verified the methylation level of some DMR genes correlates negatively with expression level. Integrated analysis between transcriptomes and methylomes identified 80 candidate genes involved in GnRH and neuroactive ligand signal pathways, of which CACNA1C and CRH were differentially expressed genes (DEGs) influenced by methylation level. The GnRH gene was the only DEG not affected by its methylation level. In summary, in this study, we identified eight genes and three pathways that are related to pubertal onset in Jining Gray goats, and the expression of CACNA1C and CRH genes of the GnRH and neuroactive ligand signal pathways were influenced by DNA methylation, while that of the GnRH gene was not affected.

Cell subpopulation deconvolution reveals breast cancer heterogeneity based on DNA methylation signature

Tumour heterogeneity describes the coexistence of divergent tumour cell clones within tumours, which is often caused by underlying epigenetic changes. DNA methylation is commonly regarded as a significant regulator that differs across cells and tissues. In this study, we comprehensively reviewed research progress on estimating of tumour heterogeneity. Bioinformatics-based analysis of DNA methylation has revealed the evolutionary relationships between breast cancer cell lines and tissues. Further analysis of the DNA methylation profiles in 33 breast cancer-related cell lines identified cell line-specific methylation patterns. Next, we reviewed the computational methods in inferring clonal evolution of tumours from different perspectives and then proposed a deconvolution strategy for modelling cell subclonal populations dynamics in breast cancer tissues based on DNA methylation. Further analysis of simulated cancer tissues and real cell lines revealed that this approach exhibits satisfactory performance and relative stability in estimating the composition and proportions of cellular subpopulations. The application of this strategy to breast cancer individuals of the Cancer Genome Atlas's identified different cellular subpopulations with distinct molecular phenotypes. Moreover, the current and potential future applications of this deconvolution strategy to clinical breast cancer research are discussed, and emphasis was placed on the DNA methylation-based recognition of intra-tumour heterogeneity. The wide use of these methods for estimating heterogeneity to further clinical cohorts will improve our understanding of neoplastic progression and the design of therapeutic interventions for treating breast cancer and other malignancies.

The effects of DNA methylation on human psychology

DNA methylation is a fundamental epigenetic modification in the human genome; pivotal in development, genomic imprinting, X inactivation, chromosome stability, gene expression and methylation aberrations are involved in an array of human diseases. Methylation at promoters is associated with transcriptional repression, whereas gene body methylation is generally associated with gene expression. Extrinsic factors such as age, diets and lifestyle affect DNA methylation which consequently alters gene expression. Stress, anxiety, depression, life satisfaction, emotion among numerous other psychological factors also modify DNA methylation patterns. This correlation is frequently investigated in four candidate genes; NR3C1, SLC6A4, BDNF and OXTR, since regulation of these genes directly impact responses to social situations, stress, threats, behaviour and neural functions. Such studies underpin the hypothesis that DNA methylation is involved in deviant human behaviour, psychological and psychiatric conditions. These candidate genes may be targeted in future to assess the correlation between methylation, social experiences and long-term behavioural phenotypes in humans; and may potentially serve as biomarkers for therapeutic intervention.

Genetic and epigenetic regulation of major histocompatibility complex class I gene expression in bovine trophoblast cells

Problem

The regulatory mechanisms governing differential expression of classical major histocompatibility complex (MHC) class I (MHC‐Ia) and non‐classical MHC class I (MHC‐Ib) genes are poorly understood.

Method of study

Quantitative reverse transcription‐ polymerase chain reaction (PCR) was used to compare the abundance of MHC‐I transcripts and related transcription factors in peripheral blood mononuclear cells (PBMC) and placental trophoblast cells (PTC). Methylation of MHC‐I CpG islands was detected by bisulfite treatment and next‐generation sequencing. Demethylation of PBMC and PTC with 5′‐aza‐deoxycytidine was used to assess the role of methylation in gene regulation.

Results

MHC‐I expression was higher in PBMC than PTC and was correlated with expression of IRF1, class II MHC transactivator (CIITA), and STAT1. The MHC‐Ia genes and BoLA‐NC1 were devoid of CpG methylation in PBMC and PTC. In contrast, CpG sites in the gene body of BoLA‐NC2, ‐NC3, and ‐NC4 were highly methylated in PBMC but largely unmethylated in normal PTC and moderately methylated in somatic cell nuclear transfer PTC. In PBMC, demethylation resulted in upregulation of MHC‐Ib by 2.8‐ to 6‐fold, whereas MHC‐Ia transcripts were elevated less than 2‐fold.

Conclusion

DNA methylation regulates bovine MHC‐Ib expression and is likely responsible for the different relative levels of MHC‐Ib to MHC‐Ia transcripts in PBMC and PTC.

Epigenetic silencing of PRSS3 provides growth and metastasis advantage for human hepatocellular carcinoma

Protease, serine, 3 (PRSS3), a member of the trypsin family of serine proteases, has been shown to be aberrantly expressed in several cancer types and to play important roles in tumor progression and metastasis. However, the expression and function of PRSS3 gene in hepatocellular carcinoma (HCC) remain unclear. Here we found that PRSS3 expression was decreased in human HCC cell lines and HCC surgical specimens. This was associated with intragenic methylation of PRSS3 gene. Treatment with DNA methyltransferase inhibitor 5-aza-2'-deoxycytidine and/or histone deacetylase inhibitor trichostatin A restored PRSS3 expression in HCC cell lines. Ectopic overexpression of PRSS3 gene in HCC cell lines significantly suppressed cell proliferation and colony formation and arrested cell cycle at G1/S phase, accompanied with downregulation of cyclin D1 (CCND1)/CDK4 and cyclin E1 (CCNE1)/CDK2 complexes. Moreover, PRSS3 overexpression in HCC cells inhibited HCC cell migration and invasion with downregulation of matrix metallopeptidase 2 (MMP2). Further study showed that PRSS3 overexpression diminished the phosphorylation of mitogen-activated protein kinase/extracellular-signal-regulated kinase signaling protein, mitogen-activated protein kinase kinase 1 (MEK1)/mitogen-activated protein kinase kinase 2 (MEK2) and extracellular-signal related kinase 1 (ERK1)/extracellular-signal related kinase 2 (ERK2), in HCC cells. In contrast, knockdown of PRSS3 by small interfering RNA resulted in opposite effects on an HCC cell line SNU-387 which constitutively expresses PRSS3. These results demonstrate that downregulation of PRSS3 by intragenic hypermethylation provides growth and metastasis advantage to HCC cells. The clinical relevance of PRSS3 to human HCC was shown by the intragenic methylation of PRSS3 in HCC specimens and its association with poor tumor differentiation in patients with HCC. Thus, PRSS3 is a potential prognostic biomarker and an epigenetic target for intervention of human HCC.

KEY MESSAGES

• PRSS3 is downregulated by intragenic hypermethylation in HCC. • Epigenetic silencing of PRSS3 facilitates growth, migration, and invasion of HCC. • PRSS3 intragenic methylation has implication in diagnosis of HCC.