Exploring the Gastrointestinal "Nemabiome": Deep Amplicon Sequencing to Quantify the Species Composition of Parasitic Nematode Communities

Molecular analysis of feces reveals gastrointestinal nematodes in reintroduced wild asses of the Negev desert

Reintroduced animals face disease risks, potentially impacting both the reintroduced and the local wildlife/domestic populations. This study focuses on the Asiatic wild asses (Equus hemionus) reintroduced to the Negev desert in southern Israel. Despite potential threats of disease spill-over to and from domesticated donkeys and horses in the area, there are no records of the gastrointestinal nematodes (GIN) of the wild ass population. We used DNA metabarcoding on fecal samples of wild asses collected across seasons and habitats, near water sources that they frequently use. Ten GIN species were detected in the feces, nine belonging to the family Strongylidae, which commonly infects and causes disease in equids worldwide, such as horses, zebras, and donkeys. Some of these Strongylidae species are also found in domesticated equids in Israel, thus raising concerns regarding potential parasite transmission between wild and domestic animals. The high prevalence of certain GIN species suggests frequent transmission, likely due to the congregation of the wild asses around water sources. While we observed statistically significant variations in some GIN species across seasons and habitats, we did not find clear overall differences between GIN communities. DNA metabarcoding proves to be a valuable tool for identifying GIN species in wild animals, with potential applications in monitoring their health and preventing disease transmission to and from domestic animals.

Efficacy of a Haemonchus contortus vaccine under field conditions in young alpacas

Anthelmintic resistance to Haemonchus contortus creates increasing management challenges with small ruminants and camelids. The commercial vaccine, Barbervax®, contains H11 and H-gal-GP antigens, derived from gut mucosal membrane enzymes of H. contortus involved in digesting blood. Antibody neutralization of these antigens causes failure of H. contortus to digest blood, resulting in parasite death. H11 and H-gal-GP are considered "hidden" antigens, meaning the host immune system does not encounter these proteins under natural infection. Therefore, repeat immunization is required to maintain protective humoral responses. One previous study evaluated the safety of Barbervax® in camelids but the efficacy could not be assessed due to lack of successful infection in the controls. The objective of the current study was to evaluate clinical parameters of anemia, fecal egg counts (FECs), and humoral immune responses of healthy alpacas after immunizing with Barbervax® compared to non-vaccinated controls, all under natural environmental exposure on parasite-laden pastures. A crossover-like study was performed where twenty alpacas (298 ± 66 days of age) were assigned to be initially vaccinated with Barbervax® (n=10) or receive no treatment (n=10). Three doses of Barbervax® were administered at three-week intervals. Feces and blood were collected on Day -10, 0, 21, 43, 64, 85, 106, and 135 to evaluate FECs, packed cell volume (PCV), and antibody titers. Each group was kept on separate adjacent pastures. Tracer sheep (n=2 per study group) were introduced on Day 43 for a three-week period to ensure parasite acquisition. For the crossover-like component on Day 85, the initial non-vaccinated group was administered Barbervax® with dosing repeated on Day 106 and 135. Results indicated all initially vaccinated alpacas produced antibody titers to vaccine antigen that corresponded to lower mean FECs compared to the initially non-vaccinated group. A reduced mean FEC in the vaccinate group was observed 21 days after peak antibody titers. Similarly, when pooled vaccinate antibody titers were noted to wane on Day 106, an increase in FEC was observed at the following time point (Day 135). Conclusions from our study support the use of Barbervax® to reduce H. contortus burdens in alpacas. Furthermore, a less than 30-day lag time between antibody titer and resultant effect in FECs was observed. Additional studies assessing the ability of Barbervax® to reduce H. contortus burdens during subsequent grazing seasons would provide even greater information regarding the use of Barbervax® within alpaca herds to modulate H. contortus infections, refugia, and anthelmintic use.

Towards precision parasite management for livestock gastrointestinal nematodes in 2030

The management of parasitic nematodes calls for a shift from conventional, indiscriminate, anthelmintic use to a more precise approach, directed by diagnostics. We should accept those parasite infection intensities that have minimal impact on production and welfare rather than attempt to eliminate them. The diagnostic toolbox for gastrointestinal nematodes (GINs) faces challenges due to anthelmintic resistance (AR), which is species-specific, drug-class-specific, and varies by region. We discuss which traditional tools may become obsolete and which tools need development to gain widespread use. Social science research highlights the need for dialogue between farmers and veterinarians that emphasises effective parasite management and upskilling the veterinary workforce for more sustainable practices centred on diagnostics to be adopted in practice by 2030.

Comparison of commercially available DNA and RNA extraction kits for wildlife feces collected from the environment

Wildlife fecal samples were collected across two Air Force installations to evaluate the effectiveness of commercially available DNA and RNA extraction kits. Four DNA kits, two DNA/RNA kits and one RNA only kit were used. Sample extracts were evaluated on nucleic acid concentration, TapeStation DNA or RNA Integrity Number values and presence of PCR inhibitors. For the DNA kits, PFP produced higher concentrations compared with PLM and RPM, while MWFM gave higher DNA Integrity Number values when compared with RPM. No PCR inhibition was detected. For the RNA kits, RPM gave higher concentrations compared with MWTV and no differences were seen in RNA Integrity Number values. PCR inhibition was observed in all RNA samples, with MWTV exhibiting higher inhibition compared with RPM.

Gastrointestinal nematodes in German outdoor-reared pigs based on faecal egg count and next-generation sequencing nemabiome data

BACKGROUND

There is a higher risk for nematode infections associated with outdoor-reared pigs. Next to Ascaris suum, Oesophagostomum dentatum and Trichuris suis, there is the potential of infections with other nodular worm species, Hyostrongylus rubidus, Stongyloides ransomi and Metastrongylus spp. lungworms. Next-generation sequencing methods describing the nemabiome have not yet been established for porcine nematodes.

METHODS

FLOTAC was used for faecal egg counts of porcine gastrointestinal nematodes and lungworms in piglets, fatteners and adults individually. A nemabiome analyses based on ITS-2 gene region metabarcoding was used to differentiate strongyle species. Additionally, questionnaire data was analysed using mixed-effect regression to identify potential risk factors associated with parasite occurrences and egg shedding intensity.

RESULTS

On 15 of 17 farms nematode eggs were detected. Ascaris suum, strongyles and T. suis were detected on 82%, 70% and 35% of the 17 farms, respectively. Lungworms were detected on one out of four farms with access to pasture. Strongyloides ransomi was not detected. 32% (CI 28-36%), 27% (24-31%), 5% (4-7%) and 3% (0.9-8%) of the samples where tested positive for strongyles, A. suum, T. suis and lungworms, respectively. The nemabiome analysis revealed three different strongyle species, with O. dentatum being the most common (mean 93.9%), followed by O. quadrispinulatum (5.9%) and the hookworm Globocephalus urosubulatus (0.1%). The bivariate and multivariate risk factor analyses showed among others that cleaning once a week compared to twice a week increased the odds significantly for being infected with A. suum (OR 78.60) and strongyles (2077.59). Access to pasture was associated with higher odds for A. suum (43.83) and strongyles (14.21). Compared to shallow litter systems, deep litter and free range systems resulted in significant higher odds for strongyles (85.74, 215.59, respectively) and T. suis (200.33, 623.08).

CONCLUSIONS

Infections with A. suum, O. dentatum, O. quadrispinulatum, T. suis, Metastrongylus spp. and G. urosubulatus are present in German outdoor-reared pigs. This is the first report of G. urosubulatus in domestic pigs in Europe. Metabarcoding based on the ITS-2 region is a suitable tool to analyse the porcine nemabiome. Furthermore, management practices have the potential of reducing the risk of parasite infections.

Faecal egg count reduction tests and nemabiome analysis reveal high frequency of multi-resistant parasites on sheep farms in north-east Germany involving multiple strongyle parasite species

