Search

Abstract Detecting aberrant cell-free DNA (cfDNA) methylation is a promising strategy for lung cancer diagnosis. In this study, our aim is to identify methylation markers to distinguish patients with lung cancer from healthy individuals. Additionally, we sought to develop a deep learning model incorporating cfDNA methylation and fragment size profiles. To achieve this, we utilized methylation data collected from The Cancer Genome Atlas and Gene Expression Omnibus databases. Then we generated methylated DNA immunoprecipitation sequencing and genome-wide Enzymatic Methyl-seq (EM-seq) form lung cancer tissue and plasma. Using these data, we selected 366 methylation markers. A targeted EM-seq panel was designed using the selected markers, and 142 lung cancer and 56 healthy samples were produced with the panel. Additionally, cfDNA samples from healthy individuals and lung cancer patients were diluted to evaluate sensitivity. Its lung cancer detection performance reached an accuracy of 81.5% and an area under the receiver operating characteristic curve of 0.87. In the serial dilution experiment, we achieved tumor fraction detection of 1% at 98% specificity and 0.1% at 80% specificity. In conclusion, we successfully developed and validated a combination of methylation panel and a deep learning model that can distinguish between patients with lung cancer and healthy individuals.

Abstract Ovarian cancer, among all gynecologic malignancies, exhibits the highest incidence and mortality rate, primarily because it is often presents with non-specific or no symptoms during its early stages. For the advancement of Ovarian Cancer Diagnosis, it is crucial to identify the potential molecular signatures that could significantly differentiate between healthy and ovarian cancerous tissues and can be used further as a diagnostic biomarker for detecting ovarian cancer. In this study, we investigated the genome-wide methylation patterns in ovarian cancer patients using Methylated DNA Immunoprecipitation (MeDIP-Seq) followed by NGS. Identified differentially methylated regions (DMRs) were further validated by targeted bisulfite sequencing for CpG site-specific methylation profiles. Furthermore, expression validation of six genes by Quantitative Reverse Transcriptase-PCR was also performed. Out of total 120 differentially methylated genes (DMGs), 68 genes were hypermethylated, and 52 were hypomethylated in their promoter region. After analysis, we identified the top 6 hub genes, namely POLR3B, PLXND1, GIGYF2, STK4, BMP2 and CRKL. Interestingly we observed Non-CpG site methylation in the case of POLR3B and CRKL which was statistically significant in discriminating ovarian cancer samples from normal controls. The most significant pathways identified were focal adhesion, the MAPK signaling pathway, and the Ras signaling pathway. Expression analysis of hypermethylated genes was correlated with the downregulation of the genes. POLR3B and GIGYF2 turned out to be the novel genes associated with the carcinogenesis of EOC. Our study demonstrated that methylation profiling through MeDIP-sequencing has effectively identified six potential hub genes and pathways that might exacerbate our understanding of underlying molecular mechanisms of ovarian carcinogenesis.

Ovarian cancer, among all gynecologic malignancies, exhibits the highest incidence and mortality rate, primarily because it is often presents with non-specific or no symptoms during its early stages. For the advancement of Ovarian Cancer Diagnosis, it is crucial to identify the potential molecular signatures that could significantly differentiate between healthy and ovarian cancerous tissues and can be used further as a diagnostic biomarker for detecting ovarian cancer. In this study, we investigated the genome-wide methylation patterns in ovarian cancer patients using Methylated DNA Immunoprecipitation (MeDIP-Seq) followed by NGS. Identified differentially methylated regions (DMRs) were further validated by targeted bisulfite sequencing for CpG site-specific methylation profiles. Furthermore, expression validation of six genes by Quantitative Reverse Transcriptase-PCR was also performed. Out of total 120 differentially methylated genes (DMGs), 68 genes were hypermethylated, and 52 were hypomethylated in their promoter region. After analysis, we identified the top 6 hub genes, namely POLR3B, PLXND1, GIGYF2, STK4, BMP2 and CRKL. Interestingly we observed Non-CpG site methylation in the case of POLR3B and CRKL which was statistically significant in discriminating ovarian cancer samples from normal controls. The most significant pathways identified were focal adhesion, the MAPK signaling pathway, and the Ras signaling pathway. Expression analysis of hypermethylated genes was correlated with the downregulation of the genes. POLR3B and GIGYF2 turned out to be the novel genes associated with the carcinogenesis of EOC. Our study demonstrated that methylation profiling through MeDIP-sequencing has effectively identified six potential hub genes and pathways that might exacerbate our understanding of underlying molecular mechanisms of ovarian carcinogenesis. [ABSTRACT FROM AUTHOR]

Heavy metal toxicity is a significant global health concern, with particular attention given to lead (Pb) exposure due to its adverse effects on cognitive development, especially in children exposed to low concentrations. While Pb neurotoxicity has been extensively studied, the analysis and molecular mechanisms underlying the transgenerational effects of Pb exposure-induced neurotoxicity remain poorly understood. In this study, we utilized Drosophila, a powerful developmental animal model, to investigate this phenomenon. Our findings demonstrated that Pb exposure during the developmental stage had a profound effect on the neurodevelopment of F0 fruit flies. Specifically, we observed a loss of correlation between the terminal motor area and muscle fiber area, along with an increased frequency of the β-lobe midline crossing phenotype in mushroom bodies. Western blot analysis indicated altered expression levels of synaptic vesicle proteins, with a decrease in Synapsin (SYN) and an increase in Bruchpilot (BRP) expression, suggesting changes in synaptic vesicle release sites. These findings were corroborated by electrophysiological data, showing an increase in the amplitude of evoked excitatory junctional potential (EJP) and an increase in the frequency of spontaneous excitatory junctional potential (mEJP) following Pb exposure. Importantly, our results further confirmed that the developmental neurotoxicity resulting from grandparental Pb exposure exhibited a transgenerational effect. The F3 offspring displayed neurodevelopmental defects, synaptic function abnormalities, and repetitive behavior despite lacking direct Pb exposure. Our MeDIP-seq analysis further revealed significant alterations in DNA methylation levels in several neurodevelopmental associated genes (eagle, happyhour, neuroglian, bazooka, and spinophilin) in the F3 offspring exposed to Pb. These findings suggest that DNA methylation modifications may underlie the inheritance of acquired phenotypic traits resulting from environmental Pb exposure.

Background: Tardive dyskinesia (TD) develops in 20–30% of schizophrenia patients and up to 50% in patients > 50 years old. DNA methylation may play an important role in the development of TD. Objective: DNA methylation analyses in schizophrenia with TD. Methods: We conducted a genome-wide DNA methylation analysis in schizophrenia with TD using methylated DNA immunoprecipitation coupled with next-generation sequencing (MeDIP-Seq) in a Chinese sample including five schizophrenia patients with TD and five without TD (NTD), and five healthy controls. The results were expressed as the log2FC, fold change of normalized tags between two groups within the differentially methylated region (DMR). For validation, the pyrosequencing was used to quantify DNA methylation levels of several methylated genes in an independent sample (n = 30). Results: Through genome-wide MeDIP-Seq analysis, we identified 116 genes that were significantly differentially methylated in promotor regions in comparison of TD group with NTD group including 66 hypermethylated genes (top 4 genes are GABRR1, VANGL2, ZNF534, and ZNF746) and 50 hypomethylated genes (top 4 genes are DERL3, GSTA4, KNCN, and LRRK1). Part of these genes (such as DERL3, DLGAP2, GABRR1, KLRG2, LRRK1, VANGL2, and ZP3) were previously reported to be associated with methylation in schizophrenia. Gene Ontology enrichment and KEGG pathway analyses identified several pathways. So far, we have confirmed the methylation of 3 genes (ARMC6, WDR75, and ZP3) in schizophrenia with TD using pyrosequencing. Conclusions: This study identified number of methylated genes and pathways for TD and will provide potential biomarkers for TD and serve as a resource for replication in other populations. [ABSTRACT FROM AUTHOR]

