Faecal egg counts and nemabiome metabarcoding highlight the genomic complexity of equine cyathostomin communities and provide insight into their dynamics in a Scottish native pony herd

A set of nematode rRNA cistron databases and a primer assessment tool to enable more flexible and comprehensive metabarcoding

The ITS-2-rRNA has been particularly useful for nematode metabarcoding but does not resolve all phylogenetic relationships, and reference sequences are not available for many nematode species. This is a particular issue when metabarcoding complex communities such as wildlife parasites or terrestrial and aquatic free-living nematode communities. We have used markerDB to produce four databases of distinct regions of the rRNA cistron: the 18S rRNA gene, the 28S rRNA gene, the ITS-1 intergenic spacer and the region spanning ITS-1_5.8S_ITS-2. These databases comprise 2645, 254, 13,461 and 10,107 unique full-length sequences representing 1391, 204, 1837 and 1322 nematode species, respectively. The comparative analysis illustrates the complementary value but also reveals a better representation of Clade III, IV and V than Clade I and Clade II nematodes in each case. Although the ITS-1 database includes the largest number of unique full-length sequences, the 18S rRNA database provides the widest taxonomic coverage. We also developed PrimerTC, a tool to assess primer sequence conservation across any reference sequence database, and have applied it to evaluate a large number of previously published rRNA cistron primers. We identified sets of primers that currently provide the broadest taxonomic coverage for each rRNA marker across the nematode phylum. These new resources will facilitate more comprehensive metabarcoding of nematode communities using either short-read or long-read sequencing platforms. Further, PrimerTC is available as a simple WebApp to guide or assess PCR primer design for any genetic marker and/or taxonomic group beyond the nematode phylum.

Molecular species determination of cyathostomins from horses in Ireland

Cyathostomins are globally important equine parasites, responsible for both chronic and acute pathogenic effects. The occurrence of mixed infections with numerous cyathostomin species hinders our understanding of parasite epidemiology, host-parasite dynamics, and species pathogenicity. There have been few studies of cyathostomin species occurring in horses in Ireland, where temperate climatic conditions with year-round rainfall provide suitable conditions for infection of grazing animals with bursate nematodes. Here, we amplified and sequenced the ITS-2 region of adult worms harvested at post-mortem from eleven adult horses between August 2018 and June 2020, and recorded species prevalence and abundance of worms recovered from the caecum, right ventral colon and left dorsal colon, using both BLAST and IDTAXA for taxonomic attribution. Phylogenetic relationships and community composition were also recorded and compared with other relevant studies, including a global meta-analysis. Overall, our results agree with previous studies that there does not seem to be a major difference in cyathostomin species occurrence in equids in different geographical regions. We confirmed the results of other workers in relation to the difficulties in discriminating between Cylicostephanus calicatus and Coronocyclus coronatus on the basis of ITS-2 sequences.

Assessing periparturient ewe characteristics and nemabiome composition to guide targeted selective treatment for sustainable gastrointestinal nematodes control in sheep

Gastrointestinal nematodes (GINs) are a significant threat to the sustainability of global sheep production. Periparturient ewes play a key role in GIN epidemiology, with increased GIN faecal egg counts (FECs) in these ewes resulting in heavy pasture contamination that facilitates parasitic gastroenteritis in immunologically naïve lambs later during the grazing period. Traditionally, blanket anthelmintic treatment would suppress GIN egg outputs in these ewes and subsequent pasture contamination. However, farmers are now advised to implement targeted selective treatment (TST) to reduce anthelmintic use and subsequent anthelmintic resistance development, yet, there is currently limited evidence to determine optimal TST strategies in ewes. In this study, the characteristics of 226 ewes on seven Welsh farms were assessed postlambing to identify factors associated with their individual strongyle FECs using negative binomial mixed model analysis. Nemabiome analysis was conducted on 34 ewes across two study farms using the Oxford Nanopore MinIon platform with an aim of identifying factors associated with variations in ewe nemabiome composition within flocks. The best-fitted model of ewe FEC incorporated ewe body condition score, dag score, breed, and an interaction effect between ewe age and litter size as fixed factors. The addition of a mean FEC value for ewes of a specific litter size on each farm further improved model fit and reduced between-farm variance in the model. Nemabiome analysis revealed significant variation in within flock nemabiome diversity on individual farms, with significantly reduced nemabiome diversity recorded in ewes exhibiting dags and in twin-bearing ewes on respective farms, whilst T. circumcincta was present as a significantly higher proportion of the nemabiome in Suffolk ewes and twin bearing ewes (P < 0.05) in respective flocks. Our data demonstrate that commonly recorded ewe characteristics can be exploited to predict individual periparturient ewe FEC and subsequently may be used as a guide for TST strategies on sheep farms once specific TST thresholds are identified to deliver the optimal balance between minimal pasture contamination and maximal GIN refugia. This study is the first to utilise Oxford Nanopore MinIon sequencing to evaluate the nemabiome of sheep, and to molecularly assess the nemabiome of individual ruminants within a flock/herd, with results indicating that significant within flock variations in nemabiome composition which may have implications for TST and flock management strategies.

