Variant Linkage Analysis Using de Novo Transcriptome Sequencing Identifies a Conserved Phosphine Resistance Gene in Insects

Abstract

Next-generation sequencing methods enable identification of the genetic basis of traits in species that have no prior genomic information available. The combination of next-generation sequencing, variant analysis, and linkage is a powerful way of identifying candidate genes for a trait of interest. Here, we used a comparative transcriptomics [RNA sequencing (RNAseq)] and genetic linkage analysis approach to identify the rph1 gene. rph1 variants are responsible for resistance to the fumigant phosphine (PH₃) that is used to control insect pests of stored grain. In each of the four major species of pest insect of grain we have investigated, there are two major resistance genes, rph1 and rph2, which interact synergistically to produce strongly phosphine-resistant insects. Using RNAseq and genetic linkage analyses, we identified candidate resistance (rph1) genes in phosphine-resistant strains of three species: Rhyzopertha dominica (129 candidates), Sitophilus oryzae (206 candidates), and Cryptolestes ferrugineus (645 candidates). We then compared these candidate genes to 17 candidate resistance genes previously mapped in Tribolium castaneum and found only one orthologous gene, a cytochrome b5 fatty acid desaturase (Cyt-b5-r), to be associated with the rph1 locus in all four species. This gene had either missense amino acid substitutions and/or insertion/deletions/frameshift variants in each of 18 phosphine-resistant strains that were not observed in the susceptible strains of the four species. We propose a model of phosphine action and resistance in which phosphine induces lipid peroxidation through reactive oxygen species generated by dihydrolipoamide dehydrogenase, whereas disruption of Cyt-b5-r in resistant insects decreases the polyunsaturated fatty acid content of membranes, thereby limiting the potential for lipid peroxidation.

Keywords: insecticide resistance; linkage; pesticide; transcriptome.

Figure 1

Venn diagram showing the number of shared homologs within the candidate gene sets for the four species. Homology was determined by best similarity by BLAST (Basic Local Alignment Search Tool) with D. melanogaster and T. castaneum.

Figure 2

Aligned predicted protein models for phosphine-resistant strains of the four species studied. The susceptible reference strains used for the analysis are denoted by an asterisk. SNPs that changed amino acids compared to the respective susceptible reference strains are denoted by triangles, with those that are likely (P < 0.05) to cause gene disruption by SIFT analysis highlighted in red. Dotted lines indicate a translation frameshift.