Brown algae have convergently evolved plant-like body plans and reproductive cycles, which in plants are controlled by differential DNA methylation. This contribution provides the first single-base methylome profiles of haploid gametophytes and diploid sporophytes of a multicellular alga. Although only c. 1.4% of cytosines in Saccharina japonica were methylated mainly at CHH sites and characterized by 5-methylcytosine (5mC), there were significant differences between life-cycle stages. DNA methyltransferase 2 (DNMT2), known to efficiently catalyze tRNA methylation, is assumed to methylate the genome of S. japonica in the structural context of tRNAs as the genome does not encode any other DNA methyltransferases. Circular and long noncoding RNA genes were the most strongly methylated regulatory elements in S. japonica. Differential expression of genes was negatively correlated with DNA methylation with the highest methylation levels measured in both haploid gametophytes. Hypomethylated and highly expressed genes in diploid sporophytes included genes involved in morphogenesis and halogen metabolism. The data herein provide evidence that cytosine methylation, although occurring at a low level, is significantly contributing to the formation of different life-cycle stages, tissue differentiation and metabolism in brown algae.

Autologous mesenchymal stem/stromal cells (MSCs) have demonstrated important therapeutic effects in several diseases. Cardiovascular risk factors may impair MSC mitochondrial structure and function, but the underlying mechanisms remain unknown. We hypothesized that metabolic syndrome (MetS) induces epigenetic alterations in mitochondria-related genes in swine MSCs. Pigs were fed a Lean or MetS diet (n = 6 each) for 16 weeks. MSCs were collected from subcutaneous abdominal fat, and DNA hydroxymethylation (5 hmC) profiles of mitochondria-related genes (MitoCarta-2.0) were analyzed by hydroxymethylated DNA immunoprecipitation and next-generation sequencing (hMeDIP-seq) in Lean- and MetS-MSCs untreated or treated with the epigenetic modulator vitamin (Vit)-C (n = 3 each). Functional analysis of genes with differential 5 hmC regions was performed using DAVID6.8. Mitochondrial structure (electron microscopy), oxidative stress, and membrane potential were assessed. hMeDIP-seq identified 172 peaks (associated with 103 mitochondrial genes) with higher and 416 peaks (associated with 165 mitochondrial genes) with lower 5 hmC levels in MetS-MSCs versus Lean-MSCs (≥2-fold, p < 0.05). Genes with higher 5 hmC levels in MetS + MSCs were primarily implicated in fatty acid metabolism, whereas those with lower 5 hmC levels were associated with electron transport chain activity. Vit-C increased 5 hmC levels in mitochondrial antioxidant genes, improved mitochondrial structure and membrane potential, and decreased oxidative stress. MetS alters 5 hmC levels of mitochondria-related genes in swine MSCs. Vit-C modulated 5 hmC levels in these genes and preserved mitochondrial structure and function in MetS-MSCs. These observations may contribute to development of strategies to overcome the deleterious effects of MetS on MSCs. [ABSTRACT FROM AUTHOR]

Background: Diabetic cardiomyopathy (DCM) is a leading cause of death in diabetic patients. Hyperglycemic myocardial microenvironment significantly alters chromatin architecture and the transcriptome, resulting in aberrant activation of signaling pathways in a diabetic heart. Epigenetic marks play vital roles in transcriptional reprogramming during the development of DCM. The current study is aimed to profile genome-wide DNA (hydroxy)methylation patterns in the hearts of control and streptozotocin (STZ)-induced diabetic rats and decipher the effect of modulation of DNA methylation by alpha-ketoglutarate (AKG), a TET enzyme cofactor, on the progression of DCM. Methods: Diabetes was induced in male adult Wistar rats with an intraperitoneal injection of STZ. Diabetic and vehicle control animals were randomly divided into groups with/without AKG treatment. Cardiac function was monitored by performing cardiac catheterization. Global methylation (5mC) and hydroxymethylation (5hmC) patterns were mapped in the Left ventricular tissue of control and diabetic rats with the help of an enrichment-based (h)MEDIP-sequencing technique by using antibodies specific for 5mC and 5hmC. Sequencing data were validated by performing (h)MEDIP-qPCR analysis at the gene-specific level, and gene expression was analyzed by qPCR. The mRNA and protein expression of enzymes involved in the DNA methylation and demethylation cycle were analyzed by qPCR and western blotting. Global 5mC and 5hmC levels were also assessed in high glucose-treated DNMT3B knockdown H9c2 cells. Results: We found the increased expression of DNMT3B, MBD2, and MeCP2 with a concomitant accumulation of 5mC and 5hmC, specifically in gene body regions of diabetic rat hearts compared to the control. Calcium signaling was the most significantly affected pathway by cytosine modifications in the diabetic heart. Additionally, hypermethylated gene body regions were associated with Rap1, apelin, and phosphatidyl inositol signaling, while metabolic pathways were most affected by hyperhydroxymethylation. AKG supplementation in diabetic rats reversed aberrant methylation patterns and restored cardiac function. Hyperglycemia also increased 5mC and 5hmC levels in H9c2 cells, which was normalized by DNMT3B knockdown or AKG supplementation. Conclusion: This study demonstrates that reverting hyperglycemic damage to cardiac tissue might be possible by erasing adverse epigenetic signatures by supplementing epigenetic modulators such as AKG along with an existing antidiabetic treatment regimen. [ABSTRACT FROM AUTHOR]

Entamoeba histolytica , an intestinal parasite of global significance, poses substantial health risks with its associated high morbidity and mortality rates. Despite the current repertoire of molecular tools for the study of gene function in, the regulatory mechanisms governing its pathogenicity remain largely unexplored. This knowledge gap underscores the need to elucidate key genetic determinants orchestrating cellular functions critical to its virulence. Previously, our group generated an avirulent strain, termed UG10, with the same genetic background as the HM1:IMSS strain. UG10 strain, despite showing normal expression levels of well-known virulence factors, was unable to perform in-vitro and in-vivo activities related to amoebic virulence. In this study, we aimed to uncover the genome-wide modifications that rendered the avirulent phenotype of the UG10 strain through whole-genome sequencing. As a complementary approach, we conducted Methylated DNA Immunoprecipitation coupled with sequencing (MeDIP-seq) analysis on both the highly virulent HM1:IMSS strain and the low-virulence UG10 strain to uncover the genome-wide methylation profile. These dual methodologies revealed two aspects of the UG10 avirulent strain. One is the random integration of fragments from the ribosomal gene cluster and tRNA genes, ranging from 120 to 400 bp; and secondly, a clear, enriched methylation profile in the coding and non-coding strand relative to the start codon sequence in genes encoding small GTPases, which is associated with the previously described avirulent phenotype. This study provides the foundation to explore other genetic and epigenetic regulatory circuitries in E. histolytica and novel targets to understand the pathogenic mechanism of this parasite. • Entamoeba histolytica has a genomic structure highly variable. • Using strains with different virulence phenotype can uncover virulence mechanisms. • UG10 strain exhibits low virulence phenotype in comparison with HM1:IMSS strain. • Whole genome sequence of UG10 strain revealed discreate modifications. • MeDIP-seq analysis revealed a differentially methylated upstream regions. [ABSTRACT FROM AUTHOR]

Abstract Background Colorectal cancer is the most common malignancy and the third leading cause of cancer-related death worldwide. This study aimed to identify potential diagnostic biomarkers for colorectal cancer by genome-wide plasma cell-free DNA (cfDNA) methylation analysis. Methods Peripheral blood from colorectal cancer patients and healthy controls was collected for cfDNA extraction. Genome-wide cfDNA methylation profiling, especially differential methylation profiling between colorectal cancer patients and healthy controls, was performed by methylated DNA immunoprecipitation coupled with high-throughput sequencing (MeDIP-seq). Logistic regression models were established, and the accuracy of this diagnostic model for colorectal cancer was verified using tissue-sourced data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) due to the lack of cfDNA methylation data in public datasets. Results Compared with the control group, 939 differentially methylated regions (DMRs) located in promoter regions were found in colorectal cancer patients; 16 of these DMRs were hypermethylated, and the remaining 923 were hypomethylated. In addition, these hypermethylated genes, mainly PRDM14, RALYL, ELMOD1, and TMEM132E, were validated and confirmed in colorectal cancer by using publicly available DNA methylation data. Conclusions MeDIP-seq can be used as an optimal approach for analyzing cfDNA methylomes, and 12 probes of four differentially methylated genes identified by MeDIP-seq (PRDM14, RALYL, ELMOD1, and TMEM132E) could serve as potential biomarkers for clinical application in patients with colorectal cancer.

