Phosphine is a small redox-active gas that is used to protect global grain reserves, which are threatened by the emergence of phosphine resistance in pest insects. We find that polymorphisms responsible for genetic resistance cluster around the redox-active catalytic disulfide or the dimerization interface of dihydrolipoamide dehydrogenase (DLD) in insects (Rhyzopertha dominica and Tribolium castaneum) and nematodes (Caenorhabditis elegans). DLD is a core metabolic enzyme representing a new class of resistance factor for a redox-active metabolic toxin. It participates in four key steps of core metabolism, and metabolite profiles indicate that phosphine exposure in mutant and wild-type animals affects these steps differently. Mutation of DLD in C. elegans increases arsenite sensitivity. This specific vulnerability may be exploited to control phosphine-resistant insects and safeguard food security.

The recent emergence of heritable high level resistance to phosphine in stored grain pests is a serious concern among major grain growing countries around the world. Here we describe the genetics of phosphine resistance in the rust red flour beetle Tribolium castaneum (Herbst), a pest of stored grain as well as a genetic model organism. We investigated three field collected strains of T. castaneum viz., susceptible (QTC4), weakly resistant (QTC1012) and strongly resistant (QTC931) to phosphine. The dose-mortality responses of their test- and inter-cross progeny revealed that most resistance was conferred by a single major resistance gene in the weakly (3.2×) resistant strain. This gene was also found in the strongly resistant (431×) strain, together with a second major resistance gene and additional minor factors. The second major gene by itself confers only 12–20× resistance, suggesting that a strong synergistic epistatic interaction between the genes is responsible for the high level of resistance (431×) observed in the strongly resistant strain. Phosphine resistance is not sex linked and is inherited as an incompletely recessive, autosomal trait. The analysis of the phenotypic fitness response of a population derived from a single pair inter-strain cross between the susceptible and strongly resistant strains indicated the changes in the level of response in the strong resistance phenotype; however this effect was not consistent and apparently masked by the genetic background of the weakly resistant strain. The results from this work will inform phosphine resistance management strategies and provide a basis for the identification of the resistance genes.

Phosphine (PH 3 ) fumigation is the primary method worldwide for controlling insect pests of stored commodities. Over-reliance on phosphine, however, has led to the emergence of strong resistance. Detailed genetic studies previously identified two loci, rph1 and rph2, that interact synergistically to create a strong resistance phenotype. We compared the genetics of phosphine resistance in strains of Rhyzopertha dominica and Tribolium castaneum from India and Australia, countries having similar pest species but widely differing in pest management practices. Sequencing analysis of the rph2 locus, dihydrolipoamide dehydrogenase (dld), identified two structurally equivalent variants, Proline494Serine (P49S) in one R. dominica strain and P45S in three strains of T. castaneum from India. These variants of the DLD protein likely affect FAD cofactor interaction with the enzyme. A survey of insects from storage facilities across southern India revealed that the P45/49S variant is distributed throughout the region at very high frequencies, in up to 94% of R. dominica and 97% of T. castaneum in the state of Tamil Nadu. The abundance of the P45/49S variant in insect populations contrasted sharply with the evolutionary record in which the variant was absent from eukaryotic DLD sequences. This suggests that the variant is unlikely to provide a strong selective advantage in the absence of phosphine fumigation.

BackgroundNext-generation sequencing technology is an important tool for the rapid, genome-wide identification of genetic variations. However, it is difficult to resolve the ‘signal’ of variations of interest and the ‘noise’ of stochastic sequencing and bioinformatic errors in the large datasets that are generated. We report a simple approach to identify regional linkage to a trait that requires only two pools of DNA to be sequenced from progeny of a defined genetic cross (i.e. bulk segregant analysis) at low coverage (<10×) and without parentage assignment of individual SNPs. The analysis relies on regional averaging of pooled SNP frequencies to rapidly scan polymorphisms across the genome for differential regional homozygosity, which is then displayed graphically.ResultsProgeny from defined genetic crosses of Tribolium castaneum (F4 and F19) segregating for the phosphine resistance trait were exposed to phosphine to select for the resistance trait while the remainders were left unexposed. Next generation sequencing was then carried out on the genomic DNA from each pool of selected and unselected insects from each generation. The reads were mapped against the annotated T. castaneum genome from NCBI (v3.0) and analysed for SNP variations. Since it is difficult to accurately call individual SNP frequencies when the depth of sequence coverage is low, variant frequencies were averaged across larger regions. Results from regional SNP frequency averaging identified two loci, tc_rph1 on chromosome 8 and tc_rph2 on chromosome 9, which together are responsible for high level resistance. Identification of the two loci was possible with only 5-7× average coverage of the genome per dataset. These loci were subsequently confirmed by direct SNP marker analysis and fine-scale mapping. Individually, homozygosity of tc_rph1 or tc_rph2 results in only weak resistance to phosphine (estimated at up to 1.5-2.5× and 3-5× respectively), whereas in combination they interact synergistically to provide a high-level resistance >200×. The tc_rph2 resistance allele resulted in a significant fitness cost relative to the wild type allele in unselected beetles over eighteen generations.ConclusionWe have validated the technique of linkage mapping by low-coverage sequencing of progeny from a simple genetic cross. The approach relied on regional averaging of SNP frequencies and was used to successfully identify candidate gene loci for phosphine resistance in T. castaneum. This is a relatively simple and rapid approach to identifying genomic regions associated with traits in defined genetic crosses that does not require any specialised statistical analysis.