Skip to main content

Advertisement

First molecular description, phylogeny and genetic variation of Taenia hydatigena from Nigerian sheep and goats based on three mitochondrial genes

Article metrics

  • 164 Accesses

Abstract

Background

Cysticercosis caused by the metacestode larval stage of Taenia hydatigena is a disease of veterinary and economic importance. A considerable level of genetic variation among isolates of different intermediate hosts and locations has been documented. Generally, data on the genetic population structure of T. hydatigena is scanty and lacking in Nigeria. Meanwhile, similar findings in other cestodes like Echinococcus spp. have been found to be of epidemiological importance. Our aim, therefore, was to characterize and compare the genetic diversity of T. hydatigena population in Nigeria based on three mitochondrial DNA markers as well as to assess the phylogenetic relationship with populations from other geographical regions.

Methods

In the present study, we described the genetic variation and diversity of T. hydatigena isolates from Nigerian sheep and goats using three full-length mitochondrial genes: the cytochrome c oxidase subunit 1 (cox1), NADH dehydrogenase subunit 1 (nad1), and NADH dehydrogenase subunit 5 (nad5).

Results

The median-joining network of concatenated cox1-nad1-nad5 sequences indicated that T. hydatigena metacestodes of sheep origin were genetically distinct from those obtained in goats and this was supported by high FST values of nad1, cox1, and concatenated cox1-nad1-nad5 sequences. Genetic variation was also found to be higher in isolates from goats than from sheep.

Conclusions

To the best of our knowledge, the present study described the genetic variation of T. hydatigena population for the first time in Nigeria using full-length mitochondrial genes and suggests the existence of host-specific variants. The population indices of the different DNA markers suggest that analysis of long mitochondrial DNA fragments may provide more information on the molecular ecology of T. hydatigena. We recommend that future studies employ long mitochondrial DNA sequence in order to provide reliable data that would explain the extent of genetic variation in different hosts/locations and the biological and epidemiological significance.

Background

Cestodes belonging to the genus Taenia infect a wide range of intermediate host species where the metacestode larval stage causes cysticercosis or coenurosis. A few members of this genus such as Taenia solium, T. saginata and T. asiatica are zoonotic and are responsible for taeniosis in humans while others are of veterinary importance.

The tapeworm T. hydatigena uses canids (primarily dogs) as definitive hosts. The metacestode larval stage of T. hydatigena infects mostly domestic animals such as goats, sheep and pigs, resulting in cysticercosis; an infection of veterinary importance that causes huge economic loses especially in livestock production as a result of mortality [1], condemnation of infected organs and carcasses, and the financial cost of diagnosis and inspection [2,3,4,5].

Taenia hydatigena has a global distribution with a prevalence range of 0.1–32.0%, differing between countries and hosts [4,5,6]. The epidemiology is such that prevalence is usually higher in sheep and goats in most African and European countries compared to countries in Asia [1, 4, 5], and higher in pigs in Asian and South American countries than in countries of other regions [4, 6, 7].

In Nigeria, cysticercosis due to T. hydatigena is poorly documented. However, existing reports have demonstrated high infection rates in livestock especially in sheep and goats [7,8,9,10,11,12,13,14]. Although the prevalence of cysticercosis seems to be quite high, knowledge of the molecular ecology and intraspecific variation of T. hydatigena species is lacking. Meanwhile, genetic variation in other cestodes, e.g. Echinococcus granulosus [15,16,17,18,19] has been found to affect host infectivity, epidemiology, as well as control strategies [20]. Similarly, studies on T. saginata and T. solium mitochondrial genome have also demonstrated intraspecific variation [21,22,23,24] suggested to influence the pathological presentations exhibited in different hosts [25].

Globally, there exists a dearth of data on the genetic variation of T. hydatigena. Although a few studies conducted mostly by using partial cox1 and nad1 mitochondrial gene sequences have reported considerable levels of genetic variation among T. hydatigena populations from different geographical regions and hosts [4, 26,27,28,29], the epidemiological implications if any are far from clear and more studies in this regard may further help us understand the genetic variation or genetic population structure of T. hydatigena that could provide better insight on the disease epidemiology.

Further, the limitation of partial gene sequences in inferring the phylogenetic status of taeniids has been largely emphasized for Echinococcus [30, 31]. Therefore, the available data based on partial gene sequences may have limitation in understanding the extent and significance of the genetic variation observed among T. hydatigena. Also, some studies have suggested that complete mitochondrial genes such as nad6, nad5, atp6, nad3 and nad2 may have some advantage over the conventional mitochondrial markers like cox1 and nad1 in investigating the molecular ecology of Taenia spp. [28, 32]. Therefore, in the present study, our aim was to identify and investigate the genetic variation of T. hydatigena population in Nigeria based on the cox1, nad1 and nad5 mitochondrial genes, compare the genetic variation between the genes, and infer the phylogenetic relationships with T. hydatigena populations from other geographical regions.

Methods

