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.

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]

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]

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]

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]

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]

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]

We conducted a comparative study of multiplexed affinity enrichment sequence methodologies (MBD-seq and MeDIP-seq) in a rodent model of schizophrenia, induced by in utero methylazoxymethanol acetate (MAM) exposure. We also examined related gene expression changes using a pooled sample approach. MBD-seq and MeDIP-seq identified 769 and 1771 differentially methylated regions (DMRs) between F2 offspring of MAM-exposed rats and saline control rats, respectively. The assays showed good concordance, with ~ 56% of MBD-seq-detected DMRs being identified by or proximal to MeDIP-seq DMRs. There was no significant overlap between DMRs and differentially expressed genes, suggesting that DNA methylation regulatory effects may act upon more distal genes, or are too subtle to detect using our approach. Methylation and gene expression gene ontology enrichment analyses identified biological processes important to schizophrenia pathophysiology, including neuron differentiation, prepulse inhibition, amphetamine response, and glutamatergic synaptic transmission regulation, reinforcing the utility of the MAM rodent model for schizophrenia research. [ABSTRACT FROM AUTHOR]

Background: DNA methylation is an important epigenetic mechanism involved in gene regulation, with alterations in DNA methylation in the nuclear genome being linked to numerous complex diseases. Mitochondrial DNA methylation is a phenomenon that is receiving ever-increasing interest, particularly in diseases characterized by mitochondrial dysfunction; however, most studies have been limited to the investigation of specific target regions. Analyses spanning the entire mitochondrial genome have been limited, potentially due to the amount of input DNA required. Further, mitochondrial genetic studies have been previously confounded by nuclear-mitochondrial pseudogenes. Methylated DNA Immunoprecipitation Sequencing is a technique widely used to profile DNA methylation across the nuclear genome; however, reads mapped to mitochondrial DNA are often discarded. Here, we have developed an approach to control for nuclear-mitochondrial pseudogenes within Methylated DNA Immunoprecipitation Sequencing data. We highlight the utility of this approach in identifying differences in mitochondrial DNA methylation across regions of the human brain and pre-mortem blood. Results: We were able to correlate mitochondrial DNA methylation patterns between the cortex, cerebellum and blood. We identified 74 nominally significant differentially methylated regions (p < 0.05) in the mitochondrial genome, between anatomically separate cortical regions and the cerebellum in matched samples (N = 3 matched donors). Further analysis identified eight significant differentially methylated regions between the total cortex and cerebellum after correcting for multiple testing. Using unsupervised hierarchical clustering analysis of the mitochondrial DNA methylome, we were able to identify tissue-specific patterns of mitochondrial DNA methylation between blood, cerebellum and cortex. Conclusions: Our study represents a comprehensive analysis of the mitochondrial methylome using pre-existing Methylated DNA Immunoprecipitation Sequencing data to identify brain region-specific patterns of mitochondrial DNA methylation. [ABSTRACT FROM AUTHOR]

Background: The association of in vitro fertilisation (IVF) and DNA methylation has been studied predominantly at regulatory regions of imprinted genes and at just thousands of the ∼28 million CpG sites in the human genome. Methods: We investigated the links between IVF and DNA methylation patterns in whole cord blood cells (n = 98) and cord blood mononuclear cells (n = 82) from newborn twins using genome-wide methylated DNA immunoprecipitation coupled with deep sequencing. Results: At a false discovery rate (FDR) of 5%, we identified one significant whole blood DNA methylation change linked to conception via IVF, which was located ∼3 kb upstream of TNP1, a gene previously linked to male infertility. The 46 most strongly associated signals (FDR of 25%) included a second region in a gene also previously linked to infertility, C9orf3, suggesting that our findings may in part capture the effect of parental subfertility. Using twin modelling, we observed that individual-specific environmental factors appear to be the main overall contributors of methylation variability at the FDR 25% IVF-associated differentially methylated regions, although evidence for methylation heritability was also obtained at several of these regions. We replicated previous findings of differential methylation associated with IVF at the H19/IGF2 region in cord blood mononuclear cells, and we validated the signal at C9orf3 in monozygotic twins. We also explored the impact of intracytoplasmic sperm injection on the FDR 25% signals for potential effects specific to male or female infertility factors. Conclusions: To our knowledge, this is the most comprehensive study of DNA methylation profiles at birth and IVF conception to date, and our results show evidence for epigenetic modifications that may in part reflect parental subfertility. [ABSTRACT FROM AUTHOR]

