Host-pathogen interactions are central components of ecological networks where the MAPK signaling pathways act as central hubs of these complex interactions. We have previously shown that an insect hormone modulated MAPK signaling cascade participates as a general switch to trans-regulate differential expression of diverse midgut genes in the diamondback moth, Plutella xylostella (L.) to cope with the insecticidal action of Cry1Ac toxin, produced by the entomopathogenic bacterium Bacillus thuringiensis (Bt). The relationship between topology and functions of this four-tiered phosphorylation signaling cascade, however, is an uncharted territory. Here, we carried out a genome-wide characterization of all the MAPK orthologs in P. xylostella to define their phylogenetic relationships and to confirm their evolutionary conserved modules. Results from quantitative phosphoproteomic analyses, combined with functional validations studies using specific inhibitors and dsRNAs lead us to establish a MAPK “road map”, where p38 and ERK MAPK signaling pathways, in large part, mount a resistance response against Bt toxins through regulating the differential expression of multiple Cry toxin receptors and their non-receptor paralogs in P. xylostella midgut. These data not only advance our understanding of host-pathogen interactions in agricultural pests, but also inform the future development of biopesticides that could suppress Cry resistance phenotypes.

ABC transporter superfamily is one of the largest and ubiquitous groups of proteins. Because of their role in detoxification, insect ABC transporters have gained more attention in recent years. In this study, we annotated ABC transporters from a newly sequenced sweetpotato whitefly genome. Bemisia tabaci Q biotype is an emerging global invasive species that has caused extensive damages to field crops as well as ornamental plants. A total of 55 ABC transporters containing all eight described subfamilies (A to H) were identified in the B. tabaci Q genome, including 8 ABCAs, 3 ABCBs, 6 ABCCs, 2 ABCDs, 1 ABCE, 3 ABCFs, 23 ABCGs and 9 ABCHs. In comparison to other species, subfamilies G and H in both phloem- and blood-sucking arthropods are expanded. The temporal expression profiles of these 55 ABC transporters throughout B. tabaci developmental stages and their responses to imidacloprid, a neonicotinoid insecticide, were investigated using RNA-seq analysis. Furthermore, the mRNA expression of 24 ABC transporters (44% of the total) representing all eight subfamilies was confirmed by the quantitative real-time PCR (RT-qPCR). Furthermore, mRNA expression levels estimated by RT-qPCR and RNA-seq analyses were significantly correlated (r = 0.684, p < 0.01). It is the first genome-wide analysis of the entire repertoire of ABC transporters in B. tabaci. The identification of these ABC transporters, their temporal expression profiles during B. tabaci development, and their response to a neonicotinoid insecticide lay the foundation for functional genomic understanding of their contribution to the invasiveness of B. tabaci.