Physio-biochemical and DNA methylation analysis of the defense response network of wheat to drought stress

In the present work, physio-biochemical and DNA methylation analysis were conducted in wheat (Triticum aestivum L.) cultivars "Bolani" (drought-tolerant) and "Sistan" (drought-sensitive) during drought treatments: well-watered (at 90% field capacity (FC)), mild stress (at 50% FC, and severe stress (at 25% FC). During severe stress, O2•- and H2O2 content in cultivar Sistan showed significant increase (by 1.3 and 2.5-fold, respectively) relative to cultivar Bolani. In Bolani, the increased levels of radical scavenging activity (by 32%), glycine betaine (GB) (by 11.44%), proline (4-fold), abscisic acid (by 63.76%), and more stability of relative water content (RWC) (2-fold) were observed against drought-induced oxidative stress. Methylation level significantly decreased from 70.26% to 60.64% in Bolani and from 69.06% to 59.85% in Sistan during stress, and higher decreased tendency was related to CG and CHG in Bolani but CG in Sistan under severe stress. Methylation patterns showed that the highest polymorphism in Bolani was mainly as CG. As the intensity of stress increased, the enhanced physio-biochemical responses of Bolani cultivar were accompanied by a more decrease in the number of unchanged bands. According to heat map analysis, the highest difference (84.38%) in methylation patterns was observed between control and severe stress. Multivariate analysis using principal component analysis (PCA) showed a cultivar-specific methylation during stress and that methylation changes between cultivars are much higher than that of within a cultivar. Higher methylation to demethylation in Bolani (30.06 vs. 22.12%) compared to that of cultivar Sistan (23.21 vs. 30.15%) indicated more demethylation did not induce tolerance responses in Sistan. Sequencing differentially methylated fragments along with qRT-PCR analysis showed the efficient role of various DNA fragments, including demethylated fragments such as phosphoenol pyruvate carboxylase (PEPC), beta-glucosidase (BGlu), glycosyltransferase (GT), glutathione S-transferase (GST) and lysine demethylase (LSD) genes and methylated fragments like ubiquitin E2 enzyme genes in the development of drought tolerance. These results suggested the specific roles of DNA methylation in development of drought tolerance in wheat landrace.

Epigenome-wide association study of DNA methylation in maternal blood leukocytes with BMI in pregnancy and gestational weight gain

OBJECTIVES

We aimed to discover CpG sites with differential DNA methylation in peripheral blood leukocytes associated with body mass index (BMI) in pregnancy and gestational weight gain (GWG) in women of European and South Asian ancestry. Furthermore, we aimed to investigate how the identified sites were associated with methylation quantitative trait loci, gene ontology, and cardiometabolic parameters.

METHODS

In the Epigenetics in pregnancy (EPIPREG) sample we quantified maternal DNA methylation in peripheral blood leukocytes in gestational week 28 with Illumina's MethylationEPIC BeadChip. In women with European (n = 303) and South Asian (n = 164) ancestry, we performed an epigenome-wide association study of BMI in gestational week 28 and GWG between gestational weeks 15 and 28 using a meta-analysis approach. Replication was performed in the Norwegian Mother, Father, and Child Cohort Study, the Study of Assisted Reproductive Technologies (MoBa-START) (n = 877, mainly European/Norwegian).

RESULTS

We identified one CpG site significantly associated with GWG (p 5.8 × 10-8) and five CpG sites associated with BMI at gestational week 28 (p from 4.0 × 10-8 to 2.1 × 10-10). Of these, we were able to replicate three in MoBa-START; cg02786370, cg19758958 and cg10472537. Two sites are located in genes previously associated with blood pressure and BMI. DNA methylation at the three replicated CpG sites were associated with levels of blood pressure, lipids and glucose in EPIPREG (p from 1.2 × 10-8 to 0.04).

CONCLUSIONS

We identified five CpG sites associated with BMI at gestational week 28, and one with GWG. Three of the sites were replicated in an independent cohort. Several genetic variants were associated with DNA methylation at cg02786379 and cg16733643 suggesting a genetic component influencing differential methylation. The identified CpG sites were associated with cardiometabolic traits.

CLINICALTRIALS

GOV REGISTRATION NO

Not applicable.

A preliminary investigation of epigenome-wide DNA methylation and temperament during infancy

This study provides preliminary evidence for an epigenetic architecture of infant temperament. At 12 months of age, blood was collected and assayed for DNA methylation and maternally reported infant temperament was assessed using the Infant Behavior Questionnaire in 67 mother-infant dyads. Epigenome-wide analyses showed that the higher order temperament dimensions Surgency and Negative Affect were associated with DNA methylation. The epigenetic signatures of Surgency and Negative Affect were situated at genes involved in synaptic signaling and plasticity. Although replication is required, these results are consistent with a biologically based model of temperament, create new avenues for hypothesis-driven research into epigenetic pathways that underlie individual differences in temperament, and demonstrate that infant temperament has a widespread epigenetic signature in the methylome.

Reshaped DNA methylation cooperating with homoeolog-divergent expression promotes improved root traits in synthesized tetraploid wheat

Polyploidization is a major event driving plant evolution and domestication. However, how reshaped epigenetic modifications coordinate gene transcription to generate phenotypic variations during wheat polyploidization is currently elusive. Here, we profiled transcriptomes and DNA methylomes of two diploid wheat accessions (SlSl and AA) and their synthetic allotetraploid wheat line (SlSlAA), which displayed elongated root hair and improved root capability for nitrate uptake and assimilation after tetraploidization. Globally decreased DNA methylation levels with a reduced difference between subgenomes were observed in the roots of SlSlAA. DNA methylation changes in first exon showed strong connections with altered transcription during tetraploidization. Homoeolog-specific transcription was associated with biased DNA methylation as shaped by homoeologous sequence variation. The hypomethylated promoters showed significantly enriched binding sites for MYB, which may affect gene transcription in response to root hair growth. Two master regulators in root hair elongation pathway, AlCPC and TuRSL4, exhibited upregulated transcription levels accompanied by hypomethylation in promoter, which may contribute to the elongated root hair. The upregulated nitrate transporter genes, including NPFs and NRTs, also are significantly associated with hypomethylation, indicating an epigenetic-incorporated regulation manner in improving nitrogen use efficiency. Collectively, these results provided new insights into epigenetic changes in response to crop polyploidization and underscored the importance of epigenetic regulation in improving crop traits.

Epigenetic and transcriptional landscapes during cerebral cortex development in a microcephaly mouse model

The cerebral cortex is a pivotal structure integral to advanced brain functions within the mammalian central nervous system. DNA methylation and hydroxymethylation play important roles in regulating cerebral cortex development. However, it remains unclear whether abnormal cerebral cortex development, such as microcephaly, could rescale the epigenetic landscape, potentially contributing to dysregulated gene expression during brain development. In this study, we characterize and compare the DNA methylome/hydroxymethylome and transcriptome profiles of the cerebral cortex across several developmental stages in wild-type (WT) mice and Mcph1 knockout (Mcph1-del) mice with severe microcephaly. Intriguingly, we discover a global reduction of 5'-hydroxymethylcytosine (5hmC) level, primarily in TET1-binding regions, in Mcph1-del mice compared to WT mice during juvenile and adult stages. Notably, genes exhibiting diminished 5hmC levels and concurrently decreased expression are essential for neurodevelopment and brain functions. Additionally, genes displaying a delayed accumulation of 5hmC in Mcph1-del mice are significantly associated with the establishment and maintenance of the nervous system during the adult stage. These findings reveal that aberrant cerebral cortex development in the early stages profoundly alters the epigenetic regulation program, which provides unique insights into the molecular mechanisms underpinning diseases related to cerebral cortex development.

Mapping parental DMRs predictive of local and distal methylome remodeling in epigenetic F1 hybrids

F1 hybrids derived from a cross between two inbred parental lines often display widespread changes in DNA methylation and gene expression patterns relative to their parents. An emerging challenge is to understand how parental epigenomic differences contribute to these events. Here, we generated a large mapping panel of F1 epigenetic hybrids, whose parents are isogenic but variable in their DNA methylation patterns. Using a combination of multi-omic profiling and epigenetic mapping strategies we show that differentially methylated regions in parental pericentromeres act as major reorganizers of hybrid methylomes and transcriptomes, even in the absence of genetic variation. These parental differentially methylated regions are associated with hybrid methylation remodeling events at thousands of target regions throughout the genome, both locally (in cis) and distally (in trans). Many of these distally-induced methylation changes lead to nonadditive expression of nearby genes and associate with phenotypic heterosis. Our study highlights the pleiotropic potential of parental pericentromeres in the functional remodeling of hybrid genomes and phenotypes.

Genome-wide identification, evolution of DNA methyltransferases and their expression under salinity stress in Larimichthys crocea

DNA methyltransferases (Dnmts) are responsible for DNA methylation which influences patterns of gene expression and plays a crucial role in response to environmental changes. In this study, 7 LcDnmt genes were identified in the genome of large yellow croaker (Larimichthys crocea). The comprehensive analysis was conducted on gene structure, protein and location site of LcDnmts. LcDnmt proteins belonged to three groups (Dnmt1, Dnmt2, and Dnmt3) according to their conserved domains and phylogenetic analysis. Although Dnmt3 can be further divided into three sub groups (Dnmt3a, Dnmt3b, and Dnmt3l), there is no Dnmnt3l member in the large yellow croaker. Phylogenetic analysis revealed that the Dnmt family was highly conserved in teleosts. Expression patterns derived from the RNA-seq, qRT-PCR and Western blot analysis revealed that 2 LcDnmt genes (LcDnmt1 and LcDnmt3a2) significantly regulated under salinity stress in the liver, which was found to be dominantly expressed in the intestine and brain, respectively. These two genes may play an important role in the salinity stress of large yellow croaker and represent candidates for future functional analysis. Our results revealed the conservation of Dnmts during evolution and indicated a potential role of Dnmts in epigenetic regulation of response to salinity stress.

DNMT1 methylation of LncRNA-ANRIL causes myocardial fibrosis pyroptosis by interfering with the NLRP3/Caspase-1 pathway

Myocardial fibrosis is a common pathological manifestation that occurs in various cardiac diseases. The present investigation aims to reveal how DNMT1/lncRNA-ANRIL/NLRP3 influences fibrosis and cardiac fibroblast pyroptosis. Here, we used ISO to induce myocardial fibrosis in mice, and LPS and ATP to induce myocardial fibroblast pyroptosis. The results showed that DNMT1, Caspase-1, and NLRP3 expression were significantly increased in fibrotic murine myocardium and pyroptotic cardiac fibroblasts, whereas LncRNA-ANRIL expression was decreased. DNMT1 overexpression decreased the level of LncRNA-ANRIL while increasing the levels of NLRP3 and Caspase-1. Contrarily, silencing DNMT1 increased the LncRNA-ANRIL and decreased the levels of NLRP3 and Caspase-1. Silencing LncRNA-ANRIL increased the levels of NLRP3 and Caspase-1. The present findings suggest that DNMT1 can methylate LncRNA-ANRIL during the development of myocardial fibrosis and CFs cell scorching, resulting in low LncRNA-ANRIL expression, thereby influencing myocardial fibrosis and cardiac fibroblast pyroptosis.

FinaleMe: Predicting DNA methylation by the fragmentation patterns of plasma cell-free DNA

Analysis of DNA methylation in cell-free DNA reveals clinically relevant biomarkers but requires specialized protocols such as whole-genome bisulfite sequencing. Meanwhile, millions of cell-free DNA samples are being profiled by whole-genome sequencing. Here, we develop FinaleMe, a non-homogeneous Hidden Markov Model, to predict DNA methylation of cell-free DNA and, therefore, tissues-of-origin, directly from plasma whole-genome sequencing. We validate the performance with 80 pairs of deep and shallow-coverage whole-genome sequencing and whole-genome bisulfite sequencing data.

Mating system is associated with seed phenotypes upon loss of RNA-directed DNA methylation in Brassicaceae

In plants, de novo DNA methylation is guided by 24-nt short interfering (si)RNAs in a process called RNA-directed DNA methylation (RdDM). Primarily targeted at transposons, RdDM causes transcriptional silencing and can indirectly influence expression of neighboring genes. During reproduction, a small number of siRNA loci are dramatically upregulated in the maternally derived seed coat, suggesting that RdDM might have a special function during reproduction. However, the developmental consequence of RdDM has been difficult to dissect because disruption of RdDM does not result in overt phenotypes in Arabidopsis (Arabidopsis thaliana), where the pathway has been most thoroughly studied. In contrast, Brassica rapa mutants lacking RdDM have a severe seed production defect, which is determined by the maternal sporophytic genotype. To explore the factors that underlie the different phenotypes of these species, we produced RdDM mutations in 3 additional members of the Brassicaceae family: Camelina sativa, Capsella rubella, and Capsella grandiflora. Among these 3 species, only mutations in the obligate outcrosser, C. grandiflora, displayed a seed production defect similar to Brassica rapa mutants, suggesting that mating system is a key determinant for reproductive phenotypes in RdDM mutants.

Differential DNA methylation landscape of miRNAs genes in mice liver fibrosis

BACKGROUND

Patients with chronic liver disease were found nearly all to have liver fibrosis, which is characterized by excess accumulation of extracellular matrix (ECM) proteins. While ECM accumulation can prevent liver infection and injury, it can destroy normal liver function and architecture. miRNA's own regulation was involved in DNA methylation change. The purpose of this study is to detect DNA methylation landscape of miRNAs genes in mice liver fibrosis tissues.

METHODS

Male mice (10-12 weeks) were injected CCl4 from abdominal cavity to induced liver fibrosis. 850 K BeadChips were used to examine DNA methylation change in whole genome. The methylation change of 16 CpG dinucleotides located in promoter regions of 4 miRNA genes were detected by bisulfite sequencing polymerase chain reaction (BSP) to verify chip data accuracy, and these 4 miRNA genes' expressions were detected by RT-qPCR methods.

RESULTS

There are 769 differential methylation sites (DMS) in total between fibrotic liver tissue and normal mice liver tissue, which were related with 148 different miRNA genes. Chips array data were confirmed by bisulfite sequencing polymerase chain reaction (R = 0.953; P < 0.01). GO analysis of the target genes of 2 miRNA revealed that protein binding, cytoplasm and chromatin binding activity were commonly enriched; KEGG pathway enrichment analysis displayed that TGF-beta signaling pathway was commonly enriched.

CONCLUSION

The DNA of 148 miRNA genes was found to have methylation change in liver fibrosis tissue. These discoveries in miRNA genes are beneficial to future miRNA function research in liver fibrosis.

Population epigenetics: DNA methylation in the plant omics era

DNA methylation plays an important role in many biological processes. The mechanisms underlying the establishment and maintenance of DNA methylation are well understood thanks to decades of research using DNA methylation mutants, primarily in Arabidopsis (Arabidopsis thaliana) accession Col-0. Recent genome-wide association studies (GWASs) using the methylomes of natural accessions have uncovered a complex and distinct genetic basis of variation in DNA methylation at the population level. Sequencing following bisulfite treatment has served as an excellent method for quantifying DNA methylation. Unlike studies focusing on specific accessions with reference genomes, population-scale methylome research often requires an additional round of sequencing beyond obtaining genome assemblies or genetic variations from whole-genome sequencing data, which can be cost prohibitive. Here, we provide an overview of recently developed bisulfite-free methods for quantifying methylation and cost-effective approaches for the simultaneous detection of genetic and epigenetic information. We also discuss the plasticity of DNA methylation in a specific Arabidopsis accession, the contribution of DNA methylation to plant adaptation, and the genetic determinants of variation in DNA methylation in natural populations. The recently developed technology and knowledge will greatly benefit future studies in population epigenomes.

Transcription factor PpNAC1 and DNA demethylase PpDML1 synergistically regulate peach fruit ripening

Fruit ripening is accompanied by dramatic changes in color, texture, and flavor and is regulated by transcription factors (TFs) and epigenetic factors. However, the detailed regulatory mechanism remains unclear. Gene expression patterns suggest that PpNAC1 (NAM/ATAF1/2/CUC) TF plays a major role in peach (Prunus persica) fruit ripening. DNA affinity purification (DAP)-seq combined with transactivation tests demonstrated that PpNAC1 can directly activate the expression of multiple ripening-related genes, including ACC synthase1 (PpACS1) and ACC oxidase1 (PpACO1) involved in ethylene biosynthesis, pectinesterase1 (PpPME1), pectate lyase1 (PpPL1), and polygalacturonase1 (PpPG1) related to cell wall modification, and lipase1 (PpLIP1), fatty acid desaturase (PpFAD3-1), and alcohol acyltransferase1 (PpAAT1) involved in volatiles synthesis. Overexpression of PpNAC1 in the tomato (Solanum lycopersicum) nor (nonripening) mutant restored fruit ripening, and its transient overexpression in peach fruit induced target gene expression, supporting a positive role of PpNAC1 in fruit ripening. The enhanced transcript levels of PpNAC1 and its target genes were associated with decreases in their promoter mCG methylation during ripening. Declining DNA methylation was negatively associated with increased transcripts of DNA demethylase1 (PpDML1), whose promoter is recognized and activated by PpNAC1. We propose that decreased methylation of the promoter region of PpNAC1 leads to a subsequent decrease in DNA methylation levels and enhanced transcription of ripening-related genes. These results indicate that positive feedback between PpNAC1 and PpDML1 plays an important role in directly regulating expression of multiple genes required for peach ripening and quality formation.

Heritable changes of epialleles near genes in maize can be triggered in the absence of CHH methylation

Trans-chromosomal interactions resulting in changes in DNA methylation during hybridization have been observed in several plant species. However, little is known about the causes or consequences of these interactions. Here, we compared DNA methylomes of F1 hybrids that are mutant for a small RNA biogenesis gene, Mop1 (Mediator of paramutation1), with that of their parents, wild-type siblings, and backcrossed progeny in maize (Zea mays). Our data show that hybridization triggers global changes in both trans-chromosomal methylation (TCM) and trans-chromosomal demethylation (TCdM), most of which involved changes in CHH methylation. In more than 60% of these TCM differentially methylated regions (DMRs) in which small RNAs are available, no significant changes in the quantity of small RNAs were observed. Methylation at the CHH TCM DMRs was largely lost in the mop1 mutant, although the effects of this mutant varied depending on the location of these DMRs. Interestingly, an increase in CHH at TCM DMRs was associated with enhanced expression of a subset of highly expressed genes and suppressed expression of a small number of lowly expressed genes. Examination of the methylation levels in backcrossed plants demonstrates that both TCM and TCdM can be maintained in the subsequent generation, but that TCdM is more stable than TCM. Surprisingly, although increased CHH methylation in most TCM DMRs in F1 plants required Mop1, initiation of a new epigenetic state of these DMRs did not require a functional copy of this gene, suggesting that initiation of these changes is independent of RNA-directed DNA methylation.

Rice requires a chromatin remodeler for Polymerase IV-small interfering RNA production and genomic immunity

Transgenes are often spontaneously silenced, which hinders the application of genetic modifications to crop breeding. While gene silencing has been extensively studied in Arabidopsis (Arabidopsis thaliana), the molecular mechanism of transgene silencing remains elusive in crop plants. We used rice (Oryza sativa) plants silenced for a 35S::OsGA2ox1 (Gibberellin 2-oxidase 1) transgene to isolate five elements mountain (fem) mutants showing restoration of transgene expression. In this study, we isolated multiple fem2 mutants defective in a homolog of Required to Maintain Repression 1 (RMR1) of maize (Zea mays) and CLASSY (CLSY) of Arabidopsis. In addition to failing to maintain transgene silencing, as occurs in fem3, in which mutation occurs in NUCLEAR RNA POLYMERASE E1 (OsNRPE1), the fem2 mutant failed to establish transgene silencing of 35S::OsGA2ox1. Mutation in FEM2 eliminated all RNA POLYMERASE IV (Pol-IV)-FEM1/OsRDR2 (RNA-DEPENDENT RNA POLYMERASE 2)-dependent small interfering RNAs (siRNAs), reduced DNA methylation on genome-wide scale in rice seedlings, caused pleiotropic developmental defects, and increased disease resistance. Simultaneous mutation in 2 FEM2 homologous genes, FEM2-Like 1 (FEL1) and FEL2, however, did not affect DNA methylation and rice development and disease resistance. The predominant expression of FEM2 over FEL1 and FEL2 in various tissues was likely caused by epigenetic states. Overexpression of FEL1 but not FEL2 partially rescued hypomethylation of fem2, indicating that FEL1 maintains the cryptic function. In summary, FEM2 is essential for establishing and maintaining gene silencing; moreover, FEM2 is solely required for Pol IV-FEM1 siRNA biosynthesis and de novo DNA methylation.

The MBD-ACD DNA methylation reader complex recruits MICRORCHIDIA6 to regulate ribosomal RNA gene expression in Arabidopsis

DNA methylation is an important epigenetic mark implicated in selective rRNA gene expression, but the DNA methylation readers and effectors remain largely unknown. Here, we report a protein complex that reads DNA methylation to regulate variant-specific 45S ribosomal RNA (rRNA) gene expression in Arabidopsis (Arabidopsis thaliana). The complex, consisting of METHYL-CpG-BINDING DOMAIN PROTEIN5 (MBD5), MBD6, ALPHA-CRYSTALLIN DOMAIN PROTEIN15.5 (ACD15.5), and ACD21.4, directly binds to 45S rDNA. While MBD5 and MBD6 function redundantly, ACD15.5 and ACD21.4 are indispensable for variant-specific rRNA gene expression. These 4 proteins undergo phase separation in vitro and in vivo and are interdependent for their phase separation. The α-crystallin domain of ACD15.5 and ACD21.4, which is essential for their function, enables phase separation of the complex, likely by mediating multivalent protein interactions. The effector MICRORCHIDIA6 directly interacts with ACD15.5 and ACD21.4, but not with MBD5 and MBD6, and is recruited to 45S rDNA by the MBD-ACD complex to regulate variant-specific 45S rRNA expression. Our study reveals a pathway in Arabidopsis through which certain 45S rRNA gene variants are silenced, while others are activated.

Detection of hepatocellular carcinoma methylation markers in salivary DNA

BACKGROUND

Alterations to DNA methylation have been identified in both hepatocellular carcinoma (HCC) tumor and circulating DNA from affected individuals. These markers have potential utility in HCC screening. Adherence to HCC screening is poor and acceptable HCC screening tests are needed.

METHODS

A feasibility study was performed on a subset of case patients and control subjects from a prior study of risk factors for HCC. Case patients (n=12) included adults aged 47-85 years with a first diagnosis of HCC between 2011 and 2016 and without viral hepatitis. Control subjects (n=12) were matched on age, sex, and state of residence. Participants provided saliva samples for DNA genotyping. Log fold change in salivary DNA methylation at 1359 CpG sites representing 25 candidate genes previously associated with HCC was compared across case patients and control subjects.

RESULTS

The quantity of DNA ranged from 9.65 to 257.79 μg. The purity of DNA isolates was good, with mean OD260/280 ratio of 1.78 (SD: 0.14). Of 25 candidate genes, 16 had at ≥1 CpG site with detectable differences in methylation across HCC case patients and control subjects. Sites differentially methylated in HCC case patients included genes encoding tumor suppressors (PRDM2, RUNX3, p15/16, and RASSF1/5), regulators of cell cycle progression (DAPK1 and TP73), and DNA repair (MGMT and GSTP1). No associations met the significance threshold 3.7 × 10-5 required for multiple comparisons.

CONCLUSIONS

Salivary DNA may be a feasible alternative to blood samples in the era of novel DNA-based screening tests for HCC. The ease of saliva-based testing supports further investigation of its potential.

RNA m6A modification facilitates DNA methylation during maize kernel development

N6-methyladenosine (m6A) in mRNA and 5-methylcytosine (5mC) in DNA have critical functions for regulating gene expression and modulating plant growth and development. However, the interplay between m6A and 5mC is an elusive territory and remains unclear mechanistically in plants. We reported an occurrence of crosstalk between m6A and 5mC in maize (Zea mays) via the interaction between mRNA adenosine methylase (ZmMTA), the core component of the m6A methyltransferase complex, and decrease in DNA methylation 1 (ZmDDM1), a key chromatin-remodeling factor that regulates DNA methylation. Genes with m6A modification were coordinated with a much higher level of DNA methylation than genes without m6A modification. Dysfunction of ZmMTA caused severe arrest during maize embryogenesis and endosperm development, leading to a significant decrease in CHH methylation in the 5' region of m6A-modified genes. Instead, loss of function of ZmDDM1 had no noteworthy effects on ZmMTA-related activity. This study establishes a direct link between m6A and 5mC during maize kernel development and provides insights into the interplay between RNA modification and DNA methylation.

Genetic control of DNA methylation is largely shared across European and East Asian populations

DNA methylation is an ideal trait to study the extent of the shared genetic control across ancestries, effectively providing hundreds of thousands of model molecular traits with large QTL effect sizes. We investigate cis DNAm QTLs in three European (n = 3701) and two East Asian (n = 2099) cohorts to quantify the similarities and differences in the genetic architecture across populations. We observe 80,394 associated mQTLs (62.2% of DNAm probes with significant mQTL) to be significant in both ancestries, while 28,925 mQTLs (22.4%) are identified in only a single ancestry. mQTL effect sizes are highly conserved across populations, with differences in mQTL discovery likely due to differences in allele frequency of associated variants and differing linkage disequilibrium between causal variants and assayed SNPs. This study highlights the overall similarity of genetic control across ancestries and the value of ancestral diversity in increasing the power to detect associations and enhancing fine mapping resolution.

A novel imprinted locus on bovine chromosome 18 homologous with human chromosome 16q24.1

Genomic imprinting is an epigenetic regulation mechanism in mammals resulting in the parentally dependent monoallelic expression of genes. Imprinting disorders in humans are associated with several congenital syndromes and cancers and remain the focus of many medical studies. Cattle is a better model organism for investigating human embryo development than mice. Imprinted genes usually cluster on chromosomes and are regulated by different methylation regions (DMRs) located in imprinting control regions that control gene expression in cis. There is an imprinted locus on human chromosome 16q24.1 associated with congenital lethal developmental lung disease in newborns. However, genomic imprinting on bovine chromosome 18, which is homologous with human chromosome 16 has not been systematically studied. The aim of this study was to analyze the allelic expressions of eight genes (CDH13, ATP2C2, TLDC1, COTL1, CRISPLD2, ZDHHC7, KIAA0513, and GSE1) on bovine chromosome 18 and to search the DMRs associated gene allelic expression. Three transcript variants of the ZDHHC7 gene (X1, X2, and X5) showed maternal imprinting in bovine placentas. In addition, the monoallelic expression of X2 and X5 was tissue-specific. Five transcripts of the KIAA0513 gene showed tissue- and isoform-specific monoallelic expression. The CDH13, ATP2C2, and TLDC1 genes exhibited tissue-specific imprinting, however, COTL1, CRISLPLD2, and GSE1 escaped imprinting. Four DMRs, established after fertilization, were found in this region. Two DMRs were located between the ZDHHC7 and KIAA0513 genes, and two were in exon 1 of the CDH13 and ATP2C2 genes, respectively. The results from this study support future studies on the molecular mechanism to regulate the imprinting of candidate genes on bovine chromosome 18.

Altered DNA methylation underlies monocyte dysregulation and immune exhaustion memory in sepsis

Monocytes can develop an exhausted memory state characterized by reduced differentiation, pathogenic inflammation, and immune suppression that drives immune dysregulation during sepsis. Chromatin alterations, notably via histone modifications, underlie innate immune memory, but the contribution of DNA methylation remains poorly understood. Using an ex vivo sepsis model, we show altered DNA methylation throughout the genome of exhausted monocytes, including genes implicated in immune dysregulation during sepsis and COVID-19 infection (e.g., Plac8). These changes are recapitulated in septic mice induced by cecal slurry injection. Methylation profiles developed in septic mice are maintained during ex vivo culture, supporting the involvement of DNA methylation in stable monocyte exhaustion memory. Methylome reprogramming is driven in part by Wnt signaling inhibition in exhausted monocytes and can be reversed with DNA methyltransferase inhibitors, Wnt agonists, or immune training molecules. Our study demonstrates the significance of altered DNA methylation in the maintenance of stable monocyte exhaustion memory.

Low expression of miR-182 caused by DNA hypermethylation accelerates acute lymphocyte leukemia development by targeting PBX3 and BCL2: miR-182 promoter methylation is a predictive marker for hypomethylation agents + BCL2 inhibitor venetoclax

BACKGROUND

miR-182 promoter hypermethylation frequently occurs in various tumors, including acute myeloid leukemia, and leads to low expression of miR-182. However, whether adult acute lymphocyte leukemia (ALL) cells have high miR-182 promoter methylation has not been determined.

METHODS

To assess the methylation status of the miR-182 promoter, methylation and unmethylation-specific PCR analysis, bisulfite-sequencing analysis, and MethylTarget™ assays were performed to measure the frequency of methylation at the miR-182 promoter. Bone marrow cells were isolated from miR-182 knockout (182KO) and 182 wild type (182WT) mice to construct BCR-ABL (P190) and Notch-induced murine B-ALL and T-ALL models, respectively. Primary ALL samples were performed to investigate synergistic effects of the hypomethylation agents (HMAs) and the BCL2 inhibitor venetoclax (Ven) in vitro.

RESULTS

miR-182 (miR-182-5P) expression was substantially lower in ALL blasts than in normal controls (NCs) because of DNA hypermethylation at the miR-182 promoter in ALL blasts but not in normal controls (NCs). Knockout of miR-182 (182KO) markedly accelerated ALL development, facilitated the infiltration, and shortened the OS in a BCR-ABL (P190)-induced murine B-ALL model. Furthermore, the 182KO ALL cell population was enriched with more leukemia-initiating cells (CD43+B220+ cells, LICs) and presented higher leukemogenic activity than the 182WT ALL population. Furthermore, depletion of miR-182 reduced the OS in a Notch-induced murine T-ALL model, suggesting that miR-182 knockout accelerates ALL development. Mechanistically, overexpression of miR-182 inhibited proliferation and induced apoptosis by directly targeting PBX3 and BCL2, two well-known oncogenes, that are key targets of miR-182. Most importantly, DAC in combination with Ven had synergistic effects on ALL cells with miR-182 promoter hypermethylation, but not on ALL cells with miR-182 promoter hypomethylation.

CONCLUSIONS

Collectively, we identified miR-182 as a tumor suppressor gene in ALL cells and low expression of miR-182 because of hypermethylation facilitates the malignant phenotype of ALL cells. DAC + Ven cotreatment might has been applied in the clinical try for ALL patients with miR-182 promoter hypermethylation. Furthermore, the methylation frequency at the miR-182 promoter should be a potential biomarker for DAC + Ven treatment in ALL patients.

MOF-mediated acetylation of UHRF1 enhances UHRF1 E3 ligase activity to facilitate DNA methylation maintenance

The multi-domain protein UHRF1 (ubiquitin-like, containing PHD and RING finger domains, 1) recruits DNMT1 for DNA methylation maintenance during DNA replication. Here, we show that MOF (males absent on the first) acetylates UHRF1 at K670 in the pre-RING linker region, whereas HDAC1 deacetylates UHRF1 at the same site. We also identify that K667 and K668 can also be acetylated by MOF when K670 is mutated. The MOF/HDAC1-mediated acetylation in UHRF1 is cell-cycle regulated and peaks at G1/S phase, in line with the function of UHRF1 in recruiting DNMT1 to maintain DNA methylation. In addition, UHRF1 acetylation significantly enhances its E3 ligase activity. Abolishing UHRF1 acetylation at these sites attenuates UHRF1-mediated H3 ubiquitination, which in turn impairs DNMT1 recruitment and DNA methylation. Taken together, these findings identify MOF as an acetyltransferase for UHRF1 and define a mechanism underlying the regulation of DNA methylation maintenance through MOF-mediated UHRF1 acetylation.

A validated restriction enzyme ddPCR cg05575921 (AHRR) assay to accurately assess smoking exposure

BACKGROUND & METHODS

In this study, a novel restriction enzyme (RE) digestion-based droplet digital polymerase chain reaction (ddPCR) assay was designed for cg005575921 within the AHRR gene body and compared with matching results obtained by bisulfite conversion (BIS) ddPCR and Illumina DNA methylation array.

RESULTS

The RE ddPCR cg05575921 assay appeared concordant with BIS ddPCR (r2 = 0.94, P < 0.0001) and, when compared with the Illumina array, had significantly better smoking status classification performance for current versus never smoked (AUC 0.96 versus 0.93, P < 0.04) and current versus ex-smoker (AUC 0.88 versus 0.83, P < 0.04) comparisons.

CONCLUSIONS

The RE ddPCR cg05575921 assay accurately predicts smoking status and could be a useful component of 'precision-medicine' chronic disease risk screening tools.

A systematic review of high impact CpG sites and regions for MGMT methylation in glioblastoma [A systematic review of MGMT methylation in GBM]

BACKGROUND

MGMT (O 6 -methylguanine-DNA methyltransferase) promoter methylation is a commonly assessed prognostic marker in glioblastoma (GBM). Epigenetic silencing of the MGMT gene by promoter methylation is associated with greater overall and progression free survival with alkylating agent regimens. To date, there is marked heterogeneity in how MGMT promoter methylation is tested and which CpG sites are interrogated.

METHODS

To further elucidate which MGMT promoter CpG sites are of greatest interest, we performed comprehensive searches in PubMed, Web of Science, and Embase and reviewed 2,925 article abstracts. We followed the GRADE scoring system to assess risk of bias and the quality of the studies we included.

RESULTS

We included articles on adult glioblastoma that examined significant sites or regions within MGMT promoter for the outcomes: overall survival, progression free survival, and/or MGMT expression. We excluded systemic reviews and articles on lower grade glioma. fifteen articles met inclusion criteria with variable overlap in laboratory and statistical methods employed, as well as CpG sites interrogated. Pyrosequencing or BeadChip arrays were the most popular methods utilized, and CpG sites between CpG's 70-90 were most frequently investigated. Overall, there was moderate concordance between the CpG sites that the studies reported to be highly predictive of prognosis. Combinations or means of sites between CpG's 73-89 were associated with improved OS and PFS. Six studies identified CpG sites associated with prognosis that were closer to the transcription start site: CpG's 8, 19, 22, 25, 27, 32,38, and CpG sites 21-37, as well as low methylation level of the enhancer regions.

CONCLUSION

The following systematic review details a comprehensive investigation of the current literature and highlights several potential key CpG sites that demonstrate significant association with OS, PFS, and MGMT expression. However, the relationship between extent of MGMT promoter methylation and survival may be non-linear and could be influenced by potential CpG hotspots, the extent of methylation at each CpG site, and MGMT enhancer methylation status. There were several limitations within the studies such as smaller sample sizes, variance between methylation testing methods, and differences in the various statistical methods to test for association to outcome. Further studies of high impact CpG sites in MGMT methylation is warranted.

Predominance of MGMT promoter methylation among Pakistani glioblastoma patients

BACKGROUND

Glioblastoma multiforme (GBM), the most prevalent subgroup of neuroepithelial tumors, is characterized by dismal overall survival (OS). Several studies have linked O6-methylguanine-DNA-methyltransferase (MGMT) promoter methylation to OS in GBM patients. However, MGMT methylation frequencies vary geographically and across ethnicities, with limited data for South Asian populations, including Pakistan. This study aimed to analyze MGMT promoter methylation in Pakistani GBM patients.

METHODS

Consecutive primary GBM patients diagnosed ≥ 18 years-of-age, with no prior chemotherapy or radiotherapy history, were retrospectively selected. DNA was isolated from formalin-fixed-paraffin-embedded tissues. MGMT promoter methylation was analyzed using methylation-specific PCR. Clinical, pathological, and treatment data were assessed using Fisher's exact/Chi-squared tests. OS was calculated using Kaplan-Meier analysis in SPSS 27.0.1.

RESULTS

The study included 48 GBM patients, comprising 38 (79.2%) males and 10 (20.8%) females. The median diagnosis age was 49.5 years (range 18-70). MGMT methylation was observed in 87.5% (42/48) of all cases. Patients with MGMT methylation undergoing radiotherapy or radiotherapy plus chemotherapy exhibited significantly improved median OS of 7.2 months (95% CI, 3.7-10.7; P < 0.001) and 16.9 months (95% CI, 15.9-17.9; P < 0.001), respectively, compared to those undergoing surgical resection only (OS: 2.2 months, 95% CI, 0.8-3.6).

CONCLUSION

This is the first comprehensive study highlighting a predominance of MGMT methylation in Pakistani GBM patients. Furthermore, our findings underscore the association of MGMT methylation with improved OS across diverse treatment modalities. Larger studies are imperative to validate our findings for better management of Pakistani GBM patients.

A maternal-effect Padi6 variant causes nuclear and cytoplasmic abnormalities in oocytes, as well as failure of epigenetic reprogramming and zygotic genome activation in embryos

Maternal inactivation of genes encoding components of the subcortical maternal complex (SCMC) and its associated member, PADI6, generally results in early embryo lethality. In humans, SCMC gene variants were found in the healthy mothers of children affected by multilocus imprinting disturbances (MLID). However, how the SCMC controls the DNA methylation required to regulate imprinting remains poorly defined. We generated a mouse line carrying a Padi6 missense variant that was identified in a family with Beckwith-Wiedemann syndrome and MLID. If homozygous in female mice, this variant resulted in interruption of embryo development at the two-cell stage. Single-cell multiomic analyses demonstrated defective maturation of Padi6 mutant oocytes and incomplete DNA demethylation, down-regulation of zygotic genome activation (ZGA) genes, up-regulation of maternal decay genes, and developmental delay in two-cell embryos developing from Padi6 mutant oocytes but little effect on genomic imprinting. Western blotting and immunofluorescence analyses showed reduced levels of UHRF1 in oocytes and abnormal localization of DNMT1 and UHRF1 in both oocytes and zygotes. Treatment with 5-azacytidine reverted DNA hypermethylation but did not rescue the developmental arrest of mutant embryos. Taken together, this study demonstrates that PADI6 controls both nuclear and cytoplasmic oocyte processes that are necessary for preimplantation epigenetic reprogramming and ZGA.

Uncovering the roles of DNA hemi-methylation in transcriptional regulation using MspJI-assisted hemi-methylation sequencing

Hemi-methylated cytosine dyads widely occur on mammalian genomic DNA, and can be stably inherited across cell divisions, serving as potential epigenetic marks. Previous identification of hemi-methylation relied on harsh bisulfite treatment, leading to extensive DNA degradation and loss of methylation information. Here we introduce Mhemi-seq, a bisulfite-free strategy, to efficiently resolve methylation status of cytosine dyads into unmethylation, strand-specific hemi-methylation, or full-methylation. Mhemi-seq reproduces methylomes from bisulfite-based sequencing (BS-seq & hpBS-seq), including the asymmetric hemi-methylation enrichment flanking CTCF motifs. By avoiding base conversion, Mhemi-seq resolves allele-specific methylation and associated imprinted gene expression more efficiently than BS-seq. Furthermore, we reveal an inhibitory role of hemi-methylation in gene expression and transcription factor (TF)-DNA binding, and some displays a similar extent of inhibition as full-methylation. Finally, we uncover new hemi-methylation patterns within Alu retrotransposon elements. Collectively, Mhemi-seq can accelerate the identification of DNA hemi-methylation and facilitate its integration into the chromatin environment for future studies.

Functional characteristics of DNA N6-methyladenine modification based on long-read sequencing in pancreatic cancer

Abnormalities of DNA modifications are closely related to the pathogenesis and prognosis of pancreatic cancer. The development of third-generation sequencing technology has brought opportunities for the study of new epigenetic modification in cancer. Here, we screened the N6-methyladenine (6mA) and 5-methylcytosine (5mC) modification in pancreatic cancer based on Oxford Nanopore Technologies sequencing. The 6mA levels were lower compared with 5mC and upregulated in pancreatic cancer. We developed a novel method to define differentially methylated deficient region (DMDR), which overlapped 1319 protein-coding genes in pancreatic cancer. Genes screened by DMDRs were more significantly enriched in the cancer genes compared with the traditional differential methylation method (P < 0.001 versus P = 0.21, hypergeometric test). We then identified a survival-related signature based on DMDRs (DMDRSig) that stratified patients into high- and low-risk groups. Functional enrichment analysis indicated that 891 genes were closely related to alternative splicing. Multi-omics data from the cancer genome atlas showed that these genes were frequently altered in cancer samples. Survival analysis indicated that seven genes with high expression (ADAM9, ADAM10, EPS8, FAM83A, FAM111B, LAMA3 and TES) were significantly associated with poor prognosis. In addition, the distinction for pancreatic cancer subtypes was determined using 46 subtype-specific genes and unsupervised clustering. Overall, our study is the first to explore the molecular characteristics of 6mA modifications in pancreatic cancer, indicating that 6mA has the potential to be a target for future clinical treatment.

Co-methylation analyses identify CpGs associated with lipid traits in Chinese discordant monozygotic twins

To control genetic background and early life milieu in genome-wide DNA methylation analysis for blood lipids, we recruited Chinese discordant monozygotic twins to explore the relationships between DNA methylations and total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). 132 monozygotic (MZ) twins were included with discordant lipid levels and completed data. A linear mixed model was conducted in Epigenome-wide association study (EWAS). Generalized estimating equation model was for gene expression analysis. We conducted Weighted correlation network analysis (WGCNA) to build co-methylated interconnected network. Additional Qingdao citizens were recruited for validation. Inference about Causation through Examination of Familial Confounding (ICE FALCON) was used to infer the possible direction of these relationships. A total of 476 top CpGs reached suggestively significant level (P < 10-4), of which, 192 CpGs were significantly associated with TG (FDR < 0.05). They were used to build interconnected network and highlight crucial genes from WGCNA. Finally, four CpGs in GATA4 were validated as risk factors for TC; six CpGs at ITFG2-AS1 were negatively associated with TG; two CpGs in PLXND1 played protective roles in HDL-C. ICE FALCON indicated abnormal TC was regarded as the consequence of DNA methylation in CpGs at GATA4, rather than vice versa. Four CpGs in ITFG2-AS1 were both causes and consequences of modified TG levels. Our results indicated that DNA methylation levels of 12 CpGs in GATA4, ITFG2-AS1, and PLXND1 were relevant to TC, TG, and HDL-C, respectively, which might provide new epigenetic insights into potential clinical treatment of dyslipidemia.