Differences in drought stress tolerance within diverse rice genotypes have been attributed to genetic diversity and epigenetic alterations. DNA methylation is an important epigenetic modification that influences diverse biological processes, but its effects on rice drought stress tolerance are poorly understood. In this study, methylated DNA immunoprecipitation sequencing and an Affymetrix GeneChip rice genome array were used to profile the DNA methylation patterns and transcriptomes of the drought-tolerant introgression line DK151 and its drought-sensitive recurrent parent IR64 under drought and control conditions. The introgression of donor genomic DNA induced genome-wide DNA methylation changes in DK151 plants. A total of 1190 differentially methylated regions (DMRs) were detected between the two genotypes under normal growth conditions, and the DMR-associated genes in DK151 plants were mainly related to stress response, programmed cell death, and nutrient reservoir activity, which are implicated to constitutive drought stress tolerance. A comparison of the DNA methylation changes in the two genotypes under drought conditions indicated that DK151 plants have a more stable methylome, with only 92 drought-induced DMRs, than IR64 plants with 506 DMRs. Gene ontology analyses of the DMR-associated genes in drought-stressed plants revealed that changes to the DNA methylation status of genotype-specific genes are associated with the epigenetic regulation of drought stress responses. Transcriptome analysis further helped to identify a set of 12 and 23 DMR-associated genes that were differentially expressed in DK151 and IR64, respectively, under drought stress compared with respective controls. Correlation analysis indicated that DNA methylation has various effects on gene expression, implying that it affects gene expression directly or indirectly through diverse regulatory pathways. Our results indicate that drought-induced alterations to DNA methylation may influence an epigenetic mechanism that regulates the expression of unique genes responsible for drought stress tolerance.

Common bean (Phaseolus vulgaris L.) is economically important for its high protein, fiber, and micronutrient contents, with a relatively small genome size of ~587 Mb. Common bean is genetically diverse with two major gene pools, Meso-American and Andean. The phenotypic variability within the genotypes is partly attributed to the genetic diversity and epigenetic changes that are largely influenced by environmental factors. It is well established that an important epigenetic regulator of gene expression is DNA methylation. Here, we present results generated from two high-throughput sequencing technologies, methylated DNA immunoprecipitation-sequencing (MeDIP-seq) and whole genome bisulfite-sequencing (BS-Seq). Our analyses revealed that common bean displays similar methylation patterns as other previously published plant methylomes, with CG ~50%, CHG ~30%, and CHH ~2.7% methylation, however, these differ from the common bean reference methylome of Andean origin. We identified higher CG methylation levels in both promoter and genic regions than CHG and CHH contexts. Moreover, we found relatively higher CG methylation levels in genes than in promoters. Conversely, the CHG and CHH methylation levels were highest in promoters than in genes. This is the first genome-wide DNA methylation profiling study in a Meso-American common bean cultivar (Sierra) using NGS approaches. Our long-term goal is to generate genome-wide epigenomic maps in common bean focusing on chromatin accessibility, histone modifications, and DNA methylation.

DNA methylation plays major roles in many biological processes, including aging, carcinogenesis, and development. Analyses of DNA methylation using next-generation sequencing offer a new way to profile and compare methylomes across the genome in the context of aging. We explored genomewide DNA methylation and the effects of short-term calorie restriction ( CR) on the methylome of aged rat kidney. Whole-genome methylation of kidney in young (6 months old), old (25 months old), and OCR (old with 4-week, short-term CR) rats was analyzed by methylated DNA immunoprecipitation and next-generation sequencing (Me DIP-Seq). CpG islands and repetitive regions were hypomethylated, but 5′- UTR, exon, and 3′- UTR hypermethylated in old and OCR rats. The methylation in the promoter and intron regions was decreased in old rats, but increased in OCR rats. Pathway enrichment analysis showed that the hypermethylated promoters in old rats were associated with degenerative phenotypes such as cancer and diabetes. The hypomethylated promoters in old rats related significantly to the chemokine signaling pathway. However, the pathways significantly enriched in old rats were not observed from the differentially methylated promoters in OCR rats. Thus, these findings suggest that short-term CR could partially ameliorate age-related methylation changes in promoters in old rats. From the epigenomic data, we propose that the hypermethylation found in the promoter regions of disease-related genes during aging may indicate increases in susceptibility to age-related diseases. Therefore, the CR-induced epigenetic changes that ameliorate age-dependent aberrant methylation may be important to CR's health- and life-prolonging effects. [ABSTRACT FROM AUTHOR]

