Transcriptomic and proteomic response of Manihot esculenta to Tetranychus urticae infestation at different densities

Tetranychus urticae (Acari: Tetranychidae) is an extremely serious cassava (Manihot esculenta) pest. Building a genomic resource to investigate the molecular mechanisms of cassava responses to T. urticae is vital for characterizing cassava resistance to mites. Based on the tolerance of cassava varieties to mite infestation (focusing on mite development rate, fecundity and physiology), cassava variety SC8 was selected to analyze transcriptomic and proteomic changes after 5 days of T. urticae feeding. Transcriptomic analysis revealed 698 and 2140 genes with significant expression changes under low and high mite infestation, respectively. More defense-related genes were found in the enrichment pathways at high mite density than at low density. In addition, iTRAQ-labeled proteomic analysis revealed 191 proteins with significant expression changes under low mite infestation. Differentially expressed genes and proteins were mainly found in the following defense-related pathways: flavonoid biosynthesis, phenylpropanoid biosynthesis, and glutathione metabolism under low-density mite feeding and plant hormone signal transduction and plant-pathogen interaction pathways under high-density mite feeding. The plant hormone signal transduction network, involving ethylene, jasmonic acid, and salicylic acid transduction pathways, was explored in relation to the M. esculenta response to T. urticae. Correlation analysis of the transcriptome and proteome generated a Pearson correlation coefficients of R = 0.2953 (P < 0.01), which might have been due to post-transcriptional or post-translational regulation resulting in many genes being inconsistently expressed at both the transcript and protein levels. In summary, the M. esculenta transcriptome and proteome changed in response to T. urticae, providing insight into the general activation of plant defense pathways in response to mite infestation.