Population structure of Haemonchus contortus from seven geographical regions in China, determined on the basis of microsatellite markers

Global scenario of genetic diversity in cox1 and nad1 genes of Moniezia expansa

Monieziasis is a parasite-borne production-limiting disease of livestock. Moniezia expansa is the most important species having cosmopolitan distribution. Despite of numerous prevalence reports, very little information is available about the evolutionary biology and population genetics of M. expansa. To close this research gap, this study was undertaken to recognize and inspect the genetic variation of M. expansa populations around the world using the cox1 and nad1 genes and deduce phylogenetic relationships with M. expansa populations. The cox1 and nad1 gene sequences were downloaded from the NCBI GenBank database. Followed by sequence alignment, median-joining networks were constructed using PopArt software. Diversity and neutrality indices were computed through DnaSp software while MEGA software was used to draw the maximum-likelihood phylogenetic tree. Thirty-two cox1 sequences, from five different countries, and 9 nad1 sequences from three different countries, were among the sequences used in this study. The cox1 and nad1 gene sequences had mutations in 97 and 36 different places, respectively. Twenty and 7 unique haplotypes were discovered for the cox1 and nad1 gene sequences, respectively. Comparable haplotype diversities were observed for both the genes under study (cox1 = 0.950; nad1 = 0.944). Negative Tajima's D and Fu Fs were found for the cox1 gene while these indices were positive for the nad1 gene. Phylogenetic analysis also showed the existence of unique haplotypes for both the cox1 and nad1 genes. The results of this study indicate that there is the existence of a huge genetic diversity in M. expansa isolates. For future studies, it is recommended that longer gene sequences should be used to describe genetic variation among M. expansa isolates as the length of the gene under study affects the genetic variation. Moreover, additional mitochondrial markers should also be investigated because the assertive strength of a group of gene targets is superior to defining genetic diversity.

Characterization and population genetics of Haemonchus contortus in Merino sheep in Lesotho

Haemonchus contortus is the most pathogenic and economically restrictive gastrointestinal nematode in the small ruminant industry globally. Morbidity, poor cross-bodily state, and mortality of sheep in Lesotho suggest the presence of H. contortus. The present study investigated the morphological, molecular, and population genetics of H. contortus third-stage larvae infecting sheep in four ecological zones (EZ) of Lesotho. Coprocultures were prepared for larval morphological identification and PCR determination. Larvae were identified morphologically as 100% H. contortus. The Second Internal Transcribed Spacer (ITS-2) gene of the ribosomal DNA of H. contortus isolates in the present study revealed nucleotide homology ranging from 97 to 100% when compared with selected GenBank reference sequences. Pairwise evolutionary divergence among H. contortus isolates was low, with 0.01318 recorded as the highest in the present study. Five haplotypes resulted from 14 Lesotho sequences. Haplotype diversity and nucleotide diversity were 0.76923 and 0.00590, respectively. Genetic differentiation among isolates was low but not statistically significant. An analysis of molecular variance revealed that most molecular variation was distributed within topographic populations at 94.79% (FST = 0.05206, p > 0.05) and 5.21% among populations. There was high gene flow and no definite population genetic structure among Lesotho isolates.

Genetic variation and population structure of Fasciola hepatica: an in silico analysis

Fasciola hepatica is a trematode leading to heavy economic setbacks to the livestock sector globally. The population's genetic information and intimate kinship level are frequently assessed using analysis of mitochondrial DNA. In this analysis, we retrieved cox1 (n = 247) and nad1 (n = 357) sequences of F. hepatica from the NCBI GenBank database and aligned the sequences with the respective reference sequences using MEGA software. The median joining network was drawn using PopArt software while neutrality and diversity indices were estimated with the help of DnaSp software. Neighbor-joining phylogenetic tree was constructed using the MEGA software package. A total of 46 and 98 distinctive haplotypes were observed for cox1 and nad1 genes, respectively. Diversity indices indicated high haplotype and nucleotide diversities in both genes. Positive Tajima's D and Fu's Fs values were found for the entire population of both the genes under study. The cox1 and nad1 gene segments in this study showed high Tajima's D values, suggesting a low likelihood of future population growth. The Tajima's D value of the nad1 gene sequence is lower (2.14910) than that of the cox1 gene sequence (3.40314), which suggests that the former is growing at a slower rate. However, the region-wise analysis revealed that both the cox1 and nad1 genes showed deviation from neutrality suggesting a recent population expansion as a result of an excess of low-frequency polymorphism. Furthermore, the overall host-wise analysis showed positive and significant Tajima's D values for the cox1 and nad1 gene sequences. To the best of our knowledge, this is the first attempt to provide insights into genetic variations and population structure of F. hepatica at a global scale using cox1 and nad1 genes. Our findings suggest the existence of specific variants of F. hepatica in different parts of the world and provide information on the molecular ecology of F. hepatica. The results of this study also mark a critical development in upcoming epidemiological investigations on F. hepatica and will also contribute to understanding the global molecular epidemiology and population structure of F. hepatica.

