In vitro insecticide resistance patterns in field strains of the sheep blowfly, Lucilia cuprina.

The control of the sheep blowfly relies on the use of insecticides. There have been several reports of in vitro and in vivo resistance to the most widely-used flystrike control chemical, dicyclanil. A recent report also described in vitro resistance to imidacloprid in a strain collected from a single property over three consecutive seasons that also showed resistance to dicyclanil. The present study aimed to use in vitro assays to examine five field-collected blowfly strains to determine if this co-occurrence of resistance to dicyclanil and imidacloprid was present more widely in field strains and to also measure resistance patterns to the other currently-used flystrike control chemicals. Each of the strains showed significant levels of resistance to both dicyclanil and imidacloprid: resistance factors at the IC50 of 9.1-23.8 for dicyclanil, and 8.7-14.1 for imidacloprid. Resistance factors at the IC95 ranged from 16.5 to 53.7, and 14.6-24.3 for dicyclanil and imidacloprid, respectively. Resistance factors were up to 8.5 for cyromazine at the IC95. Resistance to dicyclanil and imidacloprid was suppressed by co-treatment with the cytochrome P450 inhibitor, aminobenzotriazole, implicating this enzyme system in the observed resistances. We discuss the implications of the co-occurrence of resistance to dicyclanil and imidacloprid on insecticide rotation strategies for blowfly control. We also discuss the roles of insecticide resistance, environmental factors (e.g. rainfall), operational factors (e.g. insecticide application technique) and other animal health issues (e.g. scouring / diarrhoea) that together will impact on the likelihood of flystrike occurring at an earlier time point than expected after insecticide application.

Resistance to dicyclanil and imidacloprid in the sheep blowfly, Lucilia cuprina, in Australia.

BACKGROUND: The sheep blowfly, Lucila cuprina, is a myiasis-causing parasite responsible for significant production losses and welfare issues for the Australian sheep industry. Control relies largely on the use of insecticides. The pyrimidine compound, dicyclanil, is the predominant control chemical, although other insecticides also are used, including imidacloprid, ivermectin, cyromazine and spinosad. We investigated in vitro resistance patterns and mechanisms in field-collected blowfly strains. RESULTS: The Walgett 2019 strain showed significant levels of resistance to both dicyclanil and imidacloprid, with resistance factors at the IC50 of 26- and 17-fold, respectively, in in vitro bioassays. Co-treatment with the cytochrome P450 inhibitor, aminobenzotriazole, resulted in significant levels of synergism for dicyclanil and imidacloprid (synergism ratios of 7.2- and 6.1-fold, respectively), implicating cytochrome P450 in resistance to both insecticides. Cyp12d1 transcription levels were increased up to 40-fold throughout the larval life stages in the resistant strain compared to a reference susceptible strain, whereas transcription levels of some other cyp genes (6g1, 4d1, 28d1) did not differ between the strains. Similar resistance levels also were observed in flies collected from the same property in two subsequent years. CONCLUSION: This study indicates that in vitro resistance to both dicyclanil and imidacloprid in this field-collected blowfly strain is likely mediated by cytochrome P450, with Cyp12d1 implicated as the enzyme responsible; however, it remains possible that another P450 also may be involved. A common resistance mechanism for the two drugs has important implications for drug rotation strategies designed to prolong the useful life of flystrike control chemicals. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

A pair of adjacent genes, cry5Ad and orf2-5Ad, encode the typical N- and C-terminal regions of a Cry5Adelta-endotoxin as two separate proteins in Bacillus thuringiensis strain L366.

A new DNA sequence cry5Ad/orf2-5Ad (GenBank accession number EF219060) was isolated from Bacillus thuringiensis strain L366. This DNA sequence contains two ORFs: cry5Ad (a previously unreported member of the cry5A gene family) and orf2-5Ad. cry5Ad is unique among cry5A genes in that it encodes only the N-terminal region of a typical Cry5Adelta-endotoxin. The cry5Ad sequence includes homology blocks 1-5, which are present in most B. thuringiensisdelta-endotoxins. The usual C-terminal region of a Cry5Adelta-endotoxin (including homology blocks 6-8) is encoded by orf2-5Ad. Both proteins encoded by cry5Ad and orf2-5Ad were found in IPTG-induced Escherichia coli, after a copy of cry5Ad/orf2-5Ad was cloned into the pQE32 expression vector and transformed into pREP4 E. coli cells. Both proteins were also found in parasporal crystal inclusions of B. thuringiensis L366. Sequencing of cDNA derived from transformed E. coli cells showed that the two ORFs are transcribed as a single mRNA. Extracts prepared from the recombinant E. coli expressing Cry5Ad and Orf2-5Ad were not toxic to nematode larvae (Haemonchus contortus), indicating that these two proteins are most likely not responsible for the nematocidal activity seen previously in the B. thuringiensis strain L366.

