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.

Integrated pest management in cotton: exploiting behaviour-modifying (semiochemical) compounds for managing cotton pests

We review here research on semiochemicals for cotton pest management carried out in successive Cotton Co-operative Research Centres from 1998 to 2012. Australian cotton is now dominated by transgenic (Bt) varieties, which provide a strong platform for integrated pest management of key pests such as Helicoverpa spp., but new technologies are required to manage the development of resistance in Helicoverpa spp. to transgenic cotton and the problems posed by emerging and secondary pests, especially sucking insects. A long-range attractant for Helicoverpa moths, based on plant volatiles, has been commercialised as Magnet®. The product has substantial area-wide impacts on moth populations, and only limited effects on beneficial insects. Potential roles are being investigated for this product in resistance management of Helicoverpa spp. on transgenic cotton. Short-range, non-volatile compounds on organ surfaces of plants that do not support development of Helicoverpa spp. have been identified; these compounds deter feeding or oviposition, or are toxic to insect pests. One such product, Sero X®, is effective on Helicoverpa spp. and sucking pests such as whiteflies (Bemisia tabaci), green mirids (Creontiades dilutus), and other hemipteran insects, and is in the advanced stages of commercialisation.

Screening Corymbia populations for resistance to Puccinia psidii

The exotic rust pathogen Puccinia psidii is now widespread along the east coast of Australia from temperate Victoria to tropical far north Queensland, with a current host range exceeding 200 species from 37 myrtaceous genera. To determine the threat P. psidii poses to plantation and native eucalypts, artificial inoculation was used to screen germplasm of spotted gum (Corymbia spp.) for resistance to the biotype of P. psidii that has become established in Australia. The objective was to characterize resistance to P. psidii within the Corymbia species complex so that management strategies for the deployment of germplasm from existing breeding programmes of these spotted gum species could be developed. Symptom development initiated 7 days after inoculation, with resistant and susceptible seedlings identified within all species, provenances and families. Inter- and intraspecific variability in rust resistance was observed among spotted gum species. There was no apparent relationship between climatic conditions at the provenance origin and disease resistance. The heritability estimates for all assessments are moderate to high and indicate a significant level of additive genetic variance for rust resistance within the populations. The results of this study clearly identify potential to select for resistance at the family level within the tested populations. While the potential for P. psidii to detrimentally impact upon Corymbia in the nursery and in young plantations was demonstrated, estimations of the heritability of resistance suggest that efforts to enhance this trait through breeding have reasonable prospects for success.

Managing glyphosate resistance in Australian cotton farming: modelling shows how to delay evolution and maintain long-term population control

Glyphosate resistance is a rapidly developing threat to profitability in Australian cotton farming. Resistance causes an immediate reduction in the effectiveness of in-crop weed control in glyphosate-resistant transgenic cotton and summer fallows. Although strategies for delaying glyphosate resistance and those for managing resistant populations are qualitatively similar, the longer resistance can be delayed, the longer cotton growers will have choice over which tactics to apply and when to apply them. Effective strategies to avoid, delay, and manage resistance are thus of substantial value. We used a model of glyphosate resistance dynamics to perform simulations of resistance evolution in Sonchus oleraceus (common sowthistle) and Echinochloa colona (awnless barnyard grass) under a range of resistance prevention, delaying, and management strategies. From these simulations, we identified several elements that could contribute to effective glyphosate resistance prevention and management strategies. (i) Controlling glyphosate survivors is the most robust approach to delaying or preventing resistance. High-efficacy, high-frequency survivor control almost doubled the useful lifespan of glyphosate from 13 to 25 years even with glyphosate alone used in summer fallows. (ii) Two non-glyphosate tactics in-crop plus two in-summer fallows is the minimum intervention required for long-term delays in resistance evolution. (iii) Pre-emergence herbicides are important, but should be backed up with non-glyphosate knockdowns and strategic tillage; replacing a late-season, pre-emergence herbicide with inter-row tillage was predicted to delay glyphosate resistance by 4 years in awnless barnyard grass. (iv) Weed species' ecological characteristics, particularly seed bank dynamics, have an impact on the effectiveness of resistance strategies; S. oleraceus, because of its propensity to emerge year-round, was less exposed to selection with glyphosate than E. colona, resulting in an extra 5 years of glyphosate usefulness (18 v. 13 years) even in the most rapid cases of resistance evolution. Delaying tactics are thus available that can provide some or many years of continued glyphosate efficacy. If glyphosate-resistant cotton cropping is to remain profitable in Australian farming systems in the long-term, however, growers must adapt to the probability that they will have to deal with summer weeds that are no longer susceptible to glyphosate. Robust resistance management systems will need to include a diversity of weed control options, used appropriately.

Herbicide resistance modelling: Past, present and future

