Comparative genome-wide DNA methylation analysis reveals epigenomic differences in response to heat-humidity stress in Bombyx mori

DNA methylation is an important epigenetic modification and has been shown to be involved in the response to abiotic stress. However, there are few studies on DNA methylation in insect response to environmental signals. In this study, we conducted a comprehensive comparative analysis of DNA methylation profiles between two silkworm strains with significantly different resistance to heat and humidity by whole-genome bisulfite sequencing (WGBS). We identified, in total, 2934 differentially methylated regions (DMRs) between RT_48h (resistant strain with high-temperature/humidity treatment for 48 h) and ST_48h (sensitive strain with high-temperature/humidity treatment for 48 h) under cytosine context (CG), which corresponded to 1230 DMR-related genes (DMGs), and the DMRs were primarily located in the gene body (exon and intron) region. Gene ontology (GO) and KEGG analysis showed that these DMGs were most significantly enriched in binding, cellular metabolic process, and RNA transport pathways. Moreover, 10 DMGs have been revealed to be involved in the heat-humidity stress response in the silkworm. The results of this study indicated that DNA methylation plays crucial roles in silkworm response to environmental stressors and provides important clues to identify key resistance genes in silkworm under high-temperature/humidity stress response.