Nigeria is a West African country with a population of over 180 million people. It is made up of 36 states and a Federal Capital Territory located in Abuja. These states are further grouped into six geopolitical zones (North-East, North-Central, North-West, South-East, South-South, and South-West) due to ethnicity and common ancestry. The vegetation cover is majorly rainforest in the south and savannah in the north. Borno State is located in the North-East zone of the country and home to about 4 million small ruminants (sheep and goats) [33] that are managed majorly by traditional method of livestock farming.

Isolates (n = 32) analyzed in this study were collected postmortem from goats (n = 24) and sheep (n = 8) in the months of November and December 2018 during routine examination at an abattoir in Maiduguri, Borno State, northeastern Nigeria (11°50′41″N, 13°8′89″E). All isolates were of liver origin.

DNA extraction, amplification and sequencing

Genomic DNA was extracted from a cut piece of each cysticercus using Qiagen Blood and Tissue Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. PCR was performed in a final 25 µl reaction mixture containing 12.5 μl Premix Ex Taq™ version 2.0 (Takara Bio, Kusatsu, Japan), 10 pmol of each primer, 0.5 μl of genomic DNA extract (c.20–200 ng), and RNAse free water up to the final 25 μl volume. The reaction was as follows: initial denaturation at 95 °C for 5 min, followed by 35 cycles of denaturation at 95 °C for 30 s, annealing at 55 °C for 40 s, and elongation at 72 °C for 60–90 s and a final extension at 72 °C for 10 min. Amplicons were visualized by electrophoresis in 1.5% (w/v) agarose gels in 1× TAE (40 mM Tris-acetate, 2 mM EDTA, pH 8.5), stained with GelRed™, viewed under UV light and the products sequenced (Beijing Tsingke Biotechnology Co., Ltd., Beijing, China). Primers used for amplifying the complete cox1, nad1 and nad5 genes were designed by Primer Premier 5 software based on the full mitochondrial genome sequence of T. hydatigena (GenBank: NC_012896) (Table 1).

Table 1 Set of primers used for Taenia hydatigena cox1, nad1 and nad5 amplification (GenBank Taenia hydatigena reference: NC_012896)

Molecular analysis

