CpG islands in vertebrate genomes

Molecular testing of rhabdomyosarcoma in clinical trials to improve risk stratification and outcome: A consensus view from European paediatric Soft tissue sarcoma Study Group, Children's Oncology Group and Cooperative Weichteilsarkom-Studiengruppe

Rhabdomyosarcomas (RMSs) are the most common soft tissue sarcomas in children/adolescents less than 18 years of age with an annual incidence of 1-2/million. Inter/intra-tumour heterogeneity raise challenges in clinical, pathological and biological research studies. Risk stratification in European and North American clinical trials previously relied on clinico-pathological features, but now, incorporates PAX3/7-FOXO1-fusion gene status in the place of alveolar histology. International working groups propose a coordinated approach through the INternational Soft Tissue SaRcoma ConsorTium to evaluate the specific genetic abnormalities and generate and integrate molecular and clinical data related to patients with RMS across different trial settings. We review relevant data and present a consensus view on what molecular features should be assessed. In particular, we recommend the assessment of the MYOD1-LR122R mutation for risk escalation, as it has been associated with poor outcomes in spindle/sclerosing RMS and rare RMS with classic embryonal histopathology. The prospective analyses of rare fusion genes beyond PAX3/7-FOXO1 will generate new data linked to outcomes and assessment of TP53 mutations and CDK4 amplification may confirm their prognostic value. Pathogenic/likely pathogenic germline variants in TP53 and other cancer predisposition genes should also be assessed. DNA/RNA profiling of tumours at diagnosis/relapse and serial analyses of plasma samples is recommended where possible to validate potential molecular biomarkers, identify new biomarkers and assess how liquid biopsy analyses can have the greatest benefit. Together with the development of new molecularly-derived therapeutic strategies that we review, a synchronised international approach is expected to enhance progress towards improved treatment assignment, management and outcomes for patients with RMS.

Impact of the Exposome on the Epigenome in Inflammatory Bowel Disease Patients and Animal Models

Inflammatory bowel diseases (IBD) are chronic inflammatory disorders of the gastrointestinal tract that encompass two main phenotypes, namely Crohn's disease and ulcerative colitis. These conditions occur in genetically predisposed individuals in response to environmental factors. Epigenetics, acting by DNA methylation, post-translational histones modifications or by non-coding RNAs, could explain how the exposome (or all environmental influences over the life course, from conception to death) could influence the gene expression to contribute to intestinal inflammation. We performed a scoping search using Medline to identify all the elements of the exposome that may play a role in intestinal inflammation through epigenetic modifications, as well as the underlying mechanisms. The environmental factors epigenetically influencing the occurrence of intestinal inflammation are the maternal lifestyle (mainly diet, the occurrence of infection during pregnancy and smoking); breastfeeding; microbiota; diet (including a low-fiber diet, high-fat diet and deficiency in micronutrients); smoking habits, vitamin D and drugs (e.g., IBD treatments, antibiotics and probiotics). Influenced by both microbiota and diet, short-chain fatty acids are gut microbiota-derived metabolites resulting from the anaerobic fermentation of non-digestible dietary fibers, playing an epigenetically mediated role in the integrity of the epithelial barrier and in the defense against invading microorganisms. Although the impact of some environmental factors has been identified, the exposome-induced epimutations in IBD remain a largely underexplored field. How these environmental exposures induce epigenetic modifications (in terms of duration, frequency and the timing at which they occur) and how other environmental factors associated with IBD modulate epigenetics deserve to be further investigated.

Pan-Cancer Analysis Reveals the Prognostic Potential of the THAP9/THAP9-AS1 Sense-Antisense Gene Pair in Human Cancers

Human THAP9, which encodes a domesticated transposase of unknown function, and lncRNA THAP9-AS1 (THAP9-antisense1) are arranged head-to-head on opposite DNA strands, forming a sense and antisense gene pair. We predict that there is a bidirectional promoter that potentially regulates the expression of THAP9 and THAP9-AS1. Although both THAP9 and THAP9-AS1 are reported to be involved in various cancers, their correlative roles on each other's expression has not been explored. We analyzed the expression levels, prognosis, and predicted biological functions of the two genes across different cancer datasets (TCGA, GTEx). We observed that although the expression levels of the two genes, THAP9 and THAP9-AS1, varied in different tumors, the expression of the gene pair was strongly correlated with patient prognosis; higher expression of the gene pair was usually linked to poor overall and disease-free survival. Thus, THAP9 and THAP9-AS1 may serve as potential clinical biomarkers of tumor prognosis. Further, we performed a gene co-expression analysis (using WGCNA) followed by a differential gene correlation analysis (DGCA) across 22 cancers to identify genes that share the expression pattern of THAP9 and THAP9-AS1. Interestingly, in both normal and cancer samples, THAP9 and THAP9-AS1 often co-express; moreover, their expression is positively correlated in each cancer type, suggesting the coordinated regulation of this H2H gene pair.

Structure and Methylation of 35S rDNA in Allopolyploids Anemone multifida (2n = 4x = 32, BBDD) and Anemone baldensis (2n = 6x = 48, AABBDD) and Their Parental Species Show Evidence of Nucleolar Dominance

Transcriptional silencing of 35S rDNA loci inherited from one parental species is occurring relatively frequently in allopolyploids. However, molecular mechanisms by which it is selected for transcriptional silencing remain unclear. We applied NGS, silver staining and bisulfite sequencing to study the structure, expression and methylation landscape of 35S rDNA in two allopolyploids of common origin, allotetraploid Anemone multifida (2n = 4x = 32, genome composition BBDD) and allohexaploid A. baldensis (2n = 6x = 48, AABBDD), and their genome donors, A. sylvestris (2n = 16, AA), A. cylindrica (2n = 16, BB) and A. parviflora (2n = 16, DD). The size of the recovered 35S rDNA units varied from 10,489 bp in A. cylindrica to 12,084 bp in A. sylvestris. Anemone showed an organization typical of most ribosomal 35S rDNA composed of NTS, ETS, rRNA genes, TTS and TIS with structural features of plant IGS sequences and all functional elements needed for rRNA gene activity. The NTS was more variable than the ETS and consisted of SRs which are highly variable among Anemone. Five to six CpG-rich islands were found within the ETS. CpG island located adjacent to the transcription initiation site (TIS) was highly variable regarding the sequence size and methylation level and exhibited in most of the species lower levels of methylation than CpG islands located adjacent to the 18S rRNA gene. Our results uncover hypomethylation of A. sylvestris- and A. parviflora-derived 35S rDNA units in allopolyploids A. multifida and A. baldensis. Hypomethylation of A. parviflora-derived 35S rDNA was more prominent in A. baldensis than in A. multifida. We showed that A. baldensis underwent coupled A. sylvestris-derived 35S rDNA array expansion and A. parviflora-derived 35S rDNA copy number decrease that was accompanied by lower methylation level of A. sylvestris-derived 35S rDNA units in comparison to A. parviflora-derived 35S rDNA units. These observations suggest that in A. baldensis nucleolar dominance is directed toward A. sylvestris-derived chromosomes. This work broadens our current knowledge of the 35S rDNA organization in Anemone and provides evidence of the progenitor-specific 35S rDNA methylation in nucleolar dominance.

Chromatin Methylation Abnormalities in Autosomal Dominant Polycystic Kidney Disease

Autosomal dominant polycystic kidney disease (ADPKD) is the most common inherited kidney disease worldwide and is one of the major causes of end-stage renal disease. PKD1 and PKD2 are two genes that mainly contribute to the development and progression of ADPKD. The precise mechanism is not fully understood. In recent years, epigenetic modification has drawn increasing attention. Chromatin methylation is a very important category of PKD epigenetic changes and mostly involves DNA, histone, and RNA methylation. Genome hypomethylation and regional gene hypermethylation coexist in ADPKD. We found that the genomic DNA of ADPKD kidney tissues showed extensive demethylation by whole-genome bisulphite sequencing, while some regional DNA methylation from body fluids, such as blood and urine, can be used as diagnostic or prognostic biomarkers to predict PKD progression. Histone modifications construct the histone code mediated by histone methyltransferases and contribute to aberrant methylation changes in PKD. Considering the complexity of methylation abnormalities occurring in different regions and genes on the PKD epigenome, more specific therapy aiming to restore to the normal genome should lead to the development of epigenetic treatment.