Selection of genome-wide SNPs for pooled allelotyping assays useful for population monitoring

Parasitic worms are serious pests of humans, livestock and crops worldwide. Multiple management strategies are employed in order to reduce their impact, and some of these may affect their genome and population allelic frequency distribution. The evolution of chemical resistance, ecological changes, and pest dispersal have allowed an increasing number of pests to become difficult to control with current management methods. Their lifestyle limits the use of ecological and individual-based management of populations. There is a need to develop rapid, affordable, and simple diagnostics to assess the efficacy of management strategies and delay the evolution of resistance to these strategies. This study presents a multi-locus, equal-representation, whole genome pooled SNPs selection approach as a monitoring tool for the ovine nematode parasite Haemonchus contortus. The SNP selection method used two reference genomes of different quality, then validated these SNPs against a high-quality recent genome assembly. From over 11 million high-quality SNPs identified, 334 SNPs were selected, of which 262 were species-specific, yielded similar allele frequencies when assessed as multiple individuals or as pools of individuals, and suitable to distinguish mixed nematode isolate pools from single isolate pools. As a proof-of-concept, 21 Australian H. contortus populations with various phenotypes and genotypes were screened. This analysis confirmed the overall low-level of genetic differentiation between populations collected from the field, but clearly identifying highly inbred populations, and populations showing genetic signatures associated with chemical resistance. The analysis showed that 66% of the SNPs were necessary for stability in assessing population genetic patterns, and SNP pairs did not show linkage according to allelic frequencies across the 21 populations. This method demonstrates that ongoing monitoring of parasite allelic frequencies and genetic changes can be achieved as a management assessment tool to identify drug-treatment failure, population incursions, and inbreeding signatures due to selection. The SNP selection method could also be applied to other parasite species.

Genetic diversity and population structure analyses based on microsatellite DNA of parthenogenetic Fasciola flukes obtained from cattle and sika deer in Japan

Understanding the population structure of Fasciola flukes in domestic and wild animals is important for determining the extent of cross-infection between them. Although the parthenogenetic Fasciola flukes in Japan have been shown to comprise five genetic types based on the ribosomal internal transcribed spacer 1 (ITS1) and mitochondrial NADH dehydrogenase subunit 1 (nad1) regions, these genetic regions are not suitable for analyzing their population structure. In the present study, the genetic diversity and population structure of the parthenogenetic Fasciola flukes in Japan were studied using microsatellite DNA, ITS1, and nad1 regions. A total of 144 parthenogenetic Fasciola flukes, obtained from cattle and sika deer in 16 localities, were individually analyzed using PCR-RFLP for ITS1, PCR-direct sequence analysis for nad1, and post-labeling PCR and capillary electrophoresis for microsatellite DNA regions. The flukes showed higher genetic diversity in the microsatellite DNA regions than ITS1 and nad1. The population structures of parthenogenetic Fasciola flukes were unclear, however, it was suggested that the flukes are more diverse populations. We hypothesized that their distribution throughout Japan is closely related to livestock movement dependent on human activity. Moreover, it is considered that cross-infection of the flukes between cattle and sika deer possibly has occurred in the past.

Ovine haemonchosis: a review