Anthelmintic resistance in sheep parasitic gastrointestinal nematodes is widespread and a severe health and economic issue but prevalence of resistance and involved parasite species are unknown in Germany. Here, the faecal egg count reduction test (FECRT) was performed on eight farms using fenbendazole, ivermectin and moxidectin and on four farms using only moxidectin. A questionnaire was used to obtain data on management practices to potentially identify risk factors for presence of resistance. All requirements of the recently revised WAAVP guideline for diagnosing anthelmintic resistance using the FECRT were applied. Nematode species composition in pre- and post-treatment samples was analysed with the nemabiome approach. Using the eggCounts statistic package, resistance against fenbendazole, ivermectin and moxidectin was found on 7/8, 8/8 and 8/12 farms, respectively. No formal risk factor analysis was conducted since resistance was present on most farms. Comparison with the bayescount R package results revealed substantial agreement between methods (Cohen's κ = 0.774). In contrast, interpretation of data comparing revised and original WAAVP guidelines resulted in moderate agreement (Cohen's κ = 0.444). The FECR for moxidectin was significantly higher than for ivermectin and fenbendazole. Nemabiome data identified 4 to 12 species in pre-treatment samples and treatments caused a small but significant decrease in species diversity (inverse Simpson index). Non-metric multidimensional scaling and k-means clustering were used to identify common patterns in pre- and post-treatment samples. However, post-treatment samples were scattered among the pre-treatment samples. Resistant parasite species differed between farms. In conclusion, the revised FECRT guideline allows robust detection of anthelmintic resistance. Resistance was widespread and involved multiple parasite species. Resistance against both drug classes on the same farm was common. Further studies including additional drugs (levamisole, monepantel, closantel) should combine sensitive FECRTs with nemabiome data to comprehensively characterise the anthelmintic susceptibility status of sheep nematodes in Germany.

Transforming gastrointestinal helminth parasite identification in vertebrate hosts with metabarcoding: a systematic review

BACKGROUND

Gastrointestinal helminths are a very widespread group of intestinal parasites that can cause major health issues in their hosts, including severe illness or death. Traditional methods of helminth parasite identification using microscopy are time-consuming and poor in terms of taxonomic resolution, and require skilled observers. DNA metabarcoding has emerged as a powerful alternative for assessing community composition in a variety of sample types over the last few decades. While metabarcoding approaches have been reviewed for use in other research areas, the use of metabarcoding for parasites has only recently become widespread. As such, there is a need to synthesize parasite metabarcoding methodology and highlight the considerations to be taken into account when developing a protocol.

METHODS

We reviewed published literature that utilized DNA metabarcoding to identify gastrointestinal helminth parasites in vertebrate hosts. We extracted information from 62 peer-reviewed papers published between 2014 and 2023 and created a stepwise guide to the metabarcoding process.

RESULTS

We found that studies in our review varied in technique and methodology, such as the sample type utilized, genetic marker regions targeted and bioinformatic databases used. The main limitations of metabarcoding are that parasite abundance data may not be reliably attained from sequence read numbers, metabarcoding data may not be representative of the species present in the host and the cost and bioinformatic expertise required to utilize this method may be prohibitive to some groups.

CONCLUSIONS

Overall, using metabarcoding to assess gastrointestinal parasite communities is preferable to traditional methods, yielding higher taxonomic resolution, higher throughput and increased versatility due to its utility in any geographical location, with a variety of sample types, and with virtually any vertebrate host species. Additionally, metabarcoding has the potential for exciting new discoveries regarding host and parasite evolution.

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.

Surveillance of Ancylostoma caninum in naturally infected dogs in Quebec, Canada, and assessment of benzimidazole anthelmintics reveal a variable efficacy with the presence of a resistant isolate in imported dogs

Ancylostoma caninum is a widely prevalent parasitic nematode in dogs across the world. There has been a notable increase in reports of anthelmintic resistance in A. caninum within the United States of America in recent years, which has led us to investigate the potential of this scenario in Canada. The study objectives were to assess the prevalence of A. caninum in two different groups, including a colony of rescued dogs in Canada and three imported Greyhound dogs from USA, and to evaluate the efficacy of two benzimidazole (BZ) anthelmintics against A. caninum, complemented with a molecular genetic analysis adapted to low prevalence. Fecal samples were collected at pre- and post-treatment with fenbendazole for the native shelters-origin group, and a combination of anthelmintic formulations, including the pro-BZ febantel for the USA-origin group. The coprology analyses found several genera of internal parasites. Canine ancylostomiasis was the most prevalent parasitosis with 30.77% in the native group and 100% in the USA group, but with overall low average of A. caninum eggs per gram. Through the fecal egg count reduction test (FECRT), applying a cut-off at 90% as baseline of egg reduction for successful efficacy, BZ showed variable efficacy. Furthermore, molecular analysis confirmed the presence of A. caninum in both groups of dogs and found differences in the genetics linked to BZ resistance on the A. caninum β-tubulin isotype 1 gene. In the isolate from the native group, both codons 167 and 200 were homozygous without the presence of single nucleotide polymorphism (SNP). In contrast, the selected isolate from the USA group, showed a homozygous allele at position 200 and a heterozygous SNP at position 167. The latter was congruent with the low efficacy in FECRT and agrees with the recent findings of USA A. caninum isolate resistant phenotype to the BZ anthelmintics. The limitations of the study include an overall low eggs-per-gram in both canine groups, and the shortage of additional fecal samples from the USA group, restraining the molecular analysis only to one out of the three Greyhounds. This study provided some insights on the efficacy of BZs against A. caninum and revealed the presence of BZ resistant isolates in imported dogs in Quebec, Canada. All this information should be considered, for choosing the best strategy in the control of A. caninum using anthelmintic drugs.

World Association for the Advancement of Veterinary Parasitology (W.A.A.V.P.): Third edition of the guideline for evaluating efficacy of anthelmintics in ruminants (bovine, ovine, caprine)

This guideline is aimed at those who are involved in the assessment of anthelmintic efficacy in ruminant livestock species (bovine, ovine and caprine). The intent is to provide a framework that can be adopted worldwide for the testing of anthelmintics in ruminants, such that studies carried out in different countries can be compared and thereby unnecessary duplication can be reduced. Recommendations are made for the selection, housing and feeding of study animals, the type of studies required, the method used to conduct those studies, the assessment of results and the standards for defining anthelmintic efficacy.

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.

Next-generation sequencing technologies for helminth diagnostics and surveillance in ruminants: shifting diagnostic barriers

Helminth infections in grazing ruminants are a major issue for livestock farming globally, but are unavoidable in outdoor grazing systems and must be effectively managed to avoid deleterious effects to animal health, and productivity. Next-generation sequencing (NGS) technologies are transforming our understanding of the genetic basis of anthelmintic resistance (AR) and epidemiological studies of ruminant gastrointestinal parasites. They also have the potential to not only help develop and validate molecular diagnostic tests but to be directly used in routine diagnostics integrating species-specific identification and AR into a single test. Here, we review how these developments have opened the pathway for the development of multi-AR and multispecies identification in a single test, with widespread implications for sustainable livestock farming for the future.

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.

Occurrence of gastrointestinal nematodes in lambs in Norway, as assessed by copromicroscopy and droplet digital polymerase chain reaction

BACKGROUND

Gastrointestinal nematodes (GINs) have a major impact on sheep production, health, and welfare worldwide. Norway is no exception, but there are only a few studies on the prevalence of GINs in Norwegian sheep. The aim of this study was to investigate the current occurrence of the most important nematodes in sheep flocks in Norway. Faecal samples were collected from flocks in 2021/2022, mainly from three geographical regions in Norway, i.e., northern, eastern, and western. In each of 134 flocks included, individual samples from 10 lambs (autumn) were pooled. Third stage larvae (L3) were cultivated and harvested (Baermann method) from the pooled samples. The DNA was then extracted and further analysed using droplet digital PCR (ddPCR). This enables assessment of the proportions of the three most important nematode species/genera, i.e., H. contortus, T. circumcincta, and Trichostrongylus. The fractional abundance/relative proportion of each species/genus was assessed by performing duplex assays with universal strongyle and species/genus-specific primers and probe sets. In addition, the occurrence of Nematodirus eggs was assessed by standard faecal egg counts (i.e., McMaster method).