Epigenetic MLH1 silencing concurs with mismatch repair deficiency in sporadic, naturally occurring colorectal cancer in rhesus macaques

BACKGROUND

Naturally occurring colorectal cancers (CRC) in rhesus macaques share many features with their human counterparts and are useful models for cancer immunotherapy; but mechanistic data are lacking regarding the comparative molecular pathogenesis of these cancers.

METHODS

We conducted state-of-the-art imaging including CT and PET, clinical assessments, and pathological review of 24 rhesus macaques with naturally occurring CRC. Additionally, we molecularly characterized these tumors utilizing immunohistochemistry (IHC), microsatellite instability assays, DNAseq, transcriptomics, and developed a DNA methylation-specific qPCR assay for MLH1, CACNA1G, CDKN2A, CRABP1, and NEUROG1, human markers for CpG island methylator phenotype (CIMP). We furthermore employed Monte-Carlo simulations to in-silico model alterations in DNA topology in transcription-factor binding site-rich promoter regions upon experimentally demonstrated DNA methylation.

RESULTS

Similar cancer histology, progression patterns, and co-morbidities could be observed in rhesus as reported for human CRC patients. IHC identified loss of MLH1 and PMS2 in all cases, with functional microsatellite instability. DNA sequencing revealed the close genetic relatedness to human CRCs, including a similar mutational signature, chromosomal instability, and functionally-relevant mutations affecting KRAS (G12D), TP53 (R175H, R273*), APC, AMER1, ALK, and ARID1A. Interestingly, MLH1 mutations were rarely identified on a somatic or germline level. Transcriptomics not only corroborated the similarities of rhesus and human CRCs, but also demonstrated the significant downregulation of MLH1 but not MSH2, MSH6, or PMS2 in rhesus CRCs. Methylation-specific qPCR suggested CIMP-positivity in 9/16 rhesus CRCs, but all 16/16 exhibited significant MLH1 promoter hypermethylation. DNA hypermethylation was modelled to affect DNA topology, particularly propeller twist and roll profiles. Modelling the DNA topology of a transcription factor binding motif (TFAP2A) in the MLH1 promoter that overlapped with a methylation-specific probe, we observed significant differences in DNA topology upon experimentally shown DNA methylation. This suggests a role of transcription factor binding interference in epigenetic silencing of MLH1 in rhesus CRCs.