Nemabiome metabarcoding shows a high prevalence of Haemonchus contortus and predominance of Camelostrongylus mentulatus in alpaca herds in the northern UK

Gastrointestinal nematodes (GINs) are a common threat faced by pastoral livestock. Since their major introduction to the UK in the early 1990s, South American camelids have been cograzed with sheep, horses, and other livestock, allowing exposure to a range of GIN species. However, there have been no molecular-based studies to investigate the GIN populations present in these camelids. In the current study, we sampled nine alpaca herds from northern England and southern Scotland and used high-throughput metabarcoded sequencing to describe their GIN species composition. A total of 71 amplicon sequence variants (ASVs) were identified representing eight known GIN species. Haemonchus contortus was the most prevalent species found in almost all herds in significant proportions. The identification of H. contortus in other livestock species is unusual in the northern UK, implying that alpacas may be suitable hosts and potential reservoirs for infection in other hosts. In addition, the camelid-adapted GIN species Camelostrongylus mentulatus was identified predominantly in herds with higher faecal egg counts. These findings highlight the value of applying advanced molecular methods, such as nemabiome metabarcoding to describe the dynamics of gastrointestinal nematode infections in novel situations. The results provide a strong base for further studies involving cograzing animals to confirm the potential role of alpacas in transmitting GIN species between hosts.

A national survey of anthelmintic resistance in ascarid and strongylid nematodes in Australian Thoroughbred horses

This study quantified the extent of anthelmintic resistance (AR) in ascarid and strongylid nematodes against commonly used anthelmintics in Australian Thoroughbred horses. Faecal egg count reduction tests (FECRTs, n = 86) and egg reappearance period (ERP) tests were conducted on 22 farms across Australia. Faecal egg counts (FECs) were determined using the modified McMaster technique, and percent faecal egg count reduction (%FECR) was calculated using the Bayesian hierarchical model and hybrid Frequentist/Bayesian analysis method. The results were interpreted using old (published in 1992) and new (2023) research guidelines of the World Association for the Advancement of Veterinary Parasitology (WAAVP). The species composition of strongylid nematodes was detected utilising a DNA-metabarcoding method using pre- and post-treatment samples. Resistance was observed in strongylid nematodes to commonly used single-active and combination anthelmintics, including ivermectin (IVM %FECR range: 82%-92%; 95% lower credible interval (LCI) range: 80%-90%), abamectin (ABM: 73%-92%; 65%-88%), moxidectin (MOX: 89%-91%; 84%-89%), oxfendazole (OFZ: 0%-56%; 0%-31%) and its combination with pyrantel (OFZ + PYR: 0%-82%; 0%-78%). Resistance in Parascaris spp. was observed to IVM (10%-43%; 0%-36%), ABM (0%; 0%) and MOX (0%; 0%). When the new thresholds recommended by the WAAVP were used, AR was detected in six additional FECRTs for strongylids and three more tests for Parascaris spp., introducing resistance to OFZ and OFZ + PYR in the latter. Shortened ERPs (4-6 weeks) of strongylids were observed in 31 FECRTs in which AR was not detected at 2 weeks post-treatment for all the anthelmintics tested. Among cyathostomins, Cylicocyclus nassatus, Cylicostephanus longibursatus and Coronocyclus coronatus were the most prevalent species at 2 weeks post-treatment, whereas the main species appearing at five weeks following treatments with macrocyclic lactones were Cylicocyclus nassatus, Cylicostephanus longibursatus and Cylicocyclus ashworthi. After treatment with OFZ + PYR, the latter three, plus Coronocyclus coronatus and Cyathostomum catinatum, were detected at 5 weeks post-treatment. Overall, the study highlights the prevalence of AR in both ascarids and strongylid nematodes against commonly used anthelmintic products to control worms in Australian horses. The results indicate that ML combination products provided acceptable efficacy at 2 weeks. However, ERP calculations suggest that products work less effectively than previously measured. It is suggested to regularly monitor the efficacy of the anthelmintics and consider changing the worm control practices to better manage worms and AR in Australian horses.