Computer simulation modelling is an essential aid in building an integrated understanding of how different factors interact to affect the evolutionary and population dynamics of herbicide resistance, and thus in helping to predict and manage how agricultural systems will be affected. In this review, we first discuss why computer simulation modelling is such an important tool and framework for dealing with herbicide resistance. We then explain what questions related to herbicide resistance have been addressed to date using simulation modelling, and discuss the modelling approaches that have been used, focusing first on the earlier, more general approaches, and then on some newer, more innovative approaches. We then consider how these approaches could be further developed in the future, by drawing on modelling techniques that are already employed in other areas, such as individual-based and spatially explicit modelling approaches, as well as the possibility of better representing genetics, competition and economics, and finally the questions and issues of importance to herbicide resistance research and management that could be addressed using these new approaches are discussed. We conclude that it is necessary to proceed with caution when increasing the complexity of models by adding new details, but, with appropriate care, more detailed models will make it possible to integrate more current knowledge in order better to understand, predict and ultimately manage the evolution of herbicide resistance. © 2014 Society of Chemical Industry.

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.

Evaluation of the efficacy and economics of irrigation management, plant resistance and Brassica(spot)(TM) models for management of white blister on Brassica crops

Options for the integrated management of white blister (caused by Albugo candida) of Brassica crops include the use of well timed overhead irrigation, resistant cultivars, programs of weekly fungicide sprays or strategic fungicide applications based on the disease risk prediction model, Brassica(spot)(TM). Initial systematic surveys of radish producers near Melbourne, Victoria, indicated that crops irrigated overhead in the morning (0800-1200 h) had a lower incidence of white blister than those irrigated overhead in the evening (2000-2400 h). A field trial was conducted from July to November 2008 on a broccoli crop located west of Melbourne to determine the efficacy and economics of different practices used for white blister control, modifying irrigation timing, growing a resistant cultivar and timing spray applications based on Brassica(spot)(TM). Growing the resistant cultivar, 'Tyson', instead of the susceptible cultivar, 'Ironman', reduced disease incidence on broccoli heads by 99 %. Overhead irrigation at 0400 h instead of 2000 h reduced disease incidence by 58 %. A weekly spray program or a spray regime based on either of two versions of the Brassica(spot)(TM) model provided similar disease control and reduced disease incidence by 72 to 83 %. However, use of the Brassica(spot)(TM) models greatly reduced the number of sprays required for control from 14 to one or two. An economic analysis showed that growing the more resistant cultivar increased farm profit per ha by 12 %, choosing morning irrigation by 3 % and using the disease risk predictive models compared with weekly sprays by 15 %. The disease risk predictive models were 4 % more profitable than the unsprayed control.

The plasticity of NBS resistance genes in sorghum is driven by multiple evolutionary processes

Background Increased disease resistance is a key target of cereal breeding programs, with disease outbreaks continuing to threaten global food production, particularly in Africa. Of the disease resistance gene families, the nucleotide-binding site plus leucine-rich repeat (NBS-LRR) family is the most prevalent and ancient and is also one of the largest gene families known in plants. The sequence diversity in NBS-encoding genes was explored in sorghum, a critical food staple in Africa, with comparisons to rice and maize and with comparisons to fungal pathogen resistance QTL. Results In sorghum, NBS-encoding genes had significantly higher diversity in comparison to non NBS-encoding genes and were significantly enriched in regions of the genome under purifying and balancing selection, both through domestication and improvement. Ancestral genes, pre-dating species divergence, were more abundant in regions with signatures of selection than in regions not under selection. Sorghum NBS-encoding genes were also significantly enriched in the regions of the genome containing fungal pathogen disease resistance QTL; with the diversity of the NBS-encoding genes influenced by the type of co-locating biotic stress resistance QTL. Conclusions NBS-encoding genes are under strong selection pressure in sorghum, through the contrasting evolutionary processes of purifying and balancing selection. Such contrasting evolutionary processes have impacted ancestral genes more than species-specific genes. Fungal disease resistance hot-spots in the genome, with resistance against multiple pathogens, provides further insight into the mechanisms that cereals use in the “arms race” with rapidly evolving pathogens in addition to providing plant breeders with selection targets for fast-tracking the development of high performing varieties with more durable pathogen resistance.

Predictive modelling to aid the regional-scale management of a vertebrate pest

Extensive resources are allocated to managing vertebrate pests, yet spatial understanding of pest threats, and how they respond to management, is limited at the regional scale where much decision-making is undertaken. We provide regional-scale spatial models and management guidance for European rabbits (Oryctolagus cuniculus) in a 260,791 km(2) region in Australia by determining habitat suitability, habitat susceptibility and the effects of the primary rabbit management options (barrier fence, shooting and baiting and warren ripping) or changing predation or disease control levels. A participatory modelling approach was used to develop a Bayesian network which captured the main drivers of suitability and spread, which in turn was linked spatially to develop high resolution risk maps. Policy-makers, rabbit managers and technical experts were responsible for defining the questions the model needed to address, and for subsequently developing and parameterising the model. Habitat suitability was determined by conditions required for warren-building and by above-ground requirements, such as food and harbour, and habitat susceptibility by the distance from current distributions, habitat suitability, and the costs of traversing habitats of different quality. At least one-third of the region had a high probability of being highly suitable (support high rabbit densities), with the model supported by validation. Habitat susceptibility was largely restricted by the current known rabbit distribution. Warren ripping was the most effective control option as warrens were considered essential for rabbit persistence. The anticipated increase in disease resistance was predicted to increase the probability of moderately suitable habitat becoming highly suitable, but not increase the at-risk area. We demonstrate that it is possible to build spatial models to guide regional-level management of vertebrate pests which use the best available knowledge and capture fine spatial-scale processes.

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.

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.

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.