Sheep farming is the backbone of a rural economy in developing countries, and haemonchosis is a major impediment in the way of its progress. Haemonchus contortus (H. contortus) infection persists all over the world particularly in the tropical and sub-tropical regions. Various review articles have been published to substantially cover one or more aspects of its morphology, prevalence, pathogenesis, symptoms, diagnosis, immune response, drug resistance, treatment, and control measure. The objective of this paper is to briefly review past and present information available in the aforementioned areas in one place to enable the readers to fully understand the problem from a broader perspective. H. contortus parasite harbours in abomasum of affected animal and feeds on its blood, producing mild to severe symptoms and even death in acute form. The parasite thus inflicts heavy production losses and is of economic importance. H. contortus has developed diverse characters over the years leading to limited success in the production of vaccines. Indiscriminate use of the anthelmintics has produced drug resistance against almost all conventional products. Efficacy of medicinal plants and non-conventional chemicals has been reported under controlled experiments; however, research on their adverse effects on growth and fertility is yet to be studied. Research on molecular tools for identification and introduction of resistant genes into the flock is also underway but still a long journey to find its field application. Crossbreeding may compromise the production traits of the existing flock. In given circumstances, a targeted selective treatment approach along with selective breeding, culling of more susceptible animals, and maintaining a good body condition score through the provision of a balanced diet remains a workable strategy to control haemonchosis in sheep.

A Whole Genome Re-Sequencing Based GWA Analysis Reveals Candidate Genes Associated with Ivermectin Resistance in Haemonchus contortus

The most important and broad-spectrum drug used to control the parasitic worms to date is ivermectin (IVM). Resistance against IVM has emerged in parasites, and preserving its efficacy is now becoming a serious issue. The parasitic nematode Haemonchus contortus (Rudolphi, 1803) is economically an important parasite of small ruminants across the globe, which has a successful track record in IVM resistance. There are growing evidences regarding the multigenic nature of IVM resistance, and although some genes have been proposed as candidates of IVM resistance using lower magnification of genome, the genetic basis of IVM resistance still remains poorly resolved. Using the full magnification of genome, we herein applied a population genomics approach to characterize genome-wide signatures of selection among pooled worms from two susceptible and six ivermectin-resistant isolates of H. contortus, and revealed candidate genes under selection in relation to IVM resistance. These candidates also included a previously known IVM-resistance-associated candidate gene HCON_00148840, glc-3. Finally, an RNA-interference-based functional validation assay revealed the HCON_00143950 as IVM-tolerance-associated gene in H. contortus. The possible role of this gene in IVM resistance could be detoxification of xenobiotic in phase I of xenobiotic metabolism. The results of this study further enhance our understanding on the IVM resistance and continue to provide further evidence in favor of multigenic nature of IVM resistance.

Does host socio-spatial behavior lead to a fine-scale spatial genetic structure in its associated parasites?

Gastro-intestinal nematodes, especially Haemonchus contortus, are widespread pathogenic parasites of small ruminants. Studying their spatial genetic structure is as important as studying host genetic structure to fully understand host-parasite interactions and transmission patterns. For parasites having a simple life cycle (e.g., monoxenous parasites), gene flow and spatial genetic structure are expected to strongly rely on the socio-spatial behavior of their hosts. Based on five microsatellite loci, we tested this hypothesis for H. contortus sampled in a wild Mediterranean mouflon population (Ovis gmelini musimon × Ovis sp.) in which species- and environment-related characteristics have been found to generate socio-spatial units. We nevertheless found that their parasites had no spatial genetic structure, suggesting that mouflon behavior was not enough to limit parasite dispersal in this study area and/or that other ecological and biological factors were involved in this process, for example other hosts, the parasite life cycle, or the study area history.

Analysis of genome-wide SNPs based on 2b-RAD sequencing of pooled samples reveals signature of selection in different populations of Haemonchus contortus