Autologous mesenchymal stem/stromal cells (MSCs) have demonstrated important therapeutic effects in several diseases. Cardiovascular risk factors may impair MSC mitochondrial structure and function, but the underlying mechanisms remain unknown. We hypothesized that metabolic syndrome (MetS) induces epigenetic alterations in mitochondria-related genes in swine MSCs. Pigs were fed a Lean or MetS diet (n = 6 each) for 16 weeks. MSCs were collected from subcutaneous abdominal fat, and DNA hydroxymethylation (5 hmC) profiles of mitochondria-related genes (MitoCarta-2.0) were analyzed by hydroxymethylated DNA immunoprecipitation and next-generation sequencing (hMeDIP-seq) in Lean- and MetS-MSCs untreated or treated with the epigenetic modulator vitamin (Vit)-C (n = 3 each). Functional analysis of genes with differential 5 hmC regions was performed using DAVID6.8. Mitochondrial structure (electron microscopy), oxidative stress, and membrane potential were assessed. hMeDIP-seq identified 172 peaks (associated with 103 mitochondrial genes) with higher and 416 peaks (associated with 165 mitochondrial genes) with lower 5 hmC levels in MetS-MSCs versus Lean-MSCs (≥2-fold, p < 0.05). Genes with higher 5 hmC levels in MetS + MSCs were primarily implicated in fatty acid metabolism, whereas those with lower 5 hmC levels were associated with electron transport chain activity. Vit-C increased 5 hmC levels in mitochondrial antioxidant genes, improved mitochondrial structure and membrane potential, and decreased oxidative stress. MetS alters 5 hmC levels of mitochondria-related genes in swine MSCs. Vit-C modulated 5 hmC levels in these genes and preserved mitochondrial structure and function in MetS-MSCs. These observations may contribute to development of strategies to overcome the deleterious effects of MetS on MSCs.

DNA adenine N6-methylation (6mA) plays a critical role in various biological functions, but its occurrence and functions in filamentous plant pathogens are largely unexplored. Botrytis cinerea is an important pathogenic fungus worldwide. A systematic analysis of 6mA in B. cinerea was performed in this study, revealing that 6mA is widely distributed in the genome of this fungus. The 2 kb regions flanking many genes, particularly the upstream promoter regions, were susceptible to methylation. The role of BcMettl4, a 6mA methyltransferase, in the virulence of B. cinerea was investigated. BcMETTL4 disruption and point mutations of its catalytic motif "DPPW" both resulted in significant 6mA reduction in the genomic DNA and in reduced virulence of B. cinerea. RNA-Seq analysis revealed a total of 13 downregulated genes in the disruption mutant ΔBcMettl4 in which methylation occurred at the promoter sites. These were involved in oxidoreduction, secretory pathways, autophagy and carbohydrate metabolism. Two of these genes, BcFDH and BcMFS2 , were independently disrupted. Knockout of BcFDH led to reduced sclerotium formation, while disruption of BcMFS2 resulted in dramatically decreased conidium formation and pathogenicity. These observations indicated that 6mA provides potential epigenetic markers in B. cinerea and that BcMettl4 regulates virulence in this important plant pathogen. [ABSTRACT FROM AUTHOR]

Myocardial infarction (MI) is one of the deadliest diseases in the world, and the changes at the molecular level after MI and the DNA methylation features are not clear. Understanding the molecular characteristics of the early stages of MI is of significance for the treatment of the disease. In this study, RNA-seq and MeDIP-seq were performed on heart tissue from mouse models at multiple time points (0 h, 10 min, 1, 6, 24, and 72 h) to explore genetic and epigenetic features that influence MI progression. Analysis based on a single point in time, the number of differentially expressed genes (DEGs) and differentially methylated regions (DMRs) increased with the time of myocardial infarction, using 0 h as a control group. Moreover, within 10 min of MI onset, the cells are mainly in immune response, and as the duration of MI increases, apoptosis begins to occur. Analysis based on time series data, the expression of 1012 genes was specifically downregulated, and these genes were associated with energy metabolism. The expression of 5806 genes was specifically upregulated, and these genes were associated with immune regulation, inflammation and apoptosis. Fourteen transcription factors were identified in the genes involved in apoptosis and inflammation, which may be potential drug targets. Analysis based on MeDIP-seq combined with RNA-seq methodology, focused on methylation at the promoter region. GO revealed that the downregulated genes with hypermethylation at 72 h were enriched in biological processes such as cardiac muscle contraction. In addition, the upregulated genes with hypomethylation at 72 h were enriched in biological processes, such as cell-cell adhesion, regulation of the apoptotic signaling pathway and regulation of angiogenesis. Among these genes, the Tnni3 gene was also present in the downregulated model. Hypermethylation of Tnni3 at 72 h after MI may be an important cause of exacerbation of MI.

DNA adenine N6-methylation (6mA) plays a critical role in various biological functions, but its occurrence and functions in filamentous plant pathogens are largely unexplored. Botrytis cinerea is an important pathogenic fungus worldwide. A systematic analysis of 6mA in B. cinerea was performed in this study, revealing that 6mA is widely distributed in the genome of this fungus. The 2 kb regions flanking many genes, particularly the upstream promoter regions, were susceptible to methylation. The role of BcMettl4, a 6mA methyltransferase, in the virulence of B. cinerea was investigated. BcMETTL4 disruption and point mutations of its catalytic motif “DPPW” both resulted in significant 6mA reduction in the genomic DNA and in reduced virulence of B. cinerea. RNA-Seq analysis revealed a total of 13 downregulated genes in the disruption mutant ΔBcMettl4 in which methylation occurred at the promoter sites. These were involved in oxidoreduction, secretory pathways, autophagy and carbohydrate metabolism. Two of these genes, BcFDH and BcMFS2, were independently disrupted. Knockout of BcFDH led to reduced sclerotium formation, while disruption of BcMFS2 resulted in dramatically decreased conidium formation and pathogenicity. These observations indicated that 6mA provides potential epigenetic markers in B. cinerea and that BcMettl4 regulates virulence in this important plant pathogen.

Backgrounds: Acute myocardial infarction (AMI) has been one of the most fatal diseases among all types of heart diseases due to its rapid onset and high rates of fatality. Understanding accurately how multi-omics molecular features change at the early stage of AMI is crucial for its treatment. Currently, the changes involved in DNA methylation modification and gene expression of multiple genes have remained unexplored. Results: We used the RNA-seq and MeDIP-seq on heart tissues from AMI mouse models at series of time points (Sham, AMI 10-min, 1-h, 6-h, 24-h and 72-h), to comprehensively describe the transcriptome and genome-wide DNA methylation changes at above time points. We identified 18814, 18614, 23587, 26018 and 33788 differential methylation positions (DMPs) and 123, 135, 731, 1419 and 2779 differentially expressed genes (DEGs) at 10-min, 1-h, 6-h, 24-h and 72-h AMI, respectively, compared with the sham group. Remarkably, the 6-h AMI with the drastic changes of DEGs and a large number of enriched functional pathways in KEGG may be the most critical stage of AMI process. The 4, 9, 40, 26, and 183 genes were further identified at each time point, based on the negative correlation (P < 0.05) between the differential mRNA expression and the differential DNA methylation. The mRNA and the promoter methylation expressions of five genes (Ptpn6, Csf1r, Col6a1, Cyba, and Map3k14) were validated by qRT-PCR and BSP methods, and the mRNA expressions were further confirmed to be regulated by DNA methylation in cardiomyocytes in vitro. Conclusions: Our findings profiled the molecular variations from the perspective of DNA methylation in the early stage of AMI and provided promising epigenetic-based biomarkers for the early clinical diagnosis and therapeutic targets of AMI. [ABSTRACT FROM AUTHOR]