Gene expression in the skin of Bos taurus and Bos indicus cattle infested with the cattle tick, Rhipicephalus (Boophilus) microplus.

The cattle tick Rhipicephalus microplus (formerly Boophilus microplus) is responsible for severe production losses to the cattle industry worldwide. It has long been known that different breeds of cattle can resist tick infestation to varying degrees; however, the mechanisms by which resistant cattle prevent heavy infestation are largely unknown. The aim of this study was to determine whether gene expression varied significantly between skin sampling sites (neck, chest and tail region), and whether changes in gene expression could be detected in samples taken at tick attachment sites (tick attached to skin sample) compared with samples taken from non-attachment sites (no tick attachment). We present here the results of an experiment examining the expression of a panel of forty-four genes in skin sections taken from Bos indicus (Brahman) cattle of known high resistance, and Bos taurus (Holstein-Friesian) cattle of known low resistance to the cattle tick. The forty-four genes chosen for this study included genes known to be involved in several immune processes, some structural genes, and some genes previously suggested to be of importance in tick resistance by other researchers. The expression of fifteen gene transcripts increased significantly in Holstein-Friesian skin samples at tick attachment sites. The higher expression of many genes involved in innate inflammatory processes in the Holstein-Friesian animals at tick attachment sites suggests this breed is exhibiting a non-directed pathological response to infestation. Of the forty-four genes analysed, no transcripts were detected in higher abundance at tick attachment sites in the Brahman cattle compared with similar samples from the Holstein-Friesian group, nor difference between attachment site and non-attachment site samples within the Brahman group. The results presented here suggest that the means by which these two cattle breeds respond to tick infestation differ and warrant further investigation.

Mutations in the Hco-mptl-1 gene in a field-derived monepantel-resistant isolate of Haemonchus contortus.

Resistance to the anthelmintic drug monepantel (Zolvix®) has emerged in parasitic worms infecting sheep and goats. The mechanism of resistance in these cases is unknown. The drug targets nicotinic acetylcholine receptors belonging to the nematode-specific DEG-3 subfamily. We examined the receptor gene, Hco-mptl-1, in a highly Zolvix®-resistant and a -susceptible isolate of the parasitic nematode Haemonchus contortus. cDNA coding for the full length receptor protein (Hco-MPTL-1) was present in all clones prepared from a pool of susceptible larvae (21/21 clones) and approximately 50% of those from the resistant isolate (17/33). On the other hand, the remaining clones from the resistant isolate showed various mutations that resulted in truncated predicted proteins, missing at least one transmembrane domain. The most common mutation (11/33 clones) resulted in the retention of intron 15, a premature stop codon, and a truncated protein. Sequencing of intron 15 genomic DNA showed very few SNPs in susceptible larvae and in 12/18 clones from resistant larvae, alongside the presence of at least 17 SNPs in the remaining resistant clones. The present study shows that the highly resistant isolate has a number of mutations in the drug target gene that would most-likely result in a non-functional receptor, thus rendering the larvae insensitive to the drug. The presence of many wild-type sequences in this highly-resistant population suggests that there was a significant presence of heterozygotes in the survivors of the field drench treatment from which the isolate was derived, and hence that at least some of the mutations may be dominant. Alternatively, their presence may be due to the additional influence of mutations at another locus contributing to the resistance phenotype. The presence of multiple separate mutations in the Hco-mptl-1 gene in this viable field-derived worm isolate may at least partly explain why resistance to Zolvix® has arisen rapidly in the field.

Insecticidal activities of histone deacetylase inhibitors against a dipteran parasite of sheep, Lucilia cuprina.