Background: DNA methylation is an important regulator of gene expression and chromatin structure. Methylated DNA immunoprecipitation sequencing (MeDIP-Seq) is commonly used to identify regions of DNA methylation in eukaryotic genomes. Within MeDIP-Seq libraries, methylated cytosines can be found in both double-stranded (symmetric) and single-stranded (asymmetric) genomic contexts. While symmetric CG methylation has been relatively well-studied, asymmetric methylation in any dinucleotide context has received less attention. Importantly, no currently available software for processing MeDIP-Seq reads is able to resolve these strand-specific DNA methylation signals. Here we introduce DISMISS, a new software package that detects strand-associated DNA methylation from existing MeDIP-Seq analyses. Results: Using MeDIP-Seq datasets derived from Apis mellifera (honeybee), an invertebrate species that contains more asymmetric- than symmetric- DNA methylation, we demonstrate that DISMISS can identify strand-specific DNA methylation signals with similar accuracy as bisulfite sequencing (BS-Seq; single nucleotide resolution methodology). Specifically, DISMISS is able to confidently predict where DNA methylation predominates (plus or minus DNA strands - asymmetric DNA methylation; plus and minus DNA stands - symmetric DNA methylation) in MeDIP-Seq datasets derived from A. mellifera samples. When compared to DNA methylation data derived from BS-Seq analysis of A. mellifera worker larva, DISMISS-mediated identification of strand-specific methylated cytosines is 80 % accurate. Furthermore, DISMISS can correctly (p <0.0001) detect the origin (sense vs antisense DNA strands) of DNA methylation at splice site junctions in A. mellifera MeDIP-Seq datasets with a precision close to BS-Seq analysis. Finally, DISMISS-mediated identification of DNA methylation signals associated with upstream, exonic, intronic and downstream genomic loci from A. mellifera MeDIP-Seq datasets outperforms MACS2 (Model-based Analysis of ChIP-Seq2; a commonly used MeDIP-Seq analysis software) and closely approaches the results achieved by BS-Seq. Conclusions: While asymmetric DNA methylation is increasingly being found in growing numbers of eukaryotic species and is the predominant pattern observed in some invertebrate genomes, it has been difficult to detect in MeDIP-Seq datasets using existing software. DISMISS now enables more sensitive examinations of MeDIP-Seq datasets and will be especially useful for the study of genomes containing either low levels of DNA methylation or for genomes containing relatively high amounts of asymmetric methylation. [ABSTRACT FROM AUTHOR]

Background: In utero and early-life experienced environmental exposures are suggested to play an important role in many multifactorial diseases potentially mediated through lasting effects on the epigenome. As the epigenome in addition remains modifiable throughout life, identifying specific disease-relevant biomarkers may prove challenging. This has led to an increased interest in epigenome-wide association studies using dried blood spots (DBS) routinely collected in perinatal screening programs. Such programs are in place in numerous countries around the world producing large and unique biobanks. However, availability of this biological material is highly limited as each DBS is made only from a few droplets of blood and storage conditions may be suboptimal for epigenetic studies. Furthermore, as relevant markers may reside outside gene bodies, epigenome-wide interrogation is needed. Results: Here we demonstrate, as a proof of principle, that genome-wide interrogation of the methylome based on methylated DNA immunoprecipitation coupled with next-generation sequencing (MeDIP-seq) is feasible using a single 3.2 mm DBS punch (60 ng DNA) from filter cards archived for up to 16 years. The enrichment profile, sequence quality and distribution of reads across genetic regions were comparable between samples archived 16 years, 4 years and a freshly prepared control sample. Conclusions: In summary, we show that high-quality MeDIP-seq data is achievable from neonatal screening filter cards stored at room temperature, thereby providing information on annotated as well as on non-RefSeq genes and repetitive elements. Moreover, the quantity of DNA from one DBS punch proved sufficient allowing for multiple epigenome studies using one single DBS. [ABSTRACT FROM AUTHOR]