Background: Colorectal cancer is the most common malignancy and the third leading cause of cancer-related death worldwide. This study aimed to identify potential diagnostic biomarkers for colorectal cancer by genome-wide plasma cell-free DNA (cfDNA) methylation analysis. Methods: Peripheral blood from colorectal cancer patients and healthy controls was collected for cfDNA extraction. Genome-wide cfDNA methylation profiling, especially differential methylation profiling between colorectal cancer patients and healthy controls, was performed by methylated DNA immunoprecipitation coupled with high-throughput sequencing (MeDIP-seq). Logistic regression models were established, and the accuracy of this diagnostic model for colorectal cancer was verified using tissue-sourced data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) due to the lack of cfDNA methylation data in public datasets. Results: Compared with the control group, 939 differentially methylated regions (DMRs) located in promoter regions were found in colorectal cancer patients; 16 of these DMRs were hypermethylated, and the remaining 923 were hypomethylated. In addition, these hypermethylated genes, mainly PRDM14, RALYL, ELMOD1, and TMEM132E, were validated and confirmed in colorectal cancer by using publicly available DNA methylation data. Conclusions: MeDIP-seq can be used as an optimal approach for analyzing cfDNA methylomes, and 12 probes of four differentially methylated genes identified by MeDIP-seq (PRDM14, RALYL, ELMOD1, and TMEM132E) could serve as potential biomarkers for clinical application in patients with colorectal cancer. [ABSTRACT FROM AUTHOR]

Hair follicle growth and development are a complex and long-term physiological process, which is regulated by a variety of physical factors and signal pathways. Increasing the understanding of the epigenetic regulation and function of candidate genes related to hair follicle development will help to better understand the molecular regulatory mechanisms of hair follicle development. In this study, the methylated DNA immunoprecipitation sequencing (MeDIP-seq) was used to obtain the genome-wide methylation map of the hair follicular development of Super Merino sheep in six stages (fetal skin tissue at 65d, 85d, 105d, 135d, 7d, and 30d after birth). Combined with the results of previous RNA-sequencing, 65 genes were screened out that were both differential methylation and differential expression, including EDN1, LAMC2, NR1D1, RORB, MyOZ3, and WNT2 gene. Differential methylation genes were enriched in Wnt, TNF, TGF-beta, and other signaling pathways related to hair follicle development. The bisulfite sequencing PCR results and MeDIP-seq were basically consistent, indicating that the sequencing results were accurate. As a key gene in the Wnt signaling pathway, both differential methylation and expression gene identified by MeDIP-seq and RNA-seq, further exploration of the function of WNT2 gene revealed that the DNA methylation of exon 5 (CpG11 site) promoted the expression of WNT2 gene. The overexpression vector of lentivirus pLEX-MCS-WNT2 was constructed, and WNT2 gene effectively promoted the proliferation of sheep skin fibroblasts. The results showed that WNT2 gene could promote the growth and development of skin and hair follicles. The results of this study will provide a theoretical basis for further research on sheep hair follicle development and gene regulation mechanisms.

Hair follicle growth and development are a complex and long-term physiological process, which is regulated by a variety of physical factors and signal pathways. Increasing the understanding of the epigenetic regulation and function of candidate genes related to hair follicle development will help to better understand the molecular regulatory mechanisms of hair follicle development. In this study, the methylated DNA immunoprecipitation sequencing (MeDIP-seq) was used to obtain the genome-wide methylation map of the hair follicular development of Super Merino sheep in six stages (fetal skin tissue at 65d, 85d, 105d, 135d, 7d, and 30d after birth). Combined with the results of previous RNA-sequencing, 65 genes were screened out that were both differential methylation and differential expression, including EDN1, LAMC2, NR1D1, RORB, MyOZ3, and WNT2 gene. Differential methylation genes were enriched in Wnt, TNF, TGF-beta, and other signaling pathways related to hair follicle development. The bisulfite sequencing PCR results and MeDIP-seq were basically consistent, indicating that the sequencing results were accurate. As a key gene in the Wnt signaling pathway, both differential methylation and expression gene identified by MeDIP-seq and RNA-seq, further exploration of the function of WNT2 gene revealed that the DNA methylation of exon 5 (CpG11 site) promoted the expression of WNT2 gene. The overexpression vector of lentivirus pLEX-MCS-WNT2 was constructed, and WNT2 gene effectively promoted the proliferation of sheep skin fibroblasts. The results showed that WNT2 gene could promote the growth and development of skin and hair follicles. The results of this study will provide a theoretical basis for further research on sheep hair follicle development and gene regulation mechanisms. [ABSTRACT FROM AUTHOR]

Silica nanoparticles (SiNPs), which are the main inorganic components of atmospheric particulate matter, have been proved to have certain male reproductive toxicity in previous studies. Spermatogenesis involves complex epigenetic regulation, but it is still unclear if SiNPs exposure will interfere with the DNA methylation patterns in mouse spermatocytes. The present study was designed to investigate the effects of SiNPs on DNA methylation in the mouse spermatocyte GC-2spd(ts). GC-2 cells were treated with 0 and 20 μg/mL SiNPs for 24 h. MeDIP-seq assay was then performed to analyze the differentially methylated genes related to spermatogenesis. The results showed that SiNPs induced extensive methylation changes in the genome of GC-2 cells, and 24a total of 428 hyper-methylated genes and 398 hypo-methylated genes were identified. Gene Ontology and pathway analysis showed that differential DNA methylation induced by SiNPs was probably involved with abnormal transcription and translation, mitochondrial damage, and cell apoptosis. Results from qRT-PCR verification showed that the expression of spermatogenesis-related genes Akap1, Crem, Spz1, and Tex11 were dysregulated by SiNPs exposure, which was consistent with the MeDIP-seq assay. In general, this study suggested that SiNPs caused genome-wide DNA methylation changes in GC-2 cells, providing valuable reference for the future epigenetic studies in SiNPs-induced male reproductive toxicity. [ABSTRACT FROM AUTHOR]

Abstract Background Tibetan pigs, which inhabit the Tibetan Plateau, exhibit distinct phenotypic and physiological characteristics from those of lowland pigs and have adapted well to the extreme conditions at high altitude. However, the genetic and epigenetic mechanisms of hypoxic adaptation in animals remain unclear. Methods Whole-genome DNA methylation data were generated for heart tissues of Tibetan pigs grown in the highland (TH, n = 4) and lowland (TL, n = 4), as well as Yorkshire pigs grown in the highland (YH, n = 4) and lowland (YL, n = 4), using methylated DNA immunoprecipitation sequencing. Results We obtained 480 million reads and detected 280679, 287224, 259066, and 332078 methylation enrichment peaks in TH, YH, TL, and YL, respectively. Pairwise TH vs. YH, TL vs. YL, TH vs. TL, and YH vs. YL comparisons revealed 6829, 11997, 2828, and 1286 differentially methylated regions (DMRs), respectively. These DMRs contained 384, 619, 192, and 92 differentially methylated genes (DMGs), respectively. DMGs that were enriched in the hypoxia-inducible factor 1 signaling pathway and pathways involved in cancer and hypoxia-related processes were considered to be important candidate genes for high-altitude adaptation in Tibetan pigs. Conclusions This study elucidates the molecular and epigenetic mechanisms involved in hypoxic adaptation in pigs and may help further understand human hypoxia-related diseases.

Background: Avian leukosis virus subgroup J (ALV-J) is an oncogenic virus that causes serious economic losses in the poultry industry; unfortunately, there is no effective vaccine against ALV-J. DNA methylation plays a crucial role in several biological processes, and an increasing number of diseases have been proven to be related to alterations in DNA methylation. In this study, we screened ALV-J-positive and -negative chickens. Subsequently, we generated and provided the genome-wide gene expression and DNA methylation profiles by MeDIP-seq and RNA-seq of ALV-J-positive and -negative chicken samples; 8,304 differentially methylated regions (DMRs) were identified by MeDIP-seq analysis (p ≤ 0.005) and 515 differentially expressed genes were identified by RNA-seq analysis (p ≤ 0.05). As a result of an integration analysis, we screened six candidate genes to identify ALV-J-negative chickens that possessed differential methylation in the promoter region. Furthermore, TGFB2 played an important role in tumorigenesis and cancer progression, which suggested TGFB2 may be an indicator for identifying ALV-J infections.

