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.

Bioinformatics Analysis on 5' Regulatory Region Sequence of DNA Demethylatase Gene IDM1 in Soybean

Innovation of a Regulatory Mechanism Modulating Semi-determinate Stem Growth Through Artificial Selection in Soybean

[Phylogenetic Analysis of DNA Demethylase Genes in Angiosperm].

Analysis of the Molecular Motif for Inducing Response to Jasmonic Acid and Ethylene in Pib Promoter Via Rice Transfor-Mation

A new demethylase gene OsDML4 involved in high temperature induced floury endosperm formation in rice (Oryza sativa L.)

Analysis of the molecular motif for inducing response to jasmonic acid and ethylene in Pib promoter via rice transformation]

Locus- and Site-Specific DNA Methylation of 19 Kda Zein Genes in Maize.

Genome-Wide Analysis of DNA Methylation 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

A DIRIGENT Gene GmDIR26 Regulates Pod Dehiscence in Soybean

Genome-wide analysis of soybean hypoxia inducible gene domain containing genes: a functional investigation of GmHIGD3

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

Characterization of 40 soybean (Glycine max) promoters, isolated from across 5 thematic gene groups

Global Analysis Reveals The Crucial Roles Of Dna Methylation During Rice Seed Development

Methylome and transcriptome analyses of soybean response to bean pyralid larvae

Genetic dissection of the soybean dwarf mutant dm with integrated genomic, transcriptomic and methylomic analyses

Genome-wide characterization of the soybean DOMAIN OF UNKNOWN FUNCTION 679 membrane protein gene family highlights their potential involvement in growth and stress response

Over-expression of GmMYB39 leads to an inhibition of the isoflavonoid biosynthesis in soybean (Glycine max. L)

Characteristics and phylogeny of DREB gene subfamily in soybeans [Glycine max (L.) Meril]

Natural Variation in the Promoter of GmSPL9d Affects Branch Number in Soybean

DNA methylation impacts soybean early development by modulating hormones and metabolic pathways