Metabarcoding study to reveal the structural community of strongylid nematodes in domesticated horses in Thailand

BACKGROUND

Mixed strongylid infections significantly impact equine health and performance. Traditional microscopy-based methods exhibit limitations in accurately identifying strongylid species. Nemabiome deep amplicon sequencing approach previously succeeded in describing the strongylid communities in livestock including equids. However, there are no available studies that describe the structural communities of strongylid parasites in horses in Thailand. Therefore, this study was undertaken encompassing the ITS-2 rDNA metabarcoding assay to characterize strongylid species within horse fecal samples collected from a cohort of yearlings at the largest domesticated stud farm in Thailand. In addition, to investigate the capability of ITS-2 rDNA in assessing the phylogenetic relationships among the identified strongylid species.

RESULTS

The study identified 14 strongylid species in the examined equine populations, each with varying prevalence. Notably, Cylicocyclus nassatus and Cylicostephanus longibursatus were identified as the predominant species, with Strongylus spp. conspicuously absent. The phylogenetic analysis of 207 amplicon sequence variants (ASVs) displayed a complex relationship among the investigated cyathostomin species, with some species are positioned across multiple clades, demonstrating close associations with various species and genera.

CONCLUSION

The ITS-2 nemabiome sequencing technique provided a detailed picture of horse strongylid parasite species in the studied population. This establishes a foundation for future investigations into the resistance status of these parasites and enables efforts to mitigate their impact.

Patterns of Equine Small Strongyle Species Infection after Ivermectin Intervention in Thailand: Egg Reappearance Period and Nemabiome Metabarcoding Approach

The indiscriminate use of anthelmintics to control parasitic nematodes in horses has led to the emergence of anthelmintic resistance worldwide. However, there are no data available on using ivermectin for treating strongyle infections within domesticated horses in Thailand. Therefore, this study aimed to use the fecal egg count reduction (FECR) test to determine the strongylid egg reappearance period (ERP). Additionally, the nemabiome metabarcoding approach is incorporated to study patterns of strongyle species infection following ivermectin treatment. The study results indicate that, although ivermectin effectively eliminated adult strongyle parasites within two weeks post-treatment, the ERP was shortened to 6 weeks post-treatment with a mean FECR of 70.4% (95% CI 46.1-84.0). This potentially indicates a recent change in drug performance. In addition, nemabiome metabarcoding revealed that strongyle species have different levels of susceptibility in response to anthelmintic drugs. The reduction in ERP was associated with the early reappearance of specific species, dominated by Cylicostephanus longibursatus and Cylicocyclus nassatus, indicating the lower susceptibility of these species. In contrast, Poteriostomum imparidentatum, Triodontophorus nipponicus, and Triodontophorus serratus were not found post-treatment, indicating the high level of susceptibility of these species. This information is vital for comprehending the factors contributing to the emergence of resistance and for devising strategies to manage and control strongyle infections in horses.

Gastrointestinal parasite community structure in horses after the introduction of selective anthelmintic treatment strategies