Background: Epirubicin (EPI) is an important anticancer drug that is well-known for its cardiotoxic side effect. Studying epigenetic modification such as DNA methylation can help to understand the EPI-related toxic mechanisms in cardiac tissue. In this study, we analyzed the DNA methylation profile in a relevant human cell model and inspected the expression of differentially methylated genes at the transcriptome level to understand how changes in DNA methylation could affect gene expression in relation to EPI-induced cardiotoxicity. Methods: Human cardiac microtissues were exposed to either therapeutic or toxic (IC20) EPI doses during 2 weeks. The DNA and RNA were collected from microtissues in triplicates at 2, 8, 24, 72, 168, 240, and 336 hours of exposure. Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) analysis was used to detect DNA methylation levels in EPI-treated and control samples. The MeDIP-seq data were analyzed and processed using the QSEA package with a recently published workflow. RNA sequencing (RNA-seq) was used to measure global gene expression in the same samples. Results: After processing the MeDIP-seq data, we detected 35, 37, 15 candidate genes which show strong methylated alterations between all EPI-treated, EPI therapeutic and EPI toxic dose-treated samples compared to control, respectively. For several genes, gene expressions changed compatibly reflecting the DNA methylation regulation. Conclusions: The observed DNA methylation modifications provide further insights into the EPI-induced cardiotoxicity. Multiple differentially methylated genes under EPI treatment, such as SMARCA4, PKN1, RGS12, DPP9, NCOR2, SDHA, POLR2A, and AGPAT3, have been implicated in different cardiac dysfunction mechanisms. Together with other differentially methylated genes, these genes can be candidates for further investigations of EPI-related toxic mechanisms. Data Repository: The data has been generated by the HeCaToS project (http://www.ebi.ac.uk/biostudies) under accession numbers S-HECA433 and S-HECA434 for the MeDIP-seq data and S-HECA11 for the RNA-seq data. The R code is available on Github (https://github.com/NhanNguyen000/MeDIP).

As one of the most malicious cancers, pancreatic cancer is difficult to treat due to the lack of effective early diagnosis. Therefore, it is urgent to find reliable diagnostic and predictive markers for the early detection of pancreatic cancer. In recent years, the detection of circulating cell-free DNA (cfDNA) methylation in plasma has attracted global attention for non-invasive and early cancer diagnosis. Here, we carried out a genome-wide cfDNA methylation profiling study of pancreatic ductal adenocarcinoma (PDAC) patients by methylated DNA immunoprecipitation coupled with high-throughput sequencing (MeDIP-seq). Compared with healthy individuals, 775 differentially methylated regions (DMRs) located in promoter regions were identified in PDAC patients with 761 hypermethylated and 14 hypomethylated regions; meanwhile, 761 DMRs in CpG islands (CGIs) were identified in PDAC patients with 734 hypermethylated and 27 hypomethylated regions (p-value < 0.0001). Then, 143 hypermethylated DMRs were further selected which were located in promoter regions and completely overlapped with CGIs. After performing the least absolute shrinkage and selection operator (LASSO) method, a total of eight markers were found to fairly distinguish PDAC patients from healthy individuals, including TRIM73, FAM150A, EPB41L3, SIX3, MIR663, MAPT, LOC100128977, and LOC100130148. In conclusion, this work identified a set of eight differentially methylated markers that may be potentially applied in non-invasive diagnosis of pancreatic cancer.

DNA N6-methyladenine (6mA) modification has been discovered as the most prevalent DNA modification in prokaryotes and eukaryotes, involving gene expression, DNA replication and repair, and host-pathogen interactions. Single-molecule real-time sequencing (SMRT-seq) can detect 6mA events in prokaryotic and eukaryotic genomes at the single-nucleotide level. However, there are no strict and economical quality control methods for high false-positive 6mA events in eukaryotic genomes. Therefore, by analyzing the distribution of 6mA in eukaryotic and prokaryotes, we proposed a method named MASQC (MeDIP-seq assists SMRT-seq for quality control in 6mA identification), which can identify 6mA events without doing the whole genome amplification (WGA) sequencing. The proposed MASQC method was assessed on two eukaryotic genomes and six bacterial genomes, our results demonstrate that MASQC performs well in quality control of false positive 6mA identification for both eukaryotic and prokaryotic genomes.

Abstract Background DNA methylation of CpG dinucleotides is an essential epigenetic modification that plays a key role in transcription. Widely used DNA enrichment-based methods offer high coverage for measuring methylated CpG dinucleotides, with the lowest cost per CpG covered genome-wide. However, these methods measure the DNA enrichment of methyl-CpG binding, and thus do not provide information on absolute methylation levels. Further, the enrichment is influenced by various confounding factors in addition to methylation status, for example, CpG density. Computational models that can accurately derive absolute methylation levels from DNA enrichment data are needed. Results We developed “MeDEStrand,” a method that uses a sigmoid function to estimate and correct the CpG bias from enrichment results to infer absolute DNA methylation levels. Unlike previous methods, which estimate CpG bias based on reads mapped at the same genomic loci, MeDEStrand processes the reads for the positive and negative DNA strands separately. We compared the performance of MeDEStrand to that of three other state-of-the-art methods “MEDIPS,” “BayMeth,” and “QSEA” on four independent datasets generated using immortalized cell lines (GM12878 and K562) and human primary cells (foreskin fibroblasts and mammary epithelial cells). Based on the comparison of the inferred absolute methylation levels from MeDIP-seq data and the corresponding reduced-representation bisulfite sequencing data from each method, MeDEStrand showed the best performance at high resolution of 25, 50, and 100 base pairs. Conclusions The MeDEStrand tool can be used to infer whole-genome absolute DNA methylation levels at the same cost of enrichment-based methods with adequate accuracy and resolution. R package MeDEStrand and its tutorial is freely available for download at https://github.com/jxu1234/MeDEStrand.git.

Abstract Background Milk yield for Holstein cows has doubled over five decades due to genetic selection and changes to management, but the molecular mechanisms that facilitated this increase are mostly unknown. Epigenetic modifications to the cattle genome are a plausible molecular mechanism to cause variation in milk yield and our objective was to describe genome-wide DNA methylation patterns in peripheral blood mononuclear cells (PBMC) from mature Holstein dairy cows with variable milk yield. Results Whole genome MeDIP-seq was performed following DNA extraction from PBMC of 6 lactating dairy cows from 4 different herds that varied in milk yield from 13,556 kg to 23,105 kg per 305 day lactation. We describe methylation across the genome and for 13,677 protein coding genes. Repetitive element reads were primarily mapped to satellite (36.4%), SINE (29.1%), and LINE (23.7%) regions and the majority (78.4%) of CpG sites were sequenced at least once. DNA methylation was generally low upstream of genes with the nadir occurring 95 bp prior to the transcription start site (TSS). Methylation was lower in the first exon than in later exons, was highest for introns near the intron-exon junctions, and declined downstream as the distance from the gene increased. We identified 72 differentially methylated regions (DMR) between high milk yield cows and their control, and 252 DMR across herd environments. Conclusions This reference methylome for cattle with extreme variation in milk yield phenotype provides a resource to more fully evaluate relationships between DNA methylation and phenotype in populations subject to selection. The detection of DMR in cows of varying milk yield suggests potential to exploit epigenetic variation in cattle improvement programs.

Abstract Purpose We investigated the genomic DNA methylation profile of prostate cancer in transgenic adenocarcinoma of the mouse prostate (TRAMP) cancer model and to analyze the crosstalk among targeted genes and the related functional pathways. Methods Prostate DNA samples from 24-week-old TRAMP and C57BL/6 male mice were isolated. The DNA methylation profiles were analyzed by methylated DNA immunoprecipitation (MeDIP) followed by next-generation sequencing (MeDIP-seq). Canonical pathways, diseases and function and network analyses of the different samples were then performed using the Ingenuity® Pathway Analysis (IPA) software. Some target genes with significant difference in methylation were selected for validation using methylation specific primers (MSP) and qPCR. Results TRAMP mice undergo extensive aberrant CpG hyper- and hypo-methylation. There were 2147 genes with a significant (log2-change ≥ 2) change in CpG methylation between the two groups, as mapped by the IPA software. Among these genes, the methylation of 1105 and 1042 genes was significantly decreased and increased, respectively, in TRAMP prostate tumors. The top associated disease identified by IPA was adenocarcinoma; however, the cAMP response element-binding protein (CREB)-, histone deacetylase 2 (HDAC2)-, glutathione S-transferase pi (GSTP1)- and polyubiquitin-C (UBC)-related pathways showed significantly altered methylation profiles based on the canonical pathway and network analyses. MSP and qPCR results of genes of interests corroborated with MeDIP-seq findings. Conclusions This is the first MeDIP-seq with IPA analysis of the TRAMP model to provide novel insight into the genome-wide methylation profile of prostate cancer. Studies on epigenetics, such as DNA methylation, will potentially provide novel avenues and strategies for further development of biomarkers targeted for treatment and prevention approaches for prostate cancer.

As one of the most malicious cancers, pancreatic cancer is difficult to treat due to the lack of effective early diagnosis. Therefore, it is urgent to find reliable diagnostic and predictive markers for the early detection of pancreatic cancer. In recent years, the detection of circulating cell-free DNA (cfDNA) methylation in plasma has attracted global attention for non-invasive and early cancer diagnosis. Here, we carried out a genome-wide cfDNA methylation profiling study of pancreatic ductal adenocarcinoma (PDAC) patients by methylated DNA immunoprecipitation coupled with high-throughput sequencing (MeDIP-seq). Compared with healthy individuals, 775 differentially methylated regions (DMRs) located in promoter regions were identified in PDAC patients with 761 hypermethylated and 14 hypomethylated regions; meanwhile, 761 DMRs in CpG islands (CGIs) were identified in PDAC patients with 734 hypermethylated and 27 hypomethylated regions (p -value < 0.0001). Then, 143 hypermethylated DMRs were further selected which were located in promoter regions and completely overlapped with CGIs. After performing the least absolute shrinkage and selection operator (LASSO) method, a total of eight markers were found to fairly distinguish PDAC patients from healthy individuals, including TRIM73 , FAM150A , EPB41L3 , SIX3 , MIR663 , MAPT , LOC100128977 , and LOC100130148. In conclusion, this work identified a set of eight differentially methylated markers that may be potentially applied in non-invasive diagnosis of pancreatic cancer. [ABSTRACT FROM AUTHOR]

DNA methylation plays an important role in the regulation of some biological processes. Up to now, with the development of machine learning models, there are several sequence-based deep learning models designed to predict DNA methylation states, which gain better performance than traditional methods like random forest and SVM. However, convolutional network based deep learning models that use one-hot encoding DNA sequence as input may discover limited information and cause unsatisfactory prediction performance, so more data and model structures of diverse angles should be considered. In this work, we proposed a hybrid sequence-based deep learning model with both MeDIP-seq data and Histone information to predict DNA methylated CpG states (MHCpG). We combined both MeDIP-seq data and histone modification data with sequence information and implemented convolutional network to discover sequence patterns. In addition, we used statistical data gained from previous three input data and adopted a 3-layer feedforward neuron network to extract more high-level features. We compared our method with traditional predicting methods using random forest and other previous methods like CpGenie and DeepCpG, the result showed that MHCpG exceeded the other approaches and gained more satisfactory performance. [ABSTRACT FROM AUTHOR]

In our previous study, we found that subchronic exposure of chlorpyrifos (CPF) can cause reproductive damage in male rats. However, the mechanisms underlying the reproductive effects of CPF are not well understood. DNA methylation is essential for epigenetic gene regulation in development and disease. Therefore, we aim to compare DNA methylation profiles between controls and CPF-treated rats in order to identify the epigenetic mechanism of male reproductive toxicity induced by CPF. Methylated DNA immunoprecipitation with high-throughput sequencing (MeDIP-seq) was used to investigate the genome-wide DNA methylation pattern in testes of control and CPF-treated rats for 90 days.We identified 27 019 differentially methylated regions (DMRs) (14 150 upmethylated and 12 869 downmethylated) between CPF-exposed and control groups. The DMR-related genes are mainly involved in 113 pathways predicted by Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The result showed that high methylation gene PIK3CD may play a key role in epigenetic regulation of multiple pathways, such as Ras signaling pathway, AGE-RAGE signaling pathway in diabetic complications, HIF-1 signaling pathway, VEGF signaling pathway, and glioma and Fc epsilon RI signaling pathway in rats exposed to CPF. Our study provides significant explanations for the epigenetic mechanism of male reproductive toxicology induced by CPF. [ABSTRACT FROM AUTHOR]

DNA methylation is a type of epigenetic modification that alters gene expression without changing the DNA sequence and mediates some cases of phenotypic plasticity. In this study, we identified six DNA methyltransferase (DNMT) genes and two methyl-CpG binding domain protein2 (MBD2) gene from Pinctada fucata martensii. We also analyzed the genome-wide DNA methylation levels of mantle edge (ME) and mantle central (MC) from P. f. martensii via methylated immunoprecipitation sequencing (MeDIP-Seq). Results revealed that both ME and MC had 122 million reads, and had 58,702 and 55,721 peaks, respectively. The obtained methylation patterns of gene elements and repeats showed that the methylation of the protein-coding genes, particularly intron and coding exons (CDSs), was more frequent than that of other genomic elements in the pearl oyster genome. We combined the methylation data with the RNA-seq data of the ME and MC of P. f. martensii and found that promoter, CDS, and intron methylation levels were positively correlated with gene expression levels except the highest gene expression level. We also identified 313 differential methylation genes (DMGs) and annotated 212 of them. These DMGs were significantly enriched in 30 pathways, such as amino acid and protein metabolism, energy metabolism, terpenoid synthesis, and immune-related pathways. This study comprehensively analyzed the methylomes of biomineralization-related tissues and helped enhance our understanding of the regulatory mechanism underlying shell formation. [ABSTRACT FROM AUTHOR]

Heavy metal toxicity is a significant global health concern, with particular attention given to lead (Pb) exposure due to its adverse effects on cognitive development, especially in children exposed to low concentrations. While Pb neurotoxicity has been extensively studied, the analysis and molecular mechanisms underlying the transgenerational effects of Pb exposure-induced neurotoxicity remain poorly understood. In this study, we utilized Drosophila , a powerful developmental animal model, to investigate this phenomenon. Our findings demonstrated that Pb exposure during the developmental stage had a profound effect on the neurodevelopment of F0 fruit flies. Specifically, we observed a loss of correlation between the terminal motor area and muscle fiber area, along with an increased frequency of the β-lobe midline crossing phenotype in mushroom bodies. Western blot analysis indicated altered expression levels of synaptic vesicle proteins, with a decrease in Synapsin (SYN) and an increase in Bruchpilot (BRP) expression, suggesting changes in synaptic vesicle release sites. These findings were corroborated by electrophysiological data, showing an increase in the amplitude of evoked excitatory junctional potential (EJP) and an increase in the frequency of spontaneous excitatory junctional potential (mEJP) following Pb exposure. Importantly, our results further confirmed that the developmental neurotoxicity resulting from grandparental Pb exposure exhibited a transgenerational effect. The F3 offspring displayed neurodevelopmental defects, synaptic function abnormalities, and repetitive behavior despite lacking direct Pb exposure. Our MeDIP-seq analysis further revealed significant alterations in DNA methylation levels in several neurodevelopmental associated genes (eagle , happyhour , neuroglian , bazooka , and spinophilin) in the F3 offspring exposed to Pb. These findings suggest that DNA methylation modifications may underlie the inheritance of acquired phenotypic traits resulting from environmental Pb exposure. [Display omitted] • Lead exposure induces neurodevelopmental abnormalities. • Lead exposure induces abnormal synaptic function in Neuromuscular junctions. • Lead exposure has transgenerational developmental neurotoxicity. • Lead exposure can alter genomic methylation in Drosophila melanogaster. [ABSTRACT FROM AUTHOR]

Abstract Background Switchgrass is a crop with many desirable traits for bioenergy production. Plant genomes have high DNA methylation levels throughout genes and transposable elements and DNA methylation is known to play a role in silencing transposable elements. Here we analyzed methylomes in two switchgrass genotypes AP13 and VS16. AP13 is derived from a lowland ecotype and VS16, typically considered drought-tolerant, is derived from an upland ecotype, both genotypes are tetraploid (2n = 4× = 36). Results Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) and bisulfite-sequencing (BS-seq) were used to profile DNA methylation in genomic features of AP13 and VS16. The methylation patterns in genes and transposable elements were similar to other plants, however, overall CHH methylation levels were comparatively low. Differentially methylated regions (DMRs) were assessed and a total of 1777 CG-DMRs, 573 CHG-DMRs, and 3 CHH-DMRs were detected between the two genotypes. TEs and their flanking regions were higher than that of genic regions. Different types of TEs had different methylation patterns, but the two LTRs (Copia and Gypsy) were similarly methylated, while LINEs and DNA transposons typically had different methylation patterns. MeDIP-seq data was compared to BS-seq data and most of the peaks generated by MeDIP-seq were confirmed to be highly methylated by BS-seq. Conclusions DNA methylation in switchgrass genotypes obtained from the two ecotypes were found similar. Collinear gene pairs in two subgenomes (A and B) were not significantly differentially methylated. Both BS-seq and MeDIP-seq methodologies were found effective. Methylation levels were highest at CG and least in CHH. Increased DNA methylation was seen in TEs compared to genic regions. Exploitation of TE methylations can be a viable option in future crop improvement.

