BioM2: biologically informed multi-stage machine learning for phenotype prediction using omics data.

Navigating the complex landscape of high-dimensional omics data with machine learning models presents a significant challenge. The integration of biological domain knowledge into these models has shown promise in creating more meaningful stratifications of predictor variables, leading to algorithms that are both more accurate and generalizable. However, the wider availability of machine learning tools capable of incorporating such biological knowledge remains limited. Addressing this gap, we introduce BioM2, a novel R package designed for biologically informed multistage machine learning. BioM2 uniquely leverages biological information to effectively stratify and aggregate high-dimensional biological data in the context of machine learning. Demonstrating its utility with genome-wide DNA methylation and transcriptome-wide gene expression data, BioM2 has shown to enhance predictive performance, surpassing traditional machine learning models that operate without the integration of biological knowledge. A key feature of BioM2 is its ability to rank predictor variables within biological categories, specifically Gene Ontology pathways. This functionality not only aids in the interpretability of the results but also enables a subsequent modular network analysis of these variables, shedding light on the intricate systems-level biology underpinning the predictive outcome. We have proposed a biologically informed multistage machine learning framework termed BioM2 for phenotype prediction based on omics data. BioM2 has been incorporated into the BioM2 CRAN package (https://cran.r-project.org/web/packages/BioM2/index.html).

Genetic and epigenetic alterations in precursor lesions of endometrial endometrioid carcinoma.

The hyperplasia-carcinoma sequence is a stepwise tumourigenic programme towards endometrial cancer in which normal endometrial epithelium becomes neoplastic through non-atypical endometrial hyperplasia (NAEH) and atypical endometrial hyperplasia (AEH), under the influence of unopposed oestrogen. NAEH and AEH are known to exhibit polyclonal and monoclonal cell growth, respectively; yet, aside from focal PTEN protein loss, the genetic and epigenetic alterations that occur during the cellular transition remain largely unknown. We sought to explore the potential molecular mechanisms that promote the NAEH-AEH transition and identify molecular markers that could help to differentiate between these two states. We conducted target-panel sequencing on the coding exons of 596 genes, including 96 endometrial cancer driver genes, and DNA methylome microarrays for 48 NAEH and 44 AEH lesions that were separately collected via macro- or micro-dissection from the endometrial tissues of 30 cases. Sequencing analyses revealed acquisition of the PTEN mutation and the clonal expansion of tumour cells in AEH samples. Further, across the transition, alterations to the DNA methylome were characterised by hypermethylation of promoter/enhancer regions and CpG islands, as well as hypo- and hyper-methylation of DNA-binding regions for transcription factors relevant to endometrial cell differentiation and/or tumourigenesis, including FOXA2, SOX17, and HAND2. The identified DNA methylation signature distinguishing NAEH and AEH lesions was reproducible in a validation cohort with modest discriminative capability. These findings not only support the concept that the transition from NAEH to AEH is an essential step within neoplastic cell transformation of endometrial epithelium but also provide deep insight into the molecular mechanism of the tumourigenic programme. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Integrated analysis of the DNA methylome and RNA transcriptome during the development of skeletal muscle in Duroc pigs.

BACKGROUND: Skeletal muscle development plays a crucial role in yield and quality of pork; however, this process is influenced by various factors. In this study, we employed whole-genome bisulfite sequencing (WGBS) and transcriptome sequencing to comprehensively investigate the longissimus dorsi muscle (LDM), aiming to identify key genes that impact the growth and development of Duroc pigs with different average daily gains (ADGs). RESULTS: Eight pigs were selected and divided into two groups based on ADGs: H (774.89 g) group and L (658.77 g) group. Each pair of the H and L groups were half-siblings. The results of methylation sequencing revealed 2631 differentially methylated genes (DMGs) involved in metabolic processes, signalling, insulin secretion, and other biological activities. Furthermore, a joint analysis was conducted on these DMGs and the differentially expressed genes (DEGs) obtained from transcriptome sequencing of the same individual. This analysis identified 316 differentially methylated and differentially expressed genes (DMEGs), including 18 DMEGs in promoter regions and 294 DMEGs in gene body regions. Finally, LPAR1 and MEF2C were selected as candidate genes associated with muscle development. Bisulfite sequencing PCR (BSP) and quantitative real-time PCR (qRT-PCR) revealed that the promoter region of LPAR1 exhibited significantly lower methylation levels (P < 0.05) and greater expression levels (P < 0.05) in the H group than in the L group. Additionally, hypermethylation was observed in the gene body region of MEF2C, as was a low expression level, in the H group (P < 0.05). CONCLUSIONS: These results suggest that the differences in the ADGs of Duroc pigs fed the same diet may be influenced by the methylation levels and expression levels of genes related to skeletal muscle development.

SPEN haploinsufficiency causes a neurodevelopmental disorder overlapping proximal 1p36 deletion syndrome with an episignature of X chromosomes in females.

Deletion 1p36 (del1p36) syndrome is the most common human disorder resulting from a terminal autosomal deletion. This condition is molecularly and clinically heterogeneous. Deletions involving two non-overlapping regions, known as the distal (telomeric) and proximal (centromeric) critical regions, are sufficient to cause the majority of the recurrent clinical features, although with different facial features and dysmorphisms. SPEN encodes a transcriptional repressor commonly deleted in proximal del1p36 syndrome and is located centromeric to the proximal 1p36 critical region. Here, we used clinical data from 34 individuals with truncating variants in SPEN to define a neurodevelopmental disorder presenting with features that overlap considerably with those of proximal del1p36 syndrome. The clinical profile of this disease includes developmental delay/intellectual disability, autism spectrum disorder, anxiety, aggressive behavior, attention deficit disorder, hypotonia, brain and spine anomalies, congenital heart defects, high/narrow palate, facial dysmorphisms, and obesity/increased BMI, especially in females. SPEN also emerges as a relevant gene for del1p36 syndrome by co-expression analyses. Finally, we show that haploinsufficiency of SPEN is associated with a distinctive DNA methylation episignature of the X chromosome in affected females, providing further evidence of a specific contribution of the protein to the epigenetic control of this chromosome, and a paradigm of an X chromosome-specific episignature that classifies syndromic traits. We conclude that SPEN is required for multiple developmental processes and SPEN haploinsufficiency is a major contributor to a disorder associated with deletions centromeric to the previously established 1p36 critical regions.

Mdivi-1: Effective but complex mitochondrial fission inhibitor.

Mdivi-1, Mitochondrial DIVIsion inhibitor 1, has been widely employed in research under the assumption that it exclusively influences mitochondrial fusion, but effects other than mitochondrial dynamics have been underinvestigated. This paper provides transcriptome and DNA methylome-wide analysis for Mdivi-1 treated SH-SY5Y human neuroblastoma cells using RNA sequencing (RNA-seq) and methyl capture sequencing (MC-seq) methods. Gene ontology analysis of RNA sequences revealed that p53 transcriptional gene network and DNA replication initiation-related genes were significantly up and down-regulated, respectively, showing the correlation with the arrest cell cycle in the G1 phase. MC-seq, a powerful sequencing method for capturing DNA methylation status in CpG sites, revealed that although Mdivi-1 does not induce dramatic DNA methylation change, the subtle alterations were concentrated within the CpG island. Integrative analysis of both sequencing data disclosed that the p53 transcriptional network was activated while the Parkinson's disease pathway was halted. Next, we investigated several changes in mitochondria in response to Mdivi-1. Copy number and transcription of mitochondrial DNA were suppressed. ROS levels increased, and elevated ROS triggered mitochondrial retrograde signaling rather than inducing direct DNA damage. In this study, we could better understand the molecular network of Mdivi-1 by analyzing DNA methylation and mRNA transcription in the nucleus and further investigating various changes in mitochondria, providing inspiration for studying nuclear-mitochondrial communications.

CDK5R1, GSE1, HSPG2 and WDFY3 as indirect epigenetic-sensitive genes in atrial fibrillation.

BACKGROUND: Although mounting evidence supports that aberrant DNA methylation occurs in the hearts of patients with atrial fibrillation (AF), noninvasive epigenetic characterization of AF has not yet been defined. METHODS: We investigated DNA methylome changes in peripheral blood CD4+ T cells isolated from 10 patients with AF relative to 11 healthy subjects (HS) who were enrolled in the DIANA clinical trial (NCT04371809) via reduced-representation bisulfite sequencing (RRBS). RESULTS: An atrial-specific PPI network revealed 18 hub differentially methylated genes (DMGs), wherein ROC curve analysis revealed reasonable diagnostic performance of DNA methylation levels found within CDK5R1 (AUC = 0.76; p = 0.049), HSPG2 (AUC = 0.77; p = 0.038), WDFY3 (AUC = 0.78; p = 0.029), USP49 (AUC = 0.76; p = 0.049), GSE1 (AUC = 0.76; p = 0.049), AIFM1 (AUC = 0.76; p = 0.041), CDK5RAP2 (AUC = 0.81; p = 0.017), COL4A1 (AUC = 0.86; p < 0.001), SEPT8 (AUC = 0.90; p < 0.001), PFDN1 (AUC = 0.90; p < 0.01) and ACOT7 (AUC = 0.78; p = 0.032). Transcriptional profiling of the hub DMGs provided a significant overexpression of PSDM6 (p = 0.004), TFRC (p = 0.01), CDK5R1 (p < 0.001), HSPG2 (p = 0.01), WDFY3 (p < 0.001), USP49 (p = 0.004) and GSE1 (p = 0.021) in AF patients vs HS. CONCLUSIONS: CDK5R1, GSE1, HSPG2 and WDFY3 resulted the best discriminatory genes both at methylation and gene expression level. Our results provide several candidate diagnostic biomarkers with the potential to advance precision medicine in AF.

Decoding the DNA methylome of central nervous system tumors: An emerging modality for integrated diagnosis.

The definitive diagnosis and classification of individual cancers are crucial for patient care and cancer research. To achieve a robust diagnosis of central nervous system (CNS) tumors, a genotype-phenotype integrated diagnostic approach was introduced in recent versions of the World Health Organization classification, followed by the incorporation of a genome-wide DNA methylome-based classification. Microarray-based platforms are widely used to obtain DNA methylome data, and the German Cancer Research Center (Deutsches Krebsforschungszentrum [DKFZ]) has a webtool for a DNA methylation-based classifier (DKFZ classifier). Integration of DNA methylome will further enhance the precision of CNS tumor classification, especially in diagnostically challenging cases. However, in the clinical application of DNA methylome-based classification, challenges related to data interpretation persist, in addition to technical caveats, regulations, and limited accessibility. Dimensionality reduction (DMR) can complement integrated diagnosis by visualizing a profile and comparing it with other known samples. Therefore, DNA methylome-based classification is a highly useful research tool for auxiliary analysis in challenging diagnostic and rare disease cases, and for establishing novel tumor concepts. Decoding the DNA methylome, especially by DMR in addition to DKFZ classifier, emphasizes the capability of grasping the fundamental biological principles that provide new perspectives on CNS tumors.

The DNA methylation status of the serotonin metabolic pathway associated with reproductive inactivation induced by long-light exposure in Magang geese.

BACKGROUND: Domestic geese are seasonal breeders and have the lowest reproductive capacity among all poultry species. Magang geese is a topical short-day breeder, short photoperiod exposure stimulates its reproductive activity while long photoperiod inhibits. To explore epigenetic change that could influence reproductive activity, we performed whole genome bisulfite sequencing and transcriptome sequencing in the hypothalamus at three reproductive stages during long-light exposure in male Magang geese. RESULTS: A total number of 10,602 differentially methylated regions (DMRs) were identified among three comparison groups. We observed that the vast majority of DMRs were enriched in intron regions. By integrating the BS-sequencing and RNA-seq data, the correlation between methylation changes of CG DMRs and expression changes of their associated genes was significant only for genes containing CG DMRs in their intron. A total of 278 DMR-associated DEGs were obtained among the three stages. KEGG analysis revealed that the DMR-associated DEGs were mainly involved in 11 pathways. Among them, the neuroactive ligand-receptor interaction pathway was significantly enriched in both two comparisons (RA vs.RD and RD vs.RI); the Wnt signaling pathway, apelin signaling pathway, melanogenesis, calcium signaling pathway, focal adhesion, and adherens junction were significantly enriched in the RA vs. RI comparison. In addition, the expression level of two serotonin-metabolic genes was significantly altered during reproductive axis inactivation by the methylation status of their promoter region (TPH2) and intron region (SLC18A2), respectively. These results were confirmed by Bisulfite sequencing PCR (BSP), pyrosequencing, and real-time qPCR, indicating that serotonin metabolic signaling may play a key role in decreasing the reproductive activity of Magang geese induced by long-light exposure. Furthermore, we performed a metabolomics approach to investigate the concentration of neurotransmitters among the three stages, and found that 5-HIAA, the last product of the serotonin metabolic pathway, was significantly decreased in the hypothalamus during RI. CONCLUSIONS: Our study reveals that the methylation status of the serotonin metabolic pathway in the hypothalamus is associated with reproductive inactivation, and provided new insight into the effect of DNA methylation on the reproductive regulation of the hypothalamus in Magang geese.

DNA methylation trajectories during innate and adaptive immune responses of human B lymphocytes.

B lymphocytes can engage in either a rapid T cell-independent pathway (TI) or a delayed long-lasting T cell-independent (TD) response through highly ordered transcriptional programs, yet the detailed underlying mechanisms are unclear. Since DNA methylation plays a key role in controlling gene expression and lineage specification, we explored the dynamics of whole-genome cytosine modifications during the ex vivo response of human B cells isolated from normal individuals by negative selection. We found that B cell differentiation following TI and TD signalling is accompanied by extensive remodelling of the epigenome, including global gain and loss of DNA methylation. The epigenetic changes map to different regions of the B cell genome, including non- C-phosphate-G CpG islands, indicating that modifications of distal regulatory elements likely regulate specific gene transcription in B cells. Non-CpG methylation also occurs in differentiating human B cells, suggesting that this DNA modification is involved in transcriptional regulation of B cell genes with promoters exhibiting a low-density methylation, possibly by changing the chromatin shape that could have an impact on gene expression. Most strikingly, compared to TD activation, stimulation of B cells through an innate pathway induced higher levels of DNA methylation modifications at CpG, CHG and CGG contexts, supporting the view that DNA methylation modifications are used in distinct trajectories to specify the TI and TD B lymphocyte responses.

Epigenome alterations in food allergy: A systematic review of candidate gene and epigenome-wide association studies.

OBJECTIVE: The aim of this study was to systematically review the evidence across studies that assessed DNA methylome variations in association with food allergy (FA). DESIGN: A systematic review of the literature and meta-analysis were carried out within several databases. However, the risk of bias in the included articles was not evaluated. DATA SOURCES: PubMed, Cochrane Database of Systematic Reviews, and Web of Science were used to search up to July 2022. ELIGIBILITY CRITERIA: We included targeted and epigenome-wide association studies (EWASs) that assessed DNA methylome alterations in association with FA in adult or paediatric populations. RESULTS: Among 366 publications, only 16 were retained, which were mainly focused on FA in children. Seven candidate gene-targeted studies found associations in Th1/Th2 imbalance (IL4, IL5, IL10, INFG, IL2 and IL12B genes), regulatory T cell function (FOXP3 gene), Toll-like receptors pathway (TLR2, CD14 genes) and digestive barrier integrity (FLG gene). Nine EWAS assessed the association with peanut allergy (n = 3), cow's milk allergy (n = 2) or various food allergens (n = 4). They highlighted 11 differentially methylated loci in at least two studies (RPS6KA2, CAMTA1, CTBP2, RYR2, TRAPPC9, DOCK1, GALNTL4, HDAC4, UMODL1, ZAK and TNS3 genes). Among them, CAMTA1 and RPS6KA2, and CTBP2 are involved in regulatory T cell function and Th2 cell differentiation, respectively. Gene-functional analysis revealed two enriched gene clusters involved in immune responses and protein phosphorylation. ChIP-X Enrichment Analysis 3 showed eight significant transcription factors (RXRA, ZBTB7A, ESR1, TCF3, MYOD1, CTCF, GATA3 and CBX2). Ingenuity Pathway Analysis identified canonical pathways involved, among other, in B cell development, pathogen-induced cytokine storm signalling pathway and dendritic cell maturation. CONCLUSION: This review highlights the involvement of epigenomic alterations of loci in Th1/Th2 and regulatory T cell differentiation in both candidate gene studies and EWAS. These alterations provide a better insight into the mechanistic aspects in FA pathogenesis and may guide the development of epigenome-based biomarkers for FA.

Doubled-haploid induction generates extensive differential DNA methylation in Arabidopsis.

Doubled haploids (DHs) fix traits from hybrids in one generation. DH induction includes two changes in ploidy levels typically associated with variation in DNA methylation. However, DNA methylation patterns in DH plants and their biological significance are largely unknown. We generated three DH lines in Arabidopsis thaliana by crossing a haploid inducer with the accession Col-0, thus removing tissue culture and hybridization as a variable. DH induction produced thousands of differentially DNA methylated regions (DMRs), most of which were stochastic. Both haploidization and colchicine-induced genome duplication produced DMRs; the former mainly yielded DMRs at non-CG contexts, whereas the latter affected differential gene body methylation. Spontaneous genome doubling of haploid plants also induced DMRs in greater numbers than self-propagation. Our results provide the first evidence that haploid induction and genome doubling result in differential DNA methylation, offering a novel approach to induce epialleles.

Knockout of TRDMT1 methyltransferase affects DNA methylome in glioblastoma cells.

PURPOSE: We have previously shown that TRDMT1 methyltransferase is a regulator of chemotherapy-associated responses in glioblastoma cells. Despite the fact that glioblastoma, a common and malignant brain tumor, is widely characterized in terms of genetic and epigenetic markers, there are no data on TRDMT1-related changes in 5-methylcytosine pools in the genome. In the present study, the effect of TRDMT1 gene knockout (KO) on DNA methylome was analyzed. METHODS: CRISPR-based approach was used to obtain TRDMT1 KO glioblastoma cells. Total 5-methylcytosine levels in DNA, DNMT1 pools and DNMT activity were studied using ELISA. Reduced representation bisulfite sequencing (RRBS) was considered to comprehensively evaluate DNA methylome in glioblastoma cells with TRDMT1 KO. RESULTS: TRDMT1 KO cells were characterized by decreased levels of total 5-methylcytosine in DNA and DNMT1, and DNMT activity. RRBS-based methylome analysis revealed statistically significant differences in methylation-relevant DMS-linked genes in control cells compared to TRDMT1 KO cells. TRDMT1 KO-associated changes in DNA methylome may affect the activity of several processes and pathways such as telomere maintenance, cell cycle and longevity regulating pathway, proteostasis, DNA and RNA biology. CONCLUSIONS: TRDMT1 may be suggested as a novel modulator of gene expression by changes in DNA methylome that may affect cancer cell fates during chemotherapy. We postulate that the levels and mutation status of TRDMT1 should be considered as a prognostic marker and carefully monitored during glioblastoma progression.

MeDIP-sequencing for profiling global DNA methylation in buffalo embryos produced by in vitro fertilization.

Assisted reproductive technique like in vitro fertilization has contributed immensely in producing genetically improved livestock. Production of embryos under in vitro conditions can affect global DNA methylation pattern during the course of embryonic development. The present study is aimed at the generation and comparison of global DNA methylome of embryos at 2-cell, 8-cell and blastocyst stage of buffalo embryos produced by in vitro fertilization using MeDIP-Sequencing. It is observed that there is a profound difference in the global DNA methylation profile of IVF embryos at different developmental stages. These differences are manifested throughout the course of embryonic development. Pathways like Wnt signaling pathway, gonadotropin-releasing hormone receptor pathway and integrin signaling were found to be majorly affected by hypermethylation of DNA in IVF embryos throughout the development.

Integrated Analysis of DNA Methylome and Transcriptome Reveals Epigenetic Regulation of Cold Tolerance in Litopenaeus vannamei.

DNA methylation is an important epigenetic modification that has been shown to be associated with responses to non-biological stressors. However, there is currently no research on DNA methylation in response to environmental signals in shrimp. In this study, we conducted a comprehensive comparative analysis of DNA methylation profiles and differentially expressed genes between two strains of Litopenaeus vannamei with significantly different cold tolerance through whole genome bisulfite sequencing (WGBS) and transcriptome sequencing. Between Lv-C and Lv-T (constant temperature of 28 °C and low temperatures of 18 °C and 10 °C) under cytosine-guanine (CG) environments, 39,100 differentially methylated regions (DMRs) were identified, corresponding to 9302 DMR-related genes (DMRGs). The DMRs were mainly located in the gene body (exons and introns). Gene Ontology (GO) analysis showed that these DMRGs were significantly enriched in cell parts, catalytic activity, and metabolic processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed significant enrichment of these DMRGs in pathways such as proteasome (ko03050), oxidative phosphorylation (ko00190), mTOR signaling pathway (ko04150), fatty acid metabolism (ko01212), and fatty acid degradation (ko00071). The comprehensive results suggested that L. vannamei mainly regulates gene expression in response to low temperatures through hypermethylation or demethylation of some genes involved in thermogenesis, glycolysis, the autophagy pathway, the peroxisome, and drug metabolism pathways. These results provide important clues for studying DNA methylation patterns and identifying cold tolerance genes in shrimp.

Exploring crucial molecular events in pearl oyster after pre-grafting conditioning by genome-wide bisulfite sequencing for DNA methylation analysis.

Pre-grafting condition is an important method to promote recovery from transplant surgery during pearl production. In the present study, we constructed two DNA methylomes from pearl oysters with and without conditioning to investigate the molecular mechanism of the pearl oyster Pinctada fucata martensii underlying the pre-grafting condition. A total of 4,594,997 and 4,930,813 methyl CG in the control (Con) and pre-grafting group (PT) were detected, resulting in the whole genome methylation profile and methylation pattern in P. f. martensii. Results reveal that the promoter, especially the CpG island-rich region, was more infrequently methylated than the gene function elements in P. f. martensii. A total of 51,957 differently methylated regions (DMRs) between Con and PT were obtained, including 3789 DMR in the promoter and 16,021 in the gene body. Based on gene ontology and pathway enrichment analyses, these DMRs were mainly related to "cellular process", "metabolic process", "Epstein-Barr virus infection", and "Fanconi anemia pathway". The methylation site in the promoter region may be associated with the promoter activity and transcription factor binding. These results help our understanding of the mechanism of pre-grafting condition, thereby providing key information in guiding to improve the conditioning methods for enhanced pearl oyster survival rate after transplantation.

A role of age-dependent DNA methylation reprogramming in regulating the regeneration capacity of Boea hygrometrica leaves.

Plants can regenerate new individuals under appropriate culture conditions. Although the molecular basis of shoot regeneration has steadily been unraveled, the role of age-dependent DNA methylation status in the regulation of explant regeneration remains practically unknown. Here, we established an effective auxin/cytokinin-induced shoot regeneration system for the resurrection plant Boea hygrometrica via direct organogenesis and observed that regeneration was postponed with increasing age of donor plants. Global transcriptome analysis revealed significant upregulation of genes required for hormone signaling and phenylpropanoid biosynthesis and downregulation of photosynthetic genes during regeneration. Transcriptional changes in the positive/negative regulators and cell wall-related proteins involved in plant regeneration, such as ELONGATED HYPOCOTYL5 (HY5), LATERAL ORGAN BOUNDARIES DOMAIN, SHOOT-MERISTEMLESS, and WUSCHEL, were associated with the regeneration process. Comparison of DNA methylation profiling between leaves from young seedlings (YL) and mature plants (ML) revealed increased asymmetrical methylation in ML, which was predominantly distributed in promoter regions of genes, such as HY5 and a member of ABA-responsive element (ABRE) binding protein/ABRE binding factor, as well as genes encoding glycine-rich cell wall structural protein, CENTRORADIALIS-like protein, and beta-glucosidase 40-like essential for shoot meristem and cell wall architecture. Their opposite transcription response in ML explants during regeneration compared with those from YL demonstrated the putative involvement of DNA methylation in regeneration. Moreover, a significant lower expression of DNA glycosylase-lyase required for DNA demethylation in ML was coincident with its postponed regeneration compared with those in YL. Taken together, our results suggest a role of promoter demethylation in B. hygrometrica regeneration.

Identification of DNA methylation regulated novel host genes relevant to inhibition of virus replication in porcine PK15 cell using double stranded RNA mimics and DNA methyltransferase inhibitor.

During RNA viruses's replication, double-stranded RNA (dsRNA) is normally produced and induce host innate immune response. Most of gene activation due cytokine mediated but which are due to methylation mediated is still unknown. In the study, DNA methylome was integrated with our previous transcriptome data to investigate the differentially methylated regions and genes using MeDIP-chip technology. We found that the transcriptional expressions of 15, 37 and 18 genes were negatively related with their promoter DNA methylation levels in the cells treated by PolyI:C, Aza-CdR, as well as PolyI:C plus Aza-CdR, respectively, compared with the untreated cells. GO analysis revealed hypo-methylated genes (BNIP3L and CDK9) and a hyper-methylated gene (ZC3HAV1) involved in the host response to viral replication. Our results suggest that these novel genes targeted by DNA methylation can be potential markers relevant to virus replication and host innate immune response to set up a medical model of infectious diseases.

Insights from multidimensional analyses of the pan-cancer DNA methylome heterogeneity and the uncanonical CpG-gene associations.

Although the DNA methylome profiles have been available in large cancer cohorts such as The Cancer Genome Atlas (TCGA), integrative analysis of the DNA methylome architectures in a pan-cancer manner remains limited. In the present study, we aimed to systematically dissect the insightful features related to the inter-tumoral DNA methylome heterogeneity in a pan-cancer context of 21 cancers in TCGA. First, pan-cancer clustering of the DNA methylomes revealed convergence of cancers and, meanwhile, new classifications of cancer subtypes, which are often associated to prognostic differences. Next, within each type of cancer, we showed that the transcription factor (TF) genes tend to bear more dynamic promoter DNA methylation profiles than the other genes, which serves as a potential source of the transcriptome heterogeneity in cancers. Finally, we found unanticipated significant numbers of the non-canonical promoter CpG sites that are positively correlated with the gene expression. Distribution patterns of these CpG sites in the CpG islands, ChIP-seq, DNaseI-seq, PMD regions and histone modification landscapes suggested against a pervasive mechanism of transcriptional activation due to mCpG-dependent binding of TFs, which is not in complete agreement with previous hypothesis. In summary, our deep mining of the highly heterogeneous DNA methylome data in a pan-cancer context generated novel insights into the architecture of cancer epigenetics and provided a series of resources for further investigations in the related fields of cancer genomics and epigenetics.

Plasticity of DNA methylation and gene expression under zinc deficiency in Arabidopsis roots.

DNA methylation is a heritable chromatin modification that maintains chromosome stability, regulates transposon silencing and appears to be involved in gene expression in response to environmental conditions. Environmental stress alters DNA methylation patterns that are correlated with gene expression differences. Here, genome-wide differential DNA methylation was identified upon prolonged zinc (Zn) deficiency, leading to hypo- and hypermethylated chromosomal regions. Preferential CpG methylation changes occurred in gene promoters and gene bodies, but did not overlap with transcriptional start sites. Methylation changes were also prominent in transposable elements. In contrast, non-CpG methylation differences were exclusively found in promoters of protein-coding genes and in transposable elements. Strongly Zn deficiency-induced genes and their promoters were mostly non-methylated, irrespective of Zn supply. Differential DNA methylation in the CpG and CHG, but not in the CHH context, was found close to a few up-regulated Zn deficiency genes. However, the transcriptional Zn deficiency response in roots appeared little correlated with associated DNA methylation changes in promoters or gene bodies. Furthermore, under Zn deficiency, developmental defects were identified in an Arabidopsis mutant lacking non-CpG methylation. The root methylome thus responds specifically to a micronutrient deficiency and is important for efficient Zn utilization at low availability, but the relationship of differential methylation and differentially expressed genes is surprisingly poor.

DNA methylome profiling at single-base resolution through bisulfite sequencing of 5mC-immunoprecipitated DNA.

BACKGROUND: Detection of DNA methylome at single-base resolution is a significant challenge but promises to shed considerable light on human disease etiology. Current technologies could not detect DNA methylation genome-wide at single-base resolution with small amount of sequencing data and could not avoid detecting the methylation of repetitive elements which are considered as "junk DNA". METHODS: In this study, we have developed a novel DNA methylome profiling technology named MB-seq with its ability to identify genome-wide 5mC and quantify DNA methylation levels by introduced an assistant adapter AluI-linker This linker can be ligated to sonicated DNA and then be digested after the bisulfite treatment and amplification, which has no effect of MeDIP enrichment. Because many researchers are interested in investigating the methylation of functional regions such as promoters and gene bodies, we have also developed a novel alternative method named MRB-seq, which can be used to investigate the DNA methylation of functional regions by removing the repeats with Cot-1 DNA. RESULTS: In this study, we have developed MB-seq, a novel DNA methylome profiling technology combining MeDIP-seq with bisulfite conversion, which can precisely detect the 5mC sites and determine their DNA methylation level at single-base resolution in a cost-effective way. In addition, we have developed a new alternative method, MRB-seq (MeDIP-repetitive elements removal-bisulfite sequencing), which interrogates 5mCs in functional regions by depleting nearly half of repeat fragments enriched by MeDIP. Comparing MB-seq and MRB-seq to whole-genome BS-seq using the same batch of DNA from YH peripheral blood mononuclear cells. We found that the sequencing data of MB-seq and MRB-seq almost reaches saturation after generating 7-8 Gbp data, whereas BS-seq requires about 100 Gbp data to achieve the same effect. In comparison to MeDIP-seq and BS-seq, MB-seq offers several key advantages, including single-base resolution, discriminating the methylated sites within a CpG and non-CpG pattern and overcoming the false positive of MeDIP-seq due to the non-specific binding of 5-methylcytidine antibody to genomic fragments. CONCLUSION: Our novel developed method MB-seq can accelerate the decoding process of DNA methylation mechanism in human diseases because it requires 7-8 Gbp data to measure human methylome with enough coverage and sequencing depth, affording it a direct and practical application in the study of multiple samples. In addition, we have also provided a novel alternative MRB-seq method, which removes most repetitive sequences and allows researchers to genome-wide characterize DNA methylation of functional regions.