Comparison of DNA methylation measurements from EPIC BeadChip and SeqCap targeted bisulphite sequencing in PON1 and nine additional candidate genes

Epigenome-wide association studies (EWAS) are widely implemented in epidemiology, and the Illumina HumanMethylationEPIC BeadChip (EPIC) DNA microarray is the most-used technology. Recently, next-generation sequencing (NGS)-based methods, which assess DNA methylation at single-base resolution, have become more affordable and technically feasible. While the content of microarray technology is fixed, NGS-based approaches, such as the Roche Nimblegen, SeqCap Epi Enrichment System (SeqCap), offer the flexibility of targeting most CpGs in a gene. With the current usage of microarrays and emerging NGS-based technologies, it is important to establish whether data generated from the two platforms are comparable. We harnessed 112 samples from the Center for the Health Assessment of Mothers and Children of Salinas (CHAMACOS) birth cohort study and compared DNA methylation between the EPIC microarray and SeqCap for PON1 and nine additional candidate genes, by evaluating epigenomic coverage and correlations. We conducted multivariable linear regression and principal component analyses to assess the ability of the EPIC array and SeqCap to detect biological differences in gene methylation by the PON1_-108 single nucleotide polymorphism. We found an overall high concordance (r = 0.84) between SeqCap and EPIC DNA methylation, among highly methylated and minimally methylated regions. However, substantial disagreement was present between the two methods in moderately methylated regions, with SeqCap measurements exhibiting greater within-site variation. Additionally, SeqCap did not capture PON1 SNP associated differences in DNA methylation that were evident with the EPIC array. Our findings indicate that microarrays perform well for analysing DNA methylation in large cohort studies but with limited coverage.

Missense variants causing Wiedemann-Steiner syndrome preferentially occur in the KMT2A-CXXC domain and are accurately classified using AlphaFold2

Wiedemann-Steiner syndrome (WDSTS) is a neurodevelopmental disorder caused by de novo variants in KMT2A, which encodes a multi-domain histone methyltransferase. To gain insight into the currently unknown pathogenesis of WDSTS, we examined the spatial distribution of likely WDSTS-causing variants across the 15 different domains of KMT2A. Compared to variants in healthy controls, WDSTS variants exhibit a 61.9-fold overrepresentation within the CXXC domain–which mediates binding to unmethylated CpGs–suggesting a major role for this domain in mediating the phenotype. In contrast, we find no significant overrepresentation within the catalytic SET domain. Corroborating these results, we find that hippocampal neurons from Kmt2a-deficient mice demonstrate disrupted histone methylation (H3K4me1 and H3K4me3) preferentially at CpG-rich regions, but this has no systematic impact on gene expression. Motivated by these results, we combine accurate prediction of the CXXC domain structure by AlphaFold2 with prior biological knowledge to develop a classification scheme for missense variants in the CXXC domain. Our classifier achieved 92.6% positive and 92.9% negative predictive value on a hold-out test set. This classification performance enabled us to subsequently perform an in silico saturation mutagenesis and classify a total of 445 variants according to their functional effects. Our results yield a novel insight into the mechanistic basis of WDSTS and provide an example of how AlphaFold2 can contribute to the in silico characterization of variant effects with very high accuracy, suggesting a paradigm potentially applicable to many other Mendelian disorders.

Wiedemann-Steiner syndrome (WDSTS) is a neurodevelopmental pediatric disorder caused by the genetic disruption of the histone methyltransferase KMT2A. Since KMT2A has many different domains that perform different functions, we reasoned that by identifying the domains most enriched for WDSTS-causing genetic variants we would gain insights into the incompletely understood molecular pathogenesis of WDSTS. We discovered that the CXXC domain—which binds unmethylated CpGs—shows by far the greatest enrichment, suggesting that loss of the CpG-binding ability of KMT2A plays a central role in WDSTS. Next, to understand specific rules underlying the genetic disruption of the CXXC domain, we combined prior knowledge about the function/structure of the domain with 3D structure prediction by AlphaFold2 to develop an effect classifier for CXXC missense variants. We found that this classifier exhibits accurate performance, and we therefore applied it to provide classifications for any such variant that can possibly arise, in order to aid in the interpretation of such variants in the clinic. Our work provides novel insights into WDSTS and suggests a strategy for missense variant classification that can potentially be applied to many other pediatric genetic disorders.

Dynamic changes in hepatic DNA methylation during the development of nonalcoholic fatty liver disease induced by a high-sugar diet

DNA methylation is an important epigenetic mechanism of gene expression control. The present study aimed to evaluate the temporal effect of isocaloric high-sugar diet (HSD) intake on the development of nonalcoholic fatty liver disease (NAFLD) and the role of DNA methylation in this event. Newly weaned Wistar rats were divided into eight groups and fed a standard chow diet or an HSD ad libitum for 4 weeks, 8 weeks, 15 weeks, and 18 weeks. After the experimental periods, the animals were euthanized and their livers were removed for histological analysis, gene expression of maintenance methylase (Dnmt1), de novo methylases (Dnmt3a and Dnmt3b), demethylases (Tet2 and Tet3) of DNA, and global DNA methylation. HSD intake led to the gradual development of NAFLD. HSD intake for 18 weeks was associated with downregulation of Dnmt1 expression and global DNA hypomethylation; these results were negatively correlated with more severe steatosis scores observed in these animals. The HSD consumption for 18 weeks was also associated with a decrease in Dnmt3a and Tet2 expression. Interestingly, the expression of de novo methyltransferase Dnmt3b was reduced by HSD during all experimental periods. Together, these results indicate that the downregulation of de novo DNA methylation, Dnmt3b, induced by HSD is the primary factor in the development of NAFLD. On the other hand, disease progression is associated with downregulation of maintenance DNA methylation and global DNA hypomethylation. These results suggest a link between the dynamic changes in hepatic DNA methylation and the development of NAFLD induced by an HSD intake.

An Integrative Analysis of DNA Methylation Pattern in Myotonic Dystrophy Type 1 Samples Reveals a Distinct DNA Methylation Profile between Tissues and a Novel Muscle-Associated Epigenetic Dysregulation

Myotonic dystrophy type 1 (DM1) is a progressive, non-treatable, multi-systemic disorder. To investigate the contribution of epigenetics to the complexity of DM1, we compared DNA methylation profiles of four annotated CpG islands (CpGis) in the DMPK locus and neighbouring genes, in distinct DM1 tissues and derived cells, representing six DM1 subtypes, by bisulphite sequencing. In blood, we found no differences in CpGi 74, 43 and 36 in DNA methylation profile. In contrast, a CTCF1 DNA methylation gradient was found with 100% methylation in congenital cases, 50% in childhood cases and 13% in juvenile cases. CTCF1 methylation correlated to disease severity and CTG expansion size. Notably, 50% of CTCF1 methylated cases showed methylation in the CTCF2 regions. Additionally, methylation was associated with maternal transmission. Interestingly, the evaluation of seven families showed that unmethylated mothers passed on an expansion of the CTG repeat, whereas the methylated mothers transmitted a contraction. The analysis of patient-derived cells showed that DNA methylation profiles were highly preserved, validating their use as faithful DM1 cellular models. Importantly, the comparison of DNA methylation levels of distinct DM1 tissues revealed a novel muscle-specific epigenetic signature with methylation of the CTCF1 region accompanied by demethylation of CpGi 43, a region containing an alternative DMPK promoter, which may decrease the canonical promoter activity. Altogether, our results showed a distinct DNA methylation profile across DM1 tissues and uncovered a novel and dual epigenetic signature in DM1 muscle samples, providing novel insights into the epigenetic changes associated with DM1.

DNA methylation-independent long-term epigenetic silencing with dCRISPR/Cas9 fusion proteins

