Rapid genome wide mapping of phosphine resistance loci by a simple regional averaging analysis in the red flour beetle, Tribolium castaneum

Background Next-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. Results Progeny 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. Conclusion We 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.

Do Phosphine Resistance Genes Influence Movement and Dispersal Under Starvation?

Phosphine resistance alleles might be expected to negatively affect energy demanding activities such as walking and flying, because of the inverse relationship between phosphine resistance and respiration. We used an activity monitoring system to quantify walking of Rhyzopertha dominica (F.) and a flight chamber to estimate their propensity for flight initiation. No significant difference in the duration of walking was observed between the strongly resistant, weakly resistant, and susceptible strains of R. dominica we tested, and females walked significantly more than males regardless of genotype. The walking activity monitor revealed no pattern of movement across the day and no particular time of peak activity despite reports of peak activity of R. dominica and Tribolium castaneum (Herbst) under field conditions during dawn and dusk. Flight initiation was significantly higher for all strains at 28 degrees C and 55% relative humidity than at 25, 30, 32, and 35 degrees C in the first 24 h of placing beetles in the flight chamber. Food deprivation and genotype had no significant effect on flight initiation. Our results suggest that known resistance alleles in R. dominica do not affect insect mobility and should therefore not inhibit the dispersal of resistant insects in the field.

Developing effective fumigation protocols to manage strongly phosphine-resistant Cryptolestes ferrugineus (Stephens) (Coleoptera: Laemophloeidae)

BACKGROUND The emergence of high levels of resistance in Cryptolestes ferrugineus (Stephens) in recent years threatens the sustainability of phosphine, a key fumigant used worldwide to disinfest stored grain. We aimed at developing robust fumigation protocols that could be used in a range of practical situations to control this resistant pest. RESULTS Values of the lethal time to kill 99.9% (LT99.9, in days) of mixed-age populations, containing all life stages, of a susceptible and a strongly resistant C. ferrugineus population were established at three phosphine concentrations (1.0, 1.5 and 2.0 mg L−1) and three temperatures (25, 30 and 35 °C). Multiple linear regression analysis revealed that phosphine concentration and temperature both contributed significantly to the LT99.9 of a population (P < 0.003, R2 = 0.92), with concentration being the dominant variable, accounting for 75.9% of the variation. Across all concentrations, LT99.9 of the strongly resistant C. ferrugineus population was longest at the lowest temperature and shortest at the highest temperature. For example, 1.0 mg L−1 of phosphine is required for 20, 15 and 15 days, 1.5 mg L−1 for 12, 11 and 9 days and 2.0 mg L−1 for 10, 7 and 6 days at 25, 30 and 35 °C, respectively, to achieve 99.9% mortality of the strongly resistant C. ferrugineus population. We also observed that phosphine concentration is inversely proportional to fumigation period in regard to the population extinction of this pest. CONCLUSION The fumigation protocols developed in this study will be used in recommending changes to the currently registered rates of phosphine in Australia towards management of strongly resistant C. ferrugineus populations, and can be repeated in any country where this type of resistance appears.

Diagnostic Molecular Markers for Phosphine Resistance in U.S. Populations of Tribolium castaneum and Rhyzopertha dominica

Stored product beetles that are resistant to the fumigant pesticide phosphine (hydrogen phosphide) gas have been reported for more than 40 years in many places worldwide. Traditionally, determination of phosphine resistance in stored product beetles is based on a discriminating dose bioassay that can take up to two weeks to evaluate. We developed a diagnostic cleaved amplified polymorphic sequence method, CAPS, to detect individuals with alleles for strong resistance to phosphine in populations of the red flour beetle, Tribolium castaneum, and the lesser grain borer, Rhyzopertha dominica, according to a single nucleotide mutation in the dihydrolipoamide dehydrogenase (DLD) gene. We initially isolated and sequenced the DLD genes from susceptible and strongly resistant populations of both species. The corresponding amino acid sequences were then deduced. A single amino acid mutation in DLD in populations of T.castaneum and R.dominica with strong resistance was identified as P45S in T.castaneum and P49S in R.dominica, both collected from northern Oklahoma, USA. PCR products containing these mutations were digested by the restriction enzymes MboI and BstNI, which revealed presence or absence, respectively of the resistant (R) allele and allowed inference of genotypes with that allele. Seven populations of T.castaneum from Kansas were subjected to discriminating dose bioassays for the weak and strong resistance phenotypes. Application of CAPS to these seven populations confirmed the R allele was in high frequency in the strongly resistant populations, and was absent or at a lower frequency in populations with weak resistance, which suggests that these populations with a low frequency of the R allele have the potential for selection of the strong resistance phenotype. CAPS markers for strong phosphine resistance will help to detect and confirm resistant beetles and can facilitate resistance management actions against a given pest population.