CONCLUSIONS

These data indicate that epigenetic silencing suppresses MLH1 transcription, induces the loss of MLH1 protein, abrogates mismatch repair, and drives genomic instability in naturally occurring CRC in rhesus macaques. We consider this spontaneous, uninduced CRC in immunocompetent, treatment-naïve rhesus macaques to be a uniquely informative model for human CRC.

Single-cell transcriptomic sequencing data reveal aberrant DNA methylation in SMAD3 promoter region in tumor-associated fibroblasts affecting molecular mechanism of radiosensitivity in non-small cell lung cancer

OBJECTIVE

Non-small cell lung cancer (NSCLC) often exhibits resistance to radiotherapy, posing significant treatment challenges. This study investigates the role of SMAD3 in NSCLC, focusing on its potential in influencing radiosensitivity via the ITGA6/PI3K/Akt pathway.

METHODS

The study utilized gene expression data from the GEO database to identify differentially expressed genes related to radiotherapy resistance in NSCLC. Using the GSE37745 dataset, prognostic genes were identified through Cox regression and survival analysis. Functional roles of target genes were explored using Gene Set Enrichment Analysis (GSEA) and co-expression analyses. Gene promoter methylation levels were assessed using databases like UALCAN, DNMIVD, and UCSC Xena, while the TISCH database provided insights into the correlation between target genes and CAFs. Experiments included RT-qPCR, Western blot, and immunohistochemistry on NSCLC patient samples, in vitro studies on isolated CAFs cells, and in vivo nude mouse tumor models.