RESULTS

Of the 134 flocks sampled, 24 were from the northern region, 31 from eastern, and 71 from western Norway. In addition, some flocks from central (n = 7), and southern (n = 1) Norway were included. Among the sampled flocks, T. circumcincta occurred most commonly (94%), followed by H. contortus (60%) and Trichostrongylus (55%), and Nematodirus (51%). In general, mixed infections were observed, with 38% and 18% of flocks infected with three or all four genera, respectively.

CONCLUSIONS

The results of this study indicate that GINs are widespread in Norway. Teladorsagia circumcincta seems to be present in most flocks based on this screening. Moreover, the results show that Nematodirus spp. infect lambs throughout the country, predominantly N. battus, and indicate that this nematode has become more abundant, which could lead to an increase in nematodirosis.

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.

New insight into genetic diversity of zoonotic-potential Ancylostoma ceylanicum in stray cats living in Bangkok, Thailand, based on deep amplicon sequencing

AIMS

This study aimed to characterize feline hookworms from stray cats living in Bangkok.

METHODS AND RESULTS

A total of 56 hookworm-positive faecal samples were identified for hookworm species by using PCR targeting the ITS1, 5.8S, and ITS2 fragment and qPCR targeting ITS2. Of 56 samples, 96.4% (54/56) were identified as Ancylostoma ceylanicum and 1.8% (1/56) as Ancylostoma caninum. With qPCR, 89.3% (50/56) were identified as single A. ceylanicum infection and 5.4% (3/56) as coinfection of A. ceylanicum and A. caninum. For genetic characterization of A. ceylanicum, 10 samples were pooled, and the partial COI gene was amplified, followed by deep amplicon sequencing. Five pooled samples were analysed, and 99.73% were identified with A. ceylanicum sequences, which were allocated into 19 haplotypes (AC01-AC19). Genetic diversity findings for A. ceylanicum in Asia revealed that three of eight haplotypes considered of zoonotic significance occurred in humans, dogs, and cats, including haplotypes H01, H20, and H21. The predominant haplotype in this study, AC01, was clustered with H01-a zoonotic haplotype.

CONCLUSIONS

The diversity obtained by deep amplicon sequencing supported that the A. ceylanicum community had high genetic variation. Deep amplicon sequencing was a useful method to determine source, zoonotic potential, and host-parasite relationship.

The paradox of spring: Thyroid and glucocorticoid responses to cold temperatures and food availability in free living Carneddau ponies

In seasonal environments, maintaining a constant body temperature poses challenges for endotherms. Cold winters at high latitudes, with limited food availability, create opposing demands on metabolism: upregulation preserves body temperature but depletes energy reserves. Examining endocrine profiles, such as thyroid hormone triiodothyronine (T3) and glucocorticoids (GCs), proxies for changes in metabolic rate and acute stressors, offer insights into physiological trade-offs. We evaluated how environmental conditions and gestation impact on faecal hormone metabolites (fT3Ms and fGCMs) from late winter to spring in a free-living population of Carneddau ponies. Faecal T3Ms were highest in late February and March, when temperatures were lowest. Then, fT3Ms concentrations decreased throughout April and were at the lowest in May before increasing towards the end of the study. The decline in fT3M levels in April and May was associated with warmer weather but poor food availability, diet diversity and diet composition. On the other hand, fGCM levels did not display a clear temporal pattern but were associated with reproductive status, where pregnant and lactating females had higher fGCM levels as compared to adult males and non-reproductive females. The temporal profile of fT3Ms levels highlights metabolic trade-offs in a changing environment. In contrast, the ephemeral but synchronous increase in fGCM concentrations across the population suggest a shared experience of acute stressors (i.e., weather, disturbance or social). This multi-biomarker approach can evaluate the role of acute stressors versus energy budgets in the context of interventions, reproduction, seasonality and environmental change, or across multiple scales from individuals to populations.

Validation of deep amplicon sequencing of Dicrocoelium in small ruminants from Northern regions of Pakistan

Dicrocoelium lancet flukes cause significant production loss in ruminant livestock. Although co-infection with multiple Dicrocoelium species within a host is common, techniques for studying the composition of these complex parasite communities are lacking. The pathogenicity, epidemiology, and therapeutic susceptibility of different helminth species vary, and little is known about the interactions that take place between co-infecting species and their hosts. Here, we describe the first applicationof metabarcoding deep amplicon sequencing method to studythe Dicrocoelium species in sheep and goats. First, rDNA ITS-2 sequences of four Dicrocoelium species (Dicrocoelium dendriticum, Dicrocoelium hospes, Dicrocoelium orientalis, and Dicrocoelium chinensis) were extracted from the NCBI public database. Phylogenetic analysis revealed separate clades of Dicrocoelium species; hence, molecular differentiation between each species is possible in co-infections. Second, 202 flukes belonging to seventeen host populations (morphologically verified as belonging to the Dicrocoelium genus) were evaluated to determine the deep amplicon sequencing read threshold of an individual fluke for each of the four species. The accuracy of the method in proportional quantification of samples collected from single hosts was further assessed. Overall, 198 (98.01%) flukes were confirmed as D. dendriticum and 1.98% produced no reads. The comparison of genetic distances between rDNA ITS-2 revealed 86% to 98% identity between the Dicrocoelium species. Phylogenetic analysis demonstrated a distinct clustering of species, apart from D. orientalis and D. chinensis, which sit very close to each other in a single large clade whereas D. hospes and D. dendriticum are separated into their own clade. In conclusion each sample was identified as D. dendriticum based on the proportion of MiSeq reads and validated the presence of this group of parasites in the Gilgit Baltistan and Khyber Pakhtunkhwa provinces of Pakistan. The metabarcoding deep amplicon sequencing technology and bioinformatics pathway have several potential applications, including species interactions during co-infections, identifying the host and geographical distribution of Dicrocoelium in livestock, drug therapy response evaluation and understanding of the emergence and spread of drug resistance.

Faecal egg count reduction test in goats: Zooming in on the genus level

The faecal egg count reduction test (FECRT) is the most widely used method to assess treatment efficacy against gastrointestinal nematodes (GIN). Information on genera composition of the GIN community is not available with this test and it is commonly obtained by identifying cultured third-stage larvae (L3) or through molecular assays in the post-treatment survey, but results provided are usually only qualitative or semi-quantitative. The updated WAAVP guidelines now recommend assessing anthelmintic efficacy for each GIN genus/species separately (genus-specific FECRT), but this approach is poorly employed in Europe and in goats especially. For this reason, four FECRT trials were conducted using oxfendazole and eprinomectin in two Italian goat farms. Samples were processed individually using the McMaster technique and then pooled to create two samples from faeces of 5 animals each. Pooled samples were analysed using the McMaster and cultured for seven days at 26°C to obtain L3s. The genus-specific FECRT was based on larval identification, integrating coproculture and FEC results. Larvae were identified as Haemonchus, Trichostrongylus, Teladorsagia, Oesophagostomum / Chabertia and Bunostomum. Molecular assays (a multiplex real-time PCR and two end-point PCRs) were also implemented on pooled samples to support the morphological identification. The Spearmann Rho test confirmed a high correlation between the two approaches (Rho = 0.941 and Rho = 0.914 respectively for Haemonchus and Trichostrongylus, the two most common genera). Both oxfendazole and eprinomectin were effective in one farm, while none in the other farm (FECR = 75.9% and 73.3% respectively). In the second farm, the genus-specific FECRT highlighted a different response to treatment among genera: oxfendazole lacked efficacy against both Haemonchus and Trichostrongylus spp., eprinomectin only against Haemonchus, while all other genera were susceptible to both drugs. This study brings new attention on the importance of adopting a genus-specific approach to identify and quantify differences in susceptibility to anthelmintics among genera in goats, providing support for FECRT interpretation, anthelmintic resistance evaluation and evidence-based GIN control.