Prevalence and potential fitness cost of weak phosphine resistance in Tribolium castaneum (Herbst) in eastern Australia

The prevalence of resistance to phosphine in the rust-red flour beetle, Tribolium castaneum, from eastern Australia was investigated, as well as the potential fitness cost of this type of resistance. Discriminating dose tests on 115 population samples collected from farms from 2006 to 2010 showed that populations containing insects with the weakly resistant phenotype are common in eastern Australia (65.2 of samples), although the frequency of resistant phenotypes within samples was typically low (median of 2.3). The population cage approach was used to investigate the possibility that carrying the alleles for weak resistance incurs a fitness cost. Hybridized populations were initiated using a resistant strain and either of two different susceptible strains. There was no evidence of a fitness cost based on the frequency of susceptible phenotypes in hybridized populations that were reared for seven generations without exposure to phosphine. This suggests that resistant alleles will tend to persist in field populations that have undergone selection even if selection pressure is removed. The prevalence of resistance is a warning that this species has been subject to considerable selection pressure and that effective resistance management practices are needed to address this problem. The resistance prevalence data also provide a basis against which to measure management success.

Inheritance and Characterization of Strong Resistance to Phosphine in Sitophilus oryzae(L.)

Sitophilus oryzae (Linnaeus) is a major pest of stored grain across Southeast Asia and is of increasing concern in other regions due to the advent of strong resistance to phosphine, the fumigant used to protect stored grain from pest insects. We investigated the inheritance of genes controlling resistance to phosphine in a strongly resistant S. oryzae strain (NNSO7525) collected in Australia and find that the trait is autosomally inherited and incompletely recessive with a degree of dominance of -0.66. The strongly resistant strain has an LC50 52 times greater than a susceptible reference strain (LS2) and 9 times greater than a weakly resistant strain (QSO335). Analysis of F2 and backcross progeny indicates that two or more genes are responsible for strong resistance, and that one of these genes, designated Sorph1, not only contributes to strong resistance, but is also responsible for the weak resistance phenotype of strain QSO335. These results demonstrate that the genetic mechanism of phosphine resistance in Soryzae is similar to that of other stored product insect pests. A unique observation is that a subset of the progeny of an F1 backcross generation are more strongly resistant to phosphine than the parental strongly resistant strain, which may be caused by multiple alleles of one of the resistance genes.

Genetic characterization of field-evolved resistance to phosphine in the rusty grain beetle, Cryptolestes ferrugineus (Laemophloeidae: Coleoptera)

Inheritance of resistance to phosphine fumigant was investigated in three field-collected strains of rusty grain beetle, Cryptolestes ferrugineus, Susceptible (S-strain), Weakly Resistant (Weak-R) and Strongly Resistant (Strong-R). The strains were purified for susceptibility, weak resistance and strong resistance to phosphine, respectively, to ensure homozygosity of resistance genotype. Crosses were established between S-strain × Weak-R, S-strain × Strong-R and Weak-R × Strong-R, and the dose mortality responses to phosphine of these strains and their F1, F2 and F1-backcross progeny were obtained. The fumigations were undertaken at 25 °C and 55% RH for 72 h. Weak-R and Strong-R showed resistance factors of 6.3 × and 505 × compared with S-strain at the LC50. Both weak and strong resistances were expressed as incompletely recessive with degrees of dominance of − 0.48 and − 0.43 at the LC50, respectively. Responses of F2 and F1-backcross progeny indicated the existence of one major gene in Weak-R, and at least two major genes in Strong-R, one of which was allelic with the major factor in Weak-R. Phenotypic variance analyses also estimated that the number of independently segregating genes conferring weak resistance was 1 (nE = 0.89) whereas there were two genes controlling strong resistance (nE = 1.2). The second gene, unique to Strong-R, interacted synergistically with the first gene to confer a very high level of resistance (~ 80 ×). Neither of the two major resistance genes was sex linked. Despite the similarity of the genetics of resistance to that previously observed in other pest species, a significant proportion (~ 15 to 30%) of F1 individuals survived at phosphine concentrations higher than predicted. Thus it is likely that additional dominant heritable factors, present in some individuals in the population, also influenced the resistance phenotype. Our results will help in understanding the process of selection for phosphine resistance in the field which will inform resistance management strategies. In addition, this information will provide a basis for the identification of the resistance genes.