DNA methylation is an important form of epigenetic regulation that can regulate the expression of genes and the development of tissues. Muscle satellite cells play an important role in skeletal muscle development and regeneration. Therefore, the DNA methylation status of genes in satellite cells is important in the regulation of the development of skeletal muscle. This study systematically investigated the changes of genome-wide DNA methylation in satellite cells during skeletal muscle development. According to the MeDIP-Seq data, 52,809–123,317 peaks were obtained for each sample, covering 0.70–1.79% of the genome. The number of reads and peaks was highest in the intron regions followed by the CDS regions. A total of 96,609 DMRs were identified between any two time points. Among them 6198 DMRs were annotated into the gene promoter regions, corresponding to 4726 DMGs. By combining the MeDIP-Seq and RNA-Seq data, a total of 202 overlap genes were obtained between DMGs and DEGs. GO and Pathway analysis revealed that the overlap genes were mainly involved in 128 biological processes and 23 pathways. Among the biological processes, terms related to regulation of cell proliferation and Wnt signaling pathway were significantly different. Gene–gene interaction analysis showed that Wnt5a, Wnt9a, and Tgfβ1 were the key nodes in the network. Furthermore, the expression level of Wnt5a, Wnt9a, and Tgfβ1 genes could be influenced by the methylation status of promoter region during skeletal muscle development. These results indicated that the Wnt and Tgfβ signaling pathways may play an important role in functional regulation of satellite cells, and the DNA methylation status of Wnt and Tgfβ signals is a key regulatory factor during skeletal muscle development. This study provided new insights into the effects of genome-wide methylation on the function of satellite cells.