A relatively new method to study the species richness and diversity of nematode parasites in grazing animals is to perform deep sequencing on composite samples containing a mixture of parasites. In this work, we compared species composition of strongyles in two groups of horses as a function of egg count and age, based on a DNA barcoding approach. Faecal egg counts and larval cultures were obtained from nearly 300 horses, i.e., domestic horses (n = 167) and trotters (n = 130) sampled nationwide. The second internal transcribed spacer region (ITS2) of strongyle nematodes in the larval cultures was first amplified using barcoded universal primers and then sequenced on the PacBio platform. Subsequently, bioinformatic sequence analysis was performed using SCATA to assign operational taxonomic units (OTU). Finally, species occurrence and composition were assessed using R. ITS2 sequences were found in the majority (89%) of larval samples. Sequencing yielded an average of 140 (26 to 503) reads per sample. The OTUs were assigned to 28 different taxa, of which all but three could be identified as species. The average relative abundance of the seven most abundant species (all Cyathostominae) accounted for 87% of the combined data set. The three species with the highest prevalence in both horse groups were Cyathostomum catinatum, Cylicocyclus nassatus and Cylicostephanus calicatus, and they were frequently found in different combinations with other species regardless of horse group. Interestingly, this result is largely consistent with a previous Swedish study based on morphological analysis of adult worms. In addition, two migratory strongylids (Strongylus vulgaris and S. edentatus) occurred in few domestic horses and trotters. Except for C. minutus and C. nassatus, which decreased with age, and C. catinatum and S. vulgaris, which increased, no specific trends were observed with respect to horse age. Taken together, these results are broadly consistent with data obtained before the introduction of selective targeted treatment in Sweden in 2007. All in all, our results suggest that this treatment strategy has not led to a significant change in strongyle nematode community structure in Swedish horses. The study also confirms that nemabiome analysis in combination with diversity index analysis is an objective method to study strongyle communities in horses.

Advances in diagnostic approaches to Fasciola infection in animals and humans: An overviews

Fasciolosis, caused by Fasciola hepatica and F. gigantica, is an impediment to the livestock industry’s expansion and has a massively negative socio-economic impact due to its widespread prevalence in livestock. It is a waterborne zoonosis affecting human populations in the countries where rural economies are associated with livestock rearing. Conventional diagnosis of Fasciola infection is done by detecting parasite eggs in the faeces of infected animals or by immunological methods. Accurate and quick immunodiagnosis of Fasciola infection in animals and humans is based on the detection of antibodies and specific antigens expressed in the prepatent stage of the parasite. Both molecular and serodiagnostic tests developed thus far have enhanced the reliability of Fasciola diagnosis in both man and animals but are not widely available in resource-poor nations. A pen-side diagnostic test based on a lateral flow assay or a DNA test like loop-mediated isothermal amplification (LAMP) would be simple, fast, and cost-effective, enabling clinicians to treat animals in a targeted manner and avoid the development of drug resistance to the limited flukicides. This review focuses on the recent advances made in the diagnosis of this parasite infection in animals and humans.

Understanding temporal and spatial distribution of intestinal nematodes of horses using faecal egg counts and DNA metabarcoding

This study reports the spatial and temporal distribution of ascarid and strongylid nematodes in Thoroughbred horses by age category across different climatic zones in Australia over an 18-month period. Faecal samples (n = 2046) from individual horses were analysed using the modified McMaster technique for faecal egg counts (FECs). Strongylids were identified using PCR-directed next-generation sequencing of the second internal transcribed spacer (ITS-2) of the nuclear ribosomal DNA. Yearlings had the highest prevalence (82%) of strongyle eggs followed by weanlings (79%), foals (58%), wet mares (49%) and dry mares (46%). For Parascaris spp., foals had the highest prevalence (35%) followed by weanlings (21%) and yearlings (10%). The highest mean FECs for Parascaris spp. were observed in foals (525 eggs per gram [EPG] of faeces) while those for strongyles were in yearlings (962 EPG). Among horses that were classified as adults at the time of sampling, 77% (860 of 1119) of mares were low (i.e., <250 EPG) strongyle egg-shedders. Mean strongyle FEC counts were highest in the Mediterranean (818 EPG) followed by summer (599 EPG), winter (442 EPG), and non-seasonal (413 EPG) rainfall zones. Twenty-six nematode species were detected, with Cylicostephanus longibursatus (26.5%), Cylicocyclus nassatus (23.7%) and Coronocyclus coronatus (20.5%) being the most frequently detected species. Their richness and relative abundance varied with horse age, season and climatic zone. In addition, Strongylus equinus and Triodontophorus spp. (T. brevicauda and T. serratus) were also detected. This comprehensive study elucidates spatial (climatic zone) and temporal (i.e., seasonal) trends in prevalence and burdens of intestinal nematodes in Australian horses using non-invasive conventional and molecular methods. The information presented in this study is crucial for developing integrated management strategies to control horse parasites in farmed horses.