The programmable CRISPR/Cas9 DNA nuclease is a versatile genome editing tool, but it requires the host cell DNA repair machinery to alter genomic sequences. This fact leads to unpredictable changes of the genome at the cut sites. Genome editing tools that can alter the genome without causing DNA double-strand breaks are therefore in high demand. Here, we show that expression of promoter-associated short guide (sg)RNAs together with dead Cas9 (dCas9) fused to a Krüppel-associated box domains (KRABd) in combination with the transcription repression domain of methyl CpG-binding protein 2 (MeCP2) can lead to persistent gene silencing in mouse embryonic stem cells and in human embryonic kidney (HEK) 293 cells. Surprisingly, this effect is achievable and even enhanced in DNA (cytosine-5)-methyltransferase 3A and 3B (Dnmt3A-/-, Dnmt3b-/-) depleted cells. Our results suggest that dCas9-KRABd-MeCP2 fusions are useful for long-term epigenetic gene silencing with utility in cell biology and potentially in therapeutical settings.

LARGE-SCALE INTEGRATION OF DNA METHYLATION AND GENE EXPRESSION ARRAY PLATFORMS IDENTIFIES BOTH cis AND trans RELATIONSHIPS

Although epigenome-wide association studies (EWAS) have been successful in identifying DNA methylation (DNAm) patterns associated with disease states, any further characterization of etiologic mechanisms underlying disease remains elusive. This knowledge gap does not originate from a lack of DNAm–trait associations, but rather stems from study design issues that affect the interpretability of EWAS results. Despite known limitations in predicting the function of a particular CpG site, most EWAS maintain the broad assumption that altered DNAm results in a concomitant change of transcription at the most proximal gene. This study integrated DNAm and gene expression (GE) measurements in two cohorts, the Adolescent and Young Adult Twin Study (AYATS) and the Pregnancy, Race, Environment, Genes (PREG) study, to improve the understanding of epigenomic regulatory mechanisms. CpG sites associated with GE in cis were enriched in areas of transcription factor binding and areas of intermediate-to-low CpG density. CpG sites associated with trans GE were also enriched in areas of known regulatory significance, including enhancer regions. These results highlight issues with restricting DNAm-transcript annotations to small genomic intervals and question the validity of assuming a cis DNAm–GE pathway. Based on these findings, the interpretation of EWAS results is limited in studies without multi-omic support and further research should identify genomic regions in which GE-associated DNAm is overrepresented. An in-depth characterization of GE-associated CpG sites could improve predictions of the downstream functional impact of altered DNAm and inform best practices for interpreting DNAm–trait associations generated by EWAS.

Thorough statistical analyses of breast cancer co-methylation patterns

Background

Breast cancer is one of the most commonly diagnosed cancers. It is associated with DNA methylation, an epigenetic event with a methyl group added to a cytosine paired with a guanine, i.e., a CG site. The methylation levels of different genes in a genome are correlated in certain ways that affect gene functions. This correlation pattern is known as co-methylation. It is still not clear how different genes co-methylate in the whole genome of breast cancer samples. Previous studies are conducted using relatively small datasets (Illumina 27K data). In this study, we analyze much larger datasets (Illumina 450K data).

Results

Our key findings are summarized below. First, normal samples have more highly correlated, or co-methylated, CG pairs than tumor samples. Both tumor and normal samples have more than 93% positive co-methylation, but normal samples have significantly more negatively correlated CG sites than tumor samples (6.6% vs. 2.8%). Second, both tumor and normal samples have about 94% of co-methylated CG pairs on different chromosomes, but normal samples have 470 million more CG pairs. Highly co-methylated pairs on the same chromosome tend to be close to each other. Third, a small proportion of CG sites’ co-methylation patterns change dramatically from normal to tumor. The percentage of differentially methylated (DM) sites among them is larger than the overall DM rate. Fourth, certain CG sites are highly correlated with many CG sites. The top 100 of such super-connector CG sites in tumor and normal samples have no overlaps. Fifth, both highly changing sites and super-connector sites’ locations are significantly different from the genome-wide CG sites’ locations. Sixth, chromosome X co-methylation patterns are very different from other chromosomes. Finally, the network analyses of genes associated with several sets of co-methylated CG sites identified above show that tumor and normal samples have different patterns.

Conclusions

Our findings will provide researchers with a new understanding of co-methylation patterns in breast cancer. Our ability to thoroughly analyze co-methylation of large datasets will allow researchers to study relationships and associations between different genes in breast cancer.

Superstructure Detection in Nucleosome Distribution Shows Common Pattern within a Chromosome and within the Genome

Nucleosome positioning plays an important role in crucial biological processes such as replication, transcription, and gene regulation. It has been widely used to predict the genome’s function and chromatin organisation. So far, the studies of patterns in nucleosome positioning have been limited to transcription start sites, CTCFs binding sites, and some promoter and loci regions. The genome-wide organisational pattern remains unknown. We have developed a theoretical model to coarse-grain nucleosome positioning data in order to obtain patterns in their distribution. Using hierarchical clustering on the auto-correlation function of this coarse-grained nucleosome positioning data, a genome-wide clustering is obtained for Candida albicans. The clustering shows the existence beyond hetero- and eu-chromatin inside the chromosomes. These non-trivial clusterings correspond to different nucleosome distributions and gene densities governing differential gene expression patterns. Moreover, these distribution patterns inside the chromosome appeared to be conserved throughout the genome and within species. The pipeline of the coarse grain nucleosome positioning sequence to identify underlying genomic organisation used in our study is novel, and the classifications obtained are unique and consistent.

Gene Promoter-Methylation Signature as Biomarker to Predict Cisplatin-Radiotherapy Sensitivity in Locally Advanced Cervical Cancer

Despite efforts to promote health policies focused on screening and early detection, cervical cancer continues to be one of the leading causes of mortality in women; in 2020, estimated 30,000 deaths in Latin America were reported for this type of tumor. While the therapies used to treat cervical cancer have excellent results in tumors identified in early stages, those women who are diagnosed in locally advanced and advanced stages show survival rates at 5 years of <50%. Molecular patterns associated with clinical response have been studied in patients who present resistance to treatment; none of them have reached clinical practice. It is therefore necessary to continue analyzing molecular patterns that allow us to identify patients at risk of developing resistance to conventional therapy. In this study, we analyzed the global methylation profile of 22 patients diagnosed with locally advanced cervical cancer and validated the genomic results in an independent cohort of 70 patients. We showed that BRD9 promoter region methylation and CTU1 demethylation were associated with a higher overall survival (p = 0.06) and progression-free survival (p = 0.0001), whereas DOCK8 demethylation was associated with therapy-resistant patients and a lower overall survival and progression-free survival (p = 0.025 and p = 0.0001, respectively). Our results suggest that methylation of promoter regions in specific genes may provide molecular markers associated with response to treatment in cancer; further investigation is needed.

Optical Imaging of Epigenetic Modifications in Cancer: A Systematic Review

Increasing evidence has demonstrated that abnormal epigenetic modifications are strongly related to cancer initiation. Thus, sensitive and specific detection of epigenetic modifications could markedly improve biological investigations and cancer precision medicine. A rapid development of molecular imaging approaches for the diagnosis and prognosis of cancer has been observed during the past few years. Various biomarkers unique to epigenetic modifications and targeted imaging probes have been characterized and used to discriminate cancer from healthy tissues, as well as evaluate therapeutic responses. In this study, we summarize the latest studies associated with optical molecular imaging of epigenetic modification targets, such as those involving DNA methylation, histone modification, noncoding RNA regulation, and chromosome remodeling, and further review their clinical application on cancer diagnosis and treatment. Lastly, we further propose the future directions for precision imaging of epigenetic modification in cancer. Supported by promising clinical and preclinical studies associated with optical molecular imaging technology and epigenetic drugs, the central role of epigenetics in cancer should be increasingly recognized and accepted.

WGCNA-Based DNA Methylation Profiling Analysis on Allopurinol-Induced Severe Cutaneous Adverse Reactions: A DNA Methylation Signature for Predisposing Drug Hypersensitivity