RESULTS

Fifteen key genes associated with radiotherapy resistance in NSCLC cells were identified. SMAD3 was recognized as an independent prognostic factor for NSCLC, linked to poor patient outcomes. High expression of SMAD3 was correlated with low DNA methylation in its promoter region and was enriched in CAFs. In vitro and in vivo experiments confirmed that SMAD3 promotes radiotherapy resistance by activating the ITGA6/PI3K/Akt signaling pathway.

CONCLUSION

High expression of SMAD3 in NSCLC tissues, cells, and CAFs is closely associated with poor prognosis and increased radiotherapy resistance. SMAD3 is likely to enhance radiotherapy resistance in NSCLC cells by activating the ITGA6/PI3K/Akt signaling pathway.

ICBP90, an epigenetic regulator, induces DKK3 promoter methylation, promotes glioma progression, and reduces sensitivity to cis-platinum

Glioma is the most common brain malignancy, characterized by high morbidity, high mortality, and treatment-resistance. Inverted CCAAT box Binding Protein of 90 kDa (ICBP90) has been reported to be involved in tumor progression and the maintenance of DNA methylation. Herein, we constructed ICBP90 over-expression and knockdown glioma cell lines, and found that ICBP90 knockdown inhibited glioma cell proliferation, migration, and invasion. ICBP90 silencing potentially enhanced cellular sensitivity to cis-platinum (DDP) and exacerbated DDP-induced pyroptosis, manifested by the elevated levels of gasdermin D-N-terminal and cleaved caspase 1; whereas, ICBP90 over-expression exhibited the opposite effects. Consistently, ICBP90 knockdown inhibited tumor growth in an in vivo mouse xenograft study using U251 cells stably expressing sh-ICBP90 and oe-ICBP90. Further experiments found that ICBP90 reduced the expression of Dickkopf 3 homolog (DKK3), a negative regulator of β-catenin, by binding its promoter and inducing DNA methylation. ICBP90 knockdown prevented the nuclear translocation of β-catenin and suppressed the expression of c-Myc and cyclin D1. Besides, DKK3 over-expression restored the effects of ICBP90 over-expression on cell proliferation, migration, invasion, and DDP sensitivity. Our findings suggest that ICBP90 inhibits the expression of DKK3 in glioma by maintaining DKK3 promoter methylation, thereby conducing to ICBP90-mediated carcinogenesis and drug insensitivity.