Demonstration of reduced efficacy against cyathostomins without change in species composition after pyrantel embonate treatment in Swedish equine establishments

Consisting of approximately 50 different species, the cyathostomin parasites are ubiquitous in grazing horses. Co-infection with several species is common, and large burdens can cause the fatal disease of larval cyathostominosis. Due to intense anthelmintic drug use, cyathostomin resistance has developed to all available anthelmintic drug groups. Resistance to the anthelmintic drug pyrantel (PYR) has been documented in over 90% of studies published over the past two decades. In Sweden, a study performed in the early 2000s only confirmed resistance in 4.5% of farms. Further, prescription-only administration of equine anthelmintic drugs was enforced in Sweden in 2007. However, it is unknown if this conservative drug use has maintained PYR efficacy in cyathostomins. The aim of the present study was to investigate the effect of PYR on cyathostomin infection in Sweden using fecal egg count reduction tests (FECRTs). Further, the effect of PYR treatment on cyathostomin species composition was studied using metabarcoding. Sixteen farms with at least six horses excreting a minimum of 100 eggs per gram feces were included. Using the current World Association for the Advancement of Veterinary Parasitology (WAAVP) guidelines, PYR resistance was demonstrated in nine of farms, with seven farms showing full susceptibility. Farms with low biosecurity measures had significantly lower efficacy of PYR treatment. The most common cyathostomin species were Cylicocyclus nassatus, Cyathostomum catinatum, Cylicostephanus longibursatus, Cys. calicatus, Cys. goldi, Cys. minutus, Coronocyclus coronatus and Cya. pateratum, accounting for 97% of all sequence reads prior to treatment. Of these, Cyc. nassatus and Cya. catinatum had the highest occurrence, accounting for 68% of all sequence reads prior to PYR treatment. Treatment did not significantly affect the species composition. The results highlight the importance of drug efficacy testing when using PYR to treat cyathostomin infection, even when selective anthelmintic treatment and thus low treatment intensity, is used on the farm.

Prevalence and diversity of ascarid and strongylid nematodes in Australian Thoroughbred horses using next-generation sequencing and bioinformatic tools

The study presents the results of a cross-sectional survey to describe the epidemiology of ascarid and strongylid nematodes in horses, the impact of diverse climatic conditions on parasite diversity and the levels of faecal egg shedding in different age groups of managed Thoroughbred horses. Individual faecal samples (n = 1377) collected from 62 Thoroughbred farms across four climatic zones in Australia were analysed using the modified McMaster technique for faecal egg counts (FECs) and strongylid nematodes were identified utilising PCR-directed next-generation sequencing (NGS) of the second internal transcribed spacer of the nuclear ribosomal DNA (ITS-2). Across all age groups, the prevalence of ascarid and strongylid nematodes was 12% (95% confidence interval 10-14%) and 72% (70-74%), respectively. Based on strongylid FECs, yearlings had the highest prevalence (89%) followed by weanlings (83%), foals (79%), wet mares (61%), dry mares (59%) and stallions (54%). However, for Parascaris spp., foals had the highest prevalence (46%) followed by weanlings (32%) and yearlings (13%). The highest mean FECs for Parascaris spp. were observed in foals (418 eggs per gram [EPG] of faeces) while those for strongylids were in yearlings (1002 EPG). Of the adult horses (mares and stallions), 67% (489 of 729) and 11% (77 of 729) were low (i.e., ≤250 EPG) and moderate (i.e., 251-500 EPG) strongylid egg-shedders, respectively. Strongylid egg shedding varied across climatic zones, with the highest mean FECs in the summer rainfall (723 EPG) followed by non-seasonal rainfall (629 EPG), winter rainfall (613 EPG), and Mediterranean (606 EPG) rainfall zones. Twenty-three nematode species were detected using NGS, with Cylicostephanus longibursatus (28%), Cylicocyclus nassatus (23%) and Coronocyclus coronatus (23%), being the most abundant species. Three species of Strongylus (i.e., S. vulgaris, S. equinus and S. edentatus) were also detected. The nemabiome composition, species richness and relative abundance varied within horse age and between climatic zones. These empirical findings provide a comprehensive understanding of the prevalence of parasites within horse populations and the multifaceted factors that influence their occurrence, thereby allowing for the formulation of tailored strategies aimed at parasite control in domestic horses.