DNA methylation is an important form of epigenetic regulation that can regulate the expression of genes and the development of tissues. Muscle satellite cells play an important role in skeletal muscle development and regeneration. Therefore, the DNA methylation status of genes in satellite cells is important in the regulation of the development of skeletal muscle. This study systematically investigated the changes of genome-wide DNA methylation in satellite cells during skeletal muscle development. According to the MeDIP-Seq data, 52,809–123,317 peaks were obtained for each sample, covering 0.70–1.79% of the genome. The number of reads and peaks was highest in the intron regions followed by the CDS regions. A total of 96,609 DMRs were identified between any two time points. Among them 6198 DMRs were annotated into the gene promoter regions, corresponding to 4726 DMGs. By combining the MeDIP-Seq and RNA-Seq data, a total of 202 overlap genes were obtained between DMGs and DEGs. GO and Pathway analysis revealed that the overlap genes were mainly involved in 128 biological processes and 23 pathways. Among the biological processes, terms related to regulation of cell proliferation and Wnt signaling pathway were significantly different. Gene–gene interaction analysis showed that Wnt5a , Wnt9a , and Tgf β 1 were the key nodes in the network. Furthermore, the expression level of Wnt5a , Wnt9a , and Tgf β 1 genes could be influenced by the methylation status of promoter region during skeletal muscle development. These results indicated that the Wnt and Tgfβ signaling pathways may play an important role in functional regulation of satellite cells, and the DNA methylation status of Wnt and Tgfβ signals is a key regulatory factor during skeletal muscle development. This study provided new insights into the effects of genome-wide methylation on the function of satellite cells. [ABSTRACT FROM AUTHOR]

Background: DNA methylation of CpG dinucleotides is an essential epigenetic modification that plays a key role in transcription. Widely used DNA enrichment-based methods offer high coverage for measuring methylated CpG dinucleotides, with the lowest cost per CpG covered genome-wide. However, these methods measure the DNA enrichment of methyl-CpG binding, and thus do not provide information on absolute methylation levels. Further, the enrichment is influenced by various confounding factors in addition to methylation status, for example, CpG density. Computational models that can accurately derive absolute methylation levels from DNA enrichment data are needed. Results: We developed "MeDEStrand," a method that uses a sigmoid function to estimate and correct the CpG bias from enrichment results to infer absolute DNA methylation levels. Unlike previous methods, which estimate CpG bias based on reads mapped at the same genomic loci, MeDEStrand processes the reads for the positive and negative DNA strands separately. We compared the performance of MeDEStrand to that of three other state-of-the-art methods "MEDIPS," "BayMeth," and "QSEA" on four independent datasets generated using immortalized cell lines (GM12878 and K562) and human primary cells (foreskin fibroblasts and mammary epithelial cells). Based on the comparison of the inferred absolute methylation levels from MeDIP-seq data and the corresponding reduced-representation bisulfite sequencing data from each method, MeDEStrand showed the best performance at high resolution of 25, 50, and 100 base pairs. Conclusions: The MeDEStrand tool can be used to infer whole-genome absolute DNA methylation levels at the same cost of enrichment-based methods with adequate accuracy and resolution. R package MeDEStrand and its tutorial is freely available for download at https://github.com/jxu1234/MeDEStrand.git. [ABSTRACT FROM AUTHOR]

Background: Milk yield for Holstein cows has doubled over five decades due to genetic selection and changes to management, but the molecular mechanisms that facilitated this increase are mostly unknown. Epigenetic modifications to the cattle genome are a plausible molecular mechanism to cause variation in milk yield and our objective was to describe genome-wide DNA methylation patterns in peripheral blood mononuclear cells (PBMC) from mature Holstein dairy cows with variable milk yield. Results: Whole genome MeDIP-seq was performed following DNA extraction from PBMC of 6 lactating dairy cows from 4 different herds that varied in milk yield from 13,556 kg to 23,105 kg per 305 day lactation. We describe methylation across the genome and for 13,677 protein coding genes. Repetitive element reads were primarily mapped to satellite (36.4%), SINE (29.1%), and LINE (23.7%) regions and the majority (78.4%) of CpG sites were sequenced at least once. DNA methylation was generally low upstream of genes with the nadir occurring 95 bp prior to the transcription start site (TSS). Methylation was lower in the first exon than in later exons, was highest for introns near the intron-exon junctions, and declined downstream as the distance from the gene increased. We identified 72 differentially methylated regions (DMR) between high milk yield cows and their control, and 252 DMR across herd environments. Conclusions: This reference methylome for cattle with extreme variation in milk yield phenotype provides a resource to more fully evaluate relationships between DNA methylation and phenotype in populations subject to selection. The detection of DMR in cows of varying milk yield suggests potential to exploit epigenetic variation in cattle improvement programs. [ABSTRACT FROM AUTHOR]

DNA methylation is an important epigenetic modification involved in the estrous cycle and the regulation of reproduction. Here, we investigated the genome-wide profiles of DNA methylation in porcine ovaries in proestrus and estrus using methylated DNA immunoprecipitation sequencing. The results showed that DNA methylation was enriched in intergenic and intron regions. The methylation levels of coding regions were higher than those of the 5=- and 3=-flanking regions of genes. There were 4,813 differentially methylated regions (DMRs) of CpG islands in the estrus vs. proestrus ovarian genomes. Additionally, 3,651 differentially methylated genes (DMGs) were identified in pigs in estrus and proestrus. The DMGs were significantly enriched in biological processes and pathways related to reproduction and hormone regulation. We identified 90 DMGs associated with regulating reproduction in pigs. Our findings can serve as resources for DNA methylome research focused on porcine ovaries and further our understanding of epigenetically regulated reproduction in mammals. [ABSTRACT FROM AUTHOR]