Histone deacetylase inhibitors (HDACi) are being investigated for the control of various human parasites. Here we investigate their potential as insecticides for the control of a major ecto-parasite of sheep, the Australian sheep blowfly, Lucilia cuprina. We assessed the ability of HDACi from various chemical classes to inhibit the development of blowfly larvae in vitro, and to inhibit HDAC activity in nuclear protein extracts prepared from blowfly eggs. The HDACi prodrug romidepsin, a cyclic depsipeptide that forms a thiolate, was the most potent inhibitor of larval growth, with equivalent or greater potency than three commercial blowfly insecticides. Other HDACi with potent activity were hydroxamic acids (trichostatin, CUDC-907, AR-42), a thioester (KD5170), a disulphide (Psammaplin A), and a cyclic tetrapeptide bearing a ketone (apicidin). On the other hand, no insecticidal activity was observed for certain other hydroxamic acids, fatty acids, and the sesquiterpene lactone parthenolide. The structural diversity of the 31 hydroxamic acids examined here revealed some structural requirements for insecticidal activity; for example, among compounds with flexible linear zinc-binding extensions, greater potency was observed in the presence of branched capping groups that likely make multiple interactions with the blowfly HDAC enzymes. The insecticidal activity correlated with inhibition of HDAC activity in blowfly nuclear protein extracts, indicating that the toxicity was most likely due to inhibition of HDAC enzymes in the blowfly larvae. The inhibitor potencies against blowfly larvae are different from inhibition of human HDACs, suggesting some selectivity for human over blowfly HDACs, and a potential for developing compounds with the inverse selectivity. In summary, these novel findings support blowfly HDAC enzymes as new targets for blowfly control, and point to development of HDAC inhibitors as a promising new class of insecticides.

RNA interference in Haemonchus contortus: suppression of beta-tubulin gene expression in L3, L4 and adult worms in vitro.

We have used RNAi to target two beta-tubulin genes in the parasitic stages of Haemonchus contortus in vitro. Soaking exsheathed-L3, L4 and adult worms in medium containing dsRNA resulted in a significant decrease (greater than 1000-fold in some cases) in the expression of the specific beta-tubulin transcript, as measured by quantitative PCR. During the initial 24h exposure to the dsRNA, the gene suppression effect was quite specific to the targeted gene. Six days after initial exposure to dsRNA, treated L3 worms were less able to migrate through a filter mesh, indicating decreased motility, and showed less development to the L4 stage than control larvae. The gene suppression effect occurs in exsheathed L3 larvae despite the fact that this life stage does not have functioning mouthparts, indicating that uptake of dsRNA does not depend on its ingestion. Suppression occurred with dsRNA presented in 'naked' or liposome-encapsulated forms, indicating that a liposome formulation was not necessary for uptake to occur. Adult worms also showed significant gene suppression, however, they did not show any reduced motility compared to controls over a 3-day period. Adult worms treated with ivermectin to paralyse their pharynx, still showed significant gene suppression, again suggesting that uptake of dsRNA does not require ingestion. We have shown that soaking in dsRNA is an effective method for RNAi with the parasitic stages of H. contortus, and, hence, may offer significant potential as a tool for studying gene function in this parasite species.

Effects of in vitro exposure to ivermectin and levamisole on the expression patterns of ABC transporters in Haemonchus contortus larvae.

This study investigated the interaction of ATP binding cassette (ABC) transport proteins with ivermectin (IVM) and levamisole (LEV) in larvae of susceptible and resistant isolates of Haemonchus contortus in vitro by measuring transcription patterns following exposure to these anthelmintics. Furthermore, we studied the consequences of drug exposure by measuring the sensitivity of L3 to subsequent exposure to higher drug concentrations using larval migration assays. The most highly transcribed transporter genes in both susceptible and resistant L3 were pgp-9.3, abcf-1, mrp-5, abcf-2, pgp-3, and pgp-10. The resistant isolate showed significantly higher transcription of pgp-1, pgp-9.1 and pgp-9.2 compared to the susceptible isolate. Five P-gp genes and the haf-6 gene showed significantly higher transcription (up to 12.6-fold) after 3 h exposure to IVM in the resistant isolate. Similarly, five P-gp genes, haf-6 and abcf-1 were transcribed at significantly higher levels (up to 10.3-fold) following 3 h exposure to LEV in this isolate. On the other hand, there were no significant changes in transcriptional patterns of all transporter genes in the susceptible isolate following 3 and 6 h exposure to IVM or LEV. In contrast to these isolate-specific transcription changes, both isolates showed an increase in R-123 efflux following exposure to the drugs, suggesting that the drugs stimulated activity of existing transporter proteins in both isolates. Exposure of resistant larvae to IVM or LEV resulted, in some instances, in an increase in the proportion of the population able to migrate at the highest IVM concentrations in subsequent migration assays. The significant increase in transcription of some ABC transporter genes following 3 h exposure to both IVM and LEV in the resistant isolate only, suggests that an ability to rapidly upregulate protective pathways in response to drugs may be a component of the resistance displayed by this isolate.

