A Method for Identification of Biotype-Specific Salivary Effector Candidates of Aphid

Polyphagous aphids often consist of host-specialized biotypes that perform poorly in non-native hosts. The underlying mechanisms remain unknown. Host-specialized biotypes may express biotype-specific salivary effectors or elicitors that determine aphid hosts. Here, we tried three strategies to identify possible effectors in Malvaceae- (MA) and Cucurbitaceae-specialized (CU) biotypes of the cotton-melon aphid Aphis gossypii Glover. The whole-aphid RNA-seq identified 765 differentially expressed genes (DEGs), and 139 of them were possible effectors; aphid-head RNA-seq identified 523 DEGs were identified, and 98 of them were possible effectors. The homologous genes of published aphid effectors were not differentially expressed between CU and MA. Next, quantitative proteomic analyses of saliva identified 177 possible proteins, and 44 of them were different proteins. However, none of the genes of the 44 proteins were differentially expressed, reflecting the discrepancy between transcriptome and proteome data. Finally, we searched for DEGs of the 177 salivary proteins in the aphid-head transcriptomes, and the salivary proteins with expression differences were regarded as effector candidates. Through this strategy, 11 effector candidates were identified, and their expression differences were all confirmed by RT-qPCR. The combinatorial analysis has great potential to identify biotype-specific effector candidates in aphids and other sap-sucking insects.