Investigating the perceived versus actual gastrointestinal nematode challenge on extensive sheep farms

Extensive farming systems form an integral part of sheep production systems across Europe. However, with innate production handicaps, declining sheep numbers and narrow economic margins, production is becoming increasingly challenging threatening the future sustainability of the industry. Gastrointestinal nematodes (GINs) are a significant cause of production losses to the global sheep industry, with well-established resistance to the major anthelmintic groups. Traditionally, extensive farming systems are not thought to have a significant parasite challenge compared with intensive farms, but there is a need to identify the scale and importance of GINs on extensive farms to inform the need for sustainable control strategies. In this study, a questionnaire of extensive farmers (n=34) was conducted and parasitological data were collected from nine study farms to investigate the perceived versus actual GIN and anthelmintic resistance challenge faced by extensive farms. The results showed a production-limiting challenge on most farms, with a higher GIN challenge observed on improved pastures. Furthermore, over half of the extensive farmers perceived anthelmintic resistance to be a greater problem for intensive farmers, with only 20% of respondents reporting known anthelmintic resistance. However, all study farms had evidence of resistance to at least one group of anthelmintics. Consequently, this study has demonstrated that despite the traditional perception of parasitism on extensive farms, there is a need to increasingly consider its impact and take a proactive approach to sustainable control, with solutions tailored to their unique management.

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 combined amplicon approach to nematode polyparasitism occurring in captive wild animals in southern China

BACKGROUND

Gastrointestinal tract (GIT) nematodes prefer to live in the intestines of wild animals, causing damage and even death, and posing a zoonotic risk. The polyparasitism of GIT nematodes results in the complex dynamics of the nematode communities that occur naturally in wild animals. However, the nematode community in captive wild animals is poorly understood.

METHODS

We combined  microscopic examination and amplicon sequencing for community diversity.

RESULTS

We characterized GIT nematode assemblages to one order, one family, four genera, and ten species, in 512 fecal samples of 121 species from captive wild animals in southern China. The positive rate of GIT nematodes was 20.7% (106/512), including 42.3% (11/26) in reptiles, 26.5% (39/147) in herbivores, 25.0% (25/100) in non-human primates, 20.0% (5/25) in omnivores, 12.2% (9/74) in carnivores, and 12.1% (17/140) in avians. The dominant nematodes were Haemonchus contortus in herbivores and Trichuris species in primates. The nematode communities of arboreal primates differed from their terrestrial counterparts, reflecting both host phylogeny and ecological constraints. Soil-transmitted Strongyloides species were widespread throughout the herbivore, primate, avian, and carnivore communities, and tended to infect omnivorous primates and terrestrial herbivores. In addition, new Trichuris and Heterakis species were found in the nematode communities of captive porcupines and peafowls.

CONCLUSION

This study highlights the variation in the composition of the GIT nematode community and strengthens the attention to the harms induced by zoonotic nematodes and co-infective nematodes with low species richness.

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.

Assessing morphological preservation of gastrointestinal parasites from fecal samples of wild capuchin monkeys (Cebus imitator) stored in ethanol versus formalin

The copromicroscopic identification of gastrointestinal parasites is a common, cost-effective method vital to understanding host-parasite interactions. However, its efficacy depends on effective preservation of the samples. In this study, we compare the preservation of ethanol and formalin preserved gastrointestinal parasites collected from a wild population of Costa Rican capuchin monkeys (Cebus imitator). Fecal samples were collected, halved, and stored in either 10% formalin or 96% ethanol at ambient temperature, then microscopically screened for the presence of parasites. Parasites were morphologically identified and rated based on their preservation using a newly developed rubric. We identified more parasitic morphotypes in formalin-preserved samples but found no difference in the number of parasites per fecal gram (PFG) between mediums. There was no difference in the PFG of two most prevalent parasite morphotypes, Filariopsis barretoi larvae and Strongyle-type eggs, and while Filariopsis larvae were better preserved in formalin, strongyle eggs showed no preservation difference between mediums. Our results support the suitability of both ethanol and formalin for morphological parasite identification in samples stored over 1 year, describe the morphological changes and challenges associated with parasite degradation, and highlight the potential for future studies to use both morphological and molecular methods in non-invasively collected samples.

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.

Real-time single-base specific detection of the Haemonchus contortus S168T variant associated with levamisole resistance using loop-primer endonuclease cleavage loop-mediated isothermal amplification

Haemonchus contortus is a parasitic haematophagous nematode that primarily affects small ruminants and causes significant economic loss to the global livestock industry. Treatment of haemonchosis typically relies on broad-spectrum anthelmintics, resistance to which is an important cause of treatment failure. Resistance to levamisole remains less widespread than to other major anthelmintic classes, prompting the need for more effective and accurate surveillance to maintain its efficacy. Loop-primer endonuclease cleavage loop-mediated isothermal amplification (LEC-LAMP) is a recently developed diagnostic method that facilitates multiplex target detection with single nucleotide polymorphism (SNP) specificity and portable onsite testing. In this study, we designed a new LEC-LAMP assay and applied it to detect the levamisole resistance marker S168T in H. contortus. We explored multiplexing probes for both the resistant S168T and the susceptible S168 alleles in a single-tube assay. We then included a generic probe to detect the acr-8 gene in the multiplex assay, which could facilitate the quantification of both resistance markers and overall genetic material from H. contortus in a single step. Our results showed promising application of these technologies, demonstrating a proof-of-concept assay which is amenable to detection of resistance alleles within the parasite population, with the potential for multiplex detection, and point-of-care application enabled by lateral flow end-point detection. However, further optimisation and validation is necessary.

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.

Development and validation of a long-read metabarcoding platform for the detection of filarial worm pathogens of animals and humans

BACKGROUND

Filarial worms are important vector-borne pathogens of a large range of animal hosts, including humans, and are responsible for numerous debilitating neglected tropical diseases such as, lymphatic filariasis caused by Wuchereria bancrofti and Brugia spp., as well as loiasis caused by Loa loa. Moreover, some emerging or difficult-to-eliminate filarioid pathogens are zoonotic using animals like canines as reservoir hosts, for example Dirofilaria sp. 'hongkongensis'. Diagnosis of filariasis through commonly available methods, like microscopy, can be challenging as microfilaremia may wane below the limit of detection. In contrast, conventional PCR methods are more sensitive and specific but may show limited ability to detect coinfections as well as emerging and/or novel pathogens. Use of deep-sequencing technologies obviate these challenges, providing sensitive detection of entire parasite communities, whilst also being better suited for the characterisation of rare or novel pathogens. Therefore, we developed a novel long-read metabarcoding assay for deep-sequencing the filarial nematode cytochrome c oxidase subunit I gene on Oxford Nanopore Technologies' (ONT) MinION™ sequencer. We assessed the overall performance of our assay using kappa statistics to compare it to commonly used diagnostic methods for filarial worm detection, such as conventional PCR (cPCR) with Sanger sequencing and the microscopy-based modified Knott's test (MKT).

RESULTS

We confirmed our metabarcoding assay can characterise filarial parasites from a diverse range of genera, including, Breinlia, Brugia, Cercopithifilaria, Dipetalonema, Dirofilaria, Onchocerca, Setaria, Stephanofilaria and Wuchereria. We demonstrated proof-of-concept for this assay by using blood samples from Sri Lankan dogs, whereby we identified infections with the filarioids Acanthocheilonema reconditum, Brugia sp. Sri Lanka genotype and zoonotic Dirofilaria sp. 'hongkongensis'. When compared to traditionally used diagnostics, such as the MKT and cPCR with Sanger sequencing, we identified an additional filarioid species and over 15% more mono- and coinfections.

