Genetic diversity and multiplicity of infection in Fasciola gigantica isolates of Pakistani livestock

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.

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.

Dispersal direction of Malaysian Fasciola gigantica from neighboring southeast Asian countries inferred using mitochondrial DNA analysis

Fasciola gigantica and hybrid Fasciola flukes, responsible for the disease fasciolosis, are found in Southeast Asian countries. In the present study, we performed molecular species identification of Fasciola flukes distributed in Terengganu, Malaysia using multiplex PCR for phosphoenolpyruvate carboxykinase (pepck) and PCR-restriction fragment length polymorphism (RFLP) for DNA polymerase delta (pold). Simultaneously, phylogenetic analysis based on mitochondrial NADH dehydrogenase subunit 1 (nad1) was performed for the first time on Malaysian Fasciola flukes to infer the dispersal direction among neighboring countries. A total of 40 flukes used in this study were identified as F. gigantica. Eight nad1 haplotypes were identified in the F. gigantica population of Terengganu. Median-joining network analysis revealed that the Malaysian population was related to those obtained from bordering countries such as Thailand and Indonesia. However, genetic differentiation was detected using population genetics analyses. Nevertheless, the nucleotide diversity (π) value suggested that F. gigantica with the predominant haplotypes was introduced into Malaysia from Thailand and Indonesia. The dispersal direction suggested by population genetics in the present study may not be fully reliable since Fasciola flukes were collected from a single location in one state of Malaysia. Further studies analyzing more samples from many locations are required to validate the dispersal direction proposed herein.

A novel metabarcoded deep amplicon sequencing tool for disease surveillance and determining the species composition of Trypanosoma in cattle and other farm animals

The World Health Organization (WHO) and the Food and Agriculture Organization (FAO) have developed strategies to control trypanosomiasis in humans and livestock in endemic areas. These require a better understanding of the distribution of different Trypanosoma species and improved predictions of where they might appear in the future, based on accurate diagnosis and robust surveillance systems. Here, we describe a metabarcoding deep amplicon sequencing method to identify and determine the Trypanosoma species in co-infecting communities. First, four morphological verified Trypanosoma species (T. brucei, T. congolense, T. vivax and T. theileri) were used to prepare test DNA pools derived from different numbers of parasites to evaluate the method's detection threshold for each of the four species and to assess the accuracy of their proportional quantification. Having demonstrated the accurate determination of species composition in Trypanosoma communities, the method was applied to determine its detection threshold using blood samples collected from cattle with confirmed Trypanosoma infections based on a PCR assay. Each sample showed a different Trypanosoma species composition based on the proportion of MiSeq reads. Finally, we applied the assay to field samples to develop new insight into the species composition of Trypanosoma communities in cattle, camels, buffalo, horses, sheep, and goat in endemically infected regions of Pakistan. We confirmed that Trypanosoma evansi is the major species in Pakistan and for the first time showed the presence of Trypanosoma theileri. The metabarcoding deep amplicon sequencing method and bioinformatics pathway have several potential applications in animal and human research, including evaluation of drug treatment responses, understanding of the emergence and spread of drug resistance, and description of species interactions during co-infections and determination of host and geographic distribution of trypanosomiasis in humans and livestock.

A cross-sectional survey of hard ticks and molecular characterization of Rhipicephalus microplus parasitizing domestic animals of Khyber Pakhtunkhwa, Pakistan

Background

In tropical and subtropical countries, tick infestation causes major public health problems and considerable financial losses to the livestock industry. This study was aimed to assess the species composition of richness and analyze the phylogeny of Rhipicephalus microplus in the District Bannu of Khyber Pakhtunkhwa, Pakistan.

Methods

Collected ticks were identified morphologically and DNA extracted from R. microplus was amplified and subjected to sequencing.

Results