The parasitic nematode Haemonchus contortus is one of the world's most important parasites of small ruminants that causes significant economic losses to the livestock sector. The population structure and selection in its various strains are poorly understood. No study so far compared its different populations using genome-wide data. Here, we focused on different geographic populations of H. contours from China (Tibet, TB; Hubei, HB; Inner Mongolia, IM; Sichuan, SC), UK and Australia (AS), using genome-wide population-genomic approaches, to explore genetic diversity, population structure and selection. We first performed next-generation high-throughput 2b RAD pool sequencing using Illumina technology, and identified single-nucleotide polymorphisms (SNPs) in all the strains. We identified 75,187 SNPs for TB, 82,271 for HB, 82,420 for IM, 79,803 for SC, 83,504 for AS and 78,747 for UK strain. The SNPs revealed low-nucleotide diversity (pi= 0.0092-0.0133) within each strain, and a significant differentiation level (average F_st = 0.34264) among them. Chinese populations TB and SC, along with the UK strain, were more divergent populations. Chinese populations IM and HB showed affinities to the Australian strain. We then analysed signature of selection and detected 44 (UK) and 03 (AS) private selective sweeps containing 49 and 05 genes, respectively. Finally, we performed the functional annotation of selective sweeps and proposed biological significance to signature of selection. Our data suggest that 2b-RAD pool sequencing can be used to assess the signature of selection in H. contortus.

Microsatellite analysis reveals extensive gene flow, and lack of population structure in the farm populations of Haemonchus contortus in northern China

The parasitic nematode Haemonchus contortus is economically an important parasite of small ruminants across the globe. China is the world's largest producer, consumer, and importer of mutton. With ubiquitous distribution across the country H. contortus is one of the potential candidates to cause huge economic losses to small ruminant farming industry in China. We herein investigated genetic diversity and population structure of six farm populations of H. contortus in northern China, and also compared them to H. contortus isolates from UK and Australia. We first prepared individual DNA samples from 240 adult worms, and generated genotyping data using eight microsatellite markers. Obtained data was then subjected to allelic frequency and population genetic analyses. The overall allelic richness (mean/locus/pop = 7.375 ± 0.844-10.125 ± 1.109), and expected heterozygosity (mean/locus/pop = 0.646 ± 0.040-0.735 ± 0.025) indicated high degree of population genetic variation across the Chinese isolates. Low level of genetic differentiation (Fst = 0.010-0.066) was observed across all the populations. AMOVA results showed high level of variation (93%) within the populations. PCA analysis revealed mixed clustering of all the populations with no visible geographical sub-structuring. Finally the population admixture analysis resulted in extensive admixing of genotypes across all the populations. With these findings we conclude that there is no obvious population genetic structure with extensive gene flow across all the farm populations of H. contortus in northern China.

Population genetics of benzimidazole-resistant Haemonchus contortus and Haemonchus placei from buffalo and cattle: implications for the emergence and spread of resistance mutations

The population genetics of nematode parasites are poorly understood with practical reference to the selection and spread of anthelmintic resistance mutations. Haemonchus species are important to study the nematode population genetics due to their clinical importance in ruminant livestock, and the availability of genomic resources. In the present study, it has been examined that Haemonchus contortus and Haemonchus placei populations from three buffalo and nine cattle hosts. Seventy-three individual adult worms of H. contortus and 148 of H. placei were analysed using a panel of seven microsatellite markers. The number of alleles per locus in H. contortus and H. placei indicated that all populations were polymorphic for the microsatellites used in the present study. Genetic diversity parameters included high levels of allelic richness and heterozygosity, indicating effective population sizes, high mutation rates and high transmission frequencies in the area. Genetic structure parameters revealed low genetic differentiation between and high levels of genetic variation within H. contortus and H. placei populations. Population dynamic analyses showed an absence of heterozygosity excess in both species, suggesting that there was no deviation from genetic drift equilibrium. Our results provide a proof of concept for better understanding of the consequences of specific control strategies, climatic change or management strategies on the population genetics of anthelmintic resistance alleles in Haemonchus spp. infecting co-managed buffalo and cattle.

Recent Research Progress in China on Haemonchus contortus

Haemonchus contortus is one of the most important parasites of ruminants with worldwide distribution that can bring huge economic losses to the breeding industry of cattle, sheep, and goats. In recent 20 years, studies on H. contortus in China mainly focused on the epidemiology, population genetics, anthelmintic resistance, structural and functional studies of important genes regulating the development of this parasite, interaction between parasite molecules and host cells and vaccine development against haemonchosis, and achieved good progress. However, there is no systematic review about the studies by Chinese researchers on H. contortus in China. The purpose of this review is to bring together the findings from the studies on H. contortus in China in order to obtain the knowledge gained from the recent studies in China and provide foundation for identifying future research directions to establish novel diagnostic methods, discover new drug targets and vaccine candidates for use in preventing and controlling H. contortus in China.