Abstract:Abstract Background Bean pyralid is one of the major leaf-feeding insects that affect soybean crops. DNA methylation can control the networks of gene expressions, and it plays an important role in responses to biotic stress. However, at present the genome-wide DNA methylation profile of the soybean resistance to bean pyralid has not been reported so far. Results Using whole-genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq), we analyzed the highly resistant material (Gantai-2-2, HRK) and highly susceptible material (Wan82–178, HSK), under bean pyralid larvae feeding 0 h and 48 h, to clarify the molecular mechanism of the soybean resistance and explore its insect-resistant genes. We identified 2194, 6872, 39,704 and 40,018 differentially methylated regions (DMRs), as well as 497, 1594, 9596 and 9554 differentially methylated genes (DMGs) in the HRK0/HRK48, HSK0/HSK48, HSK0/HRK0 and HSK48/HRK48 comparisons, respectively. Through the analysis of global methylation and transcription, 265 differentially expressed genes (DEGs) were negatively correlated with DMGs, there were 34, 49, 141 and 116 negatively correlated genes in the HRK0/HRK48, HSK0/HSK48, HSK0/HRK0 and HSK48/HRK48, respectively. The MapMan cluster analysis showed that 114 negatively correlated genes were clustered in 24 pathways, such as protein biosynthesis and modification; primary metabolism; secondary metabolism; cell cycle, cell structure and component; RNA biosynthesis and processing, and so on. Moreover, CRK40; CRK62; STK; MAPK9; L-type lectin-domain containing receptor kinase VIII.2; CesA; CSI1; fimbrin-1; KIN-14B; KIN-14 N; KIN-4A; cytochrome P450 81E8; BEE1; ERF; bHLH25; bHLH79; GATA26, were likely regulatory genes involved in the soybean responses to bean pyralid larvae. Finally, 5 DMRs were further validated that the genome-wide DNA data were reliable through PS-PCR and 5 DEGs were confirmed the relationship between DNA methylation and gene expression by qRT-PCR. The results showed an excellent agreement with deep sequencing. Conclusions Genome-wide DNA methylation profile of soybean response to bean pyralid was obtained for the first time. Several specific DMGs which participated in protein kinase, cell and organelle, flavonoid biosynthesis and transcription factor were further identified to be likely associated with soybean response to bean pyralid. Our data will provide better understanding of DNA methylation alteration and their potential role in soybean insect resistance.

Functional Characterization of the Lysine-Specific Histone Demethylases Family in Soybean

Expression and Functional Analysis of a Novel Group of Legume-specific WRKY and Exo70 Protein Variants from Soybean

[Phylogenetic Analysis of DNA Demethylase Genes in Angiosperm].

Genome-Wide Identification of the Soybean AlkB Homologue Gene Family and Functional Characterization of GmALKBH10Bs as RNA m6A Demethylases and Expression Patterns under Abiotic Stress

Characterization and Subcellular Localization of Histone Deacetylases and Their Roles in Response to Abiotic Stresses in Soybean

Genome-Wide Identification and Characterization of the Abiotic-Stress-Responsive LACS Gene Family in Soybean (Glycine max)

Genome-Wide Identification, Evolution, and Expression of GDSL-Type Esterase/Lipase Gene Family in Soybean

Genome-Wide DNA Methylation Analysis of Soybean Curled-Cotyledons Mutant and Functional Evaluation of a Homeodomain-Leucine Zipper (HD-Zip) I Gene GmHDZ20

Structural Basis of a Histone H3 Lysine 4 Demethylase Required for Stem Elongation in Rice.

Genome-Wide Analysis of Soybean Lateral Organ Boundaries Domain Gene Family Reveals the Role in Phytophthora Root and Stem Rot

Genome-Wide Identification and Characterization of Soybean GmLOR Gene Family and Expression Analysis in Response to Abiotic Stresses

Genome-wide identification and expression analysis of DNA demethylase family in cotton

Identification and Molecular Characterization of FKF1 and GI Homologous Genes in Soybean

Identification and Expression Analysis of Soybean (Glycine max L.) Dynamin Genes Reveal Their Involvements in Plant Development and Stress Response

Genome-wide analysis of the GLP gene family and overexpression of GLP1-5–1 to promote lignin accumulation during early somatic embryo development in Dimocarpus longan

Genome-wide identification of genes encoding SWI/SNF components in soybean and the functional characterization of GmLFR1 in drought-stressed plants

Genome-wide identification of actin-depolymerizing factor gene family and their expression patterns under various abiotic stresses in soybean (Glycine max)

Genome-wide identification and characterization of DEAD-box helicase family associated with early somatic embryogenesis in Dimocarpus longan Lour

GmLecRlk, a Lectin Receptor-like Protein Kinase, Contributes to Salt Stress Tolerance by Regulating Salt-Responsive Genes in Soybean

Genome-wide investigation of CCCH zinc finger family in longan (Dimocarpus longan Lour): characteristic identification and expression profiles in longan somatic embryo

Methylome and transcriptome analyses of soybean response to bean pyralid larvae