Genetic characterization of field-evolved resistance to phosphine in the rusty grain beetle, Cryptolestes ferrugineus (Laemophloeidae: Coleoptera)

Inheritance of resistance to phosphine fumigant was investigated in three field-collected strains of rusty grain beetle, Cryptolestes ferrugineus, Susceptible (S-strain), Weakly Resistant (Weak-R) and Strongly Resistant (Strong-R). The strains were purified for susceptibility, weak resistance and strong resistance to phosphine, respectively, to ensure homozygosity of resistance genotype. Crosses were established between S-strain × Weak-R, S-strain × Strong-R and Weak-R × Strong-R, and the dose mortality responses to phosphine of these strains and their F1, F2 and F1-backcross progeny were obtained. The fumigations were undertaken at 25 °C and 55% RH for 72 h. Weak-R and Strong-R showed resistance factors of 6.3 × and 505 × compared with S-strain at the LC50. Both weak and strong resistances were expressed as incompletely recessive with degrees of dominance of − 0.48 and − 0.43 at the LC50, respectively. Responses of F2 and F1-backcross progeny indicated the existence of one major gene in Weak-R, and at least two major genes in Strong-R, one of which was allelic with the major factor in Weak-R. Phenotypic variance analyses also estimated that the number of independently segregating genes conferring weak resistance was 1 (nE = 0.89) whereas there were two genes controlling strong resistance (nE = 1.2). The second gene, unique to Strong-R, interacted synergistically with the first gene to confer a very high level of resistance (~ 80 ×). Neither of the two major resistance genes was sex linked. Despite the similarity of the genetics of resistance to that previously observed in other pest species, a significant proportion (~ 15 to 30%) of F1 individuals survived at phosphine concentrations higher than predicted. Thus it is likely that additional dominant heritable factors, present in some individuals in the population, also influenced the resistance phenotype. Our results will help in understanding the process of selection for phosphine resistance in the field which will inform resistance management strategies. In addition, this information will provide a basis for the identification of the resistance genes.

Developing effective fumigation protocols to manage strongly phosphine-resistant Cryptolestes ferrugineus (Stephens) (Coleoptera: Laemophloeidae)

BACKGROUND The emergence of high levels of resistance in Cryptolestes ferrugineus (Stephens) in recent years threatens the sustainability of phosphine, a key fumigant used worldwide to disinfest stored grain. We aimed at developing robust fumigation protocols that could be used in a range of practical situations to control this resistant pest. RESULTS Values of the lethal time to kill 99.9% (LT99.9, in days) of mixed-age populations, containing all life stages, of a susceptible and a strongly resistant C. ferrugineus population were established at three phosphine concentrations (1.0, 1.5 and 2.0 mg L−1) and three temperatures (25, 30 and 35 °C). Multiple linear regression analysis revealed that phosphine concentration and temperature both contributed significantly to the LT99.9 of a population (P < 0.003, R2 = 0.92), with concentration being the dominant variable, accounting for 75.9% of the variation. Across all concentrations, LT99.9 of the strongly resistant C. ferrugineus population was longest at the lowest temperature and shortest at the highest temperature. For example, 1.0 mg L−1 of phosphine is required for 20, 15 and 15 days, 1.5 mg L−1 for 12, 11 and 9 days and 2.0 mg L−1 for 10, 7 and 6 days at 25, 30 and 35 °C, respectively, to achieve 99.9% mortality of the strongly resistant C. ferrugineus population. We also observed that phosphine concentration is inversely proportional to fumigation period in regard to the population extinction of this pest. CONCLUSION The fumigation protocols developed in this study will be used in recommending changes to the currently registered rates of phosphine in Australia towards management of strongly resistant C. ferrugineus populations, and can be repeated in any country where this type of resistance appears. © 2014 Commonwealth of Australia. Pest Management Science © 2014 Society of Chemical Industry