Background: The role of aberrant DNA methylation in allopurinol-induced severe cutaneous adverse reactions (SCARs) is incompletely understood. To fill the gap, we analyze the DNA methylation profiling in allopurinol-induced Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN) patients and identify the DNA methylation signature for predisposing allopurinol hypersensitivity. Methods: Genome-scale methylation analysis was conducted using the Illumina® HumanMethylation450 BeadChip. Weighted Gene Co-Expression Network Analysis (WGCNA) was utilized to analyze the data. Results: A total of 21,497 annotated promoter regions were analyzed. Ten modules were identified between allopurinol hypersensitivity and tolerance, with turquoise and yellow modules being the most significant correlation. ATG13, EPM2AIP1, and SRSF11 were the top three hub genes in the turquoise module. MIR412, MIR369, and MIR409 were the top three hub genes in the yellow module. Gene Ontology (GO) analysis revealed that the turquoise module was related to the metabolic process in intracellular organelles and the binding of various compounds, proteins, or nucleotides. The yellow module, however, was related to stimulus sensory perception in cytoskeletal elements and the activity of the receptor or transducer. Conclusion: DNA methylation plays a vital role in allopurinol-induced severe cutaneous adverse reactions (SCARs). DNA methylation profiling of SJS/TEN is significantly related to autophagy and microRNAs (miRNAs).

Aryl hydrocarbon receptor as a DNA methylation reader in the stress response pathway

The aryl hydrocarbon receptor (AhR) mediates various cellular responses upon exposure to exogenous and endogenous stress factors. In these responses, AhR plays a dual role as a stress sensor for detecting various AhR ligands and as a transcription factor that upregulates the expression of downstream effector genes, such as those encoding drug-metabolizing enzymes. As a transcription factor, it selectively binds to the unmethylated form of a specific sequence called the xenobiotic responsive element (XRE). We suggest that AhR is a novel DNA methylation reader, unlike classical methylation readers, such as methyl-CpG-binding protein 2, which binds to methylated sequences. Under physiological conditions of continuous exposure to endogenous AhR ligands, such as kynurenine, methylation states of the individual target XREs must be strictly regulated to select and coordinate the expression of downstream genes responsible for maintaining homeostasis in the body. In contrast, long-term exposure to AhR ligands frequently leads to changes in the methylation patterns around the XRE sequence. These data indicate that AhR may contribute to the adaptive cellular response to various stresses by modulating DNA methylation. Thus, the DNA methylation profile of AhR target genes should be dynamically controlled through a balance between robustness and flexibility under both physiological and stress conditions. AhR is a pivotal player in the regulation of stress response as it shows versatility by functioning as a stress sensor, methylation reader, and putative methylation modulator.

In silico identification of single nucleotide variations at CpG sites regulating CpG island existence and size

Genetic and epigenetic modifications of genes involved in the key regulatory pathways play a significant role in the pathophysiology and progression of multifactorial diseases. The present study is an attempt to identify single nucleotide variations (SNVs) at CpG sites of promoters of ACAT1, APOB, APOE, CYBA, FAS, FLT1, KSR2, LDLR, MMP9, PCSK9, PHOX2A, REST, SH2B3, SORT1 and TIMP1 genes influencing CpG island (CGI) existence and size associated with the pathophysiology of Diabetes mellitus, Coronary artery disease and Cancers. Promoter sequences located between -2000 to + 2000 bp were retrieved from the EPDnew database and predicted the CpG island using MethPrimer. Further, SNVs at CpG sites were accessed from NCBI, Ensembl while transcription factor (TF) binding sites were accessed using AliBaba2.1. CGI existence and size were determined for each SNV at CpG site with respect to wild type and variant allele by MethPrimer. A total of 200 SNVs at CpG sites were analyzed from the promoters of ACAT1, APOB, APOE, CYBA, FAS, FLT1, KSR2, LDLR, MMP9, PCSK9, PHOX2A, REST, SH2B3, SORT1 and TIMP1 genes. Of these, only 17 (8.5%) SNVs were found to influence the loss of CGI while 70 (35%) SNVs were found to reduce the size of CGI. It has also been found that 59% (10) of CGI abolishing SNVs are showing differences in binding of TFs. The findings of the study suggest that the candidate SNVs at CpG sites regulating CGI existence and size might influence the DNA methylation status and expression of genes involved in molecular pathways associated with several diseases. The insights of the present study may pave the way for new experimental studies to undertake challenges in DNA methylation, gene expression and protein assays.

GC-Profile 2.0: an extended web server for the prediction and visualization of CpG islands

SUMMARY

Due to the spontaneous deamination of 5'-methylcytosine into thymine, the number of CpG dinucleotides is less than expected in vertebrate genomes. Exceptionally, there are a large number of CpG dinucleotides clustered at certain genomic loci, known as CpG islands (CGIs), where CpG dinucleotides are free from methylation. Identification of CGIs is of great significance in the field of genomics and epigenetics because they can serve as important gene markers or regulatory elements. Here, GC-Profile 2.0 has been presented as a newly extended application for CGIs detection and visualization. Based on a benchmark test of assembled sequences, GC-Profile 2.0 has shown better overall performance compared with other four popular methods. In addition, cumulative CpG profile, a visualization tool of CpG content variation, is also proposed to intuitively display the change trend of CpG content.

AVAILABILITY AND IMPLEMENTATION

GC-Profile 2.0 is freely available at http://tubic.org/GC-Profile2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

TET2 Regulates 5-Hydroxymethylcytosine Signature and CD4+ T-Cell Balance in Allergic Rhinitis

Purpose

Previous studies have shown the role of ten-eleven translocation 2 (TET2) in CD4⁺ T cells. However, its function in CD4⁺ T cells under allergic inflammation is unclear. We aimed to investigate the epigenomic distribution of DNA 5-hydroxymethylcytosine (5hmC) and the role of TET2 in CD4⁺ T cells of allergic rhinitis (AR).

Methods

The hMeDIP-seq was performed to identify sequences with 5hmC deposition in CD4⁺ T cells of AR patients. Tet2-deficient or wild type mice were stimulated with ovalbumin (OVA) to develop an AR mouse model. The histopathology in nasal mucosae, Th1/Th2/Treg/Th17 cell percentage, concentrations of Th-related cytokines, expression of Tet and differential hydroxymethylated genes (DhMG), and the global deposition of 5hmC in sorted CD4⁺ T cells were detected.

Results

Epigenome-wide 5hmC landscape and DhMG in the CD4⁺ T cells of AR patients were identified. Tet2 depletion did not led to spontaneous inflammation. However, under the stimulation of allergen, OVA, loss of Tet2 resulted in the exacerbation of allergic inflammation, which was characterized by severer allergic symptoms, more inflammatory cells infiltrating the nasal lamina propria, sharper imbalances between Th1/Th2 and Treg/Th17 cells, and excessive secretion of OVA-specific IgE and Th2-related cytokines. Moreover, altered mRNA production of several DhMG and sharp decrease in 5hmC deposition were also observed in Tet2-deficient OVA-exposed mice.

Conclusions

TET2 may regulate DNA 5hmC, DhMG expressions, and CD4⁺ T cell balance in AR.

BCOR and BCORL1 Mutations Drive Epigenetic Reprogramming and Oncogenic Signaling by Unlinking PRC1.1 from Target Genes

Polycomb repressive epigenetic complexes are recurrently dysregulated in cancer. Unlike polycomb repressive complex 2 (PRC2), the role of PRC1 in oncogenesis and therapy resistance is not well-defined. Here, we demonstrate that highly recurrent mutations of the PRC1 subunits BCOR and BCORL1 in leukemia disrupt assembly of a noncanonical PRC1.1 complex, thereby selectively unlinking the RING-PCGF enzymatic core from the chromatin-targeting auxiliary subcomplex. As a result, BCOR-mutated PRC1.1 is localized to chromatin but lacks repressive activity, leading to epigenetic reprogramming and transcriptional activation at target loci. We define a set of functional targets that drive aberrant oncogenic signaling programs in PRC1.1-mutated cells and primary patient samples. Activation of these PRC1.1 targets in BCOR-mutated cells confers acquired resistance to treatment while sensitizing to targeted kinase inhibition. Our study thus reveals a novel epigenetic mechanism that explains PRC1.1 tumor-suppressive activity and identifies a therapeutic strategy in PRC1.1-mutated cancer.