CONCLUSIONS

Our developed metabarcoding assay may show broad applicability for the metabarcoding and diagnosis of the full spectrum of filarioids from a wide range of animal hosts, including mammals and vectors, whilst the utilisation of ONT' small and portable MinION™ means that such methods could be deployed for field use.

VESPA: an optimized protocol for accurate metabarcoding-based characterization of vertebrate eukaryotic endosymbiont and parasite assemblages

Protocols for characterizing taxonomic assemblages by deep sequencing of short DNA barcode regions (metabarcoding) have revolutionized our understanding of microbial communities and are standardized for bacteria, archaea, and fungi. Unfortunately, comparable methods for host-associated eukaryotes have lagged due to technical challenges. Despite 54 published studies, issues remain with primer complementarity, off-target amplification, and lack of external validation. Here, we present VESPA (Vertebrate Eukaryotic endoSymbiont and Parasite Analysis) primers and optimized metabarcoding protocol for host-associated eukaryotic community analysis. Using in silico prediction, panel PCR, engineered mock community standards, and clinical samples, we demonstrate VESPA to be more effective at resolving host-associated eukaryotic assemblages than previously published methods and to minimize off-target amplification. When applied to human and non-human primate samples, VESPA enables reconstruction of host-associated eukaryotic endosymbiont communities more accurately and at finer taxonomic resolution than microscopy. VESPA has the potential to advance basic and translational science on vertebrate eukaryotic endosymbiont communities, similar to achievements made for bacterial, archaeal, and fungal microbiomes.

DNA metabarcoding reveals spatial and temporal variation of fish eye fluke communities in lake ecosystems

Eye flukes (Diplostomidae) are diverse and abundant trematode parasites that form multi-species communities in fish with negative effects on host fitness and survival. However, the environmental factors and host-related characteristics that determine species diversity, composition, and coexistence in such communities remain poorly understood. Here, we developed a cost-effective cox1 region-specific DNA metabarcoding approach to characterize parasitic diplostomid communities in two common fish species (Eurasian perch and common roach) collected from seven temperate lakes in Estonia. We found considerable inter- and intra-lake, as well as inter-host species, variation in diplostomid communities. Sympatric host species characterization revealed that parasite communities were typically more diverse in roach than perch. Additionally, we detected five positive and two negative diplostomid species associations in roach, whereas only a single negative association was observed in perch. These results indicate that diplostomid communities in temperate lakes are complex and dynamic systems exhibiting both spatial and temporal heterogeneity. They are influenced by various environmental factors and by host-parasite and inter-parasite interactions. We expect that the described methodology facilitates ecological and biodiversity research of diplostomid parasites. It is also adaptable to other parasite groups where it could serve to improve current understanding of diversity, distribution, and interspecies interactions of other understudied taxa.

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.

First insight into strongylid nematode diversity and anthelmintic treatment effectiveness in beef cattle in the Czech Republic explored by HTS metagenomics

Parasitic diseases and mitigation of their effects play an important role in the health management of grazing livestock worldwide, with gastrointestinal strongylid nematodes being of prominent importance. These helminths typically occur in complex communities, often composed of species from numerous strongylid genera. Detecting the full diversity of strongylid species in non-invasively collected faecal samples is nearly impossible using conventional methods. In contrast, high-throughput amplicon sequencing (HTS) can effectively identify co-occurring species. During the four-year project, we collected and analysed faecal samples from beef cattle on >120 farms throughout the Czech Republic. Strongylids were the predominant nematodes, detected in 56% of the samples, but at a low level of infection. The apparent limitations in identifying strongylid taxa prompted this pilot study on a representative group of samples testing positive for strongylids using ITS-2 metabarcoding. The most widespread genera parasitizing Czech cattle were Ostertagia (O. ostertagi) and Oesophagostomum spp., followed by Trichostrongylus and Cooperia, while Bunostomum, Nematodirus and Chabertia were present only in a minority. As comparative material, 21 samples of cattle from the Danube Delta in Romania were used, which, in contrast, were dominated by Haemonchus placei. Finally, the effect of ivermectin treatment was tested at two Czech farms. After treatment with the anthelmintic, there was a shift in the strongylid communities, with a dominance of Cooperia and Ostertagia.

Identification of Trichinella taxa by ITS-1 amplicon next-generation sequencing with an improved resolution for detecting underrepresented genotypes in mixed natural infections

BACKGROUND

Amplicon-based next-generation sequencing (NGS) has rapidly gained popularity as a powerful method for delineating taxa in complex communities, including helminths. Here, we applied this approach to identify species and genotypes of zoonotic nematodes of the Trichinella genus. A known limitation of the current multiplex PCR (mPCR) assay recommended by the International Commission on Trichinellosis is that it does not differentiate Trichinella nativa from T. chanchalensis.

METHODS

The new assay entails deep sequencing of an amplified variable fragment of the ribosomal cistron's (rDNA) internal transcribed spacer 1 using the Illumina platform. The assay was evaluated using first-stage larvae (L1) of select laboratory strains of various Trichinella taxa mixed in known proportions and then validated using archived L1 from 109 wildlife hosts. The species/genotypes of these L1 isolates from wildlife were previously determined using mPCR.

RESULTS

NGS data analysis for Trichinella laboratory strains selected as representative of North American fauna revealed a sequence representation bias. Trichinella pseudospiralis, a non-encapsulated species, was the most underrepresented when mixed with T. spiralis, T. murrelli, T. nativa and Trichinella T6 in equal quantities. However, five L1 of T. pseudospiralis were readily revealed by NGS in a mix with 2000 L1 of T. nativa (1:400 ratio). From naturally infected wildlife, all Trichinella taxa revealed by mPCR were also identified by NGS in 103 of 107 (96.3%) samples amplified on both assays. NGS identified additional taxa in 11 (10.3%) samples, whereas additional taxa were revealed by mPCR in only four (3.7%) samples. Most isolates comprised single or mixed infections of T. nativa and Trichinella T6. On NGS, T. chanchalensis (T13) was detected in combination with Trichinella T6 in a wolverine (Gulo gulo) and in combination with T. nativa and Trichinella T6 in a marten (Martes americana) from the Northwest Territories, Canada.

CONCLUSIONS

This new NGS assay demonstrates strong potential as a single assay for identifying all recognised Trichinella taxa as well as improved sensitivity for detecting under-represented and novel genotypes in mixed infections. In addition, we report a new host record for T. chanchalensis in American marten.

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.

Understanding anthelmintic resistance in livestock using "omics" approaches

Widespread and improper use of various anthelmintics, genetic, and epidemiological factors has resulted in anthelmintic-resistant (AR) helminth populations in livestock. This is currently quite common globally in different livestock animals including sheep, goats, and cattle to gastrointestinal nematode (GIN) infections. Therefore, the mechanisms underlying AR in parasitic worm species have been the subject of ample research to tackle this challenge. Current and emerging technologies in the disciplines of genomics, transcriptomics, metabolomics, and proteomics in livestock species have advanced the understanding of the intricate molecular AR mechanisms in many major parasites. The technologies have improved the identification of possible biomarkers of resistant parasites, the ability to find actual causative genes, regulatory networks, and pathways of parasites governing the AR development including the dynamics of helminth infection and host-parasite infections. In this review, various "omics"-driven technologies including genome scan, candidate gene, quantitative trait loci, transcriptomic, proteomic, and metabolomic approaches have been described to understand AR of parasites of veterinary importance. Also, challenges and future prospects of these "omics" approaches are also discussed.

Comparison of diagnostic methods for assessment of Ostertagia ostertagi exposure in Norwegian dairy herds

BACKGROUND