Susceptibility to sulfuryl fluoride and lack of cross-resistance to phosphine in developmental stages of the red flour beetle, Tribolium castaneum (Coleoptera: Tenebrionidae)

BACKGROUND Our aim was to ascertain the potential of sulfuryl fluoride (SF) as an alternative fumigant to manage phosphine-resistant pests. We tested the susceptibility of all life stages of red flour beetle, Tribolium castaneum (Herbst), to SF and assessed the presence of cross-resistance to this fumigant in phosphine-resistant strains of this species. RESULTS Analysis of dose–response data indicated that the egg was the stage most tolerant to SF under a 48 h exposure period. At LC50, eggs were 29 times more tolerant than other immature stages and adults, and required a relatively high concentration of 48.2 mg L−1 for complete mortality. No significant differences in tolerance to SF were observed among the three larval instars, pupae and adults, and all of these stages were controlled at a low concentration of 1.32 mg L−1. Phosphine-resistant strains did not show cross-resistance to SF. CONCLUSION Our research concluded that the current maximum registered rate of SF, 1500 gh m−3, is adequate to control all the post-embryonic life stages of T. castaneum over a 48 h fumigation period, but it will fail to achieve complete mortality of eggs, indicating the risk of some survival of eggs under this short exposure period. As there is no cross-resistance to SF in phosphine-resistant insects, it will play a key role in managing phosphine resistance in stored-grain insect pests. © 2014 Commonwealth of Australia. Pest Management Science © 2014 Society of Chemical Industry

Phosphine resistance in India is characterised by a dihydrolipoamide dehydrogenase variant that is otherwise unobserved in eukaryotes