In this study, the methylation of mitochondrial genome in the immature testis of Chinese mitten crab Eriocheir sinensis of the Yangtze River system was determined for the first time using MeDIP-seq. Our methylated DNA fragments covered more than 99% of the mitochondrial genome in E. sinensis loaded from GenBank. There were 8 mutated bases and 42 SNPs in the crab mitochondrial genome. The methylation presented in all genes as well as in an A + T region, but less in intergenic regions in the mitochondrial genome. However, the level of methylation of most genes coding proteins and the A + T region were high. But, the majority of genes encoding tRNAs were hypomethylated, and both the rRNA genes also showed methylation of low or median frequency. Especially, the level of methylation of the intergenic regions is the lowest. Those features indicated that the methylation of DNA may play an important role in gene expressing regulation in the mitochondrial genome of immature testis in E. sinensis. [ABSTRACT FROM AUTHOR]

The biological structure and function of the mammalian testis undergo important developmental changes during prepuberty and DNA methylation is dynamically regulated during testis development. In this study, we generated the first genome-wide DNA methylation profile of prepubertal porcine testis using methyl-DNA immunoprecipitation (MeDIP) combined with high-throughput sequencing (MeDIP-seq). Over 190 million high-quality reads were generated, containing 43 642 CpG islands. There was an overall downtrend of methylation during development, which was clear in promoter regions but less so in gene-body regions. We also identified thousands of differentially methylated regions (DMRs) among the three prepubertal time points (1 month, T1; 2 months, T2; 3 months, T3), the majority of which showed decreasing methylation levels over time. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that many genes in the DMRs were linked with cell proliferation and some important pathways in porcine testis development. Our data suggest that DNA methylation plays an important role in prepubertal development of porcine testis, with an obvious downtrend of methylation levels from T1 to T3. Overall, our study provides a foundation for future studies and gives new insights into mammalian testis development. [ABSTRACT FROM AUTHOR]

Purpose: We investigated the genomic DNA methylation profile of prostate cancer in transgenic adenocarcinoma of the mouse prostate (TRAMP) cancer model and to analyze the crosstalk among targeted genes and the related functional pathways. Methods: Prostate DNA samples from 24-week-old TRAMP and C57BL/6 male mice were isolated. The DNA methylation profiles were analyzed by methylated DNA immunoprecipitation (MeDIP) followed by next-generation sequencing (MeDIP-seq). Canonical pathways, diseases and function and network analyses of the different samples were then performed using the Ingenuity ® Pathway Analysis (IPA) software. Some target genes with significant difference in methylation were selected for validation using methylation specific primers (MSP) and qPCR. Results: TRAMP mice undergo extensive aberrant CpG hyper- and hypo-methylation. There were 2147 genes with a significant (log2-change ⩽ 2) change in CpG methylation between the two groups, as mapped by the IPA software. Among these genes, the methylation of 1105 and 1042 genes was significantly decreased and increased, respectively, in TRAMP prostate tumors. The top associated disease identified by IPA was adenocarcinoma; however, the cAMP response element-binding protein (CREB)-, histone deacetylase 2 (HDAC2)-, glutathione S-transferase pi (GSTP1)- and polyubiquitin-C (UBC)-related pathways showed significantly altered methylation profiles based on the canonical pathway and network analyses. MSP and qPCR results of genes of interests corroborated with MeDIP-seq findings. Conclusions: This is the first MeDIP-seq with IPA analysis of the TRAMP model to provide novel insight into the genome-wide methylation profile of prostate cancer. Studies on epigenetics, such as DNA methylation, will potentially provide novel avenues and strategies for further development of biomarkers targeted for treatment and prevention approaches for prostate cancer. [ABSTRACT FROM AUTHOR]

Tibetan chickens have unique adaptations to the extreme high-altitude environment that they inhabit. Epigenetic DNA methylation affects many biological processes, including hypoxic adaptation; however, the regulatory genes for DNA methylation in hypoxic adaptation remain unknown. In this study, methylated DNA immunoprecipitation with high-throughput sequencing (MeDIP-seq) was used to provide an atlas of the DNA methylomes of the heart tissue of hypoxic highland Tibetan and lowland Chahua chicken embryos. A total of 31.2 gigabases of sequence data were generated from six MeDIP-seq libraries. We identified 1,049 differentially methylated regions (DMRs) and 695 related differentially methylated genes (DMGs) between the two chicken breeds. The DMGs are involved in vascular smooth muscle contraction, VEGF signaling pathway, calcium signaling pathway, and other hypoxia-related pathways. Five candidate genes that had low methylation (EDNRA, EDNRB2, BMPR1B, BMPRII, and ITGA2) might play key regulatory roles in the adaptation to hypoxia in Tibetan chicken embryos. Our study provides significant explanations for the functions of genes and their epigenetic regulation for hypoxic adaptation in Tibetan chickens.

Chromatin structure, DNA methylation, and histone modifications act in a concerted manner to influence gene expression and therefore plant phenotypes. Environmental stresses are often associated with extensive chromatin rearrangements and modifications of epigenetic levels and patterns. Stress-tolerant plants can be a good tool to unveil potential connections between specific epigenetic modifications and stress tolerance capacity. We analyzed genome wide DNA methylation of a salt-tolerant rice variety under salinity and identified a set of differentially methylated regions (DMRs) between control and stress samples using high-throughput sequencing of DNA immunoprecipitated with the 5-methylcytosine antibody (MeDIP-Seq). The examination of DNA methylation pattern at DMRs regions revealed a general tendency for demethylation events in stress samples as compared to control. In addition, DMRs appear to influence the expression of genes located in their vicinity. We hypothesize that short regions as DMRs can shape the chromatin landscape of specific genomic regions and, therefore, may modulate the function of several genes. In this sense, the identification of DMRs represents one step towards to uncover new players in the regulation of stress-responsive genes and new target genes with potential application in enhancement of plant salinity-tolerance.

Wide hybridization is a common and efficient breeding strategy for enhancing crop yield and quality. An interesting phenomenon is that the reciprocal hybrids usually show different phenotypes, and its underlying mechanism is not well understood. Here, we reported our comparative analysis of the DNA methylation patterns in Solanum lycopersicum, Solanum pimpinellifolium and their reciprocal hybrids by methylated DNA immunoprecipitation sequencing. The reciprocal hybrids had lower levels of DNA methylation in CpG islands and LTR retroelements when compared with those of their parents. Importantly, remarkable differences in DNA methylation patterns, mainly in introns and CDS regions, were revealed between the reciprocal hybrids. These different methylated regions were mapped to 79 genes, 14 of which were selected for analysis of gene expression levels. While there was an inverse correlation between DNA methylation and gene expression in promoter regions, the relationship was complicated in gene body regions. Further association analysis revealed that there were 15 differentially methylated genes associated with siRNAs, and that the methylation levels of these genes were inversely correlated with respective siRNAs. All these data raised the possibility that the direction of hybridization induced the divergent epigenomes leading to changes in the transcription levels of reciprocal hybrids. [ABSTRACT FROM AUTHOR]

Background: Switchgrass is a crop with many desirable traits for bioenergy production. Plant genomes have high DNA methylation levels throughout genes and transposable elements and DNA methylation is known to play a role in silencing transposable elements. Here we analyzed methylomes in two switchgrass genotypes AP13 and VS16. AP13 is derived from a lowland ecotype and VS16, typically considered drought-tolerant, is derived from an upland ecotype, both genotypes are tetraploid (2n = 4× = 36). Results: Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) and bisulfite-sequencing (BS-seq) were used to profile DNA methylation in genomic features of AP13 and VS16. The methylation patterns in genes and transposable elements were similar to other plants, however, overall CHH methylation levels were comparatively low. Differentially methylated regions (DMRs) were assessed and a total of 1777 CG-DMRs, 573 CHG-DMRs, and 3 CHH-DMRs were detected between the two genotypes. TEs and their flanking regions were higher than that of genic regions. Different types of TEs had different methylation patterns, but the two LTRs (Copia and Gypsy) were similarly methylated, while LINEs and DNA transposons typically had different methylation patterns. MeDIP-seq data was compared to BS-seq data and most of the peaks generated by MeDIP-seq were confirmed to be highly methylated by BS-seq. Conclusions: DNA methylation in switchgrass genotypes obtained from the two ecotypes were found similar. Collinear gene pairs in two subgenomes (A and B) were not significantly differentially methylated. Both BS-seq and MeDIP-seq methodologies were found effective. Methylation levels were highest at CG and least in CHH. Increased DNA methylation was seen in TEs compared to genic regions. Exploitation of TE methylations can be a viable option in future crop improvement. [ABSTRACT FROM AUTHOR]