Weighted Gene Coexpression Network Analysis of Candidate Pathways and Genes in Soft Rot Resistance of Amorphophallus

Amorphophallus species are one of the main economic crops in the mountainous areas of southwest China. However, soft rot disease (Pectobacterium carotovorum ssp. carotovorum) is devastating for this crop. This study explored the Amorphophallus resistance mechanism against soft rot disease by analyzing transcriptome data using a weighted gene coexpression network analysis. The RNA sequencing of plants infected for 0, 12, 24, and 48 hours produced a total of 52.25 Gb of clean reads. A total of 29,096 genes were divided into 34 modules. Six modules of interest with the highest correlation with the target traits were selected to elucidate the resistance genes and pathways. The selected modules were enriched in the α-linolenic acid metabolism, phenylpropane biosynthesis, plant hormone signal transduction, and plant pathogen interaction pathways. Ultimately, AmBGLU, AmCAML, AmCDPK, AmLOX, and AmRBOHD were identified as genes of interest in the four significantly related metabolic pathways for real-time fluorescence quantitative polymerase chain reaction verification. The determination of salicylic acid (SA) and jasmonic acid (JA) in Amorphophallus muelleri and Amorphophallus konjac that suffered from soft rot disease showed that SA and JA were involved in the A. muelleri and A. konjac defense response against soft rot disease. Methyl jasmonate treatment delayed the onset of A. konjac soft rot disease. This study provides a reference for the interaction between Amorphophallus species and soft rot disease and the breeding of broad-spectrum and specific Amorphophallus cultivars that are resistant to soft rot disease.