The gastrointestinal nematode (GIN) Ostertagia ostertagi can cause severe disease in first season grazers (FSG) and impaired performance due to subclinical infections in adult cows. Diagnostic methods to assess exposure include faecal egg count and detection of specific antibodies using antibody-ELISAs resulting in an optical density ratio (ODR). Using the ELISA test on bulk tank milk (BTM) allows for a herd level diagnosis. Appropriate use of diagnostic methods for evaluation of O. ostertagi exposure is required to optimize herd parasite surveillance and aid in a sustainable control regime. The aim of this study was to describe the relationship between different diagnostic tests used to assess GIN exposure in Norwegian production systems. A cross-sectional field study was carried out in twenty herds in Norway in the fall of 2020. Serum and faecal samples were taken from 380 individuals, of which 181 were FSG and 199 were cows. In addition, milk was collected from every cow and one BTM sample was taken from each herd. Faecal egg counts were performed. The distribution of ODR values in individual samples within and between herds and the associations between BTM ODR and individual ODR values were described. The data were analysed using visual assessment of scatter plots, Pearson correlation coefficients and linear regression.

RESULTS

A high variability of the within-herd individual ODR values in serum and milk in every herd was detected. The ODR in BTM explained a low degree of the variation in the individual serum and milk samples. When plotting the ODR results in milk or serum according to four BTM categories, the distribution of ODR values were notably different in the highest and lowest BTM categories. The correlation between individual milk and serum samples was moderate (r = 0.68), while the highest correlation (r = 0.81) was between the BTM ODR and the group average individual milk samples.

CONCLUSIONS

A poor predictive ability for BTM ODR to assess individual ODR values in both FSG and cows was demonstrated. However, the study indicates that the evaluation by ELISA test on BTM to assess exposure to GIN could be useful in herds with a very high or low BTM ODR.

Antagonism between co-infecting gastrointestinal nematodes: A meta-analysis of experimental infections in Sheep

Gastrointestinal nematodes (GIN) have enormous global impacts in humans, wildlife and grazing livestock. Within grazing livestock, sheep are of particular global importance and the economics and sustainability of sheep production are greatly constrained by GIN infections. Natural infections are composed of co-infections with multiple species, and while some past work suggests species can interact negatively with one another within the same host, there is wide variation in reported patterns. Here, we undertook a systematic literature search and meta-analysis of experimental GIN co-infections of sheep to determine whether these experimental studies support the hypothesis of antagonistic interactions between different co-infecting GIN, and test whether aspects of parasite biology or experimental design influence the observed effects. A systematic search of the literature yielded 4848 studies, within which, we identified 19 experimental sheep studies comparing post-mortem worm counts across two co-infecting GIN species. Meta-analysis of 67 effects obtained from these studies provides strong evidence for interactions between GIN species. There was wide variation in the strength and direction of these interactions, but the global effect was significantly antagonistic. On average, there was a decrease in the number of worms of one species when a co-infecting species was also present, relative to a mono-infection with that species alone. This effect was dependent on the infectious dose and was rapidly lost after anthelmintic treatment, suggesting that live worms are required for the effect to occur. Individual parasite species varied in the extent to which they both exerted, and were subject to, these interspecies interactions, and these differences are more complex than simply co-localisation within the gastrointestinal tract. Antagonistic interactions between co-infecting GIN may feedback into their epidemiology as well as potentially affecting the clinical impacts of infection. Furthermore, the consequences of these interactions may be heightened when clinical interventions affect only one species within the co-infecting network. Whilst it was not possible to identify the causes of variation between GIN species in the impact of co-infection, these findings point to new avenues for epidemiological, clinical and mechanistic research on GIN co-infections.

Molecular identification and phylogenetic analysis of gastrointestinal nematodes in different populations of Kazakh sheep

Gastrointestinal nematode (GIN) infection in sheep has been recognized globally as a major problem challenging animal health and production. The objective of this study is to use a molecular diagnosis of the prevalence for gastrointestinal nematode (GIN) dominant species of Kazakh sheep and its hybrid (Kazakh × Texel). The internal transcribed spacer 2 (ITS-2) sequences of ribosomal DNA (rDNA) were used as the target sequence. In the study, three dominant species of nematodes, namely Haemonchus contortus, Trichostrongylus spp., and Teladorsagia (Ostertagia) circumcincta from the Kazakh sheep and the F1 and F2 generations of Texel × Kazakh sheep hybrids were subjected to molecular identification and phylogenetic analysis. The fecal and single larva genomic DNA were extracted and amplified by PCR using specific primers to determine the infection rate of the three nematode species. In addition, the PCR products were sequenced and analyzed using bioinformatics methods to construct a phylogenetic tree. The results showed that all the three species had their ITS-2 specific amplified. According to the sequence homology analysis of PCR products, the results showed a high homology (above 98.5% homology) with H. contortus, Trichostrongylus spp., T. circumcincta ITS-2 sequences in GenBank. Phylogenetic analysis showed that the ITS-2 sequences of the three species were on the same branch as the ITS-2 sequences of the same species in NCBI. And on different branches from those of the ITS-2 sequences of different families, genera and species. Sequences carried out on three species from different samples showed a close relationship and little genetic difference in phylogenetic tree. The infection rates based on fecal DNA were 35.59, 25.55, and 11.24% for H. contortus, Trichostrongylus spp., and T. circumcincta, respectively. While the infection rates based on larva DNA, were 24.07, 18.89, and 13.26% for H. contortus, Trichostrongylus spp., and T. circumcincta, respectively. The seasonal prevalence of the three dominant species in spring was significantly higher than that in autumn and winter. And there was no significant difference between Kazakh, F1 and F2 sheep considering the infection rate of the studied three species of nematodes. This study provides valuable molecular approaches for epidemiological surveillance and for assisting in the control of Nematodirus infection in sheep.

Benzimidazole-resistance associated mutation in Haemonchus contortus in Norwegian sheep, as detected by droplet digital PCR

The aim of this study was to investigate the occurrence of benzimidazole-resistant Haemonchus contortus in Norwegian sheep flocks. Screening was based on detection of one of the resistance-conferring mutations in the β tubulin isotype 1 gene (F200Y, TAC) in larvae (L3) cultivated from H. contortus eggs from naturally infected sheep. Faecal samples were collected in 2021/2022 from flocks in the northern (n = 34), central (n = 7), eastern (n = 40), southern (n = 1), and western (n = 87) areas of Norway. In total, samples were taken from 169 flocks (spring-ewes samples: 167, autumn-lambs samples: 134). Individual faecal samples were collected from 10 randomly selected ewes (spring) and 10 randomly selected lambs (autumn) in each flock. Faecal samples collected from each flock on each occasion were pooled (lamb and ewe samples pooled separately) and cultured for L3 development. After harvest of larvae (Baermann method), DNA was extracted and then analysed using droplet digital PCR with primer/probe sets targeting the BZ-associated F200Y (TAC) mutation. Haemonchus was found in 60% (80/134) of samples from lambs, and in 63% (106/167) from ewes. Among these, the F200Y mutation was detected in 73% (58/80) of larval samples from lambs and 69% (73/106) of larval samples from ewes, respectively. Although regional differences were evident, the mutation was detected in all areas indicating a widespread distribution and a strong potential for an increasing problem with treatment-resistant haemonchosis in Norwegian sheep flocks.

Deep Sequencing of Porcine Reproductive and Respiratory Syndrome Virus ORF7: A Promising Tool for Diagnostics and Epidemiologic Surveillance

Porcine reproductive and respiratory syndrome virus (PRRSV) is a major concern worldwide. Control of PRRSV is a challenging task due to various factors, including the viral diversity and variability. In this study, we evaluated an amplicon library preparation protocol targeting the ORF7 region of both PRRSV species, Betaarterivirus suid 1 and Betaarterivirus suid 2. We designed tailed primers for a two-step PCR procedure that generates ORF7-specific amplicon libraries suitable for use on Illumina sequencers. We tested the method with serum samples containing common laboratory strains and with pooled serum samples (n = 15) collected from different pig farms during 2019-2021 in Hungary. Testing spiked serum samples showed that the newly designed method is highly sensitive and detects the viral RNA even at low copy numbers (corresponding to approx. Ct 35). The ORF7 sequences were easily assembled even from clinical samples. Two different sequence variants were identified in five samples, and the Porcilis MLV vaccine strain was identified as the minor variant in four samples. An in-depth analysis of the deep sequencing results revealed numerous polymorphic sites along the ORF7 gene in a total of eight samples, and some sites (positions 12, 165, 219, 225, 315, 345, and 351) were found to be common in several clinical specimens. We conclude that amplicon deep sequencing of a highly conserved region of the PRRSV genome could support both laboratory diagnosis and epidemiologic surveillance of the disease.