Increased expression of ATP binding cassette transporter genes following exposure of Haemonchus contortus larvae to a high concentration of monepantel in vitro.

BACKGROUND: There is some evidence that ATP binding cassette (ABC) transporters play a role in resistance to anthelmintics, particularly against macrocyclic lactones. Some anthelmintics, including ivermectin (IVM), have been shown to induce transcription of multiple ABC transporters in nematodes; however, the effects of monepantel (MPL) on transcription of these transporter genes has not been studied. METHODS: Larvae of two MPL-susceptible isolates of Haemonchus contortus were exposed to MPL at two concentrations (2.5 and 250 µg/ml) for periods of 3, 6 and 24 h. Transcription levels of sixteen ABC transporter genes were measured at the end of the incubation periods. The consequences of MPL exposure were examined by measuring rhodamine-123 efflux from the larvae, and their sensitivity to subsequent treatment with IVM or levamisole. RESULTS: Multiple ABC transporter genes showed significantly higher transcription in both worm isolates following exposure to MPL at 250 µg/ml for 3, 6 or 24 h, particularly the P-glycoprotein (P-gp) genes pgp-11, pgp-12 and pgp-14. Of these, only pgp-11 maintained the elevated levels 24 h after the end of the drug exposure period. In contrast, there was only a single instance of low-level upregulation as a result of exposure to MPL at 2.5 µg/ml. Larvae exposed to MPL at 250 µg/ml showed an increased efflux of rhodamine-123 and a proportion of the larval population showed an ability to subsequently tolerate higher concentrations of IVM in migration assays. There was no increased tolerance to IVM following pre-exposure to MPL at 2.5 µg/ml. CONCLUSIONS: Exposure of H. contortus larvae to 250 µg/ml MPL results in increased transcription of multiple transporter genes and increased R-123 efflux. The subsequent ability of a proportion of the larvae to tolerate IVM suggests a protective role of ABC transporters across different chemical entities. However, these observations were only made at a concentration of MPL well above that experienced by parasitic life stages in vivo, and hence their significance remains unclear.

Histone deacetylase enzymes as drug targets for the control of the sheep blowfly, Lucilia cuprina.

The Australian sheep blowfly, Lucilia cuprina, is an ecto-parasite that causes significant economic losses in the sheep industry. Emerging resistance to insecticides used to protect sheep from this parasite is driving the search for new drugs that act via different mechanisms. Inhibitors of histone deacetylases (HDACs), enzymes essential for regulating eukaryotic gene transcription, are prospective new insecticides based on their capacity to kill human parasites. The blowfly genome was found here to contain five HDAC genes corresponding to human HDACs 1, 3, 4, 6 and 11. The catalytic domains of blowfly HDACs 1 and 3 have high sequence identity with corresponding human and other Dipteran insect HDACs (Musca domestica and Drosophila melanogaster). On the other hand, HDACs 4, 6 and 11 from the blowfly and the other Dipteran species showed up to 53% difference in catalytic domain amino acids from corresponding human sequences, suggesting the possibility of developing HDAC inhibitors specific for insects as desired for a commercial insecticide. Differences in transcription patterns for different blowfly HDACs through the life cycle, and between the sexes of adult flies, suggest different functions in regulating gene transcription within this organism and possibly different vulnerabilities. Data that supports HDACs as possible new insecticide targets is the finding that trichostatin A and suberoylanilide hydroxamic acid retarded growth of early instar blowfly larvae in vitro, and reduced the pupation rate. Trichostatin A was 8-fold less potent than the commercial insecticide cyromazine in inhibiting larval growth. Our results support further development of inhibitors of blowfly HDACs with selectivity over human and other mammalian HDACs as a new class of prospective insecticides for sheep blowfly.

Lucilia cuprina genome unlocks parasitic fly biology to underpin future interventions.

Lucilia cuprina is a parasitic fly of major economic importance worldwide. Larvae of this fly invade their animal host, feed on tissues and excretions and progressively cause severe skin disease (myiasis). Here we report the sequence and annotation of the 458-megabase draft genome of Lucilia cuprina. Analyses of this genome and the 14,544 predicted protein-encoding genes provide unique insights into the fly's molecular biology, interactions with the host animal and insecticide resistance. These insights have broad implications for designing new methods for the prevention and control of myiasis.

Relative level of thiabendazole resistance associated with the E198A and F200Y SNPs in larvae of a multi-drug resistant isolate of Haemonchus contortus.