Phosphine (PH3) 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, Proline49>Serine (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.

Genetic conservation of Phosphine Resistance in the Rice Weevil Sitophilus oryzae (L.)

High levels of resistance to phosphine in the rice weevil Sitophilus oryzae have been detected in Asian countries including China and Vietnam, however there is limited knowledge of the genetic mechanism of resistance in these strains. We find that the genetic basis of strong phosphine resistance is conserved between strains of S. oryzae from China, Vietnam and Australia. Each of four strongly resistant strains has an identical amino acid variant in the encoded dihydrolipoamide dehydrogenase (DLD) enzyme that was previously identified as a resistance factor in Rhyzopertha dominica and Tribolium castaneum. The unique amino acid substitution, Asparagine > Threonine (N505T) of all strongly resistant S. oryzae corresponds to the position of an Asparagine > Histidine variant (N506H) that was previously reported in strongly resistant R. dominica. Progeny (F16 and F18) from two independent crosses showed absolute linkage of N505T to the strong resistance phenotype, indicating that if N505T was not itself the resistance variant that it resided within 1 or 2 genes of the resistance factor. Non-complementation between the strains confirmed the shared genetic basis of strong resistance, which was supported by the very similar level of resistance between the strains, with LC50 values ranging from 0.20 to 0.36 mgL-1 for a 48 hour exposure at 25°C. Thus, the mechanism of high level resistance to phosphine is strongly conserved between R. dominica, T. castaneum and S. oryzae. A fitness cost associated with strongly resistant allele was observed in segregating populations in the absence of selection.

Use of sulfuryl fluoride in the management of strongly phosphine-resistant insect pest populations in bulk grain storages in Australia

Sulfuryl fluoride (SF), an effective structural fumigant, is registered recently as Profume™ for controlling insect pests of stored grains and processed commodities. Information on its effectiveness in disinfestation of bulk grain, however, is limited. The ongoing problem with the strong level of resistance to phosphine has been addressed recently through deployment of SF as a ‘resistance breaker’ in bulk storages in Australia. This paper discusses important results on the efficacy of SF against key phosphine- resistant insect pests, lesser grain borer, Rhyzopertha dominca, red flour beetle, Tribolium castaneum, rice weevil, Sitophilus oryzae and the rusty grain beetle, Cryptolestes ferrugineus. We have established CT (g-hm3) profiles for SF against these insect pests at two temperature regimes 25 and 30°C, that showed that both temperature and exposure period (t) has significant influence on the effectiveness of SF than the concentration. Over a seven days fumigation period, CTs of 800 and 400 g-hm3 achieved complete control of all the target pests, including the most strongly phosphine - resistant species, C. ferrugineus at 25 and 30°C, respectively. Results from four industry scale field trials involving currently registered rate of SF (1500 g-hm3) over 2–14 d exposure period, confirmed its effectiveness in achieving complete control of the target pest species. The assessment of postfumigation grain samples across all the test storages indicated that the reinfestation occurs after three months. Monitoring resistance to phosphine in C. ferrugineus over a six year period (2009–2015), showed a significant reduction in resistant populations after the introduction of SF into the fumigation strategy at problematic storage sites. Overall our research concludes that SF is a good candidate to be used as a ‘resistance breaker’ where phosphine resistance is prevalent.

Status of resistance to phosphine in insect pests of stored products in Vietnam

In response to numerous reports of failures to control insect pests of stored products with phosphine in Vietnam, a national survey for resistance to this key fumigant was undertaken in 2009–2011. Data from a more limited survey undertaken by the authors in 2002 in northern Vietnam are also presented. Samples collected in the 2002 survey (Sitophilus oryzae, n=8; Tribolium castaneum, n=8) were tested using a full dose- response assay, while for the 2009–11 survey, F1 generations were tested for resistance with two discriminating dosages of phosphine to detect frequency of weak and strong resistance phenotypes. Compared with a susceptible reference strain, in 2002, resistance to phosphine was indicated in six T. castaneum samples but only two of S. oryzae. Resistance factor, however, did not exceed 2.8-fold in T. castaneum and 1.7 in S. oryzae indicating relatively low frequency and weak expression of resistance. In 2009–11 survey, 176 samples were collected from a range of food and feed storages along the supply chain and from all major regions of Vietnam (125 sites). Rhyzopertha dominica and S. oryzae were the most common species found infesting stored commodities. Resistance was detected at high frequency in all the species. Weak and strong resistance phenotype frequencies were, respectively: Cryptolestes ferrugineus (37 and 58%, n=19), R. dominica (1.5 and 97%, n=65), S. oryzae (34 and 59%, n=82) and T. castaneum (70 and 30%, n=10). Strong resistance phenotype was detected in all the major regions and all parts of the supply chain but frequency was the highest in central storages and animal feed establishments. The increase in frequency and strength of resistance to phosphine in the eight years between the two surveys has been rapid and dramatic. The survey demonstrates the threat of resistance to grain protection in Vietnam and highlights the need for training of fumigators, and the development and adoption of phosphine resistance management tactics nationally.

Resistance to phosphine in Sitophilus oryzae in Australia: A national analysis of trends and frequencies over time and geographical spread

In Australia, along with many other parts of the world, fumigation with phosphine is a vital component in controlling stored grain insect pests. However, resistance is a factor that may limit the continued efficacy of this fumigant. While strong resistance to phosphine has been identified and characterised, very little information is available on the causes of its development and spread. Data obtained from a unique national resistance monitoring and management program were analysed, using Bayesian hurdle modelling, to determine which factors may be responsible. Fumigation in unsealed storages, combined with a high frequency of weak resistance, were found to be the main criteria that led to the development of strong resistance in Sitophilus oryzae. Independent development, rather than gene flow via migration, appears to be primarily responsible for the geographic incidence of strong resistance to phosphine in S. oryzae. This information can now be utilised to direct resources and education into those areas at high risk and to refine phosphine resistance management strategies.