Author:
Zeng Wei-Ying,Tan Yu-Rong,Long Sheng-Feng,Sun Zu-Dong,Lai Zhen-Guang,Yang Shou-Zhen,Chen Huai-Zhu,Qing Xia-Yan
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.
Publisher
Springer Science and Business Media LLC
Reference88 articles.
1. Choi KH, Hong YK, Chang YJ, Moon JS, Kim CS, Choi DC, et al. Development under constant temperatures and seasonal prevalence in soybean field of the bean pyralid, Omiodes indicatus (Lepidoptera: Crambidae). Korean J Appl Entomol. 2008;47(4):353–8. https://doi.org/10.5656/KSAE.2008.47.4.353.
2. Li GJ. Detection of QTL for photooxidation- related traits, dynamic chlorophyll content and resistance to bean pyralid (Lamprosema indicata Fabricius) in soybean (Glycine max Merrill). Nanjing: Nanjing Agricultural University. 2010; p. 38.
3. Sun ZD, Yang SZ, Chen HZ, Li CY, Long LP. Identification of soybean resistance to bean pyralid (Lamprosema indicata Fabricicus) and obiposition preference of bean pyralid on soybean varieties. Chin J Oil Crop Sci. 2005;27(4):69–71.
4. Xing GN, Zhou B, Wang YF, Zhao TJ, Yu DY, Chen SY, et al. Genetic components and major QTL confer resistance to bean pyralid (Lamprosema indicata Fabricius) under multiple environments in four RIL populations of soybean. Theor Appl Genet. 2012;125(5):859–75. https://doi.org/10.1007/s00122-012-1878-7.
5. Li GJ, Cheng LG, Zhang GZ, He XH, Zhi HJ, Zhang YM. Mixed major-gene plus polyegens inheritance analysis for resistance in soybean to bean pyralid (Lamprosema indicata Fabricius). Soybean Sci. 2008;27(1):33–6, 41.