While the F200Y SNP in the beta-tubulin gene is most commonly associated with benzimidazole resistance in trichostrongylid nematodes, other SNPs as well as drug efflux pathways have been implicated in the resistance. The relative contributions of all these mechanisms are not understood sufficiently to allow expected drug efficacy to be inferred from molecular data. As a component of developing better means to interpret molecular resistance tests, the present study utilised a drug resistant Haemonchus contortus isolate which possesses two of the principal benzimidazole resistance SNPs (E198A and F200Y) in order to assess the relative degree of resistance conferred by the two SNPs. We exposed larvae to a range of thiabendazole concentrations in in vitro development assays, and collected the surviving L3 larvae at each drug concentration to establish sub-populations showing increasing levels of resistance. We then sequenced the isotype 1 beta-tubulin gene in pooled larval samples, and measured allele frequencies at the two SNP positions. The frequency of the resistance allele at the 198 position increased as the thiabendazole concentration increased, while the frequency of the resistance allele at the 200 position decreased. Genotyping of individual larvae showed that the highest drug concentration was associated with the removal of all genotypes except for homozygous resistance at the 198 position alongside homozygous susceptible at the 200 position. This indicates that, at least for larval life stages, the E198A SNP is able to confer higher levels of resistance to benzimidazole drugs than the F200Y SNP, and that the homozygosity at 198 in the highly resistant individuals is mutually exclusive with heterozygosity or resistant homozygosity at the 200 position. This study illustrates the need to understand the relative contributions of different resistance mechanisms in order to maximise the degree to which molecular tests are able to inform on drug resistance phenotype.

Activation of several key components of the epidermal differentiation pathway in cattle following infestation with the cattle tick, Rhipicephalus (Boophilus) microplus.

The cattle tick, Rhipicephalus (Boophilus) microplus, and the diseases it transmits pose a persistent threat to tropical beef production. Genetic selection of host resistance has become the method of choice for non-chemical control of cattle tick. Previous studies have suggested that larval stages are most susceptible to host resistance mechanisms. To gain insights into the molecular basis of host resistance that occurs during R. microplus attachment, we assessed the abundance of proteins (by isobaric tag for relative and absolute quantitation (iTRAQ) and Western blot analyses) and mRNAs (by quantitative reverse transcription PCR (qRT-PCR)) in skin adjacent to tick bite sites from high tick-resistant (HR) and low tick-resistant (LR) Belmont Red cattle following challenge with cattle tick. We showed substantially higher expression of the basal epidermal keratins KRT5 and KRT14, the lipid processing protein, lipocalin 9 (LCN9), the epidermal barrier catalysing enzyme transglutaminase 1 (TGM1), and the transcriptional regulator B lymphocyte-induced maturation protein 1 (Blimp1) in HR skin. Our data reveals the essential role of the epidermal permeability barrier in conferring greater resistance of cattle to tick infestation, and suggest that the physical structure of the epidermal layers of the skin may represent the first line of defence against ectoparasite invasion.

The complexity of the secreted NPA and FAR lipid-binding protein families of Haemonchus contortus revealed by an iterative proteomics-bioinformatics approach.

Two different classes of small nematode specific lipid-binding proteins, the nematode polyprotein allergens/antigens (NPAs) and the fatty acid- and retinol-binding (FAR) proteins, are secreted by helminth parasites. Until now, there was no evidence of the expression or secretion of these two families of proteins in Haemonchus contortus. In this study, we applied proteomic and bioinformatic tools in an iterative manner to reveal the expression and complexity of these proteins in the excretory/secretory products (ESP) of adult H. contortus at the protein and gene levels. Initial examination of the mass spectra of ESP fractions against standard databases returned nine peptides mapping to Ostertagia ostertagi NPA and FAR sequences. Searches of the H. contortus EST and genomic contig databases with the O. ostertagi and Caenorhabditis elegans homologues retrieved diverse sequences encoding H. contortus NPA and FAR proteins. H. contortus sequences were then integrated into a customized database and a new search of the mass spectra achieved a 10-fold improvement in coverage of the predicted H. contortus NPAs. The final analyses of the mass spectra achieved 49-60% coverage of H. contortus NPAs and 7-47% coverage of H. contortus FARs. Moreover, the diversity in structures of the encoding genes was revealed by assembling the genomic sequence data with predicted protein sequences confirmed by the peptide evidence. We predict there are at least one Hc-NPA gene and six Hc-FAR genes in H. contortus, and life stage gene expression of Hc-FAR-1 to -6 revealed unique transcription patterns for each of these genes.