Determining changes in the distribution and abundance of a Rhyzopertha dominica phosphine resistance allele in farm grain storages using a DNA marker

Development, application and evaluation of three novel TaqMan qPCR assays for phosphine resistance monitoring in major stored product pests Tribolium castaneum and Rhyzopertha dominica

BACKROUND

Stored product protection from insect pests relies heavily on the use of phosphine. The most serious drawback of phosphine is the development of resistance in major stored product insects worldwide, including the red flour beetle, Tribolium castaneum (Herbst) and the lesser grain borer, Rhyzopertha dominica (F.). Two genetic loci are responsible for phosphine resistance: the rph1 (S349G mutation in the cyt-b5-r homolog) in T. castaneum and the rph2 (P45/49S mutation in the dihydrolipoamide dehydrogenase (dld) gene) in T. castaneum and R. dominica.

RESULTS

In this study, we have developed and applied high-throughput, practical and specific molecular diagnostics (TaqMan qPCR) for monitoring mutations S349G, P45S and P49S. In our pilot monitoring application, we have included phosphine-resistant and susceptible populations from different parts of the world (USA, Australia, Brazil) and European strains from Greece and Serbia. Our results for the resistant T. castaneum showed a P45S mutant allele frequency (MAF) of 100% and 75.0% in the populations from Serbia and Brazil, respectively. Regarding the susceptible T. castaneum, P45S was detected in Greece (MAF = 62.5%) and was absent in Australia (MAF = 0.0%). Additionally, the S349G mutation was found to be fixed in all resistant populations, while it was also detected in susceptible ones (frequencies: 65.0% and 100.0%). The only case where both mutations were fixed (100%) was a resistant population from Serbia. In R. dominica, the P49S mutation was found only in the two resistant R. dominica populations from Serbia and Greece (50.0% and 100%) and was absent from the susceptible one from Greece; thus, P49S seems to be a satisfactory indicator for monitoring phosphine resistance.

CONCLUSIONS

Our P49S detection assay in R. dominica seems to be a viable option in this direction, yet its utilization needs additional large-scale confirmatory work. The identification of additional resistance markers also should be prioritized. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Differential Activities of Antioxidant Enzymes, Superoxide Dismutase, Peroxidase, and Catalase vis-à-vis Phosphine Resistance in Field Populations of Lesser Grain Borer (Rhyzopertha dominica) from India

Susceptibility to phosphine was compared in 15 populations of lesser grain borer (Rhyzopertha dominica) collected from grain storage godowns across India. A high level of resistance to phosphine was noticed in R. dominica collected from northern India compared to those collected from northeastern regions of India. The median lethal concentration values varied from 0.024 mg/L to 1.991 mg/L, with 1.63 to 82.96-fold resistance compared to laboratory susceptible checks. Antioxidant enzymes have been reported to negate the reactive oxygen species generated upon encountering the fumigant phosphine. Distinct differences in the activity of antioxidant enzymes were noticed in the field populations exposed to phosphine. Peroxidase activity varied between 1.28 and 336.8 nmol H₂O₂ reduced/min/mg protein. The superoxide dismutase inhibition rate was between 81.29 and 99.66%, and catalase activity varied between 6.28 and 320.13 nmol H₂O₂ reduced/min/mg protein. The findings of our investigation show that the activities of peroxidase and superoxide dismutase are positively linked (p < 0.01) with an increase in resistance ratios, whereas catalase was found to have a negative association with resistance to phosphine. The reported results elucidate the differential activities of principal antioxidant enzymes in scavenging the oxyradicals (O₂^•-, H₂O_2,^•OH) associated with tolerance to phosphine in R. dominica.

Amplification refractory mutation system based real-time PCR (ARMS-qPCR) for rapid resistance characterization of Tribolium castaneum to phosphine

Resistance of Tribolium castaneum to phosphine is related to point mutations in DNA code corresponding to amino acid changes associated with a core metabolic enzyme dihydrolipoamide dehydrogenase (DLD), but the mutation patterns vary among different resistant populations. Thus, there is a great need to develop a cost-effective method to detect core mutations in T. castaneum, which would be the key factor to understand the molecular basis of phosphine resistance. Amplification refractory mutation system-based quantitative Real-Time PCR (ARMS-qPCR) is an ideal method that can rapidly detect point mutations. Here, the P45S and G131D mutations existed in the DLD of T. castaneum selected from strong Chinese resistance phenotypes, and the DLD P45S mutation, which represents a strong phosphine resistance allele, was confirmed as the most abundant mutation to determine strong resistance genotypes. Our study found that 85 out of 120 beetles carried the P45S resistance allele, including 51 homozygous and 34 heterozygous individuals. Moreover, there was a strong linear relationship (R² = 0.917) between the resistance ratio and the resistance allele frequency among the strongly resistant populations. Our data showed that the ARMS-qPCR method that we developed could rapidly determine strong resistance phenotypes of T. castaneum to phosphine by detecting the DLD P45S mutation. These results not only provide a detailed example for developing an ARMS-qPCR-based method to characterize pesticide resistance, but also support further elucidation of the molecular basis of phosphine resistance.

The Genome of Rhyzopertha dominica (Fab.) (Coleoptera: Bostrichidae): Adaptation for Success