SIGNIFICANCE

We demonstrate that BCOR and BCORL1 mutations in leukemia unlink PRC1.1 repressive function from target genes, resulting in epigenetic reprogramming and activation of aberrant cell signaling programs that mediate treatment resistance. Our study provides mechanistic insights into the pathogenesis of PRC1.1-mutated leukemia that inform novel therapeutic approaches. This article is highlighted in the In This Issue feature, p. 85.

Intragenic CpG Islands and Their Impact on Gene Regulation

The mammalian genome is depleted in CG dinucleotides, except at protected regions where they cluster as CpG islands (CGIs). CGIs are gene regulatory hubs and serve as transcription initiation sites and are as expected, associated with gene promoters. Advances in genomic annotations demonstrate that a quarter of CGIs are found within genes. Such intragenic regions are repressive environments, so it is surprising that CGIs reside here and even more surprising that some resist repression and are transcriptionally active within a gene. Hence, intragenic CGI positioning within genes is not arbitrary and is instead, selected for. As a wealth of recent studies demonstrate, intragenic CGIs are embedded within genes and consequently, influence ‘host’ gene mRNA isoform length and expand transcriptome diversity.

Laser Capture Microdissection in the Spatial Analysis of Epigenetic Modifications in Skin: A Comprehensive Review

Each cell in the body contains an intricate regulation for the expression of its relevant DNA. While every cell in a multicellular organism contains identical DNA, each tissue-specific cell expresses a different set of active genes. This organizational property exists in a paradigm that is largely controlled by forces external to the DNA sequence via epigenetic regulation. DNA methylation and chromatin modifications represent some of the classical epigenetic modifications that control gene expression. Complex tissues like skin consist of heterogeneous cell types that are spatially distributed and mixed. Furthermore, each individual skin cell has a unique response to physiological and pathological cues. As such, it is difficult to classify skin tissue as homogenous across all cell types and across different environmental exposures. Therefore, it would be prudent to isolate targeted tissue elements prior to any molecular analysis to avoid a possibility of confounding the sample with unwanted cell types. Laser capture microdissection (LCM) is a powerful technique used to isolate a targeted cell group with extreme microscopic precision. LCM presents itself as a solution to tackling the problem of tissue heterogeneity in molecular analysis. This review will cover an overview of LCM technology, the principals surrounding its application, and benefits of its application to the newly defined field of epigenomics, in particular of cutaneous pathology. This presents a comprehensive review about LCM and its use in the spatial analysis of skin epigenetics. Within the realm of skin pathology, this ability to isolate tissues under specific environmental stresses, such as oxidative stress, allows a far more focused investigation.

Potential differences in chitin synthesis ability cause different sensitivities to diflubenzuron among three strains of Daphnia magna

Ecotoxicity testing of crustaceans using Daphnia magna has been implemented in the chemical management systems of various countries. While the chemical sensitivity of D. magna varies depending on genetically different clonal lineages, the strain used in ecotoxicity tests, including the acute immobilization test (OECD TG202), has not been specified. We hypothesized that comprehensive gene expression profiles could provide useful information on phenotypic differences among strains, including chemical sensitivity. To test this hypothesis, we performed mRNA sequencing on three different strains (NIES, England, and Clone 5) of D. magna under culture conditions. The resulting expression profile of the NIES strain was clearly different compared to the profiles of the other two strains. Gene ontology (GO) enrichment analysis suggested that chitin metabolism was significantly enriched in the NIES strain compared to that in the England strain. Consistent with the GO analysis, evidence of high levels of chitin metabolism in the NIES strain were observed across multiple levels of biological organization, such as expression of chitin synthase genes, chitin content, and chitinase activity, which suggested that the different strains would exhibit different sensitivities to chemicals used to inhibit chitin synthesis. We found that among all strains, the NIES strain was more tolerant to diflubenzuron, a chitin synthesis inhibitor, with a 14-fold difference in the 48 h-EC₅₀ value for the acute immobilization test compared to the England strain. The present study demonstrates that the differences among strains in chitin metabolism may lead to sensitivity difference to diflubenzuron, and serves as a case study of the usefulness of comprehensive gene expression profiles in finding sensitivity differences.

Rlip Depletion Alters Oncogene Transcription at Multiple Distinct Regulatory Levels

Simple Summary

Rlip76 is a multifunctional membrane protein that facilitates cancer growth, and its depletion kills cancer cells. We recently found that Rlip depletion also results in broad changes to oncogene and tumor suppressor transcription. The present studies were designed to decipher the unknown downstream signaling pathways and transcriptional regulatory mechanisms driving the effect. Building on prior findings that Rlip depletion induces broad methylomic changes, we found using bioluminescence reporter assays that depletion of Rlip also exerts transcriptional control over several cancer genes through methylation-independent changes in transcription factor-mediated activation of their promoter regions and through additional as yet unidentified mechanisms. These findings have important implications for Rlip-targeted cancer therapy.

Abstract

Rlip76 (Rlip) is a multifunctional membrane protein that facilitates the high metabolic rates of cancer cells through the efflux of toxic metabolites and other functions. Rlip inhibition or depletion results in broad-spectrum anti-cancer effects in vitro and in vivo. Rlip depletion effectively suppresses malignancy and causes global reversion of characteristic CpG island methylomic and transcriptomic aberrations in the p53-null mouse model of spontaneous carcinogenesis through incompletely defined signaling and transcriptomic mechanisms. The methylome and transcriptome are normally regulated by the concerted actions of several mechanisms that include chromatin remodeling, promoter methylation, transcription factor interactions, and miRNAs. The present studies investigated the interaction of Rlip depletion or inhibition with the promoter methylation and transcription of selected cancer-related genes identified as being affected by Rlip depletion in our previous studies. We constructed novel promoter CpG island/luciferase reporter plasmids that respond only to CpG methylation and transcription factors. We found that Rlip depletion regulated expression by a transcription factor-based mechanism that functioned independently of promoter CpG methylation, lipid peroxidation, and p53 status.

Gastric cancer: An epigenetic view

Gastric cancer (GC) poses a serious threat worldwide with unfavorable prognosis mainly due to late diagnosis and limited therapies. Therefore, precise molecular classification and search for potential targets are required for diagnosis and treatment, as GC is complicated and heterogeneous in nature. Accumulating evidence indicates that epigenetics plays a vital role in gastric carcinogenesis and progression, including histone modifications, DNA methylation and non-coding RNAs. Epigenetic biomarkers and drugs are currently under intensive evaluations to ensure efficient clinical utility in GC. In this review, key epigenetic alterations and related functions and mechanisms are summarized in GC. We focus on integration of existing epigenetic findings in GC for the bench-to-bedside translation of some pivotal epigenetic alterations into clinical practice and also describe the vacant field waiting for investigation.

DNA Methylation Data-Based Classification and Identification of Prognostic Signature of Children With Wilms Tumor

Background: As an epigenetic alteration, DNA methylation plays an important role in early Wilms tumorigenesis and is possibly used as marker to improve the diagnosis and classification of tumor heterogeneity.

Methods: Methylation data, RNA-sequencing (RNA-seq) data, and corresponding clinical information were downloaded from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. The prognostic values of DNA methylation subtypes in Wilms tumor were identified.

Results: Four prognostic subtypes of Wilms tumor patients were identified by consensus cluster analysis performed on 312 independent prognostic CpG sites. Cluster one showed the best prognosis, whereas Cluster two represented the worst prognosis. Unique CpG sites identified in Cluster one that were not identified in other subtypes were assessed to construct a prognostic signature. The prognostic methylation risk score was closely related to prognosis, and the area under the curve (AUC) was 0.802. Furthermore, the risk score based on prognostic signature was identified as an independent prognostic factor for Wilms tumor in univariate and multivariate Cox regression analyses. Finally, the abundance of B cell infiltration was higher in the low-risk group than in the high-risk group, based on the methylation data.