Barley AGO4 proteins show overlapping functionality with distinct small RNA-binding properties in heterologous complementation

KEY MESSAGE

Barley AGO4 proteins complement expressional changes of epigenetically regulated genes in Arabidopsis ago4-3 mutant and show a distinct affinity for the 5' terminal nucleotide of small RNAs, demonstrating functional conservation and divergence. The function of Argonaute 4 (AGO4) in Arabidopsis thaliana has been extensively characterized; however, its role in monocots, which have large genomes abundantly supplemented with transposable elements (TEs), remains elusive. The study of barley AGO4 proteins can provide insights into the conserved aspects of RNA-directed DNA methylation (RdDM) and could also have further applications in the field of epigenetics or crop improvement. Bioinformatic analysis of RNA sequencing data identified two active AGO4 genes in barley, HvAGO4a and HvAGO4b. These genes function similar to AtAGO4 in an Arabidopsis heterologous complementation system, primarily binding to 24-nucleotide long small RNAs (sRNAs) and triggering methylation at specific target loci. Like AtAGO4, HvAGO4B exhibits a preference for binding sRNAs with 5' adenine residue, while also accepting 5' guanine, uracil, and cytosine residues. In contrast, HvAGO4A selectively binds only sRNAs with a 5' adenine residue. The diverse binding capacity of barley AGO4 proteins is reflected in TE-derived sRNAs and in their varying abundance. Both barley AGO4 proteins effectively restore the levels of extrachromosomal DNA and transcript abundancy of the heat-activated ONSEN retrotransposon to those observed in wild-type Arabidopsis plants. Our study provides insight into the distinct binding specificities and involvement in TE regulation of barley AGO4 proteins in Arabidopsis by heterologous complementation.

