Spirotetramat resistance adaption analysis of Aphis gossypii Glover by transcriptomic survey

A resistant strain of the cotton aphid (SR) developed 441.26-fold and 11.97-fold resistance to spirotetramat for adult aphids and nymphs, respectively, compared with the susceptible (SS) strain. Solexa sequencing technology was employed to identify differentially expressed genes (DEGs) in the spirotetramat-resistant cotton aphid. Respective totals of 22,430,522 and 21,317,732 clean reads were obtained from SR and SS cDNA libraries and assembled into 35,222 non-redundant (Nr) consensus sequences. A total of 14,913, 9,220, 7,922, 4,314 and 4,686 sequences were annotated using Nr, Swiss-Prot, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Clusters of Orthologous Groups (COG), respectively. Compared with the SS strain, the SR strain had 1287 significantly changed unigenes, of which 130 genes were up-regulated and 1157 genes were down-regulated (P ≤ 0.001). Among these genes, 440 unigenes were annotated, consisting of 114 up-regulated and 326 down-regulated genes. The expression levels of heat shock protein 70 (Hsp70) and UDP-glucuronosyltransferase were significantly up-regulated in the SR strain compared to the SS strain. The genes encoding cuticle proteins, salivary glue protein, fibroin heavy chain, energy ATP synthase, and cytochrome c oxidase were dramatically decreased. Among the DEGs, cytochrome P450 6A2 (c20965.graph_c0) was the only P450 gene up-regulated in the SR strain. The expression levels of 10 DEGs were confirmed by real-time qPCR, and the trends in gene expression observed by qPCR matched those of the Solexa expression profiles. The acetyl-CoA carboxylase (ACC) genes in the SR and SS libraries both contain four single nucleotide polymorphisms (SNPs), with three common SNPs: 1227 (C/T), 1811 (A/T: F/Y) and 3759 (C/T); however, 7540 (A/T) and 108 (G/A) occurred solely in the SS and SR strains, respectively.