Evaluation of genome skimming to detect and characterise human and livestock helminths

The identification of gastrointestinal helminth infections of humans and livestock almost exclusively relies on the detection of eggs or larvae in faeces, followed by manual counting and morphological characterisation to differentiate species using microscopy-based techniques. However, molecular approaches based on the detection and quantification of parasite DNA are becoming more prevalent, increasing the sensitivity, specificity and throughput of diagnostic assays. High-throughput sequencing, from single PCR targets through to the analysis of whole genomes, offers significant promise towards providing information-rich data that may add value beyond traditional and conventional molecular approaches; however, thus far, its utility has not been fully explored to detect helminths in faecal samples. In this study, low-depth whole genome sequencing, i.e. genome skimming, has been applied to detect and characterise helminth diversity in a set of helminth-infected human and livestock faecal material. The strengths and limitations of this approach are evaluated using three methods to characterise and differentiate metagenomic sequencing data based on (i) mapping to whole mitochondrial genomes, (ii) whole genome assemblies, and (iii) a comprehensive internal transcribed spacer 2 (ITS2) database, together with validation using quantitative PCR (qPCR). Our analyses suggest that genome skimming can successfully identify most single and multi-species infections reported by qPCR and can provide sufficient coverage within some samples to resolve consensus mitochondrial genomes, thus facilitating phylogenetic analyses of selected genera, e.g. Ascaris spp. Key to this approach is both the availability and integrity of helminth reference genomes, some of which are currently contaminated with bacterial and host sequences. The success of genome skimming of faecal DNA is dependent on the availability of vouchered sequences of helminths spanning both taxonomic and geographic diversity, together with methods to detect or amplify minute quantities of parasite nucleic acids in mixed samples.

Species interactions, stability, and resilience of the gut microbiota - Helminth assemblage in horses

The nature and strength of interactions entertained among helminths and their host gut microbiota remain largely unexplored. Using 40 naturally infected Welsh ponies, we tracked the gut microbiota-cyathostomin temporal dynamics and stability before and following anthelmintic treatment and the associated host blood transcriptomic response. High shedders harbored 14 species of cyathostomins, dominated by Cylicocyclus nassatus. They exhibited a highly diverse and temporal dynamic gut microbiota, with butyrate-producing Clostridia likely driving the ecosystem steadiness and host tolerance toward cyathostomins infection. However, anthelmintic administration sharply bent the microbial community. It disrupted the ecosystem stability and the time-dependent network of interactions, affecting longer term microbial resilience. These observations highlight how anthelmintic treatments alter the triangular relationship of parasite, host, and gut microbiota and open new perspectives for adding nutritional intervention to current parasite management strategies.

Comparison of two molecular barcodes for the study of equine strongylid communities with amplicon sequencing

Basic knowledge on the biology and epidemiology of equine strongylid species still needs to be improved to contribute to the design of better parasite control strategies. Nemabiome metabarcoding is a convenient tool to quantify and identify species in bulk samples that could overcome the hurdle that cyathostomin morphological identification represents. To date, this approach has relied on the internal transcribed spacer 2 (ITS-2) of the ribosomal RNA gene, with a limited investigation of its predictive performance for cyathostomin communities. Using DNA pools of single cyathostomin worms, this study aimed to provide the first elements to compare performances of the ITS-2 and a cytochrome c oxidase subunit I (COI) barcode newly developed in this study. Barcode predictive abilities were compared across various mock community compositions of two, five and 11 individuals from distinct species. The amplification bias of each barcode was estimated. Results were also compared between various types of biological samples, i.e., eggs, infective larvae or adults. Bioinformatic parameters were chosen to yield the closest representation of the cyathostomin community for each barcode, underscoring the need for communities of known composition for metabarcoding purposes. Overall, the proposed COI barcode was suboptimal relative to the ITS-2 rDNA region, because of PCR amplification biases, reduced sensitivity and higher divergence from the expected community composition. Metabarcoding yielded consistent community composition across the three sample types. However, imperfect correlations were found between relative abundances from infective larvae and other life-stages for Cylicostephanus species using the ITS-2 barcode. While the results remain limited by the considered biological material, they suggest that additional improvements are needed for both the ITS-2 and COI barcodes.