Human papillomavirus (HPV) DNA methylation changes in HPV-associated head and neck cancer

Despite the rising incidence, currently, there are no early detection methods for HPV-driven HNC (HPV-HNC). Cervical cancer studies suggest that HPV DNA methylation changes can be used as a biomarker to discriminate cancer patients from HPV-infected individuals. As such, this study was designed to establish a protocol to evaluate DNA methylation changes in HPV late genes and long control region (LCR) in saliva samples of HPV-HNC patients and HPV-positive controls. Higher methylation levels were detected in HPV late genes (L1 and L2) in both tumour and saliva samples of HPV-HNC patients compared with HPV-positive controls. Moreover, methylation patterns between tumours and corresponding saliva samples were observed to have a strong correlation (Passing-Bablok regression analysis; τ = 0.7483, P < 0.0001). Considering the differences between HNC and controls in methylation levels in late genes, and considering primer amplification efficiencies, 13 CpG sites located at L1 and L2 genes were selected for further evaluation. A total of 18 HNC saliva samples and 10 control saliva samples were assessed for the methylation levels in the selected sites. From the CpG sites evaluated statistically significant differences were identified for CpG sites at L2-CpG 6 (P = 0.0004), L1-CpG 3 (P = 0.0144), L1-CpG 2 (P = 0.0395) and L2-CpG 19 (P = 0.0455). Our pilot data indicate that higher levels of DNA methylation in HPV late genes are indicative of HPV-HNC risk, and it is a potential supplementary biomarker for salivary HPV detection-based HPV-HNC screening.

Alternative polyadenylation quantitative trait methylation mapping in human cancers provides clues into the molecular mechanisms of APA

Genetic variants are involved in the orchestration of alternative polyadenylation (APA) events, while the role of DNA methylation in regulating APA remains unclear. We generated a comprehensive atlas of APA quantitative trait methylation sites (apaQTMs) across 21 different types of cancer (1,612 to 60,219 acting in cis and 4,448 to 142,349 in trans). Potential causal apaQTMs in non-cancer samples were also identified. Mechanistically, we observed a strong enrichment of cis-apaQTMs near polyadenylation sites (PASs) and both cis- and trans-apaQTMs in proximity to transcription factor (TF) binding regions. Through the integration of ChIP-signals and RNA-seq data from cell lines, we have identified several regulators of APA events, acting either directly or indirectly, implicating novel functions of some important genes, such as TCF7L2, which is known for its involvement in type 2 diabetes and cancers. Furthermore, we have identified a vast number of QTMs that share the same putative causal CpG sites with five different cancer types, underscoring the roles of QTMs, including apaQTMs, in the process of tumorigenesis. DNA methylation is extensively involved in the regulation of APA events in human cancers. In an attempt to elucidate the potential underlying molecular mechanisms of APA by DNA methylation, our study paves the way for subsequent experimental validations into the intricate biological functions of DNA methylation in APA regulation and the pathogenesis of human cancers. To present a comprehensive catalog of apaQTM patterns, we introduce the Pancan-apaQTM database, available at https://pancan-apaqtm-zju.shinyapps.io/pancanaQTM/.

Epigenome-wide association study of total nicotine equivalents in multiethnic current smokers from three prospective cohorts

