Validation of ITS-2 rDNA nemabiome sequencing for ovine gastrointestinal nematodes and its application to a large scale survey of UK sheep farms

Abstract

We have validated ITS-2 rDNA nemabiome next-generation amplicon sequencing to determine relative species abundance of gastrointestinal nematode species in ovine fecal samples. In order to determine species representation biases, ITS-2 rDNA amplicon sequencing was applied to mock communities or field populations with known proportions of L3 for eight of the major ovine gastrointestinal nematode species: Teladorsagia circumcincta, Trichostrongylus vitrinus, Haemonchus contortus, Cooperia curticei, Trichostrongylus axei, Trichostrongylus colubriformis, Chabertia ovina and Oesophagostumum venulosum. Correction factors, calculated from this data, were shown to reduce species representation biases when applied to an independent set of field samples of known composition. We compared ITS-2 rDNA amplicon sequencing data that was generated from harvested eggs, freshly hatched L1 or L3 larvae following fecal culture and no statistically significant differences were found for the more abundant parasite species. We then applied the validated ITS-2 rDNA nemabiome amplicon sequencing assay to a set of archived L1 gastrointestinal nematode populations, collected in 2008 from fecal samples from 93 groups of 20 ewes and 61 groups of 20 lambs derived from 99 UK sheep farms. The presence of the major gastrointestinal nematode species had previously been determined on this large sample set by species-specific PCR. We show how the ITS-2rDNA amplicon sequencing data provided much more detailed information on species abundance than the previous species-specific PCR. This new data represents the most comprehensive overview of the relative abundance of the major gastrointestinal nematode species across UK sheep farms to date. Substantial variance in the relative abundance of both T. circumcincta and T. vitrinus between farms was revealed with the former species being of statistically significantly higher abundance in all three regions sampled (England, Scotland and Wales). The data also revealed that the relative abundance of T. circumcinta was statistically significantly higher in ewes than in lambs with the opposite pattern being the case for T. vitrinus. The nemabiome sequencing data also clearly illustrated the sporadic nature and skewed distribution of H. controtus across UK sheep farms as well as a higher relative abundance on farms from England compared to Wales and Scotland. The nemabiome survey also provides the first widescale data on the relative abundance of the two major large intestinal nematodes C. ovina and O. venulosum. This work validates ITS-2 rDNA nemabiome sequencing for use in sheep and illustrates the power of the approach for large scale surveillance of ovine gastrointestinal nematodes.