INTRODUCTION
Cystic echinococcosis (CE) is a parasitic zoonosis caused by Echinococcus granulosus sensu lato (s.l.) [1]. The recent taxonomic classification of E. granulosus s.l. based on genetic information from both mitochondrial and nuclear genes proposed 5 species have been demarcated from the 10 genotypes (G1 to G10): E. granulosus sensu stricto (s.s.) (genotypes G1–G3), E. equinus (genotype G4), E. ortleppi (genotype G5), E. canadensis (genotypes G6–G10), and E. felidis [2,3]. In addition, E. granulosus s.l. has intraspecific microdiversity within and between sub-populations that constructed haplotype networks associated with the hosts and geographic distribution [4,5]. The variations in the genotypes of E. granulosus s.l. may affect the location, fertility of cysts, the severity of CE infection in intermediate hosts, and response to chemotherapy [1]. E. granulosus s.s. shows a large geographical distribution and host range [6,7].
It was reported that CE is endemic in Tanzania [8] and prevalent in the intermediate hosts [9] and human [10,11]. A retrospective survey in Ngorongoro district (one of the Maasai areas) in the northern Tanzania reported CE prevalence in humans, 10 cases per 100,000 people per year [12]. In the same district, incidence of CE was estimated as 63.8% in sheep, 34.7% in goats, and 48.7% in cattle [9]. The data on the prevalence and genetic diversity of Echinococcus spp. in the African countries including Tanzania is scarce. This situation calls for an investigation on the genetic features of CE in Tanzania. Therefore, this study was, conducted to clarify the genotypes of Echinococcus spp. from cattle, an intermediate host in Tanzania.
MATERIALS AND METHODS
Study area and sample collection
The cysts were collected between December 2018 and February 2019 from the lung of indigenous cattle breed (Tanzania Shorthorn Zebu) in slaughterhouses located in 2 regions; one in the city area (Mwanza, northwestern Tanzania) (2°31′ 0″S, 32°54′ 0″E), and the other in a far remote area (Loliondo, northern Tanzania) (2° 3′ 0″S, 35° 37′ 0″E) a specialized territory of Maasai people (Fig. 1).
A single hydatid cyst was collected from each carcass of cattle during the general postmortem inspection in the slaughterhouses. The collected cysts were individually kept in 70% ethanol in a plastic lab container and transferred to a laboratory. The germinal layer and hydatid liquid were separated for further study. A total of 7 cysts from the lung were sampled, 4 in Mwanza and 3 in Loliondo. The cysts collected in Mwanza were coded as EG_TZ01, EG_TZ03, EG_TZ06, and EG_TZ07, while those in Loliondo as EG_TZ10, EG_TZ12, and EG_TZ13.
Morphological analysis on the cysts
All cysts were assessed for protoscoleces and rostellar hooks in the germinal layer, and hydatid fluid using a dissecting microscope [13]. The protoscoleces found in the fertile cysts (Fig. 2) were isolated and stained in Hoyer’s medium [14]. The rostellar hooks of E. granulosus s.l. protoscoleces were analyzed for 5 variables: number per rostellum, blade and total lengthes according to the protocol of Ponce Gordo and Cuesta Bandera [15]. The obtained measurements were compared with previous studies from Iran [16], Jordan [17], Poland [18], and Spain [15] (Table 1).
DNA extraction and PCR amplification
The CE DNA was extracted using the DNeasy Blood and Tissue Kit (QIAGEN®) according to the manufacturer’s protocol with some modifications. A small (<2 mm) piece of germinal layer kept in 70% ethanol was washed within PBS using a shaker overnight, and collected by centrifuging at 8,000 rpm. The samples were crushed with a pestle and motor, genomic DNA was extracted and kept at −20°C until use.
Polymerase chain reaction (PCR) was carried out targeting the mitochondrial cytochrome C oxidase 1 (cox1) gene. A primer pair JB3 (5′-TTT TTT GGG CAT CCT GAG GTT TAT-3′) and JB4.5 (5′-TAA AGA AAG AAC ATA ATG AAA ATG-3′) was used to amplify a 408 bp fragment of the cox1 gene [19]. The PCR reaction mix was in 40 μl volume, containing 8 μl HiPi Plus 5× PCR Master Mix, 12.5 pmol each forward and reverse primers, 26 μl distilled water, and 4–20 ng CE DNA. The PCR reaction was carried out using an automatic thermal cycler (SuperCycler SC-200, Kyratec, Mansfield/Queensland, Australia), whereby a pre-denaturation was set for 3 min at 95°C. Thirty-five cycles of DNA denaturation for 30 sec at 95°C, annealing for 30 sec at 47°C, and extention for 1 min at 72°C, with a final extension at 72°C for 10 min. The PCR products were run on 1% agarose gel, and the images were taken on a gel-imaging device (Gel DocTM XR+System, BIO-RAD, Hercules, California, USA).
Genotype and haplotype analyses
The amplicons were sequenced by a biomolecular company (Cosmogenetech Co, Daejeon, KOREA),and aligned using Geneious software Version 9.0 [20]. The sequenced samples were trimmed and assembled using De novo sequence assemblers, and subjected to homolog search In NCBI (http://blast.ncbi.nlm.nih.gov/BLAST.cgi).
The obtained sequences were compared with the reference sequences of E. granulosus s.s. and E. ortleppi from the GenBank to observe the phylogenetic relationships between individuals. Taenia asiatica were added as an outgroup [21]. MEGA v.6 software (https://www.megasoftware.net/resources) [22] and maximum likelihood (ML) method [23] were used to root the phylogenetic tree. The Bayesian information criterion (BIC) value was determined with the HKY +I (Hasegawa-Kishino-Yano+invariant sites) method, with 1,000 bootstrap replications estimated to get high confidence. According to mutation steps, the haplotype networks were generated using the NETWORK software that depends on statistical parsimony after creating the data file for the cox1 genetic locus. Finally, DNA Sequence Polymorphism analyses were performed by utilizing the DnaSP v.6 software to estimate several measures of DNA sequence variation within and between populations [24].
RESULTS
Fertility of cysts and morphology of protoscolex
All cysts (n=4) in the Mwanza were confirmed fertility, while those (n=3) in Loliondo were found sterile. The rostellar hooks (28–36) were lined in 2 rows. Total length of large hooks was 24.5 μm (22.2–29.5) and small hooks was 17.2 μm (14.8–19.7). The blade length of large and small hooks were 17.2 μm (14.8–19.7) and 8.1 μm (7.4–9.8), respectively (Table 1).
Echinococcus genotype and sequence polymorphism
The 7 DNA samples were amplified by PCR and yielded 408 bp fragments as partial cox1. The resulting gene sequences were translated into amino acids under the flatworm mitochondrial genetic code conditions and uploaded to the GenBank with accession numbers (MN540096-MN540101 and MW729426). The DNAs of EG_TZ03, EG_TZ06, and EG_TZ07 (MN540101, MN540100, and MN540099) obtained from Mwanza showed 100% sequence similarity to E. ortleppi genotype (G5) (MH428013, JX854035, and MK492625) from India and Egypt [25]. However, EG_TZ01 (MW729426) showed 99.4% similarity, one substitution (373 G/A), to E. ortleppi genotype (G5) from India and Egypt [25]. This single nucleotide substitution results in a nonsynonymous change at position 126 with Alanine to Threonine.
The EG_TZ10 and, EG_TZ13 from Loliondo showed 100% similarity with E. granulosus s.s. (G1) sequences) from Australia [26], and (MH542383) from Iran [27]. The sequence of EG_TZ12 from Loliondo showed 99.4% similarity with those of EG_TZ10 and EG_TZ13, with single nucleotide substitution (135 C/T). At the same time, the sequence of EG_TZ12 showed 100% similarity with E. granulosus s.s. (G1) sequences from Mongolia, China [28], and Jordan [6] (MK370108, AB491414, and JQ250806).
Echinococcus ortleppi and E. granulosus s.s. phylogeny
A phylogenetic tree drawn with the collected samples revealed 2 main clades: one with E. granulosus s.s. (G1), another with E. ortleppi (G5) (Fig. 3), distinguished from the previously reported Echinococcus genotypes [29–34] (Supplementary Table S1). In Loliondo, all isolated samples were identified to genotype (G1), and sequences were divided into 2 haplotypes (marked as EgTZ01 and EgTZ02). One sample (EG_TZ12) grouped with the haplotypes from Iran, China, Jordan, Peru, and Mongolia [6], the other 2 samples (EG_TZ10 and EG_TZ13) grouped with one haplotype of Iran, Turkey, Armenia, and Australia. [26].
In Mwanza, 2 haplotypes were found in the 4 collected samples under the genotype (G5). A new haplotype (EoTZ02) was found from one sample (EG_TZ01), and another haplotype (EoTZ01) was found from the other 3 samples. The haplotypes from Mwanza were grouped with all reported sequences. As a result, 10 haplotypes (Hap 1, Hap 2, Eo04, Eo06, Eo07, Eo09, Eo14, Eo15, EoTZ02, and EoTZ02) were inferred from the E. ortleppi isolates’ sequences, which included 10 polymorphic sites. The haplotype network had an orbit-shaped profile with a dominant haplotype in the center similar to previously reported common haplotypes in Africa, Asia, Europe, and South America. However, both of the haplotypes from Mwanza belonged to the minor haplotypes. A haplotype network of E. ortleppi (G5) was illustrated (Fig. 4) to understand the phylogenetic relationships between the haplotypes from this study and other Sub-Saharan countries [35] (Supplementary Table S1).
DISCUSSION
In the present study, the CE was only found in the cattle of the exanimated carcasses in Mwanza and Loliondo. However, numerous studies reported that livestock such as sheep, goats [36], and pigs [37] are at risk of CE infection in Tanzania. The small ruminants play an essential role in the life cycle of Echinococcus spp., compared to other animals [36].
In Africa, particularly in Algeria, Ethiopia, Morocco, and Tunisia, the dominant genotype is (G1); it infects the intermediate host (sheep, camel, cattle, and goat), humans, and the definitive host (dogs) [38,39]. However, other studies revealed that the E. canadensis (G6) was the dominant genotype in Sudan [39], Egypt [40], Mauritania [41], and Nigeria [42].
The present study found that the E. granulosus s.s. (G1–G3) in cattle from Loliondo has a similar predominant haplotype (G1) to the one previously described from Australia [26], China, Iran, Jordan, Mongolia, and Peru [6]. The similarity of the obtained genotype E. granulosus s.s. is related to its population expansions during the anthropogenic movement of intermediate and definitive hosts [28]. It was mentioned that the origin of E. granulosus s.s. was from Middle East [6]. Hence, it can be assumed that E. granulosus s.s. in Tanzania comes from the Middle East via the migration of the intermediate and definitive hosts.
The present study obtained the main genotype E. ortleppi (G5) infecting cattle only in Mwanza, Tanzania. In Africa, E. ortleppi CEs were reported from the livestock (cattle, goat, camel and, pig), wildlife (oryx) [35], and human [43]. E. ortleppi (G5) was the predominant genotype that causes echinococcosis in Namibia and Zambia [39]. In contrast, in other African countries, the proportion of infected animals with E. ortleppi was low [39,43].
Although it is interesting to note that different genotypes were found in small numbers of slaughtered cattle in 2 distant regions, it is somewhat difficult to guess the dominant genotype from these results. In addition, we predicted the existence of other species such as E. felidis in Loliondo district because this land uses wildlife-controlled areas inhabited predominantly by Maasai, who practice traditional pastoralism. Also, the disease was reported in wildlife in neighbor countries (Uganda) [9]. However, to our knowledge, this is the first molecular report of E. ortleppi and E. granulosus s.s. infection in the intermediate host from Tanzania. Both species of Echinococcus reported in this study are pathogens to humans with a significant public health concern [1,9], which sets the basis to explore further whether E. ortleppi and E. granulosus s.s. are also responsible for human infection in Tanzania.
The haplotypes in this study differed in low intraspecific diversity compared with the predominant haplotypes of E. ortleppi genotype reported from the neighboring countries (Kenya and Zambia) and Brazil, France, and Namibia [35]. The low intraspecific diversity was also reported [35] compared to 178 cases of E. ortleppi in Sub-Saharan Africa, Europe, and South America.
The fertility of cysts is one of the essential factors for the spread of the parasites in the environment and continuity of the life cycle of Echinococcus species. The cysts of the E. ortleppi in cattle are usually fertile [35] as obtained in Mwanza.
The morphological characteristic of protoscoleces, the number and length of rostellar hooks, have been used to distinguish species and strains of Echinococcus [44]. The rostellar hooks in Tanzanian samples had a similar size to that of Spanish ones [15] and smaller than that of Jordanes [44] and Iranian [16] (Table 1). However, with the morphology characteristic of protoscoleces only, we cannot find a way to distinguish the species of E. granulosus as did Ihsan et al. [17]. The identification of Echinococcus species should not be based only on the morphology of rostellar hooks of protoscoleces.
In conclusion, our findings obtained that cattle may play a significant part in the life cycle of CE, highlighting the possibility of potential transmission risks to the human population and other intermediate hosts in Tanzania. More importantly, we obtained the molecular data E. granulosus s.s. and E. ortleppi for the first time in Tanzania; additional epidemiological research on the prevalence of this species and its involvement in humans, stray dogs, and wildlife to understand CE in Tanzania better are required.