Sulfonamide-induced DNA hypomethylation disturbed sugar metabolism in rice (Oryza sativa L.)

[Display omitted] •Sugar metabolism and accumulation can be significantly disturbed by sulfonamides.•DNA hypomethylation promoted sugar hydrolysis to favor sulfonamide detoxification.•Sulfonamides deformed DNA and inhibited its methylation by minor groove binding.•Sulfadiazine without –CH3 induced more serious metabolic and epigenetic disorders. DNA methylation is well-accepted as a bridge to unravel the complex interplay between genome and environmental exposures, and its alteration regulated the cellular metabolic responses towards pollutants. However, the mechanism underlying site-specific aberrant DNA methylation and metabolic disorders under pollutant stresses remained elusive. Herein, the multilevel omics interferences of sulfonamides (i.e., sulfadiazine and sulfamerazine), a group of antibiotics pervasive in farmland soils, towards rice in 14 days of 1 mg/L hydroponic exposure were systematically evaluated. Metabolome and transcriptome analyses showed that 57.1–71.4 % of mono- and disaccharides were accumulated, and the differentially expressed genes were involved in the promotion of sugar hydrolysis, as well as the detoxification of sulfonamides. Most differentially methylated regions (DMRs) were hypomethylated ones (accounting for 87–95 %), and 92 % of which were located in the CHH context (H = A, C, or T base). KEGG enrichment analysis revealed that CHH-DMRs in the promoter regions were enriched in sugar metabolism. To reveal the significant hypomethylation of CHH, multi-spectroscopic and thermodynamic approaches, combined with molecular simulation were conducted to investigate the molecular interaction between sulfonamides and DNA in different sequence contexts, and the result demonstrated that sulfonamides would insert into the minor grooves of DNA, and exhibited a stronger affinity with the CHH contexts of DNA compared to CG or CHG contexts. Computational modeling of DNA 3D structures further confirmed that the binding led to a pitch increase of 0.1 Å and a 3.8° decrease in the twist angle of DNA in the CHH context. This specific interaction and the downregulation of methyltransferase CMT2 (log2FC = –4.04) inhibited the DNA methylation. These results indicated that DNA methylation-based assessment was useful for metabolic toxicity prediction and health risk assessment.