Which larvae are they? Use of single larva for the molecular confirmation of Cooperia pectinata and Cooperia punctata in Australian cattle

A molecular assessment of Ostertagia leptospicularis and Spiculopteragia asymmetrica among wild fallow deer in Northern Ireland and implications for false detection of livestock-associated species

BACKGROUND

Wild deer populations utilizing livestock grazing areas risk cross-species transmission of gastrointestinal nematode parasites (GINs), including GINs with anthelmintic resistance (AR) traits. Wild deer have been shown to carry problematic GIN species such as Haemonchus contortus and Trichostrongylus species in the UK, but the presence of livestock GINs in Northern Ireland deer populations is unknown. Also, is it not known whether AR traits exist among GINs of deer such as Ostertagia leptospicularis and Spiculopteragia asymmetrica in pastureland where anthelmintics are heavily used.

METHODS

Adult-stage GIN samples were retrieved from Northern Irish wild fallow deer abomasa. Individual specimens were subject to a species-specific PCR analysis for common sheep and cattle GIN species with ITS-2 sequence analysis to validate species identities. In addition, the beta-tubulin gene was subject to sequencing to identify benzimidazole (BZ) resistance markers.

RESULTS

ITS-2 sequencing revealed O. leptospicularis and S. asymmetrica, but species-specific PCR yielded false-positive hits for H. contortus, Teladorsagia circimcincta, Trichostrongylus axei, T. colubriformis, T. vitrinus and Ostertagia ostertagi. For beta-tubulin, O. leptospicularis and S. asymmetrica yielded species-specific sequences at the E198 codon, but no resistance markers were identified in either species at positions 167, 198 or 200 of the coding region.

DISCUSSION

From this report, no GIN species of significance in livestock were identified among Northern Ireland fallow deer. However, false-positive PCR hits for sheep and cattle-associated GINs is concerning as the presence of deer species in livestock areas could impact both deer and livestock diagnostics and lead to overestimation of both GIN burden in deer and the role as of deer as drivers of these pathogens. ITS-2 sequences from both O. leptospicularis and S. asymmetrica show minor sequence variations to geographically distinct isolates. AR has been noted among GINs of deer but molecular analyses are lacking for GINs of wildlife. In producing the first beta-tubulin sequences for both O. leptospicularis and S. asymmetrica, we report no BZ resistance in this cohort.

CONCLUSIONS

This work contributes to genetic resources for wildlife species and considers the implications of such species when performing livestock GIN diagnostics.

A mixed amplicon metabarcoding and sequencing approach for surveillance of drug resistance to levamisole and benzimidazole in Haemonchus spp

Anthelmintic-resistant parasitic nematodes present a significant threat to sustainable livestock production worldwide. The ability to detect the emergence of anthelmintic resistance at an early stage, and therefore determine which drugs remain most effective, is crucial for minimising production losses. Despite many years of research into the molecular basis of anthelmintic resistance, no molecular-based tools are commercially available for the diagnosis of resistance as it emerges in field settings. We describe a mixed deep amplicon sequencing approach to determine the frequency of the levamisole (LEV)-resistant single nucleotide polymorphism (SNP) within arc-8 exon 4 (S168T) in Haemonchus spp., coupled with benzimidazole (BZ)-resistant SNPs within β-tubulin isotype-1 and the internal transcribed spacer-2 (ITS-2) nemabiome. This constitutes the first known multi-drug and multi-species molecular diagnostic developed for helminths of veterinary importance. Of the ovine, bovine, caprine and camelid Australian field isolates we tested, S168T was detected in the majority of Haemonchus spp. populations from sheep and goats, but rarely at a frequency greater than 16%; an arbitrary threshold we set based on whole genome sequencing (WGS) of LEV-resistant Haemonchus contortus GWBII. Overall, BZ resistance was far more prevalent in Haemonchus spp. than LEV resistance, confirming that LEV is still an effective anthelmintic class for small ruminants in New South Wales, Australia. The mixed amplicon metabarcoding approach described herein paves the way towards the use of large scale sequencing as a surveillance technology in the field, the results of which can be translated into evidence-based recommendations for the livestock sector.

Refugia or reservoir? Feral goats and their role in the maintenance and circulation of benzimidazole-resistant gastrointestinal nematodes on shared pastures