Sequences were assembled stepwise with the help of DNAstar v7.1 program and Unipro UGENE v1.29.0 software while making sure that the overlap sequences were identical and then viewed manually for correction of any nucleotide misread followed by alignment in BioEdit v7.2.6 [34]. The identity of each isolate was confirmed with their nucleotide sequence in the GenBank database using the NCBI BLAST algorithm (https://blast.ncbi.nlm.nih.gov/Blast.cgi). Nucleotide and haplotype diversity indices were estimated in DnaSP v6 [35]. Median-joining network [36] was inferred based on the sequences of mitochondrial cox1, nad1, nad5, and concatenated cox1-nad1-nad5 genes using PopART (http://popart.otago.ac.nz). Pairwise nucleotide difference was calculated using MEGA 7 [37]. Population neutrality indices, Tajima’s D [38] and Fu’s Fs [39] were calculated using DnaSP v6 [35]. FST was calculated using the Arlequin 3.5.2.2 software package [40]. The Bayesian phylogenetic relationship of the Nigerian T. hydatigena isolates with other Taenia spp. was carried out based on cox1-nad1-nad5 (4083 bp) concatenated sequences with MrBayes v.3.1.2. The General Time Reversible model with a proportion of invariable sites and a gamma-shaped distribution of rates across sites (GTR + I + G) was used as the best-fit model of sequence evolution as determined by JModelTest [41]. Markov Chain Monte Carlo (MCMC) sampling was used to assess the posterior distribution of parameters with a chain length of 2,000,000 states, and 10% was discarded as ‘burn-in’. Parameters were logged every 1000 states. TreeView v.1.6.6. (http://taxonomy.zoology.gla.ac.uk/rod/treeview.html) was used to display the phylogenetic tree.

Results

All 32 isolates (24 goats and 8 sheep) were identified as T. hydatigena. Overall, 29, 24 and 27 isolates were successfully amplified for nad1, cox1 and nad5, respectively. A BLAST search of the resulting nucleotide sequences showed > 99% similarity with T. hydatigena sequences deposited in the GenBank database.

Sequence variation

On analysis of the complete cox1 (1620 bp), nad1 (894 bp) and nad5 (1569 bp) mitochondrial nucleotide sequences, we observed 24, 24 and 33 polymorphic sites, respectively for each gene, of which 62.5% (15/24), 75% (18/24) and 36.36% (12/33) were parsimony informative, respectively. The number of haplotypes (Hap) observed was 10, 10, 9 and 9 for nad1, cox1, nad5 and cox1-nad1-nad5 sequences, respectively (Table 2). According to the median-joining network, the haplotype aNIG1 appeared centralized for nad1 gene sequences and constituted 24.13% (7/29) of the total population with 1–7 mutational differences from the other haplotypes (Fig. 1a). For cox1 network, no central haplotype was observed while for nad5, bNIG9 was at the center of the network with 1–10 mutational steps from the other haplotypes constituting 22.22% (6/27) of the total population and made up mostly of sheep isolates (Fig. 1b). Both nad1 and nad5 central haplotype comprised isolates from both hosts but did not constitute the majority of the population (Fig. 1a, b). In addition, nad5 and cox1 network showed distinct haplotypes for sheep isolates (Fig. 1b, c). Interestingly, on analysis of the concatenated sequences of all three genes, all sheep isolates formed two distinct haplotypes (Fig. 1d).

Table 2 Diversity and neutrality indices for Taenia hydatigena populations from Nigeria based on cox1, nad1, nad5 and concatenated cox1-nad1-nad5 mitochondrial gene sequences
Fig. 1
figure1

Median-joining networks of Taenia hydatigena isolates from Nigeria based on complete cox1 (a) nad1 (b) nad5 (c) and concatenated cox1-nad1-nad5 gene sequences (d). Circle sizes are proportional to the haplotype frequencies. Numbers in parentheses represent the number of mutations. Small black circles are median vectors (i.e. hypothetical haplotypes)

The observed nucleotide polymorphism between haplotypes resulted in amino acid change. Comparison of haplotype 1 with other haplotypes of each respective gene showed the following changes: cox1 gene, a single amino acid substitution (cNIG6: 447A-447V); nad1, 6 amino acid substitutions (aNIG3: 10V-10G, 74G-74S, 86V-86I, and 253V-253G; aNIG4: 199F-199L, and 227C-227R) and nad5, 8 amino acid substitutions (bNIG2: 75I-75V, and 159L-159S; bNIG3: 385I-385V; bNIG4: 93I-93V and 385I-385V; bNIG5: 57L-57S; bNIG6: 353I-353V; bNIG7: 234I-234T) (see Additional file 1: Tables S1–S3).

Representative haplotype sequences of T. hydatigena nad5, nad1, and cox1 genes can be found in GenBank under accession numbers, MN175571-MN175579, MN175580-MN175589 and MN175590-MN175599, respectively.

Population indices

The pairwise nucleotide differences were highest for nad1 (0.1–2.5%), followed by nad5 (0.1–1.5%) and cox1 (0.1–0.7%). A high haplotype (Hd) and low nucleotide (π) diversity were observed for all three genes in both intermediate hosts (Table 2). Overall, Hd and π were as follows: cox1 (Hd = 0.906, π = 0.00264), nad1 (Hd = 0.867, π = 0.00628) and nad5 (Hd = 0.835, π = 0.00270). Taenia hydatigena population from sheep showed positive Tajima’s D for all genes and was significant for cox1. In contrast, isolates from goats showed negative insignificant Tajima’s D for all genes (Table 2). Fu’s Fs was positive for all genes except for cox1 isolates from goats (Table 2). Analysis of the concatenated gene sequences (cox1-nad1-nad5; 4083 bp), resolved most inconsistencies with an overall Hd and π of 0.900 and 0.00358, respectively, and a positively significant Tajima’s D (2.20374) and Fu’s Fs (5.778) for isolates of sheep origin, suggesting the possibility of a balance selection or sudden population contraction. In contrast, a negative (insignificant) Tajima’s D (− 0.85760) and a significantly positive Fu’s Fs (5.873) were observed for isolates from goats. Using the cox1, nad1 and nad5 genes, FST values for pairwise comparison between isolates from goats and sheep were as follows: cox1 (FST = 0.183, P = 0.009), nad1 (FST = 0.147, P = 0.009), nad5 (FST = 0.052, P = 0.117), and cox1-nad1-nad5 (FST = 0.149, P = 0.045) suggesting the possibility of genetically differentiated strains.

Phylogenetic analyses

Phylogenetic analysis showed low posterior probability (pp) values indicating weak nodal support (< 0.80) for the different haplotypes (Fig. 2). Also, haplotype NIG3 appeared as a sister taxon to other T. hydatigena cluster (the same haplotype according to the median-joining network also presented at least, a 15-point mutation difference from its closest neighbour, Fig. 1d). The cluster showed three sub-clades with varying pp values of which two Chinese reference isolates formed the first sub-clade with a pp value of 0.99 compared to other sub-clades (< 0.80) which were constituted by haplotypes from Nigeria (Fig. 2).

Fig. 2
figure2

Bayesian phylogenetic relationships of the Nigerian Taenia hydatigena isolates based on cox1-nad1-nad5 (4083 bp) concatenated gene sequences. Posterior probability (pp) values are depicted at the nodes. Echinococcus granulosus was used as an outgroup. Asterisks (*) indicate haplotypes representing isolates from this study, GenBank: MN175594, MN175582, MN175571 (NIG1) MN175595, MN175580, MN175572 (NIG2) MN175596, MN175584, MN175573 (NIG3) MN175592, MN175585, MN175575 (NIG4) MN175599, MN175587, MN175577 (NIG5) MN175598, MN175588, MN175576 (NIG6) MN175597, MN175589, MN175571 (NIG7) MN175590, MN175581, MN175578 (NIG8) MN175591, MN175580, MN175579 (NIG9)

Comparison of Nigerian T. hydatigena isolates with populations from other geographical regions

Due to the scarcity of full-length mitochondrial DNA (mtDNA) data from other geographical locations as in this study, we could not employ the full-length mitochondrial genes for comparison. Nonetheless, on trimming, a final dataset of 59 nad1 sequences of 435 bp (including 29 sequences from this study) (Additional file 1: Table S4) was used to compare the geographical relatedness of the Nigerian isolates and those from different regions. Results of the analysis showed 32 haplotypes, a high Hd of 0.9573, and a π of 0.00853. Overall neutrality indices showed significantly negative values (Tajima’s D= − 1.81180, Fu’s Fs = − 25.300). Furthermore, the median-joining network presented two major haplogroups (Fig. 3) that comprised haplotypes from different hosts and locations (Fig. 3). In addition, FST values of T. hydatigena population based on pairwise comparison showed significant differences only between isolates of Nigerian and European origin (FST= 0.116, P = 0.009) and European and Asian/Middle Eastern isolates (FST = 0.088, P = 0.05) (Additional file 1: Table S5).

Fig. 3
figure3

Median-joining network of Taenia hydatigena isolates from different hosts and geographical locations based on the NADH dehydrogenase subunit 1 (nad1) (435 bp) mitochondrial gene. Circle sizes are proportional to the haplotype frequencies. Numbers in parentheses represent the number of mutations. Small black circles are median vectors (i.e. hypothetical haplotypes)

Discussion

Despite the fact that cysticercosis due to T. hydatigena causes inconceivable damage in livestock production in endemic countries [1,2,3,4], the genetic population structure and the epidemiological significance of the observed genetic variation is poorly understood [4, 29]. In Nigeria, previous epidemiological surveys have shown that cysticercosis caused by T. hydatigena in sheep and goats could reach as high as 20–30% [10, 13] and below 2% in pigs [7], causing remarkable setbacks in livestock production [3]. Similarly, in other African countries, a number of epidemiological studies have also reported high prevalence of cysticercosis in ruminants [4, 42, 43].

Meanwhile, mtDNA has been widely employed in investigating the intraspecific variation of metazoans due to the absence of recombination, maternal inheritance, conserved structure, higher mutation rate and a relatively high evolutionary rate [44,45,46,47]. Another method for studying intraspecific variation or genetic diversity in cestodes is the use microsatellite DNA, which has been reported to be highly informative and also commonly used in genetic population studies [48,49,50]. However, in this study, we report the genetic diversity of T. hydatigena isolates from sheep and goats collected from a slaughterhouse in Nigeria based on complete sequences for three mitochondrial genes, cox1, nad1 and nad5, in contrast to the commonly employed partial gene sequences [4, 28, 51, 52].

Our analysis of the three mitochondrial genes revealed a considerable degree of genetic variation. Higher haplotype diversity was recorded for goat isolates and was highest for cox1 gene while nucleotide diversity was similarly higher in goats and highest for nad1 gene. Overall, the diversity and neutrality indices based on concatenation of all three genes were higher in goats and demonstrated inconsistency with population expansion which suggest the possibility of a bottleneck event in the course of evolution. However, the diversity indices of the individual genes cox1 and nad1 were comparable to estimates from Iranian and Italian (Sardinia) T. hydatigena populations from sheep and goats [28] but higher than Palestinian (West Bank) T. hydatigena populations of sheep origin [51]. In Palestine, lower prevalence and transmission rate of cysticercus tenuicollis (the larval stage of T. hydatigena), smaller area of Palestine, low small ruminant population (about 1.5 million) [53] as well as the prevailing management system have been suggested to influence the genetic diversity of T. hydatigena [28, 51]. In contrast, Iran and Sardinia (Italy) experience a higher prevalence and transmission rate of cysticercosis, higher population of small ruminants (Iran ≥ 50 million, Sardinia ≥ 3 million), with traditional farming methods still being applied in raising sheep and goats [28]. These differences could possibly have an influence on the genetic diversity of T. hydatigena, such that a higher prevalence in intermediate hosts, as seen in Italy and Iran, may result in multiple infections (cysts from different intermediate hosts) in definitive hosts favouring sexual reproduction and consequently higher genetic variation within the population.

Although the significant differences in the number of examined Nigerian T. hydatigena isolates (from both intermediate hosts) could have influenced the outcome of the neutrality and nucleotide diversity indices, the hypothetical existence of host-specific strains which has been reported cannot be completely ruled out [28, 54, 55]

Based on the individual genes, the median-joining network featured a centralized haplotype only for nad1 and nad5 genes but these did not constitute the majority of the population. On concatenation of the three genes, which is believed to be more informative than the individual genes, no centralized haplotype was found. However, two distinct haplotypes were revealed in all sheep isolates. Meanwhile, the formation of a founder haplotype (29.2% of the population) based on partial cox1 (324 bp) gene was reported in a pooled population of T. hydatigena from different countries and hosts [28]. Conversely, analysis of the nad1 gene in the latter study did not reveal a centralised haplotype. The same was true for Italian isolates from sheep, goats and wild boars [28]. In Palestine, the formation of a founder haplotype was also observed in sheep based on cox1 fragment (444 bp) and constituted about 55% of the examined population [51]. However, in a more recent investigation on the genetic diversity of T. hydatigena isolates from Turkish sheep and goats based on nad1 fragment (471 bp), no centralised haplotype was also observed [52]. Be that as it may, clarification of such discrepancies and improved phylogenetic resolution rest mostly on the analysis of longer or complete mtDNA sequences.

We observed high FST values between T. hydatigena populations from sheep and goats, suggesting genetic differentiation. Meanwhile, reports of host specificity have been previously documented between isolates from sheep and goats in Iran [54] and India [55], such that T. hydatigena metacestodes from both hosts were found to be morphologically different. Similarly, analysis of the biochemical components of cysticerci from pigs and goats also suggested genetic differences [56]. More so, pairwise comparisons in the present study between goats and sheep isolates yielded a significant FST value. This is in-line with previous FST values suggesting genetic differentiation between sheep and goat isolates from Italy, China and Greece, as well as between pig isolates from Italy, China and Poland [28], which further supports the idea of genetic distinction among T. hydatigena populations infecting different hosts. Then again, the observed pairwise fixation values and population indices could have been influenced by the small sample size analysed in the present study as sample size has been found to influence estimates of population indices [38].

Phylogenetic analysis of the concatenated gene sequences showed that all isolates from Nigeria were correctly identified as T. hydatigena as they clustered closely with other T. hydatigena isolates from other regions with relative distances from other Taenia spp. (Fig. 2). However, haplotype NIG3 which consists only of isolates from goats formed a separate clade (Fig. 2). The weak nodal support observed between Nigerian haplotypes is similar to previous observations of low posterior probability values (0.51–1.00) and bootstrap values (79%) between T. hydatigena haplotypes from Iran and Tanzania [4, 27].

As previously suggested [56], strain specificity in T. hydatigena, may present a similar feature like that of Echinococcus in which case, deep insight may be provided potentially on examination of longer mtDNA fragments and possibly nuclear genes of larger datasets comprising isolates from different hosts and geographical locations.

Furthermore, analysis of the trimmed nad1 gene sequence dataset (435 bp) from T. hydatigena isolates from other geographical regions and hosts revealed a population network that was consistent with population expansion based on the significant negative Tajima’s D and Fu’s F, indicating the presence of rare haplotypes as well as the characteristic star-shaped configuration with a centrally placed haplotype similar to previous observation by Boufana et al. [28]. In the network, some haplotypes were formed by isolates from different hosts and locations. Conversely, the distinct haplotypes of Nigerian sheep origin previously classified based on 4083 bp of mitochondrial DNA, formed a new haplotype that comprised isolates from different hosts (including Nigerian goat isolates) and locations, thus questioning the reliability of partial gene sequences in investigating the genetic population structure or resolving the phylogenetic status of T. hydatigena. Moreover, the pitfall of short gene sequences in inferring the phylogenetic status and the genetic diversity of cestodes has been emphasized for Echinococcus spp. [57]. Further, the outcome of longer DNA fragment analysed in the present study demonstrated that using either of the genes alone is insufficient to describe the intraspecific variation existing between T. hydatigena, as the 4083 bp of the mitochondrial genome provided a better resolution as indicated by the phylogenetic network and population indices.

The influence of geography or animal husbandry practices could not have plausibly explained the genetic differences observed between T. hydatigena isolates in this study. This is because small ruminants (sheep and goats) in Nigerian pastoral communities are mostly raised in the same environment and similar animal husbandry system. Moreover, each geopolitical zone of the country share similar geography (vegetation and climatic conditions), and in this study, the examined ruminants originated from the same state.

Conclusion

The present study, to the best of our knowledge, described the genetic variation of T. hydatigena in Nigeria for the first time using full-length mitochondrial gene sequences and suggests the possible existence of host-specific strains. The scarcity of complete T. hydatigena mitochondrial cox1, nad1 and nad5 sequences from different geographical regions in the GenBank repository limited a wide-range comparison and interpretation. However, the results suggest that longer DNA fragments or complete mitochondrial genome analysis may provide a better resolution in our understanding of the genetic diversity of T. hydatigena and the possible epidemiological significance.

Availability of data and materials

All data supporting the conclusions of this article are included in the article and its additional files. Representative nucleotide sequences of cox1, nad1 and nad5 genes from the present study are available in the GenBank database under the accession numbers MN175571-MN175599.

Abbreviations

PCR:

polymerase chain reaction

pp:

posterior probability

cox1:

cytochrome c oxidase subunit 1 gene

nad1:

NADH dehydrogenase subunit 1 gene

nad5:

NADH dehydrogenase subunit 5 gene

mtDNA:

mitochondrial DNA

MCMC:

Markov Chain Monte Carlo

Hd:

haplotype diversity

References

  1. 1.

    Scala A, Urrai G, Varcasia A, Nicolussi P, Mulas M, Goddi L, et al. Acute visceral cysticercosis by Taenia hydatigena in lambs and treatment with praziquantel. J Helminthol. 2016;90:113–6.

  2. 2.

    Christodoulopoulos G, Theodoropoulos G, Petrakos G. Epidemiological survey of cestode-larva disease in Greek sheep flocks. Vet Parasitol. 2008;153:368–73.

  3. 3.

    Saulawa M, Magaji A, Faleke O, Mohammed A, Kudi A, Musawa A, et al. Prevalence of Cysticercus tenuicollis cysts in sheep slaughtered at Sokoto abattoir, Sokoto state, Nigeria. Sokoto J Vet Sci. 2011;9:24–7.

  4. 4.

    Braae UC, Kabululu M, Nørmark ME, Nejsum P, Ngowi HA, Johansen MV. Taenia hydatigena cysticercosis in slaughtered pigs, goats, and sheep in Tanzania. Trop Anim Health Prod. 2015;47:1523–30.

  5. 5.

    Scala A, Pipia AP, Dore F, Sanna G, Tamponi C, Marrosu R, et al. Epidemiological updates and economic losses due to Taenia hydatigena in sheep from Sardinia, Italy. Parasitol Res. 2015;114:3137–43.

  6. 6.

    Nguyen MT, Gabriël S, Abatih EN, Dorny P. A systematic review on the global occurrence of Taenia hydatigena in pigs and cattle. Vet Parasitol. 2016;226:97–103.

  7. 7.

    Ajayi J, Arabs W, Adeleye G. Helminths and protozoa of pigs on the Jos Plateau, Nigeria: occurrence, age incidence and seasonal distribution. Bull Anim Health Prod Afr. 1988;36:47–54.

  8. 8.

    Dada B, Belino E. Prevalence of hydatidosis and cysticercosis in slaughtered livestock in Nigeria. Vet Rec. 1978;103:311–2.

  9. 9.

    Alonge DO, Fasanmi EF. A survey of abattoir data in northern Nigeria. Trop Anim Health Prod. 1979;11:57–62.

  10. 10.

    Akinboade OA, Ajiboye A. Studies on cysticercosis of small ruminants in Nigeria. Int J Zoonoses. 1983;10:164–6.

  11. 11.

    Opasina BA. Cysticercus tenuicollis of village sheep and goats in southwest Nigeria. Ann Trop Med Parasitol. 1985;79:657–8.

  12. 12.

    Okafor FC. Epizootiology of Cysticercus bovis in Imo state, Nigeria. Angew Parasitol. 1988;29:25–30.

  13. 13.

    Fakae B. The epidemiology of helminthosis in small ruminants under the traditional husbandry system in eastern Nigeria. Vet Res Com. 1990;14:381–91.

  14. 14.

    Nwosu C, Ogunrinade A, Fagbemi B. Prevalence and seasonal changes in the gastro-intestinal helminths of Nigerian goats. J Helminthol. 1996;70:329–33.

  15. 15.

    Bowles J, Blair D, McManus DP. Genetic variants within the genus Echinococcus identified by mitochondrial DNA sequencing. Mol Biochem Parasitol. 1992;54:165–73.

  16. 16.

    Kinkar L, Laurimäe T, Acosta-Jamett G, Andresiuk V, Balkaya I, Casulli A, et al. Global phylogeography and genetic diversity of the zoonotic tapeworm Echinococcus granulosus sensu stricto genotype G1. Int J Parasitol. 2018;48:729–42.

  17. 17.

    Nakao M, Li T, Han X, Ma X, Xiao N, Qiu J, et al. Genetic polymorphisms of Echinococcus tapeworms in China as determined by mitochondrial and nuclear DNA sequences. Int J Parasitol. 2010;40:379–85.

  18. 18.

    Laurimäe T, Kinkar L, Romig T, Omer RA, Casulli A, Umhang G, et al. The benefits of analysing complete mitochondrial genomes: deep insights into the phylogeny and population structure of Echinococcus granulosus sensu lato genotypes G6 and G7. Infect Genet Evol. 2018;64:85–94.

  19. 19.

    Ohiolei JA, Yan HB, Li L, Magaji AA, Luka J, Zhu GQ, et al. Cystic echinococcosis in Nigeria: first insight into the genotypes of Echinococcus granulosus in animals. Parasit Vectors. 2019;12:392.

  20. 20.

    Carmena D, Cardona GA. Echinococcosis in wild carnivorous species: epidemiology, genotypic diversity, and implications for veterinary public health. Vet Parasitol. 2014;202:69–94.

  21. 21.

    Ito A, Yamasaki H, Nakao M, Sako Y, Okamoto M, Sato MO, et al. Multiple genotypes of Taenia solium ramifications for diagnosis, treatment and control. Acta Trop. 2003;87:95–101.

  22. 22.

    Maravilla P, Gonzalez-Guzman R, Zuniga G, Peniche A, Dominguez-Alpizar JL, Reyes-Montes R, et al. Genetic polymorphism in Taenia solium cysticerci recovered from experimental infections in pigs. Infect Genet Evol. 2008;8:213–6.

  23. 23.

    Zhang Y, Zhao W, Yang D, Tian Y, Zhang W, Liu A. Genetic characterization of three mitochondrial gene sequences of goat/sheep-derived Coenurus cerebralis and Cysticercus tenuicollis isolates in Inner Mongolia, China. Parasite. 2018;25:1.

  24. 24.

    Rostami S, Salavati R, Beech RN, Babaei Z, Sharbatkhori M, Harandi MF. Genetic variability of Taenia saginata inferred from mitochondrial DNA sequences. Parasitol Res. 2015;114:1365–76.

  25. 25.

    Vega R, Pinero D, Ramanankandrasana B, Dumas M, Bouteille B, Fleury A, et al. Population genetic structure of Taenia solium from Madagascar and Mexico: implications for clinical profile diversity and immunological technology. Int J Parasitol. 2003;33:1479–85.

  26. 26.

    Kedra AH, Tkach VV, Swiderski Z, Pawłowski Z. Intraspecific variability among NADH dehydrogenase subunit 1 sequences of Taenia hydatigena. Parasitol Int. 2001;50:145–8.

  27. 27.

    Rostami S, Salavati R, Beech RN, Babaei Z, Sharbatkhori M, Baneshi MR, et al. Molecular and morphological characterization of the tapeworm Taenia hydatigena (Pallas, 1766) in sheep from Iran. J Helminthol. 2015;89:150–7.

  28. 28.

    Boufana B, Scala A, Lahmar S, Pointing S, Craig PS, Dessì G, et al. A preliminary investigation into the genetic variation and population structure of Taenia hydatigena from Sardinia, Italy. Vet Parasitol. 2015;214:67–74.

  29. 29.

    Omar MA, Elmajdoub LO, Al-Aboody MS, Elsify AM, Elkhtam AO, Hussien AA. Molecular characterization of Cysticercus tenuicollis of slaughtered livestock in Upper Egypt governorates. Asian Pac J Trop Biomed. 2016;6:706–8.

  30. 30.

    Kinkar L, Laurimäe T, Simsek S, Balkaya I, Casulli A, Manfredi MT, et al. High-resolution phylogeography of zoonotic tapeworm Echinococcus granulosus sensu stricto genotype G1 with an emphasis on its distribution in Turkey, Italy and Spain. Parasitology. 2016;143:1790–801.

  31. 31.

    Kinkar L, Laurimäe T, Acosta-Jamett G, Andresiuk V, Balkaya I, Casulli A, et al. Distinguishing Echinococcus granulosus sensu stricto genotypes G1 and G3 with confidence: a practical guide. Infect Genet Evol. 2018;64:178–84.

  32. 32.

    Jia WZ, Yan HB, Guo AJ, Zhu XQ, Wang YC, Shi WG, et al. Complete mitochondrial genomes of Taenia multiceps, T. hydatigena and T. pisiformis: additional molecular markers for a tapeworm genus of human and animal health significance. BMC Genomics. 2010;11:447.

  33. 33.

    National Agriculture Sample Survey, Nigeria. National Bureau of Statistics. Public Access Dataset; 2010. http://nigerianstat.gov.ng/download/66. Accessed 8 Aug 2019.

  34. 34.

    Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.

  35. 35.

    Rozas J, Ferrer-Mata A, Sanchez-DelBarrio JC, Guirao Librado P, Ramos-Onsins SE, Sanchez-Gracia A. DnaSP v6: DNA sequence polymorphism analysis of large datasets. Mol Biol Evol. 2017;34:3299–302.

  36. 36.

    Bandelt H, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.

  37. 37.

    Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

  38. 38.

    Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.

  39. 39.

    Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147:915–25.

  40. 40.

    Excoffier L, Laval G, Schneider S. Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005;1:47–50.

  41. 41.

    Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.

  42. 42.

    Samuel W, Zewde GG. Prevalence, risk factors, and distribution of Cysticercus tenuicollis in visceral organs of slaughtered sheep and goats in central Ethiopia. Trop Anim Health Prod. 2010;42:1049–51.

  43. 43.

    Wondimu A, Abera D, Hailu Y. A study on the prevalence, distribution and economic importance of Cysticercus tenuicollis in visceral organs of small ruminants slaughtered at an abattoir in Ethiopia. J Vet Med Anim Health. 2011;3:67–74.

  44. 44.

    Brown WM, George M Jr, Wilson AC. Rapid evolution of animal mitochondrial DNA. Proc Natl Acad Sci USA. 1979;76:1967–71.

  45. 45.

    Mueller RL, Macey JR, Jaekel M, Wake DB, Boore JL. Morphological homoplasy, life history evolution, and historical biogeography of plethodontid salamanders inferred from complete mitochondrial genomes. Proc Natl Acad Sci USA. 2004;101:13820–5.

  46. 46.

    Shen X, Wang H, Ren J, Tian M, Wang M. The mitochondrial genome of Euphausia superba (Prydz Bay) (Crustacea: Malacostraca: Euphausiacea) reveals a novel gene arrangement and potential molecular markers. Mol Biol Rep. 2010;37:771–84.

  47. 47.

    Wei SJ, Tang P, Zheng LH, Shi M, Chen XX. The complete mitochondrial genome of Evania appendigaster (Hymenoptera: Evaniidae) has low A+T content and a long intergenic spacer between atp8 and atp6. Mol Biol Rep. 2010;37:1931–42.

  48. 48.

    Bart JM, Knapp J, Gottstein B, El-Garch F, Giraudoux P, Glowatzki ML, et al. EmsB, a tandem repeated multi-loci microsatellite, new tool to investigate the genetic diversity of Echinococcus multilocularis. Infect Genet Evol. 2006;6:390–400.

  49. 49.

    Maillard S, Gottstein B, Haag KL, Ma S, Colovic I, Benchikh-Elfegoun MC, et al. The EmsB tandemly repeated multilocus microsatellite: a new tool to investigate genetic diversity of Echinococcus granulosus sensu lato. J Clin Microbiol. 2009;47:3608–16.

  50. 50.

    Umhang G, Grenouillet F, Bastid V, M’Rad S, Valot B, Oudni-M’Rad M, et al. Investigating the genetic diversity of Echinococcus granulosus sensu stricto with new microsatellites. Parasitol Res. 2018;117:2743–55.

  51. 51.

    Adwan K, Jayousi A, Abuseir S, Abbasi I, Adwan G, Jarrar N. Genetic diversity of Taenia hydatigena in the northern part of the West Bank, Palestine as determined by mitochondrial DNA sequences. Acta Parasitol. 2018;63:299–303.

  52. 52.

    Kilinc SG, Kesik HK, Simsek S. Molecular characterization and haplotypes of sheep and goat isolates of Cysticercus tenuicollis in Turkey. Parasitol. 2019;146:1047–54.

  53. 53.

    Palestinian Central Bureau of Statistics. Livestock Statistics. 2017. http://pcbs.gov.ps/site/lang__en/734/default.aspx. Accessed 8 Aug 2019.

  54. 54.

    Radfar M, Tajalli HS, Jalalzadeh M. Prevalence and morphological characterization of Cysticercus tenuicollis (Taenia hydatigena cysticerci) from sheep and goats in Iran. Vet Arch. 2005;75:469–76.

  55. 55.

    Singh BB, Sharma R, Gill JPS, Sharma JK. Prevalence and morphological characterisation of Cysticercus tenuicollis (Taenia hydatigena cysts) in sheep and goat from north India. J Parasit Dis. 2015;39:80–4.

  56. 56.

    Abidi SMA, Nazami WA, Khan P, Ahmad M, Irshadullah M. Biochemical characterization of Taenia hydatigena cysticerci from goats and pigs. J Helminthol. 1989;63:333–7.

  57. 57.

    Laurimäe T, Kinkar L, Moks E, Romig T, Omer RA, Casulli A, et al. Molecular phylogeny based on six nuclear genes suggests that Echinococcus granulosus sensu lato genotypes G6/G7 and G8/G10 can be regarded as two distinct species. Parasitology. 2018;145:1929–37.

Download references

Acknowledgments

We would like to thank Dr. Etinosa O. Igbinosa, Head of Applied Microbial Processes & Environmental Health Research Group, University of Benin, Nigeria for allowing us to access his laboratory facilities.

Funding

This study was supported by the Central Public-interest Scientific Institution Basal Research Fund (1610312017001; 1610312016012), the National Key Research and Development Plan (2017YFD0501301; 2018YFC1602504) and the National Key Basic Research Programme (973 Program) of China (2015CB150300). The PhD programme of JAO was funded by the Chinese Government Scholarship Programme, the People’s Republic of China.

Author information

JAO, HBY and WZJ conceived and designed the experiments. JAO, JL and AAM collected the samples. JAO, GQZ, LL, MAA, YTW, HBY and JQL performed the experiments and data analyses. JAO wrote the first draft of the manuscript, and HBY, BQF and WZJ provided comments and constructive suggestions for revisions. All authors read and approved the final manuscript.

Correspondence to Wan-Zhong Jia.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ohiolei, J.A., Luka, J., Zhu, G. et al. First molecular description, phylogeny and genetic variation of Taenia hydatigena from Nigerian sheep and goats based on three mitochondrial genes. Parasites Vectors 12, 520 (2019) doi:10.1186/s13071-019-3780-5

Download citation

Keywords

  • Genetic variation
  • Haplotype
  • Taenia hydatigena
  • Sheep
  • Goat
  • Cysticercosis

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.