The impact of tobacco exposure on health varies by race and ethnicity and is closely tied to internal nicotine dose, a marker of carcinogen uptake. DNA methylation is strongly responsive to smoking status and may mediate health effects, but study of associations with internal dose is limited. We performed a blood leukocyte epigenome-wide association study (EWAS) of urinary total nicotine equivalents (TNEs; a measure of nicotine uptake) and DNA methylation measured using the MethylationEPIC v1.0 BeadChip (EPIC) in six racial and ethnic groups across three cohort studies. In the Multiethnic Cohort Study (discovery, n = 1994), TNEs were associated with differential methylation at 408 CpG sites across >250 genomic regions (p < 9 × 10-8). The top significant sites were annotated to AHRR, F2RL3, RARA, GPR15, PRSS23, and 2q37.1, all of which had decreasing methylation with increasing TNEs. We identified 45 novel CpG sites, of which 42 were unique to the EPIC array and eight annotated to genes not previously linked with smoking-related DNA methylation. The most significant signal in a novel gene was cg03748458 in MIR383;SGCZ. Fifty-one of the 408 discovery sites were validated in the Singapore Chinese Health Study (n = 340) and the Southern Community Cohort Study (n = 394) (Bonferroni corrected p < 1.23 × 10-4). Significant heterogeneity by race and ethnicity was detected for CpG sites in MYO1G and CYTH1. Furthermore, TNEs significantly mediated the association between cigarettes per day and DNA methylation at 15 sites (average 22.5%-44.3% proportion mediated). Our multiethnic study highlights the transethnic and ethnic-specific methylation associations with internal nicotine dose, a strong predictor of smoking-related morbidities.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation provides a crucial epigenetic mark linking genetic variations to environmental influence. We have analyzed array-based DNA methylation profiles of 160 human retinas with co-measured RNA-seq and >8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 methylation quantitative trait loci and 12,505 expression quantitative trait loci) and 13,747 DNA methylation loci affecting gene expression, with over one-third specific to the retina. Methylation and expression quantitative trait loci show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration. Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of macular degeneration pathology by genotype-environment interaction in retina.

Unmethylated Mosaic Full Mutation Males without Fragile X Syndrome

Fragile X syndrome (FXS) is the leading inherited cause of intellectual disability (ID) and single gene cause of autism. Although most patients with FXS and the full mutation (FM) have complete methylation of the fragile X messenger ribonucleoprotein 1 (FMR1) gene, some have mosaicism in methylation and/or CGG repeat size, and few have completely unmethylated FM alleles. Those with a complete lack of methylation are rare, with little literature about the cognitive and behavioral phenotypes of these individuals. A review of past literature was conducted regarding individuals with unmethylated and mosaic FMR1 FM. We report three patients with an unmethylated FM FMR1 alleles without any behavioral or cognitive deficits. This is an unusual presentation for men with FM as most patients with an unmethylated FM and no behavioral phenotypes do not receive fragile X DNA testing or a diagnosis of FXS. Our cases showed that mosaic males with unmethylated FMR1 FM alleles may lack behavioral phenotypes due to the presence of smaller alleles producing the FMR1 protein (FMRP). However, these individuals could be at a higher risk of developing fragile X-associated tremor/ataxia syndrome (FXTAS) due to the increased expression of mRNA, similar to those who only have a premutation.

The role of promoter methylation of the genes encoding the enzymes metabolizing di- and tricarboxylic acids in the regulation of plant respiration by light

We discuss the role of epigenetic changes at the level of promoter methylation of the key enzymes of carbon metabolism in the regulation of respiration by light. While the direct regulation of enzymes via modulation of their activity and post-translational modifications is fast and readily reversible, the role of cytosine methylation is important for providing a prolonged response to environmental changes. In addition, adenine methylation can play a role in the regulation of transcription of genes. The mitochondrial and extramitochondrial forms of several enzymes participating in the tricarboxylic acid cycle and associated reactions are regulated via promoter methylation in opposite ways. The mitochondrial forms of citrate synthase, aconitase, fumarase, NAD-malate dehydrogenase are inhibited while the cytosolic forms of aconitase, fumarase, NAD-malate dehydrogenase, and the peroxisomal form of citrate synthase are activated. It is concluded that promoter methylation represents a universal mechanism of the regulation of activity of respiratory enzymes in plant cells by light. The role of the regulation of the mitochondrial and cytosolic forms of respiratory enzymes in the operation of malate and citrate valves and in controlling the redox state and balancing the energy level of photosynthesizing plant cells is discussed.

Reconstruction of cell-specific models capturing the influence of metabolism on DNA methylation in cancer

The imbalance of epigenetic regulatory mechanisms such as DNA methylation, which can promote aberrant gene expression profiles without affecting the DNA sequence, may cause the deregulation of signaling, regulatory, and metabolic processes, contributing to a cancerous phenotype. Since some metabolites are substrates and cofactors of epigenetic regulators, their availability can be affected by characteristic cancer cell metabolic shifts, feeding cancer onset and progression through epigenetic deregulation. Hence, there is a need to study the influence of cancer metabolic reprogramming in DNA methylation to design new effective treatments. In this study, a generic Genome-Scale Metabolic Model (GSMM) of a human cell, integrating DNA methylation or demethylation reactions, was obtained and used for the reconstruction of Genome-Scale Metabolic Models enhanced with Enzymatic Constraints using Kinetic and Omics data (GECKOs) of 31 cancer cell lines. Furthermore, cell-line-specific DNA methylation levels were included in the models, as coefficients of a DNA composition pseudo-reaction, to depict the influence of metabolism over global DNA methylation in each of the cancer cell lines. Flux simulations demonstrated the ability of these models to provide simulated fluxes of exchange reactions similar to the equivalent experimentally measured uptake/secretion rates and to make good functional predictions. In addition, simulations found metabolic pathways, reactions and enzymes directly or inversely associated with the gene promoter methylation. Two potential candidates for targeted cancer epigenetic therapy were identified.

High-throughput methylome analysis reveals differential methylation for early and late onset preeclampsia for mothers and their children

Preeclampsia is a hypertensive disorder of pregnancy that affects ∼2%-5% of all pregnancies, contributes to 4 of the top 10 causes of pregnancy-related deaths, and remains a long-term risk factor for cardiometabolic diseases. Yet, little is still known about the molecular mechanisms that lead to this disease. There is evidence that some cases have a genetic cause. However, it is well appreciated that harmful factors in the environment, such as poor nutrition, stress, and toxins, may lead to epigenetics changes that can contribute to this disease. DNA methylation is one of the epigenetic modifications known to be fairly stable and impact gene expression. Using DNA from buccal swabs, we analyzed global DNA methylation among three groups of individuals: mothers who experienced 1) early-stage preeclampsia (<32 wk), 2) late-stage preeclampsia (>37 wk), or 3) no complications during their pregnancies, as well as the children from these three groups. We found significant differentially methylated regions (DMRs) between mothers who experienced preeclampsia compared with those with no complications adjacent or within genes that are important for placentation, embryonic development, cell adhesion, and inflammation (e.g., the cadherin pathway). A significant portion of DMR genes showed expression in tissues relevant to preeclampsia (i.e., the brain, heart, kidney, uterus, ovaries, and placenta). As this study was performed on DNA extracted from cheek swabs, this opens the way to future studies in different tissues, aimed at identifying possible biomarkers of risk and early detection, developing targeted interventions, and reducing the progression of this life-threatening disease.NEW & NOTEWORTHY Preeclampsia is a life-threatening hypertensive disorder, affecting 2%-5% of pregnancies, that remains poorly understood. This study analyzed DNA methylation from buccal swabs from mothers who experienced early and late-stage preeclampsia and those with uncomplicated pregnancies, along with their children. Differentially methylated regions were found near and within genes crucial for placental function, embryonic development, and inflammation. Many of these genes are expressed in preeclampsia-related tissues, offering hope for future biomarker development for this condition.

Polygenic risk for schizophrenia, social dispositions, and pace of epigenetic aging: Results from the Young Finns Study

Schizophrenia is often regarded as a disorder of premature aging. We investigated (a) whether polygenic risk for schizophrenia (PRSsch ) relates to pace of epigenetic aging and (b) whether personal dispositions toward active and emotionally close relationships protect against accelerated epigenetic aging in individuals with high PRSsch . The sample came from the population-based Young Finns Study (n = 1348). Epigenetic aging was measured with DNA methylation aging algorithms such as AgeAccelHannum , EEAAHannum , IEAAHannum , IEAAHorvath , AgeAccelHorvath , AgeAccelPheno , AgeAccelGrim , and DunedinPACE. A PRSsch was calculated using summary statistics from the most comprehensive genome-wide association study of schizophrenia to date. Social dispositions were assessed in terms of extraversion, sociability, reward dependence, cooperativeness, and attachment security. We found that PRSsch did not have a statistically significant effect on any studied indicator of epigenetic aging. Instead, PRSsch had a significant interaction with reward dependence (p = 0.001-0.004), cooperation (p = 0.009-0.020), extraversion (p = 0.019-0.041), sociability (p = 0.003-0.016), and attachment security (p = 0.007-0.014) in predicting AgeAccelHannum , EEAAHannum , or IEAAHannum . Specifically, participants with high PRSsch appeared to display accelerated epigenetic aging at higher (vs. lower) levels of extraversion, sociability, attachment security, reward dependence, and cooperativeness. A rather opposite pattern was evident for those with low PRSsch . No such interactions were evident when predicting the other indicators of epigenetic aging. In conclusion, against our hypothesis, frequent social interactions may relate to accelerated epigenetic aging in individuals at risk for psychosis. We speculate that this may be explained by social-cognitive impairments (perceiving social situations as overwhelming or excessively arousing) or ending up in less supportive or deviant social groups.