Gastrointestinal nematodes threaten the productivity of grazing livestock and anthelmintic resistance has emerged globally. It is broadly understood that wild ruminants living in sympatry with livestock act as a positive source of refugia for anthelmintic-susceptible nematodes. However, they might also act as reservoirs of anthelmintic-resistant nematodes, contributing to the spread of anthelmintic resistance at a regional scale. Here, we sampled managed sheep and cattle together with feral goats within the same property in New South Wales, Australia. Internal transcribed spacer 2 (ITS-2) nemabiome metabarcoding identified 12 gastrointestinal nematodes (Cooperia oncophora, Cooperia punctata, Haemonchus contortus, Haemonchus placei, Nematodirus spathiger, Ostertagia ostertagi, Teladorsagia circumcincta, Oesophagostomum radiatum, Oesophagostomum venulosum, Trichostrongylus axei, Trichostrongylus colubriformis and Trichostrongylus rugatus). Isotype-1 β-tubulin metabarcoding targeting benzimidazole resistance polymorphisms identified 6 of these nematode species (C. oncophora, C. punctata, H. contortus, H. placei, O. ostertagi and T. circumcincta), with the remaining 3 genera unable to be identified to the species level (Nematodirus, Oesophagostomum, Trichostrongylus). Both ITS-2 and β-tubulin metabarcoding showed the presence of a cryptic species of T. circumcincta, known from domestic goats in France. Of the gastrointestinal nematodes detected via β-tubulin metabarcoding, H. contortus, T. circumcincta, Nematodirus and Trichostrongylus exhibited the presence of at least one resistance genotype. We found that generalist gastrointestinal nematodes in untreated feral goats had a similarly high frequency of the benzimidazole-resistant F200Y polymorphism as those nematodes in sheep and cattle. This suggests cross-transmission and maintenance of the resistant genotype within the wild ruminant population, affirming that wild ruminants should be considered potential reservoirs of anthelmintic resistance.

A new diagnostic approach to fast-track and increase the accessibility of gastrointestinal nematode identification from faeces: FECPAKG2 egg nemabiome metabarcoding

Effective gastrointestinal nematode (GIN) management in livestock industries is becoming increasingly difficult due to the rise of anthelmintic resistance and changes in the temporal and geographical distribution of major GINs. Underpinning the response to these challenges is the need for a fast-tracked diagnostic identification technique, making it easier for livestock producers to make informed GIN management decisions. The traditional 'gold-standard' approach, larval culture followed by morphological differentiation, is laborious and potentially inaccurate. We developed a new diagnostic approach to identify GINs that integrates a remote-location digital faecal egg count platform, FECPAK^G2, with internal transcribed spacer 2 (ITS2) nemabiome metabarcoding. The technique involves a quick and simple protocol to harvest concentrated strongyle eggs from the FECPAK^G2 cassette utilising a repurposed pipette tip, followed by DNA isolation and Illumina next generation amplicon sequencing. The GIN compositions and alpha diversity generated by our FECPAK^G2 egg nemabiome metabarcoding approach was not significantly different to traditional morphological larval differentiation and nemabiome metabarcoding of larval and faecal samples. We demonstrated that storing FECPAK^G2 harvested eggs in either DNA isolation lysis buffer or 80% ethanol (v/v) had no impact on GIN identification outcomes for at least 60 days; enabling the transport of biological samples from their remote origins to a molecular diagnostic facility for nemabiome metabarcoding, in the absence of a cold chain. We discovered that sustained GIN egg embryonation in the lysis buffer storage solution lead to higher yields of DNA compared with ethanol-stored GIN eggs or faeces with GIN eggs. Taking advantage of an already well-established platform such as FECPAK^G2, and providing the livestock producers that use it with the option to identify GINs in their samples and contribute to large-scale GIN distribution and/or anthelmintic resistance surveys, is an important future direction for the FECPAK^G2 egg nemabiome metabarcoding approach.

A 4 year observation of gastrointestinal nematode egg counts, nemabiomes and the benzimidazole resistance genotypes of Teladorsagia circumcincta on a Scottish sheep farm

Anthelmintic resistance threatens the sustainability of sheep production globally. Advice regarding strategies to reduce the development of anthelmintic resistance incorporates the outcomes of modelling exercises. Further understanding of gastrointestinal nematode species diversity, and population dynamics and genetics (which may vary between species) is required to refine these models; and field studies combining faecal egg outputs, species composition and resistance genetics are needed to calibrate them. In this study, faecal samples were taken from ewes and lambs on a commercial farm in south-eastern Scotland at approximately 3 t-4 week intervals between spring and autumn over a period of 4 years. Faecal egg counts were performed on these samples, and L₃ were collected from pooled coprocultures. Deep amplicon sequencing was used to determine both the species composition of these L₃ and the proportions of benzimidazole-resistant single nucleotide polymorphisms in the isotype-1 β-tubulin locus of the predominant species, Teladorsagia circumcincta L₃. Despite consistent management throughout the study, the results show variation in gastrointestinal nematode species composition with time and between age groups, that was potentially associated with weather conditions. The F200Y benzimidazole resistance mutation is close to genetic fixation in the T. circumcincta population on this farm. There was no evidence of variation in isotype-1 β-tubulin single nucleotide polymorphisms frequency between age groups, and no genetic evidence of reversion to benzimidazole susceptibility, despite targeted benzimidazole usage. This study highlights the need to include speciation when investigating gastrointestinal nematode epidemiology and anthelmintic resistance, and serves as an example of how genetic data may be analysed alongside species diversity and faecal egg counts, when markers for other anthelmintic classes are identified.