Conclusion: Collectively, we divided Wilms tumor cases into four prognostic subtypes, which could efficiently identify high-risk Wilms tumor patients. Prognostic methylation risk scores that were significantly associated with the adverse clinical outcomes were established, and this prognostic signature was able to predict the prognosis of Wilms tumor in children, which may be useful in guiding clinicians in therapeutic decision-making. Further independent studies are needed to validate and advance this hypothesis.

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and −0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNA^Ser(AGC) and tRNA^Lys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

Genomic studies controvert the existence of the CUX1 p75 isoform

CUX1, encoding a homeodomain-containing transcription factor, is recurrently deleted or mutated in multiple tumor types. In myeloid neoplasms, CUX1 deletion or mutation carries a poor prognosis. We have previously established that CUX1 functions as a tumor suppressor in hematopoietic cells across multiple organisms. Others, however, have described oncogenic functions of CUX1 in solid tumors, often attributed to truncated CUX1 isoforms, p75 and p110, generated by an alternative transcriptional start site or post-translational cleavage, respectively. Given the clinical relevance, it is imperative to clarify these discrepant activities. Herein, we sought to determine the CUX1 isoforms expressed in hematopoietic cells and find that they express the full-length p200 isoform. Through the course of this analysis, we found no evidence of the p75 alternative transcript in any cell type examined. Using an array of orthogonal approaches, including biochemistry, proteomics, CRISPR/Cas9 genomic editing, and analysis of functional genomics datasets across a spectrum of normal and malignant tissue types, we found no data to support the existence of the CUX1 p75 isoform as previously described. Based on these results, prior studies of p75 require reevaluation, including the interpretation of oncogenic roles attributed to CUX1.

Methylation status of hypothalamic Mkrn3 promoter across puberty

Makorin RING finger protein 3 (MKRN3) is an important factor located on chromosome 15 in the imprinting region associated with Prader-Willi syndrome. Imprinted MKRN3 is expressed in hypothalamic regions essential for the onset of puberty and mutations in the gene have been found in patients with central precocious puberty. The pubertal process is largely controlled by epigenetic mechanisms that include, among other things, DNA methylation at CpG dinucleotides of puberty-related genes. In the present study, we investigated the methylation status of the Mkrn3 promoter in the hypothalamus of the female mouse before, during and after puberty. Initially, we mapped the 32 CpG dinucleotides in the promoter, the 5'UTR and the first 50 nucleotides of the coding region of the Mkrn3 gene. Moreover, we identified a short CpG island region (CpG islet) located within the promoter. Methylation analysis using bisulfite sequencing revealed that CpG dinucleotides were methylated regardless of developmental stage, with the lowest levels of methylation being found within the CpG islet region. In addition, the CpG islet region showed significantly lower methylation levels at the pre-pubertal stage when compared with the pubertal or post-pubertal stage. Finally, in silico analysis of transcription factor binding sites on the Mkrn3 CpG islet identified the recruitment of 29 transcriptional regulators of which 14 were transcriptional repressors. Our findings demonstrate the characterization and differential methylation of the CpG dinucleotides located in the Mkrn3 promoter that could influence the transcriptional activity in pre-pubertal compared to pubertal or post-pubertal period. Further studies are needed to clarify the possible mechanisms and effects of differential methylation of the Mkrn3 promoter.

Genome-Wide DNA Methylation Protocol for Epigenetics Studies

Epigenetic modification are heritable changes in gene expression not encoded by the DNA sequence therefore playing a significant role in a broad range of biological processes and diseases.A key player of the epigenetic modifications is the DNA methylation, a process involving the covalent transfer of a methyl group to the C-5 position of the cytosine ring of DNA forming 5-methylcytosine (5mC), catalyzed by DNA methyltransferases. Altering the structure of the chromatin, DNA methylation has the potential to down-regulate gene expression.The here presented protocol shows a method to obtain DNA samples ready for NGS sequencing for genome-wide methylation analysis.

OUP accepted manuscript

Conflicts between transcription and replication machinery are a potent source of replication stress and genome instability; however, no technique currently exists to identify endogenous genomic locations prone to transcription–replication interactions. Here, we report a novel method to identify genomic loci prone to transcription–replication interactions termed transcription–replication immunoprecipitation on nascent DNA sequencing, TRIPn-Seq. TRIPn-Seq employs the sequential immunoprecipitation of RNA polymerase 2 phosphorylated at serine 5 (RNAP2s5) followed by enrichment of nascent DNA previously labeled with bromodeoxyuridine. Using TRIPn-Seq, we mapped 1009 unique transcription–replication interactions (TRIs) in mouse primary B cells characterized by a bimodal pattern of RNAP2s5, bidirectional transcription, an enrichment of RNA:DNA hybrids, and a high probability of forming G-quadruplexes. TRIs are highly enriched at transcription start sites and map to early replicating regions. TRIs exhibit enhanced Replication Protein A association and TRI-associated genes exhibit higher replication fork termination than control transcription start sites, two marks of replication stress. TRIs colocalize with double-strand DNA breaks, are enriched for deletions, and accumulate mutations in tumors. We propose that replication stress at TRIs induces mutations potentially contributing to age-related disease, as well as tumor formation and development.

Misexpression of genes lacking CpG islands drives degenerative changes during aging

Cellular aging is characterized by disruption of the nuclear lamina and its associated heterochromatin. How these structural changes within the nucleus contribute to age-related degeneration of the organism is unclear. Genes lacking CpG islands (CGI− genes) generally associate with heterochromatin when they are inactive. Here, we show that the expression of these genes is globally activated in aged cells and tissues. This CGI− gene misexpression is a common feature of normal and pathological aging in mice and humans. We report evidence that CGI− gene up-regulation is directly responsible for age-related physiological deterioration, notably for increased secretion of inflammatory mediators.

Sex-specific chromatin remodelling safeguards transcription in germ cells

Stability of the epigenetic landscape underpins maintenance of the cell-type-specific transcriptional profile. As one of the main repressive epigenetic systems, DNA methylation has been shown to be important for long-term gene silencing; its loss leads to ectopic and aberrant transcription in differentiated cells and cancer¹. The developing mouse germ line endures global changes in DNA methylation in the absence of widespread transcriptional activation. Here, using an ultra-low-input native chromatin immunoprecipitation approach, we show that following DNA demethylation the gonadal primordial germ cells undergo remodelling of repressive histone modifications, resulting in a sex-specific signature in mice. We further demonstrate that Polycomb has a central role in transcriptional control in the newly hypomethylated germline genome as the genetic loss of Ezh2 leads to aberrant transcriptional activation, retrotransposon derepression and dramatic loss of developing female germ cells. This sex-specific effect of Ezh2 deletion is explained by the distinct landscape of repressive modifications observed in male and female germ cells. Overall, our study provides insight into the dynamic interplay between repressive chromatin modifications in the context of a developmental reprogramming system.

DNA methylations of brain-derived neurotrophic factor exon VI are associated with major depressive disorder and antidepressant-induced remission in females

BACKGROUND

Brain-derived neurotrophic factor (BDNF) has been suggested to play important roles in major depressive disorder (MDD) and antidepressant treatment. The main purpose of this study was to evaluate the association of DNA methylation changes in the BDNF gene with MDD and antidepressant treatment.

METHODS

A total of 291 MDD patients and 100 healthy controls were included and followed up for 6 weeks. The Hamilton Depression Rating Scale-17 (HDRS-17) was used to measure treatment improvement. The life events scales (LES) and childhood trauma questionnaire (CTQ) were used to rate recent and early life stress. DNA methylation levels of CpG sites in the BDNF gene were measured.

RESULTS

Two CpG sites in BDNF exon VI (BDNF133 and BDNF134) were demonstrated to have significantly higher methylation in MDD patients than in controls (both FDR-adjusted P = 0.001). A logistics regression model indicated that the interaction between the hypermethylation of BDNF133 and negative subscore of LES was associated to MDD (OR=0.0075, P<0.001). Methylation of BDNF140 at baseline was significantly elevated in remitters (FDR-adjusted P = 0.046) at week 6. In subgroup analyses, these findings could be replicated in females, but not in males.

LIMITATIONS

The methylation status of BDNF after 6 weeks of antidepressant treatment was not measured and the DNA methylation were detected in peripheral blood cells.

CONCLUSIONS