The impact of parental and developmental stress on DNA methylation in the avian hypothalamic-pituitary-adrenal axis

The hypothalamic-pituitary-adrenal (HPA) axis coordinates an organism's response to environmental stress. The responsiveness and sensitivity of an offspring's stress response may be shaped not only by stressors encountered in their early post-natal environment but also by stressors in their parent's environment. Yet, few studies have considered how stressors encountered in both of these early life environments may function together to impact the developing HPA axis. Here, we manipulated stressors in the parental and post-natal environments in a population of house sparrows (Passer domesticus) to assess their impact on changes in DNA methylation (and corresponding gene expression) in a suite of genes within the HPA axis. We found that nestlings that experienced early life stress across both life-history periods had higher DNA methylation in a critical HPA axis gene, the glucocorticoid receptor (NR3C1). In addition, we found that the life-history stage when stress was encountered impacted some genes (HSD11B1, NR3C1 and NR3C2) differently. We also found evidence for the mitigation of parental stress by post-natal stress (in HSD11B1 and NR3C2). Finally, by assessing DNA methylation in both the brain and blood, we were able to evaluate cross-tissue patterns. While some differentially methylated regions were tissue-specific, we found cross-tissue changes in NR3C2 and NR3C1, suggesting that blood is a suitable tissue for assessing DNA methylation as a biomarker of early life stress. Our results provide a crucial first step in understanding the mechanisms by which early life stress in different life-history periods contributes to changes in the epigenome of the HPA axis.

Identification of DNA methylation episignature for the intellectual developmental disorder, autosomal dominant 21 syndrome, caused by variants in the CTCF gene

PURPOSE

The main objective of this study was to assess clinical features and genome-wide DNA methylation profiles in individuals affected by intellectual developmental disorder, autosomal dominant 21 (IDD21) syndrome, caused by variants in the CCCTC-binding factor (CTCF) gene.

METHODS

DNA samples were extracted from peripheral blood of 16 individuals with clinical features and genetic findings consistent with IDD21. DNA methylation analysis was performed using the Illumina Infinium Methylation EPIC Bead Chip microarrays. The methylation levels were fitted in a multivariate linear regression model to identify the differentially methylated probes. A binary support vector machine classification model was constructed to differentiate IDD21 samples from controls.

RESULTS

We identified a highly specific, reproducible, and sensitive episignature associated with CTCF variants. Six variants of uncertain significance were tested, of which 2 mapped to the IDD21 episignature and clustered alongside IDD21 cases in both heatmap and multidimensional scaling plots. Comparison of the genomic DNA methylation profile of IDD21 with that of 56 other neurodevelopmental disorders provided insights into the underlying molecular pathophysiology of this disorder.

CONCLUSION

The robust and specific CTCF/IDD21 episignature expands the growing list of neurodevelopmental disorders with distinct DNA methylation profiles, which can be applied as supporting evidence in variant classification.

Detecting DNA hydroxymethylation: exploring its role in genome regulation

DNA methylation is one of the most extensively studied epigenetic regulatory mechanisms, known to play crucial roles in various organisms. It has been implicated in the regulation of gene expression and chromatin changes, ranging from global alterations during cell state transitions to locus-specific modifications. 5-hydroxymethylcytosine (5hmC) is produced by a major oxidation, from 5-methylcytosine (5mC), catalyzed by the ten-eleven translocation (TET) enzymes, and is gradually being recognized for its significant role in genome regulation. With the development of state-of-the-art experimental techniques, it has become possible to detect and distinguish 5mC and 5hmC at base resolution. Various techniques have evolved, encompassing chemical and enzymatic approaches, as well as thirdgeneration sequencing techniques. These advancements have paved the way for a thorough exploration of the role of 5hmC across a diverse array of cell types, from embryonic stem cells (ESCs) to various differentiated cells. This review aims to comprehensively report on recent techniques and discuss the emerging roles of 5hmC. [BMB Reports 2024; 57(3): 135-142].

The regulation of methylation on the Z chromosome and the identification of multiple novel Male Hyper-Methylated regions in the chicken

DNA methylation is a key regulator of eukaryote genomes, and is of particular relevance in the regulation of gene expression on the sex chromosomes, with a key role in dosage compensation in mammalian XY systems. In the case of birds, dosage compensation is largely absent, with it being restricted to two small Male Hyper-Methylated (MHM) regions on the Z chromosome. To investigate how variation in DNA methylation is regulated on the Z chromosome we utilised a wild x domestic advanced intercross in the chicken, with both hypothalamic methylomes and transcriptomes assayed in 124 individuals. The relatively large numbers of individuals allowed us to identify additional genomic MHM regions on the Z chromosome that were significantly differentially methylated between the sexes. These regions appear to down-regulate local gene expression in males, but not remove it entirely (unlike the lncRNAs identified in the initial MHM regions). These MHM regions were further tested and the most balanced genes appear to show decreased expression in males, whilst methylation appeared to be far more correlated with gene expression in the less balanced, as compared to the most balanced genes. In addition, quantitative trait loci (QTL) that regulate variation in methylation on the Z chromosome, and those loci that regulate methylation on the autosomes that derive from the Z chromosome were mapped. Trans-effect hotspots were also identified that were based on the autosomes but affected the Z, and also one that was based on the Z chromosome but that affected both autosomal and sex chromosome DNA methylation regulation. We show that both cis and trans loci that originate from the Z chromosome never exhibit an interaction with sex, whereas trans loci originating from the autosomes but affecting the Z chromosome always display such an interaction. Our results highlight how additional MHM regions are actually present on the Z chromosome, and they appear to have smaller-scale effects on gene expression in males. Quantitative variation in methylation is also regulated both from the autosomes to the Z chromosome, and from the Z chromosome to the autosomes.

Aberrant PRDM2 methylation as an early event in serrated lesions destined to evolve into microsatellite-instable colorectal cancers

Up to 30% of colorectal cancers (CRCs) develop from sessile serrated lesions (SSLs). Within the serrated neoplasia pathway, at least two principally distinct oncogenetic routes exist generating microsatellite-stable and microsatellite-instable CRCs, respectively. Aberrant DNA methylation (DNAm) is found early in the serrated pathway and might play a role in both oncogenetic routes. We studied a cohort of 23 SSLs with a small focus (<10 mm) of dysplasia or cancer, 10 of which were MLH1 deficient and 13 MLH1 proficient. By comparing, for each SSL, the methylation status of (1) the region of dysplasia or cancer (SSL-D), (2) the nondysplastic SSL (SSL), and (3) adjacent normal mucosa, differentially methylated probes (DMPs) and regions (DMRs) were assessed both genome-wide as well as in a tumor-suppressor gene-focused approach. By comparing DNAm of MLH1-deficient SSL-Ds with their corresponding SSLs, we identified five DMRs, including those annotating for PRDM2 and, not unexpectedly, MLH1. PRDM2 gene promotor methylation was associated with MLH1 expression status, as it was largely hypermethylated in MLH1-deficient SSL-Ds and hypomethylated in MLH1-proficient SSL-Ds. Significantly increased DNAm levels of PRDM2 and MLH1, in particular at 'critical' MLH1 probe sites, were to some extent already visible in SSLs as compared to normal mucosa (p = 0.02, p = 0.01, p < 0.0001, respectively). No DMRs, nor DMPs, were identified for SSLs destined to evolve into MLH1-proficient SSL-Ds. Our data indicate that, within both arms of the serrated CRC pathway, the majority of the epigenetic alterations are introduced early during SSL formation. Promoter hypermethylation of PRDM2 and MLH1 on the other hand specifically initiates in SSLs destined to transform into MLH1-deficient CRCs suggesting that the fate of SSLs may not necessarily result from a stochastic process but possibly is already imprinted and predisposed.

DeepSF-4mC: A deep learning model for predicting DNA cytosine 4mC methylation sites leveraging sequence features

N4-methylcytosine (4mC) is a DNA modification involving the addition of a methyl group to the fourth nitrogen atom of the cytosine base. This modification may influence gene regulation, providing potential insights into gene control mechanisms. Traditional laboratory methods for detecting 4mC DNA methylation have limitations, but the rise of artificial intelligence has introduced efficient computational strategies for 4mC site prediction. Despite this progress, challenges persist in terms of model performance and interpretability. To tackle these challenges, we propose DeepSF-4mC, a deep learning model specifically designed for predicting DNA cytosine 4mC methylation sites by leveraging sequence features. Our approach incorporates multiple encoding techniques to enhance prediction accuracy, increase model stability, and reduce the computational resources needed. Leveraging transfer learning, we harness existing models to enhance performance through learned representations or fine-tuning. Ensemble learning techniques combine predictions from multiple models, boosting robustness and accuracy. This research contributes to DNA methylation analysis and lays the groundwork for understanding 4mC's multifaceted role in biological processes. The web server for DeepSF-4mC is accessible at: http://deepsf-4mc.top/and the original code can be found at: https://github.com/754131799/DeepSF-4mC.