The lesser grain borer, Rhyzopertha dominica (F.) (Coleoptera: Bostrichidae), is a major global pest of cereal grains. Infestations are difficult to control as larvae feed inside grain kernels, and many populations are resistant to both contact insecticides and fumigants. We sequenced the genome of R. dominica to identify genes responsible for important biological functions and develop more targeted and efficacious management strategies. The genome was assembled from long read sequencing and long-range scaffolding technologies. The genome assembly is 479.1 Mb, close to the predicted genome size of 480.4 Mb by flow cytometry. This assembly is among the most contiguous beetle assemblies published to date, with 139 scaffolds, an N₅₀ of 53.6 Mb, and L₅₀ of 4, indicating chromosome-scale scaffolds. Predicted genes from biologically relevant groups were manually annotated using transcriptome data from adults and different larval tissues to guide annotation. The expansion of carbohydrase and serine peptidase genes suggest that they combine to enable efficient digestion of cereal proteins. A reduction in the copy number of several detoxification gene families relative to other coleopterans may reflect the low selective pressure on these genes in an insect that spends most of its life feeding internally. Chemoreceptor genes contain elevated numbers of pseudogenes for odorant receptors that also may be related to the recent ontogenetic shift of R. dominica to a diet consisting primarily of stored grains. Analysis of repetitive sequences will further define the evolution of bostrichid beetles compared to other species. The data overall contribute significantly to coleopteran genetic research.

Evidence for the natural occurrence of Wolbachia in Aedes aegypti mosquitoes

Wolbachia, a Gram-negative bacterium, is widely known to be present in arthropods and nematodes. Of late, great impetus is given to employ this intracellular bacterium, as an alternative to conventional biocontrol agents for the control of mosquitoes because of its inherent ability to induce sperm-egg incompatibility, feminisation etc. By employing molecular tools, we have shown the presence of Wolbachia from Aedes aegypti mosquito population collected from Coimbatore, India by PCR amplifying the Ae. aegypti mosquito genome with Wolbachia specific 16S rRNA, wsp and ftsZ gene primers. The phylogenetic analysis of these gene sequences incorporating MLST and GenBank reference sequences has confirmed the occurrence of Wolbachia supergroup B in Ae. aegypti. In addition, qRT-PCR results have shown the dynamics of Wolbachia across the developmental stages of mosquito. The absence of Wolbachia in tetracycline-treated Ae. aegypti mosquitoes evidenced by transmission electron microscopy reinforced our finding conclusively. After confirming their persistence through generations, we have designated Wolbachia from Ae. aegypti as wAegB. In our considered view, wAegB could play a dynamic role in impeding mosquito multiplication and consequently impinging transmission of the dreadful dengue.

Resistance to the Fumigant Phosphine and Its Management in Insect Pests of Stored Products: A Global Perspective

Development of resistance in major grain insect pest species to the key fumigant phosphine (hydrogen phosphide) across the globe has put the viability and sustainability of phosphine in jeopardy. The resistance problem has been aggravated over the past two decades, due mostly to the lack of suitable alternatives matching the major attributes of phosphine, including its low price, ease of application, proven effectiveness against a broad pest spectrum, compatibility with most storage conditions, and international acceptance as a residue-free treatment. In this review, we critically analyze the published literature in the area of phosphine resistance with special emphasis on the methods available for detection of resistance, the genetic basis of resistance development, key management strategies, and research gaps that need to be addressed.

Development of a Quick Knockdown Test for Diagnosing Resistance to Phosphine in Sitophilus oryzae (Coleoptera: Curculionidae), a Major Pest of Stored Products

A key component in the management of resistance to fumigant phosphine in stored products pests is their early detection and implementation of control strategies. Currently, resistance testing involves exposing adults to a specific discriminating concentration over a fixed time period (20-48 h). Although it is widely adopted, this test takes significant time for assay preparation (up to 4 wk) as well as diagnosis (1-2 wk). To address these lacunae, we have established a 'quick knockdown test' using a key grain insect pest, rice weevil, Sitophilus oryzae (L.). Susceptible, weakly and strongly phosphine-resistant reference strains were exposed to a threshold concentration of phosphine over short exposure periods (min to h). The time to knockdown (KT) responses to phosphine were characterized at 2 (1,440 ppm) and 5 mg/liter (3,600 ppm). The time to 99.9% KT (KT99.9) at 2 mg/liter was 12.52 min for the susceptible adults, compared with 167.9 and 1,510 min in the case of weakly and strongly resistant phenotypes, respectively. As anticipated, increasing the concentration of phosphine to 5 mg/liter halved the KT99.9 (81.57 min) to separate weakly and strongly resistant populations than it was required at 2 mg/liter. We validated the KT99.9 value for the 5 mg/liter against field-derived populations of S. oryzae. The results were aligned with the existing Food and Agriculture Organization approach, confirming that the proposed 'quick test' is a reliable tool to rapidly diagnose resistance in this species.

A high-throughput system used to determine frequency and distribution of phosphine resistance across large geographical regions

BACKGROUND

Next-generation sequencing can enable genetic surveys of large numbers of individuals. We developed a genotyping-by-sequencing assay for detecting strong phosphine resistance alleles in the dihydrolipoamide dehydrogenase (dld) gene of Rhyzopertha dominica populations. The assay can estimate the distribution and frequency of resistance variants in thousands of individual insects in a single run.

RESULTS

We analysed 1435 individual insects collected over a 1-year period from 59 grain-storage sites including farms (n = 29) and central storages (n = 30) across eastern Australia. Resistance alleles were detected in 49% of samples, 38% of farms and 60% of central storages. Although multiple alleles were detected, only two resistance variants (P49S and K142E) were widespread and each appeared to have a distinct but overlapping geographical distribution.

CONCLUSION

The type of structure in which the grain is stored had a strong effect on resistance allele frequency. We observed higher frequencies of resistance alleles in bunker storages at central sites compared with other storage types. This contributed to the higher frequencies of resistance alleles in bulk-handling facilities relative to farms. The discovery of a storage structure that predisposes insects to resistance highlights the utility of our high-throughput assay system for improvement of phosphine resistance management practices. © 2018 Society of Chemical Industry.

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 72h. 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.

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.

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.