Unambiguous identification of Ancylostoma caninum and Uncinaria stenocephala in Australian and New Zealand dogs from faecal samples

Hookworms (Ancylostomatidae) are well-known parasites in dogs due to their health impacts and zoonotic potential. While faecal analysis is the traditional method for detection, improvements in husbandry and deworming have decreased their prevalence in urban owned dogs. Drug resistance in Ancylostoma caninum is becoming a discussion point in small animal practices across the region. This study aimed to identify hookworm species present in Australian and New Zealand dogs using molecular techniques. The ITS-2 and isotype-1 β-tubulin assays were used to identify and quantify hookworm species. Results showed absence of coinfection in Australian samples from Greater Sydney region belonging either to A. caninum or Uncinaria stenocephala, while New Zealand samples were a mixture of A. caninum and U. stenocephala. The amplified isotype-1 β-tubulin sequences exhibited susceptibility to benzimidazole drugs. Rare mutations were identified in A. caninum and U. stenocephala sequences, representing a small percentage of reads. This study highlights the importance of molecular techniques in accurately identifying and quantifying hookworm species in dog populations.

The complete ITS2 barcoding region for Strongylus vulgaris and Strongylus edentatus

Gastrointestinal nematode parasites are of major concern for horses, where Strongylus vulgaris is considered the most pathogenic among the Strongylus species. Diagnosis of S. vulgaris infections can be determined with next generation sequencing techniques, which are inherently dependent on reference sequences. The best marker for parasitic nematodes is internal transcribed spacer 2 (ITS2) and we provide the first complete ITS2 sequences from five morphologically identified S. vulgaris and additional sequences from two S. edentatus. These sequences have high similarity to already published partial sequences and amplicon sequence variants (ASV) based on next generation sequencing (NGS). The ITS2 sequences from S. vulgaris matched available partial ITS2 sequences and the full ASVs, whereas the S. edentatus sequence matched another complete sequence. We also compare Sanger sequencing and NGS methods and conclude that the ITS2 variation is better represented with NGS methods. Based on this, we recommend that further sequencing of morphologically identified specimens of various species should be performed with NGS cover the intraspecific variation in the ITS2.

A mixed-model approach for estimating drivers of microbiota community composition and differential taxonomic abundance

Next-generation sequencing (NGS) and metabarcoding approaches are increasingly applied to wild animal populations, but there is a disconnect between the widely applied generalized linear mixed model (GLMM) approaches commonly used to study phenotypic variation and the statistical toolkit from community ecology typically applied to metabarcoding data. Here, we describe the suitability of a novel GLMM-based approach for analyzing the taxon-specific sequence read counts derived from standard metabarcoding data. This approach allows decomposition of the contribution of different drivers to variation in community composition (e.g., age, season, individual) via interaction terms in the model random-effects structure. We provide guidance to implementing this approach and show how these models can identify how responsible specific taxonomic groups are for the effects attributed to different drivers. We applied this approach to two cross-sectional data sets from the Soay sheep population of St. Kilda. GLMMs showed agreement with dissimilarity-based approaches highlighting the substantial contribution of age and minimal contribution of season to microbiota community compositions, and simultaneously estimated the contribution of other technical and biological factors. We further used model predictions to show that age effects were principally due to increases in taxa of the phylum Bacteroidetes and declines in taxa of the phylum Firmicutes. This approach offers a powerful means for understanding the influence of drivers of community structure derived from metabarcoding data. We discuss how our approach could be readily adapted to allow researchers to estimate contributions of additional factors such as host or microbe phylogeny to answer emerging questions surrounding the ecological and evolutionary roles of within-host communities. IMPORTANCE NGS and fecal metabarcoding methods have provided powerful opportunities to study the wild gut microbiome. A wealth of data is, therefore, amassing across wild systems, generating the need for analytical approaches that can appropriately investigate simultaneous factors at the host and environmental scale that determine the composition of these communities. Here, we describe a generalized linear mixed-effects model (GLMM) approach to analyze read count data from metabarcoding of the gut microbiota, allowing us to quantify the contributions of multiple host and environmental factors to within-host community structure. Our approach provides outputs that are familiar to a majority of field ecologists and can be run using any standard mixed-effects modeling packages. We illustrate this approach using two metabarcoding data sets from the Soay sheep population of St. Kilda investigating age and season effects as worked examples.

Diversity of fecal parasitomes of wild carnivores inhabiting Korea, including zoonotic parasites and parasites of their prey animals, as revealed by 18S rRNA gene sequencing

Consisting of diverse groups of organisms, parasites are among the least studied pathogens despite their enormous impacts on humans, livestock, and wildlife. In particular, little is known about their host specificity and diversity in wildlife. Here, using multiple primer pairs and sequencing 18S rRNA genes of diverse groups of parasites, we aimed to investigate fecal parasitomes of carnivorous wildlife in Korea, namely, the raccoon dog (Nyctereutes procyonoides), the leopard cat (Prionailurus bengalensis), and the Eurasian otter (Lutra lutra). A total of 5 host-specific parasite species were identified, including 2 from raccoon dogs, 2 from leopard cats, and 1 from Eurasian otters. In addition, numerous parasite species of their prey animals were detected in their feces. It was found that the parasitome composition varied between host animals, and it was thought that the difference was attributed to the difference in prey animals, as numerous small mammal parasites were detected from feces of leopard cats inhabiting inland areas and fish parasites from feces of Eurasian otters and raccoon dogs inhabiting waterside areas. Furthermore, 5 zoonotic parasites known to infect humans were identified at the species level. Wildlife-associated zoonoses are expected to increase as the proximity between humans and wildlife increases due to urbanization. Vigilance may be necessary, such as by monitoring parasites in wildlife feces, as was done in this study.

Integration of ITS-2 rDNA nemabiome metabarcoding with Fecal Egg Count Reduction Testing (FECRT) reveals ivermectin resistance in multiple gastrointestinal nematode species, including hypobiotic Ostertagia ostertagi, in western Canadian beef cattle

A large-scale Fecal Egg Count Reduction Test (FECRT) was integrated with ITS-2 rDNA nemabiome metabarcoding to investigate anthelmintic resistance in gastrointestinal nematode (GIN) parasites in western Canadian beef cattle. The study was designed to detect anthelmintic resistance with the low fecal egg counts that typically occur in cattle in northern temperate regions. Two hundred and thirty-four auction market-derived, fall-weaned steer calves coming off pasture were randomized into three groups in feedlot pens: an untreated control group, an injectable ivermectin treatment group, and an injectable ivermectin/oral fenbendazole combination treatment group. Each group was divided into six replicate pens with 13 calves per pen. Individual fecal samples were taken pre-treatment, day 14 post-treatment, and at monthly intervals for six months for strongyle egg counting and metabarcoding. Ivermectin treatment resulted in an 82.4% mean strongyle-type fecal egg count reduction (95% CI 67.8-90.4) at 14 days post-treatment, while the combination treatment was 100% effective, confirming the existence of ivermectin-resistant GIN. Nemabiome metabarcoding of third-stage larvae from coprocultures revealed an increase in the relative abundance of Cooperia oncophora, Cooperia punctata, and Haemonchus placei at 14 days post-ivermectin treatment indicating ivermectin resistance in adult worms. In contrast, Ostertagia ostertagi third-stage larvae were almost completely absent from day 14 coprocultures, indicating that adult worms of this species were not ivermectin resistant. However, there was a recrudescence of O. ostertagi third stage larvae in coprocultures at three to six months post-ivermectin treatment, which indicated ivermectin resistance in hypobiotic larvae. The calves were recruited from the auction market and, therefore, derived from multiple sources in western Canada, suggesting that ivermectin-resistant parasites, including hypobiotic O. ostertagi larvae, are likely widespread in western Canadian beef herds. This work demonstrates the value of integrating ITS-2 rDNA metabarcoding with the FECRT to enhance anthelmintic resistance detection and provide GIN species- and stage-specific information.

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.