A total of 3,600 animals were examined among them 1,494 animals were found to be infested with ticks, including 669 cows, 476 buffaloes, 163 goats, and 186 sheep (p = 0.001). Tick infestation was significantly high (43.58%) in animals of age group (<1 year) (p-value = 0.027). Female animals were more (44.05%) infested with ticks than males (34.43%) (p = 0.001). The intensity of infestation was significantly higher in summer (77.49%) (p = 0.001). A total of 5,557 ticks were collected comprising three genera and six species. R. microplus was predominantly prevalent (n = 1,474; 26.52%), followed by Rhipicephalus annulatus (n = 1,215; 21.86%), Hyalomma anatolicum (n = 1,139; 20.49%), Hyalomma marginatum (n = 1,086; 19.54%), and Rhipicephalus turanicus (n = 761; 13.69%), while the least common was Haemaphysalis aciculifer (n = 80; 1.43%) (p = 0.001). Morphologically identified R. microplus species were also analyzed genetically by using two genetic markers 16S ribosomal RNA (16S rRNA) and internal transcribed spacer 2 (ITS2) genes. The phylogenetic study revealed that R. microplus is genetically diversified and clustered in clade B with R. microplus species from China, India, and Pakistan.

Conclusion

Ticks infestation was significantly correlated with various factors including age, sex, season, and animal type. R. microplus genetically resembled species reported from India and China. However, major knowledge gaps concerning various species of ticks exist and many areas are still unexplored in Pakistan. Therefore, it is necessary to explore the epidemiological and molecular aspects of various tick species in other regions of southern Khyber Pakhtunkhwa.

High-throughput sequencing of Fasciola spp. shows co-infection and intermediate forms in Balochistan, but only Fasciola gigantica in the Punjab province of Pakistan

Fasciola gigantica and Fasciola hepatica are digenetic trematodes causing fasciolosis in ruminants. The host and geographical distribution of both Fasciola species are influenced by environmental and climatic conditions favouring survival and development of free-living stages and intermediate hosts, and livestock management practices. The aim of the present study was to describe the host distribution of the two Fasciola species in buffalo, cattle, goats, and sheep in the Balochistan and Punjab provinces of Pakistan. 359 flukes were collected from a total of 32 livers from the four livestock species. Deep amplicon sequencing of the internal transcribed spacer region 2 of ribosomal DNA (rDNA ITS-2) and mitochondrial nicotinamide adenine dinucleotide dehydrogenase 1 (mtDNA ND-1) loci confirmed co-infection of F. hepatica and F. gigantica in Balochistan and single species F. gigantica infection in Punjab. In Balochistan, co-infections and hybrids of both Fasciola species were identified in cattle, with more F. hepatica detected than F. gigantica. However, F. hepatica was the only species identified in goats, and F. gigantica was the only species identified in buffalo. In Punjab, all flukes were confirmed as F. gigantica in each of the four livestock species. Overall, the results indicate differences in the host and geographical distribution of F. gigantica and F. hepatica, and provide useful knowledge for the development of control strategies for livestock and humans.

Towards the comprehension of fasciolosis (re-)emergence: an integrative overview

The increasing distribution and prevalence of fasciolosis in both human and livestock are concerning. Here, we examine the various types of factors influencing fasciolosis transmission and burden and the interrelations that may exist between them. We present the arsenal of molecules, 'adjusting' capabilities and parasitic strategies of Fasciola to infect. Such features define the high adaptability of Fasciola species for parasitism that facilitate their transmission. We discuss current environmental perturbations (increase of livestock and land use, climate change, introduction of alien species and biodiversity loss) in relation to fasciolosis dynamics. As Fasciola infection is directly and ultimately linked to livestock management, living conditions and cultural habits, which are also changing under the pressure of globalization and climate change, the social component of transmission is also discussed. Lastly, we examine the implication of increasing scientific and political awareness in highlighting the current circulation of fasciolosis and boosting epidemiological surveys and novel diagnostic techniques. From a joint perspective, it becomes clear that factors weight differently at each place and moment, depending on the biological, environmental, social and political interrelating contexts. Therefore, the analyses of a disease as complex as fasciolosis should be as integrative as possible to dissect the realities featuring each epidemiological scenario. Such a comprehensive appraisal is presented in this review and constitutes its main asset to serve as a fresh integrative understanding of fasciolosis.