These findings highlight gender-specific alteration of methylation at several CpG sites in BDNF exon VI as a promising candidate indicator of MDD and antidepressant-induced remission.

Predicting environmentally responsive transgenerational differential DNA methylated regions (epimutations) in the genome using a hybrid deep-machine learning approach

BACKGROUND

Deep learning is an active bioinformatics artificial intelligence field that is useful in solving many biological problems, including predicting altered epigenetics such as DNA methylation regions. Deep learning (DL) can learn an informative representation that addresses the need for defining relevant features. However, deep learning models are computationally expensive, and they require large training datasets to achieve good classification performance.

RESULTS

One approach to addressing these challenges is to use a less complex deep learning network for feature selection and Machine Learning (ML) for classification. In the current study, we introduce a hybrid DL-ML approach that uses a deep neural network for extracting molecular features and a non-DL classifier to predict environmentally responsive transgenerational differential DNA methylated regions (DMRs), termed epimutations, based on the extracted DL-based features. Various environmental toxicant induced epigenetic transgenerational inheritance sperm epimutations were used to train the model on the rat genome DNA sequence and use the model to predict transgenerational DMRs (epimutations) across the entire genome.

CONCLUSION

The approach was also used to predict potential DMRs in the human genome. Experimental results show that the hybrid DL-ML approach outperforms deep learning and traditional machine learning methods.

Nucleic Acid-Sensing and Interferon-Inducible Pathways Show Differential Methylation in MZ Twins Discordant for Lupus and Overexpression in Independent Lupus Samples: Implications for Pathogenic Mechanism and Drug Targeting

Systemic lupus erythematosus (SLE) is a chronic, multisystem, autoimmune inflammatory disease with genomic and non-genomic contributions to risk. We hypothesize that epigenetic factors are a significant contributor to SLE risk and may be informative for identifying pathogenic mechanisms and therapeutic targets. To test this hypothesis while controlling for genetic background, we performed an epigenome-wide analysis of DNA methylation in genomic DNA from whole blood in three pairs of female monozygotic (MZ) twins of European ancestry, discordant for SLE. Results were replicated on the same array in four cell types from a set of four Danish female MZ twin pairs discordant for SLE. Genes implicated by the epigenetic analyses were then evaluated in 10 independent SLE gene expression datasets from the Gene Expression Omnibus (GEO). There were 59 differentially methylated loci between unaffected and affected MZ twins in whole blood, including 11 novel loci. All but two of these loci were hypomethylated in the SLE twins relative to the unaffected twins. The genes harboring these hypomethylated loci exhibited increased expression in multiple independent datasets of SLE patients. This pattern was largely consistent regardless of disease activity, cell type, or renal tissue type. The genes proximal to CpGs exhibiting differential methylation (DM) in the SLE-discordant MZ twins and exhibiting differential expression (DE) in independent SLE GEO cohorts (DM-DE genes) clustered into two pathways: the nucleic acid-sensing pathway and the type I interferon pathway. The DM-DE genes were also informatically queried for potential gene–drug interactions, yielding a list of 41 drugs including a known SLE therapy. The DM-DE genes delineate two important biologic pathways that are not only reflective of the heterogeneity of SLE but may also correlate with distinct IFN responses that depend on the source, type, and location of nucleic acid molecules and the activated receptors in individual patients. Cell- and tissue-specific analyses will be critical to the understanding of genetic factors dysregulating the nucleic acid-sensing and IFN pathways and whether these factors could be appropriate targets for therapeutic intervention.

Effect of DNA modifications on the transition between canonical and non-canonical DNA structures in CpG islands during senescence

Patterns and levels of DNA modifications play important roles in senescence. Two major epigenetic modifications of DNA, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), target CpG sites. Importantly, CpG concentrated regions, known as CpG islands, contain GC-rich sequences, which have the potential to fold into non-canonical DNA structures such as i-motifs and G-quadruplexes. In this study, we investigated the effect of 5mC and 5hmC modifications on the transition between a duplex, and i-motif and G-quadruplexes. To examine the transition, we firstly investigated the stability and structure of the i-motif and G-quadruplexes, considering the molecular environment in senescent cells. Analyses of their stability showed that the modifications did not drastically affect the stability. However, noteworthily, the modification can weaken the (de)stabilisation effect on G-quadruplexes caused by cosolute(s) and cations. Circular dichroism analyses indicated that the surrounding environments, including the molecular crowding and the type of cations such as K⁺ and Na⁺, regulate the topology of G-quadruplexes, while neither 5mC nor 5hmC had a drastic effect. On the other hand, the modifications changed the transition between duplexes and quadruplexes. Unmodified DNA preferred to fold into quadruplexes, whereas DNA with 5mC and 5hmC preferred to fold into duplexes in the absence of PEG200; on the other hand, DNA with or without modifications tended to fold into i-motifs under crowded conditions. Furthermore, an investigation of quadruplexes forming sequences in CpG islands, which are hyper- or hypomethylated during senescence, followed by gene ontology enrichment analysis for each gene group classified by the presence of quadruplexes, showed a difference in function between genes with and without quadruplexes in the CpG region. These results indicate that it is important to consider the effects of patterns and levels of DNA modifications on the transition between canonical and non-canonical DNA structures to understand gene regulation by epigenetic modification during senescence.

The modification of DNA can regulate the transition between a duplex and quadruplexes during senescence responding to surrounding environments.

Quantifying the tissue-specific regulatory information within enhancer DNA sequences

Recent efforts to measure epigenetic marks across a wide variety of different cell types and tissues provide insights into the cell type-specific regulatory landscape. We use these data to study whether there exists a correlate of epigenetic signals in the DNA sequence of enhancers and explore with computational methods to what degree such sequence patterns can be used to predict cell type-specific regulatory activity. By constructing classifiers that predict in which tissues enhancers are active, we are able to identify sequence features that might be recognized by the cell in order to regulate gene expression. While classification performances vary greatly between tissues, we show examples where our classifiers correctly predict tissue-specific regulation from sequence alone. We also show that many of the informative patterns indeed harbor transcription factor footprints.

Antitumoral activity of liraglutide, a new DNMT inhibitor in breast cancer cells in vitro and in vivo

Breast cancer (BC) is the most frequently diagnosed female cancer and second leading cause of death. Despite the discovery of many antineoplastic drugs for BC, the current therapy is not totally efficient. In this study, we investigated the potential of repurposing the well-known diabetes type II drug liraglutide to modulate epigenetic modifications in BC cells lines in vitro and in vivo via Ehrlich mice tumors models. The in vitro results revealed a significant reduction on cell viability, migration, DNMT activity and displayed lower levels of global DNA methylation in BC cell lines after liraglutide treatment. The interaction between liraglutide and the DNMT enzymes resulted in a decrease profile of DNA methylation for the CDH1, ESR1 and ADAM33 gene promoter regions and, consequently, increased their gene and protein expression levels. To elucidate the possible interaction between liraglutide and the DNMT1 protein, we performed an in silico study that indicates liraglutide binding in the catalytic cleft via hydrogen bonds and salt bridges with the interdomain contacts and disturbs the overall enzyme conformation. The in vivo study was also able to reveal that liraglutide and the combined treatment of liraglutide and paclitaxel or methotrexate were effective in reducing tumor growth. Moreover, the modulation of CDH1 and ADAM33 mouse gene expression by DNA demethylation suggests a role for liraglutide in DNMT activity in vivo. Altogether, these results indicate that liraglutide may be further analysed as a new adjuvant treatment for BC.

Sequence features of retrotransposons allow for epigenetic variability