Low-cost molecular methods to characterise gastrointestinal nematode co-infections of goats in Africa

BACKGROUND

Veterinary diagnostics aid intervention strategies, track zoonoses, and direct selective breeding programs in livestock. In ruminants, gastrointestinal nematode (GIN) parasites are a major cause of production losses, but morphologically similar species limit our understanding of how specific GIN co-infections impact health in resource-limited settings. To estimate the presence and relative abundance of GINs and other helminths at the species level, we sought to develop a low-cost and low-resource molecular toolkit applied to goats from rural Malawi smallholdings.

METHODS

Goats were subjected to health scoring and faecal sampling on smallholdings in Lilongwe district, Malawi. Infection intensities were estimated by faecal nematode egg counts with a faecal subsample desiccated for DNA analysis. Two DNA extraction methods were tested (low-resource magbead kit vs high-resource spin-column kit), with resulting DNA screened by endpoint polymerase chain reaction (PCR), semi-quantitative PCR, quantitative PCR (qPCR), high-resolution melt curve analysis (HRMC), and 'nemabiome' internal transcribed spacer 2 (ITS-2) amplicon sequencing.

RESULTS

Both DNA isolation methods yielded comparable results despite poorer DNA purity and faecal contaminant carryover from the low-resource magbead method. GINs were detected in 100% of samples regardless of infection intensity. Co-infections with GINs and coccidia (Eimeria spp.) were present in most goats, with GIN populations dominated by Haemonchus contortus, Trichostrongylus colubriformis, Trichostrongylus axei, and Oesophagostomum columbianum. Both multiplex PCR and qPCR were highly predictive of GIN species proportions obtained using nemabiome amplicon sequencing; however, HRMC was less reliable than PCR in predicting the presence of particular species.

CONCLUSIONS

These data represent the first 'nemabiome' sequencing of GINs from naturally infected smallholder goats in Africa and show the variable nature of GIN co-infections between individual animals. A similar level of granularity was detected by semi-quantitative PCR methods, which provided an accurate summary of species composition. Assessing GIN co-infections is therefore possible using cost-efficient low-resource DNA extraction and PCR approaches that can increase the capacity of molecular resources in areas where sequencing platforms are not available; and also open the door to affordable molecular GIN diagnostics. Given the diverse nature of infections in livestock and wildlife, these approaches have potential for disease surveillance in other areas.

Amplicon sequencing allows differential quantification of closely related parasite species: an example from rodent Coccidia (Eimeria)

BACKGROUND

Quantifying infection intensity is a common goal in parasitological studies. We have previously shown that the amount of parasite DNA in faecal samples can be a biologically meaningful measure of infection intensity, even if it does not agree well with complementary counts of transmission stages (oocysts in the case of Coccidia). Parasite DNA can be quantified at relatively high throughput using quantitative polymerase chain reaction (qPCR), but amplification needs a high specificity and does not simultaneously distinguish between parasite species. Counting of amplified sequence variants (ASVs) from high-throughput marker gene sequencing using a relatively universal primer pair has the potential to distinguish between closely related co-infecting taxa and to uncover the community diversity, thus being both more specific and more open-ended.

METHODS

We here compare qPCR to the sequencing-based amplification using standard PCR and a microfluidics-based PCR to quantify the unicellular parasite Eimeria in experimentally infected mice. We use multiple amplicons to differentially quantify Eimeria spp. in a natural house mouse population.

RESULTS

We show that sequencing-based quantification has high accuracy. Using a combination of phylogenetic analysis and the co-occurrence network, we distinguish three Eimeria species in naturally infected mice based on multiple marker regions and genes. We investigate geographical and host-related effects on Eimeria spp. community composition and find, as expected, prevalence to be largely explained by sampling locality (farm). Controlling for this effect, the novel approach allowed us to find body condition of mice to be negatively associated with Eimeria spp. abundance.

CONCLUSIONS

We conclude that amplicon sequencing provides the underused potential for species distinction and simultaneous quantification of parasites in faecal material. The method allowed us to detect a negative effect of Eimeria infection on the body condition of mice in the natural environment.

Preliminary results of the recombinase polymerase amplification technique for the detection of Haemonchus contortus from Hungarian field samples

Haemonchus contortus is a parasitic nematode of small ruminants responsible for significant economic losses and animal health concerns globally. Detection of gastrointestinal nematode (GIN) infection in veterinary practice typically relies on microscopy-based methods such as the faecal egg count and morphological identification of larval culture. However, mixed co-infections are common and species-specific identification is typically time-consuming and expertise-intensive. Compounded by increasing anthelmintic resistance, there is an urgent need to implement the molecular diagnosis of GIN in the livestock industry, preferably in field settings. Advances in isothermal amplification techniques including recombinase polymerase amplification (RPA) assays could improve this. Yet, constraints in RPA kit availability and amplicon detection systems limit the use of this technology in point of care settings. In this study, we present an early-stage, proof-of-concept demonstration of RPA targeting the internal transcribed spacer (ITS2) region of H. contortus. Having tested against eight closely related nematodes and also against five farm isolates in Eastern Hungary, preliminary results derived from a comparative analysis of 3 primer sets showed the assay detects H. contortus DNA and has a limit of detection of 10^-5 ng/μl. We also tested an end-result naked eye detection system using various DNA binding dyes, of which EvaGreen® dye was successful for a qualitative RPA detection that could be adaptable at farm sites.

Broad-Spectrum Inhibitors for Conserved Unique Phosphoethanolamine Methyltransferases in Parasitic Nematodes Possess Anthelmintic Efficacy

In humans, nematode infections are prevalent in developing countries, causing long-term ill health, particularly in children. Worldwide, nematode infections are prevalent in livestock and pets, affecting productivity and health. Anthelmintic drugs are the primary means of controlling nematodes, but there is now high prevalence of anthelmintic resistance, requiring urgent identification of new molecular targets for anthelmintics with novel mechanisms of action. Here, we identified orthologous genes for phosphoethanolamine methyltransferases (PMTs) in nematodes within the families Trichostrongylidae, Dictyocaulidae, Chabertiidae, Ancylostomatoidea, and Ascarididae. We characterized these putative PMTs and found that they possess bona fide PMT catalytic activities. By complementing a mutant yeast strain lacking the ability to synthesize phosphatidylcholine, the PMTs were validated to catalyze the biosynthesis of phosphatidylcholine. Using an in vitro phosphoethanolamine methyltransferase assay with PMTs as enzymes, we identified compounds with cross-inhibitory effects against the PMTs. Corroboratively, treatment of PMT-complemented yeast with the PMT inhibitors blocked growth of the yeast, underscoring the essential role of the PMTs in phosphatidylcholine synthesis. Fifteen of the inhibitors with the highest activity against complemented yeast were tested against Haemonchus contortus using larval development and motility assays. Among them, four were found to possess potent anthelmintic activity against both multiple drug-resistant and susceptible isolates of H. contortus, with IC50 values (95% confidence interval) of 4.30 μM (2.15-8.28), 4.46 μM (3.22-6.16), 28.7 μM (17.3-49.5), and 0.65 μM (0.21-1.88). Taken together, we have validated a molecular target conserved in a broad range of nematodes and identified its inhibitors that possess potent in vitro anthelmintic activity.