Molecular characterization and phylogenetic analyses of Fasciola gigantica of buffaloes and goats in Punjab, Pakistan

Fasciola gigantica is considered to be a major pathogen causing fasciolosis in the Indian subcontinent, resulting in production losses of millions of dollars in the livestock industry. Understading the dispersal origin and the patterns of spread of F. gigantica is important. A total of 53 Fasciola flukes collected from buffaloes and goats in Punjab, Pakistan between 2017 and 2018 were identified as F. gigantica based on the multiplex PCR for the phosphoenolpyruvate carboxykinase (pepck) and the PCR-restriction fragment length polymorphism (RFLP) for DNA polymerase delta (pold). A significant genetic difference between F. gigantica from buffaloes and goats was indicated by the genetic analyses of mitochondrial markers, NADH dehydrogenase subunit 1 (nad1) and cytochrome C oxidase subunit 1 (cox1). Phylogenetic analysis of the seventeen nad1 haplotypes of F. gigantica from Pakistan with those in neighbouring countries of the Indian subcontinent revealed that all the haplotypes identified in Pakistan were clustered in haplogroup A. fasciola gigantica with the eight haplotypes might be expanded in Pakistan from Indian origin, along with the migration of the domestic animals, since they were related to Indian haplotypes. In contrast, the remaining nine haplotypes were not shared with any neighbouring countries, suggesting independent origin, probably from neighbouring Middle East countries. However, cautious interpretation is required due to the very limited samples size of this study. Our study provides a proof of concept for a method that could be used to investigate the epidemiology of F. gigantica.

Molecular confirmation of Dicrocoelium dendriticum in the Himalayan ranges of Pakistan

Lancet liver flukes of the genus Dicrocoelium (Trematoda: Digenea) are recognised parasites of domestic and wild herbivores. The aim of the present study was to confirm the species identity of Dicrocoeliid flukes collected from the Chitral valley in the Himalayan ranges of Pakistan. The morphology of 48 flukes belonging to eight host populations was examined; but overlapping traits prevented accurate species designation. Phylogenetic comparison of published D. dendriticum ribosomal cistron DNA, and cytochrome oxidase-1 (COX-1) mitochondrial DNA sequences with those from D. chinensis was performed to assess within and between species variation and re-affirm the use of species-specific single nucleotide polymorphism markers. PCR and sequencing of 34 corresponding fragments of ribosomal DNA and 14 corresponding fragments of mitochondrial DNA from the Chitral valley flukes, revealed 10 and 4 unique haplotypes, respectively. These confirmed for the first time the molecular species identity of Pakistani lancet liver flukes as D. dendriticum. This work provides a preliminary illustration of a phylogenetic approach that could be developed to study the ecology, biological diversity, and epidemiology of Dicrocoeliid lancet flukes when they are identified in new settings.

Contrasting population genetics of co-endemic cattle- and buffalo- derived Theileria annulata

A study was designed to improve understanding of the genetics of Theileria annulata populations in sympatric cattle and Asian buffalo (Bubalus bubalus). The study was undertaken in the Punjab province of Pakistan, where the prevalence of tropical theileriosis is high. Parasite materials were collected from infected animals in defined regions, where cattle and Asian buffalo are kept together. Six satellite DNA markers and a mitochondrial cytochrome b marker were used to explore the multiplicity of T. annulata infection and patterns of emergence and spread of different parasite genotypes. The results show differences in the numbers of unique satellite locus alleles, suggesting that T. annulata is genetically more diverse in cattle- than in buffalo-derived populations. Heterozygosity (H_e) indices based on satellite and cytochrome b loci data show high levels of genetic diversity among the cattle- and buffalo-derived T. annulata populations. When considered in the context of high parasite transmission rates and frequent animal movements between different regions, the predominance of multiple T. annulata genotypes and multiple introductions of infection may have practical implications for the spread of parasite genetic adaptations; such as those conferring vaccine cross-protection against different strains affecting cattle and Asian buffalo, or resistance to antiprotozoal drugs.