Transposable elements (TEs) are mobile genetic elements that make up a large fraction of mammalian genomes. While select TEs have been co-opted in host genomes to have function, the majority of these elements are epigenetically silenced by DNA methylation in somatic cells. However, some TEs in mice, including the Intracisternal A-particle (IAP) subfamily of retrotransposons, have been shown to display interindividual variation in DNA methylation. Recent work has revealed that IAP sequence differences and strain-specific KRAB zinc finger proteins (KZFPs) may influence the methylation state of these IAPs. However, the mechanisms underlying the establishment and maintenance of interindividual variability in DNA methylation still remain unclear. Here, we report that sequence content and genomic context influence the likelihood that IAPs become variably methylated. IAPs that differ from consensus IAP sequences have altered KZFP recruitment that can lead to decreased KAP1 recruitment when in proximity of constitutively expressed genes. These variably methylated loci have a high CpG density, similar to CpG islands, and can be bound by ZF-CxxC proteins, providing a potential mechanism to maintain this permissive chromatin environment and protect from DNA methylation. These observations indicate that variably methylated IAPs escape silencing through both attenuation of KZFP binding and recognition by ZF-CxxC proteins to maintain a hypomethylated state.

Changes in PGC-1α-Dependent Mitochondrial Biogenesis Are Associated with Inflexible Hepatic Energy Metabolism in the Offspring Born to Dexamethasone-Treated Mothers

In the present study we investigated the participation of hepatic peroxisome proliferator-activated receptor-gamma coactivator 1 alpha (PGC-1α) in the metabolic programming of newborn rats exposed in utero to dexamethasone (DEX). On the 21st day of life, fasted offspring born to DEX-treated mothers displayed increased conversion of pyruvate into glucose with simultaneous upregulation of PEPCK (phosphoenolpyruvate carboxykinase) and G6Pase (glucose-6-phosphatase). Increased oxidative phosphorylation, higher ATP/ADP ratio and mitochondrial biogenesis and lower pyruvate levels were also found in the progeny of DEX-treated mothers. On the other hand, the 21-day-old progeny of DEX-treated mothers had increased hepatic triglycerides (TAG) and lower CPT-1 activity when subjected to short-term fasting. At the mechanistic level, rats exposed in utero to DEX exhibited increased hepatic PGC-1α protein content with lower miR-29a-c expression. Increased PGC-1α content was concurrent with increased association to HNF-4α and NRF1 and reduced PPARα expression. The data presented herein reveal that changes in the transcription machinery in neonatal liver of rats born to DEX-treated mothers leads to an inflexible metabolic response to fasting. Such programming is hallmarked by increased oxidative phosphorylation of pyruvate with impaired FFA oxidation and hepatic TAG accumulation.

Low folate induces abnormal neuronal maturation and DNA hypomethylation of neuronal differentiation-related genes in cultured mouse neural stem and progenitor cells

Folate deficiency in a fetus is well known to cause neurodevelopment defects and development disorders. A low level of folate is also thought to be a risk for depression in adults. We have previously shown that post-weaning low folate induces neuronal immaturity in the dentate gyrus in mice, which suggests that low folate causes neuropsychological disorders via inhibition of neuronal maturation. In this study, we examined the effects of low folate on expression and epigenetic modification of genes involved in neuronal differentiation and maturation in primary mouse neural stem/progenitor cells (NSPCs) in vitro. An increase in Nestin (NSPC marker)-positive cells was observed in cells differentiated in a low folate medium for 3 days. An increase in βIII-tubulin (Tuj1: immature neuron marker)-positive cells and a decrease in microtubule-associated protein 2 (MAP2: mature neuron marker)-positive cells were observed in cells differentiated in a low folate medium for 7 days. In these cells, mRNA levels for genes involved in neuronal differentiation and maturation were altered. Hypomethylation of DNA, but not of histone proteins, was also observed at some promoters of these neuronal genes. The level of S-adenosylmethionine (SAM), a methyl donor, was decreased in these cells. The abnormalities in neural maturation and changes in gene expression in culture under low folate conditions were partially normalized by addition of SAM (5 μM). Based on these results, decreased SAM may induce DNA hypomethylation at genes involved in neuronal differentiation and maturation under low folate conditions, and this hypomethylation may be associated with low folate-induced neuronal immaturity.

Human Induced Pluripotent Stem Cell-Derived Mesenchymal Stem Cells Acquire Rejuvenation and Reduced Heterogeneity

Despite the uniform selection criteria for the isolation of human mesenchymal stem cells (MSCs), considerable heterogeneity exists which reflects the distinct tissue origins and differences between individuals with respect to their genetic background and age. This heterogeneity is manifested by the variabilities seen in the transcriptomes, proteomes, secretomes, and epigenomes of tissue-specific MSCs. Here, we review literature on different aspects of MSC heterogeneity including the role of epigenetics and the impact of MSC heterogeneity on therapies. We then combine this with a meta-analysis of transcriptome data from distinct MSC subpopulations derived from bone marrow, adipose tissue, cruciate, tonsil, kidney, umbilical cord, fetus, and induced pluripotent stem cells derived MSCs (iMSCs). Beyond that, we investigate transcriptome differences between tissue-specific MSCs and pluripotent stem cells. Our meta-analysis of numerous MSC-related data sets revealed markers and associated biological processes characterizing the heterogeneity and the common features of MSCs from various tissues. We found that this heterogeneity is mainly related to the origin of the MSCs and infer that microenvironment and epigenetics are key drivers. The epigenomes of MSCs alter with age and this has a profound impact on their differentiation capabilities. Epigenetic modifications of MSCs are propagated during cell divisions and manifest in differentiated cells, thus contributing to diseased or healthy phenotypes of the respective tissue. An approach used to reduce heterogeneity caused by age- and tissue-related epigenetic and microenvironmental patterns is the iMSC concept: iMSCs are MSCs generated from induced pluripotent stem cells (iPSCs). During iMSC generation epigenetic and chromatin remodeling result in a gene expression pattern associated with rejuvenation thus allowing to overcome age-related shortcomings (e.g., limited differentiation and proliferation capacity). The importance of the iMSC concept is underlined by multiple clinical trials. In conclusion, we propose the use of rejuvenated iMSCs to bypass tissue- and age-related heterogeneity which are associated with native MSCs.

Genome-wide DNA methylation analysis of dogs with high lead exposure living near a lead mining area in Kabwe, Zambia

Lead (Pb) is a heavy metal that has been proven to be toxic to both animals and humans. Genom-wide DNA methylation in domestic dogs exposed to high levels of Pb in Kabwe, Zambia was analyzed in this study. Using next-generation sequencing on samples from 20 domestic dogs (mean blood Pb concentration: 43.6 μg/dL and 7.2 μg/dL in the high and low exposure groups), a digital restriction enzyme analysis of methylation was performed to identify the genomic locations of differentially methylated CpG sites. A validation study on an additional 20 dogs followed (blood Pb concentration: 4.9-29.7 μg/dL). The cluster analysis resolved two broad clusters indicating high and low Pb exposure. The study identified 827 (1.2%) CpG sites with differences in methylation (101 CpG sites were hypermethylated in the low exposure group and 726 were hypermethylated in the high exposure group). The sites corresponded to 26 genes with differentially methylated CpG sites at their promoter regions, including the NGF gene. The methylation of four CpG sites was validated using bisulfite pyrosequencing. The results indicate that aberrant hypermethylation is prevalent in dogs exposed to Pb. The altered DNA methylation of the genes identified in this study contributes to a greater understanding of the epigenetic changes caused by Pb exposure and highlights novel biomarker discoveries across species.

Genome-wide DNA methylation profile in feline haematological tumours: A preliminary study

Although DNA methylation has been analysed in few studies for a limited number of loci in cats with diseases, genome-wide profile of DNA methylation has never been addressed. The hypothesis for this study is that next-generation sequencing with sequential digestion of genomic DNA with SmaI and XmaI enzymes could provide highly quantitative information on methylation levels in cats. Using blood from four healthy control cats and two disease cats as well as three feline lymphoma/leukemia cell lines, approximately 74-94 thousand CpG sites across the cat genome could be analysed. CpG sites in CpG island (CGI) were broadly either methylated or unmethylated in normal blood, while CpG sites in non-CpG islands (NCGI) are largely methylated. Lymphoma cell lines showed thousands of CpG sites with gain of methylation at normally unmethylated CGI sites and loss of methylation at normally methylated NCGI sites. Hypermethylated CpG sites located at promoter regions included genes annotated with 'developmental process' and 'anatomical structure morphogenesis' such as HOXD10. This highly quantitative method would be suitable for studies of DNA methylation changes not only in cancer but also in other common diseases in cats.