Skip to main content

Widespread co-endemicity of Trypanosoma species infecting cattle in the Sudano-Sahelian and Guinea Savannah zones of Cameroon

Abstract

Background

African animal trypanosomosis remains the major constraint of livestock production and livelihood of pastoral communities in Cameroon. Despite several decades of vector and parasite control efforts, it has not been eradicated. Alternative and sustainable control strategies require a sound knowledge of the local species, strains and vectors. In the Sudano-Sahelian and Guinea Savannah of Cameroon the prevalence and genetic diversity of trypanosomes infecting cattle was investigated by microscopy of cattle blood buffy coat and molecular methods using generic primers targeting parts of the internal transcribed spacer 1 (ITS-1) and encoded glycosomal glyceraldehyde 3-phosphate dehydrogenase-gene (gGAPDH).

Results

A total of 1176 randomly chosen cattle from five divisions in the Sudano-Sahelian and Guinea Savannah of Cameroon were examined. The overall prevalence of trypanosomes by microscopy was 5.9% (56/953) in contrast to 53.2% (626/1176) when molecular tools were used. This indicated a limited sensitivity of microscopy in subclinical infections with frequently low parasitemia. Three trypanosome species were identified by light microscopy: T. vivax (2.3%), T. brucei (3.7%) and T. congolense (3.0%), whereas five were identified by PCR, namely T. grayi/T. theileri (30.8%), T. vivax (17.7%), T. brucei (14.5%) and T. congolense (5.1%). Unexpected cases of T. grayi (n = 4) and T. theileri (n = 26) were confirmed by sequencing. Phylogenetic analysis of the gGAPDH revealed the presence of T. vivax, clade A and T. vivax clade C, which were co-endemic in the Faro et Deo division.

T. grayi/T. theileri were the predominant species infecting cattle in tsetse free areas. In contrast, T. vivax, T. brucei and T. congolense were more abundant in areas where the Glossina-vectors were present.

Conclusions

The abundance of pathogenic trypanosomes in tsetse infested areas is alarming and even more, the occurrence of T. vivax, T. brucei, T. congolense, T. theileri and T. grayi in tsetse-free areas implies that tsetse control alone is not sufficient to control trypanosomosis in livestock. To implement control measures that reduce the risk of spread in tsetse free areas, close monitoring using molecular tools and a thorough search for alternative vectors of trypanosomes is recommended.

Background

In tropical Africa and South America, hemoparasitic flagellates of the genus Trypanosoma cause severe to fatal diseases in wild and domestic mammals, including the human host. Trypanosomes infecting mammals are divided into two major families: Salivaria and Stercoraria [1]. Members of the Salivaria include human and veterinary medically important pathogens Trypanosoma vivax, T. congolense and T. brucei spp. They develop as mammalian infective forms in the mouthparts, e.g. proboscis and salivary glands of tsetse (Glossina spp.). Transmission to the vertebrate host occurs during the blood meal of an infective tsetse [2]. In contrast, the Stercoraria comprise the South American parasite T. cruzi and the worldwide-distributed Megatrypanum, e.g. T. theileri, where the final stages of the parasite develop in the posterior digestive tract of the arthropod vectors. These species are transmitted by contamination of the bite puncture with infectious excreta from the vector [3]. Trypanosomes can also be transmitted by mechanical vectors, like tabanid and stomoxine horse flies and by hard-ticks [4, 5].

In Cameroon, 90% of the population of the estimated six million cattle are at risk of trypanosome infection [6]. The Adamawa highland plateau in North Cameroon is the country’s main area of cattle rearing supplying animal products to all neighboring countries. This was made possible through the control of tsetse on this up to 1000 m high plateau [7], whilst Glossina morsitans, G. fuscipes fuscipes and G. tachinoides still occur in high numbers in the savannah pastures of the Eastern and Northern regions making cattle rearing problematic [7]. However, conventional operations employed over the last three decades have not eradicated the Glossina vectors so that pasture lands previously cleared and declared free of Glossina have recently been re-invaded [7, 8]. Disease control in these areas depended on continuing diagnosis and treatment of suspected cases with the few trypanocidal drugs available on the market [9]. Isometamidium, diminazene and homidium bromide are the only drugs widely used during more than four decades for trypanosome control. Furthermore, there are reports of drug resistance coming from North Cameroon [10] and elsewhere [11, 12]. Therefore, the unequivocal identification of the prevailing trypanosome species and strains has gotten more attention to prevent unnecessary treatment of non-pathogenic parasites and thereby promoting the development of resistance.

In Northern Cameroon, little is yet known about the genetic diversity of trypanosomes infecting cattle. Most epizootiological data available were based on microscopy, such as phase-contrast or dark field examination of the buffy coat, thin or thick blood smears, and to a lesser extent also serological analyses [13,14,15]. These investigations indicated T. congolense, T. brucei and T. vivax as the only prevalent species in these areas [7,8,9,10]. Microscopy, albeit easy to perform in a fieldwork setting, needs a high investment in time and training, risks to misinterpret rare, emerging or in other ways unexpected specimens and fails to detect immature infections during the first stages of infection [16]. Advances in molecular biology have expanded the limits of the traditional methods in sensitivity and specificity. Generic and specific primers have been designed to amplify the internal transcribed spacer 1 (ITS-1) region of the ribosomal RNA gene locus of trypanosomes, chosen because of its high copy number and inter-species length variation [17,18,19]. Thus, trypanosome species are recognizable by the fragment length of their PCR-amplified ITS-1 region [17]. This method has evolved to improve sensitivity and detection of trypanosomes in animal blood [18,19,20]. In addition, the glycosomal glyceraldehyde 3-phosphate dehydrogenase gene (gGAPDH), an ubiquitous and essential glycolytic enzyme, has been used for the species differentiation of trypanosomes because of its lower rate of molecular evolution [21]. Despite the fact that it has no band size separation among different Trypanosoma species, it has been a marker of choice for phylogenetic analysis [22, 23].

A recent study in two restricted areas in Northern Cameroon relying on molecular tools for parasite detection [24] revealed active foci of AAT on the Adamawa region in the Faro et Deo close to the border with Nigeria and in the North region near the town of Gamba. The results revealed the crucial need of molecular tools to monitor the diversity of trypanosomes together with their vectors in hyper-endemic foci. A higher diversity of trypanosomes was seen in cattle and tsetse vectors than previously known. Those observations were however based only on a few Glossina-infested localities and on less than four hundred cattle examined. Therefore, this study has investigated the epizootiological picture of bovine trypanosomosis in the northern regions of Cameroon comparing tsetse infested areas in the high Guinea savannah and the Sudano-Sahelian zone with areas cleared of tsetse in both agro-ecological zones. Furthermore, the different susceptibilities of the various indigenous cattle breeds found in these zones have been addressed.

Gudali, White Fulani and Red Fulani are the major local zebu cattle breeds [25]. They are claimed to be more susceptible to trypanosomiasis than the autochthonous taurine cattle breed called Namchi (Doayo) [26], which nevertheless is at high risk of becoming extinct [27]. Only few located herds of Doayo cattle remain in the Faro division. The Kapsiki, another taurine cattle breed, with a higher introgression of Zebu genes, found mainly in the Mayo Tsanaga (Rhumsiki) area of the Far North region and also being on the verge of becoming extinct, were earlier shown to be trypanosusceptible [26].

The present research used both microscopy and molecular methods to study the occurrence and genetic diversity of trypanosomes in cattle from two agro-ecological zones (AEZ), focusing on areas with and without tsetse vectors [28].

Results

Body condition and packed cell volume in relation to breed and age

A total of 1176 animals were randomly sampled. These comprised more female (907; 77.1%) than male (269; 22.9%). Examined animals were from five divisions as follows: Vina (n = 283), Faro et Deo (n = 196), Mayo Rey (n = 316), Faro (n = 176) and Mayo-Tsanaga (n = 205). In the Faro and Mayo-Tsanaga divisions, only the indigenous taurine breeds, Namchi (Doayo) and Kapsiki, respectively, were examined. Here, the mean PCV of Namchi (Doayo) was significantly higher (F = 13.88; P < 0.001) than that of Kapsiki (Fig. 1a, Additional file 3: Table S1). Overall, animals with poor body condition had PCVs (average 29.66 ± 6.68) significantly lower (F = 22.062, P < 0.001) than that of animals in good (32.82 ± 4.99) and very good (34.26 ± 5.46) condition (Fig. 1b). Young cattle aged between 0 to 2.5 years had significantly lower PCVs (31.22 ± 6.82) than the other age groups (F = 5.38, P = 0.005, Fig. 1c, Additional file 3: Table S2). 97.6% of the Kapsiki cattle had the best body condition score (4 or higher, Additional file 3: Table S1) as compared to those of other cattle breeds. Comparing the different regions, animals in the Faro et Deo division had a mean PCV (28,13 ± 5.76) that was significantly lower (F = 49.13, P < 0.001) than those found in the Faro division (34.74 ± 5.35; Additional file 3: Table S2).

Fig. 1
figure1

Effect of cattle breed on packed cell volume (a). Comparison of the mean of PCV of five indigenous cattle breeds examined. Effect of body condition score on packed cell volume (b). Animals were grouped as described in the section “Materials and Methods” without breed distinction and the PCVs were compared. Effect of age group on body condition score (c). Animal were grouped by age as described in in the section “Materials and Methods” and PCV was compared. Details of sample collections and processing are indicated in the section “Materials and Methods”

Parasitological and molecular detection of trypanosomes

Microscopic detection of motile trypanosomes showed that 56 blood samples (5.9% of 971 cattle) carried at least one trypanosome species (Table 1). The highest prevalence was recorded in Faro et Deo (15.8%), followed by Faro (4.5%), Mayo- Rey (3.5%) and Vina (2.3%). In the Mayo-Tsanaga region no microscopy was carried out due to insecure work environment.

Table 1 Distribution of trypanosome species detected by microscopy in the study area

The most frequently identified trypanosome species was T. brucei spp., followed by T. vivax and T. congolense (Table 1). However, 7.1% of trypanosomes were not clearly identified according to their motility and morphological characteristics and were recorded as T. brucei-like trypanosome species.

In contrast, out of 1176 samples examined by ITS-1 nested PCR, 626 samples showed the presence of one or more trypanosome species, giving an overall prevalence of 53.2% (Table 2). The highest prevalence was recorded in Mayo-Tsanaga (67.8%), followed by Faro et Deo (59.2%) and lowest in Faro (34.1%). From the 56 samples classed positive by microscopy, 41 (71.9%) were also detected by nested PCR (Table 3).

Table 2 Distribution of trypanosome species detected by ITS-1 PCR in the study areas
Table 3 Comparison of the diagnostic test results obtained by parasite microscopy and molecular (ITS-1 PCR) methods

ITS-1 sequences analysis

Samples were identified according to ITS-1 amplicon size as described previously [19, 24] (Table 4). Three representative samples with a product size of 426 bp considered to be T. brucei spp. were sequenced and the results aligned to sequences retrieved from databases searches. The results showed that all sequences belonged to Trypanozoon, either to T. brucei spp. or T. evansi. They differ only in their maxi-circles DNA and additional species-specific markers are needed to distinguish these species. Additionally, six PCR amplicons in the range of 645 bp and considered to be T. congolense savannah or forest types were sequenced and showed a similarity of 73 to 85% with T. congolense isolates from South Africa and Gabon, respectively [GenBank: KX870079, KX452163].

Table 4 Trypanosome ITS-1 amplicon sizes of different Trypanosoma spp.

Interestingly, the PCR products of 180 bp and 250 bp (n = 6) both corresponded to T. vivax sequences isolated from Ethiopia ([GenBank: KM391818, KM391825], 91 to 93% identical). For PCR products in the range of 320 bp, out of 30 samples analyzed, 26 (87%) corresponded to T. theileri sequences published in Genbank (98 to 100% identical). The other four sequences (13%) matched with entries of T. grayi (90 to 96% identical) with closest similarity to T. grayi ANR4 isolated from a tsetse in The Gambia [TriTrypDB: JMRU01000589] and 94 to 99% identical with sequences [NCBI Blastn: MG255201, MG255205] obtained from cattle and tsetse in North Cameroon in Gamba and Kontcha, respectively [24].

Genetic diversity of trypanosome species

In total, five different trypanosomes were identified: T. congolense, T. brucei spp., T. theileri, T. grayi and T. vivax, respectively, using ITS-1 makers and sequencing analysis (Table 4). Due to the inability to discriminate between T. theileri and T. grayi just on the basis of the ITS1 amplicon size, samples with amplicons in the range of 320 bp were considered as T. theileri/T. grayi. Molecular analysis showed these to be the most prevalent species in all five study areas (30.8%, n = 362/1176). T. theileri/T. grayi was also the species most often missed or misidentified for T. brucei or T. congolense by microscopic observation, followed by T. vivax (Tables 1, 2 and 3). The overall prevalence of mixed infections was 11.4% (n = 134/1176). Co-infections of T. brucei spp. and T. vivax were the most common (n = 91/1176), followed by triple infections with T. congolense, T. vivax and T. theileri/T. grayi (n = 23/1176). We found eight animals co-infected by T. brucei spp., T. vivax and T. congolense savannah/forest-type, and four animals co-infected by T. brucei spp., T. vivax, T. congolense savannah/forest and T. theileri/T. grayi (Table 2).

The effect of study site, breed and age on the prevalence of trypanosomosis and correlation with the body condition score

Doayo (Namchi) cattle from Faro were significantly less infected (34.6%; X2 = 51.78, p < 0.000) with any trypanosome species than the other taurine cattle Kapsiki (67.8%) and the Zebu breeds Gudali (54.1%), Red Fulani (58.1%) and White Fulani (54.1%). There was also a significant difference between the five sampled study sites. The overall trypanosome infection rate was higher in Mayo-Tsanaga (67.8%) than in Faro et Deo (59.2%). However, 56.2% of the infected animals in Mayo-Tsanaga were infected with T. theileri/T. grayi, compared to only 7.5% in Faro et Deo. In contrast, when looking only at the species classically considered to be pathogenic such as T. congolense, T. brucei spp. and T. vivax, these were most prevalent in Faro et Deo (44.9%), followed by Mayo-Rey (42.7%), Faro (15.3%) and Mayo-Tsanaga (11.7%). The area with lowest prevalence was Vina (7.8%) (Table 5), a former tsetse-cleared area.

Table 5 Effect of age, breed, study areas and body condition score on trypanosome prevalence

Comparison of areas with or without Glossina-vectors

The overall prevalence of trypanosomes was similar or even higher in the tsetse free areas (Vina 53.7% and Mayo Tsanaga 67.8%, Tables 2 and 5) than in the Glossina-infested zones (Mayo-Rey 50.3%, Faro et Deo 59.2% and Faro 34.1%). T. theileri/T. grayi were the most abundant trypanosome species in the tsetse-free zones. In contrast, in the Glossina–infested areas T. vivax, T. brucei and T. congolense were the predominant species (Table 2).

Some T. congolense, T. brucei and T. vivax cases were even detected in the areas of Vina and Mayo-Tsanaga, although these areas have been declared tsetse-free (Table 2).

Phylogenetic analysis of gGAPDH

Two main clusters were observed in the 37 gGAPDH sequences examined, comprising the stercorarian T. grayi and T. theileri on the one hand, and the salivarian T. congolense, T. brucei brucei and T. vivax on the other (Fig. 2). Interestingly, two clades of T. theileri were observed (IIB and IA/IB) as previously described [29]. Furthermore, the occurrence of two lineages was also observed in the main group of T. vivax, cluster C and cluster A [30]. Cluster C had previously been reported in various regions in Africa and America, while cluster A was described only in Tanzania [FM164789; FM164787]. T. vivax C and A were found co-infecting cattle in the Faro et Deo region.

Fig. 2
figure2

Molecular phylogenetic analysis by Maximum Likelihood method based on the gGAPDH–encoding gene sequence as detailed under “Material and Methods”. It contains an alignment of 535 bp stretches of 37 sequences obtained in this study plus reference sequences [HQ664796; FM164792; HQ664805; HQ664784, HQ664792; HF545654; FM164789; XM_840453; FN400713] retrieved from Garcia et al. [29] and Hamilton et al. [23]. The bootstrap support values (> 70% in 1000 replications) are shown for the nodes

Correlation of packed cell volume with infection status

Animals with single or mixed infections had lower PCV values when compared to those with no infection (Fig. 3). When comparing the mean PCV with the type of infection, animals with single-infections of T. vivax (31.68 ± 5.40) and T. congolense (31.29 ± 6.92) showed no significant differences from uninfected. Animals carrying T. theileri had a mean PCV of 31.9 ± 4.5 (n = 16) for clade IIB while that for clade IA and IB it was 35.8 ± 3.4 (n = 8) (Additional file 3: Table S3). The observed difference was close to significance (F = 2.043, p = 0.056). Animal infected with T. grayi had the lowest PCV (29 ± 5.5, Ν = 4) of all the groups. However, because of the small sample size of the T. grayi group we could not test for statistical significance.

Fig. 3
figure3

Effect of the species of trypanosomes detected by PCR on the Packed Cell Volume (PCV). Mixed infection is defined as the combination of at least two trypanosome species identified in the same animal. Details of sample collections and processing are indicated in the section “Materials and Methods”

Discussion

The present study was carried out to determine the prevailing species and genetic diversity of trypanosomes infecting cattle in five divisions located in two agro-ecological zones of northern Cameroon, using both microscopy and molecular methods. The overall prevalence using microscopy is in agreement with previously reported prevalences of 3.7 to 20%, which were also determined by microscopy only [10]. However, infection rates determined by molecular analysis with ITS-1 nested PCR (53.2%) were much higher. This underpins the difficulty of microscopy to detect parasites at low levels of parasitemia in subclinical infections.

On the other hand, out of 56 trypanosome-positive cases by microscopy, only 41were detected by nested PCR giving the concordance rate of 73.2% between both techniques. This discrepancy has already been reported by Takeet et al. [31] and Adams et al. [19], the latter developing the primers used in our study. They also failed to amplify 56% of samples previously detected positive by microscopy and attributed this failure to the quality and quantity of the extracted parasite DNA. It is also possible that the primers do not amplify all trypanosome parasites [32, 33] or that Borrelia bacteria present in the blood are misinterpreted as trypanosomes, since based on their shape, size and movement, under the microscope they appear similar to T. brucei in buffy coat slide preparations [34]. Actually, recent molecular studies showed that 17.7% of cattle in the northern Cameroon are infected with Borrelia theileri (B. Abanda, A. Paguem, M. Abdoulmoumini, TK. Manchang, A. Renz and A. Eisenbarth. personal communications).

We distinguished only three species of trypanosomes by microscopy, namely T. congolense, T. vivax and T. brucei spp., while others, which we named T. brucei-like, could not be identified beyond doubt. By using PCR, we were able to identify five species of trypanosomes in the study area. This can be explained by the high sensitivity of the generic primers (ITS-1), which can detect traces of DNA up to one parasite per mL of blood of both pathogenic and non-pathogenic species [17, 18]. In contrast, microscopy of the Buffy-coat extracted from a microcapillary tube can reliably detect motile parasites only at a concentration being higher than 1.25 × 103 parasites/mL of blood [15,16,17]. Such high parasite titers in blood are more typical for trypanosomes causing pathology, like T. brucei spp., T. congolense and T. vivax at the acute clinical stage, and chronic infections are likely to be missed.

Surprisingly, the stercorarian parasites T. theileri/T. grayi were the most prevalent species (30.5%) in our study. These two parasites cannot be distinguished by ITS-1 size estimation, but only by sequence analysis, because they have a similar band size of 320 bp on the gel. Four out of 30 samples analyzed by sequencing were identified as T. grayi whereas the other 26 were T. theileri.

Trypanosoma grayi was found in two of 7 cattle from Mayo-Tsanaga and in one of 6 from Vina and yet another one from 12 cattle examined at Mayo-Rey. Previously, this species was known only to be a parasite of reptiles [1]. However, recently this parasite has been detected in a White Fulani cattle in Faro et Deo [24] and has now also been found in Kapsiki and Gudali cattle. By extrapolation on our 358 T. theileri/T. grayi cases we could expect almost 50 animals to be infected with T. grayi. This observation raises concerns whether these parasites may represent a strain undergoing a change of host range [24]. Further investigations are essential to characterize those T. grayi strains and evaluate their pathogenic potential for cattle and/or other livestock. In our study areas animals infected with this parasite correlated with lower PCV which may be an indicator of potential pathogenic effects on animal health. In this context, it is noteworthy that a recent study in Nigeria has observed a high frequency of tsetse colonised with T. grayi-like parasites (J. Weber. personal communication). Furthermore, these parasites revealed a high genetic diversity suggesting a dynamic evolution in this region. The 320 bp amplicon representative for the stercorarian parasites of T. theileri/T. grayi was most prevalent in the tsetse-free Vina (47.3%) and Mayo-Tsanaga (60.0%) regions and much less frequent in the tsetse-infested areas Faro (21.0%), Faro et Deo (19.4%) and Mayo-Rey (8.2%) (Fig. 4). This observation suggests that abundant mechanical vectors are the drivers of transmission of Stercoraria in the presumably tsetse-free areas [10, 35]. The entomological survey by Lendzele et al. [36] in the Vina and Mayo-Rey division identified seven species of tabanids as potential mechanical vectors: Tabanus gratus, Ta. par, Ta. taeniola, Ta. biguttatus, Ta. sufis and Chrysops distinctipennis. Furthermore, four prevailing species of tabanids were found in the Far North region: Atylotus agrestis, Ta. taeniola, Ta. par and Ancala spec [37].. Desquesnes and Dia [38, 39] have proved experimentally the mechanical transmission of T. vivax and T. congolense by tabanids (Atylotus agrestis). In addition, Ta. par and Ta. taeniola were tested PCR positive for the presence of T. congolense, T. theileri, T. evansi and T. brucei in South Africa and the Gambia [40]. Additionally, ixodid ticks were identified as vectors of T. theileri in Germany and in Sudan [4, 5]. However, to our knowledge no detailed studies on mechanical vectors have been performed in the study areas until now.

Fig. 4
figure4

Distribution of Salivaria (T. brucei, T. vivax and T. congolense) and Stercoraria (T. theileri/T. grayi) in tsetse free and tsetse infested areas in Northern Cameroon. Details of sample collections and processing are indicated in the section “Materials and Methods”.(map depicted in Fig. 4 is from our own)

Infections with T. brucei spp. (5.0%) and T. vivax (6.7%) were the most prevalent classical pathogenic trypanosomes found in our study areas. They were significantly more prevalent in Faro et Deo and Mayo-Rey (Tables 1 and 2) compared to the other locations. This observation was expected, because Faro et Deo is situated between the tsetse-infested Gashaka Forest Reserve on the Nigerian border and the Faro Game Reserve, and Mayo-Rey between the hunting zones and the Bouba Ndjida National Park, which both harbor a large population of known reservoir species for trypanosomes (antelopes, buffalos, etc.) with particularly abundant tsetse populations [9, 41]. The high infection rate observed in Faro et Deo is in agreement with prevalences of 10 to 41%, obtained in earlier studies [7, 9, 10, 24]. In this area Glossina morsitans submorsitans and G. palpalis palpalis are the main prevailing vector species [24, 41]. In Mayo-Rey, G. tachinoides was also abundant, together with G. m. submorsitans [9].

Lower prevalences of T. vivax, T. congolense and T. brucei spp. were observed in Mayo-Tsanaga (11.7%) and Vina (7.8%), both considered tsetse free. However, the presence of these trypanosomes may indicate presence of tsetse in these areas, perhaps due to re-infestation of tsetse from the nearby wildlife reserves which had not been subject to tsetse control. It can be as well due to the introduction of infected animals from tsetse infested zones.

Trypanosoma congolense was detected in the Vina (n = 4/283; 1.4%) and in Mayo-Tsanaga (n = 12/205; 5.8%) only by molecular methods, a status which does not exclude the activities of tsetse in these areas. In the Adamawa plateau prevalence of 3% by microscopy and 21% by serological tests were previously reported [42]. For the Mayo-Tsanaga division this is the first report of T. congolense in cattle. However, since no molecular confirmation was done before, it is possible that these infections were misidentified previously. Or they may have been recently introduced by infected tsetse or infected Fulani animals coming from transhumance through tsetse infested areas of neighboring countries like Nigeria.

Trypanosoma vivax sequence analysis revealed the occurrence of two phylogenetically distinct strains: T. vivax type C [30], previously described to be distributed across Africa and America, and T. vivax type A, which was isolated so far only in Tanzania [FM164789; FM164787]. In our study areas, we found both strains sympatric with other trypanosomes in the Faro et Deo division. The type A has been reported to be responsible for several outbreaks of bovine trypanosomosis in East Africa [30]. This raises the concern for potential outbreaks in the Faro et Deo region, and the potential to spread further throughout the country.

When looking at the PCV values, animals carrying mixed infections had significantly lower values than the non-infected and single-species infected animals. Furthermore, when comparing the sampling areas, Faro et Deo had the lowest PCV values both in infected and uninfected cattle. It has also to be considered that anaemia may be the result of other hemoprotozoan and/or helminths infections [43]. Infected Kapsiki cattle showed the lowest PCVs when compared to the other indigenous Bos taurus breed Doayo (Namchi). It has been previously reported [44] that the Doayo cattle were trypanotolerant while the Kapsiki were trypano-susceptible and this was associated with higher introgression of zebu alleles in the Kapsiki [26]. In a previous study, it was observed that T. theileri clade IIB, though considered non-pathogenic in cattle, correlated with low PCV in infected animals [24]. This tendency to become pathogenic was attributed to the genetic association to a previously described clade [29]. Comparing the PCV values of all animals in this study in which DNA of T. theileri was found, the mean PCV of cattle infected with clade IIB (31.9 ± 4.5) was slightly lower than those of animals with clade IA and IB (35.8 ± 3.4, p < 0.057). This implies that infections with clade IIB may be pathogenic to cattle and should be further investigated and considered during clinical control operations for cattle kept under local husbandry conditions. Once more, this underlines the importance to further investigate the development and evolution of trypanosome species, especially as these two clades of T. theileri are found worldwide. The prevalence of T. brucei spp, T. vivax and T. congolense in the tsetse-free areas of Mayo-Tsanaga and Vina raised questions whether the areas are still free and if tsetse control is sufficient enough to eradicate bovine trypanosomosis. Therefore, an entomological survey is urgently needed to check whether these previously tsetse-free areas have been re-infested by Glossina or whether these parasites are transmitted by non-tsetse vectors. Both scenarios call for close monitoring of the situation including molecular tools as used in this study as well as a thorough search for alternate vectors.

Conclusions

Bovine trypanosomosis is more prevalent in the two ecological zones of northern Cameroon than previously thought. Five trypanosome species and subtypes were identified. Unexpectedly several cases of T. grayi were detected in cattle. Therefore, it may not be excluded that this parasite is already adapted to the cattle host. Trypanosoma vivax, clade A, which had previously only been identified in Tanzania was found to be co-endemic with T. vivax clade A and T. vivax clade C in the Faro et Deo region. Furthermore, the presence of two strains of T. theileri, clades IIB and IA/IB, was confirmed. This high diversity of Trypanosoma species makes monitoring and local control more complex than previously thought. Finally, the abundance of pathogenic trypanosomes in tsetse infested areas is alarming and even more, the occurrence of T. vivax, T. brucei, T. congolense, T. theileri and T. grayi in tsetse-free areas implies that tsetse control alone is not sufficient to control trypanosomosis in livestock.

Methods

Study areas

This study was carried out in the Far North, North and Adamawa region of Cameroon (Fig. 5: Additional file 3: Table S4). These three regions are localized in two large Agro-Ecological Zones: the Sudano-Sahelian (Far North region and a larger part of the North region) and the Guinea savannah of the Adamawa plateau (Adamawa region with a little part of the North region). Cattle rearing is most abundant in the Guinea savanna of the Adamawa plateau with suitable climate and pasturelands for extensive cattle rearing. Overall, this plateau contributes to about 38% of beef production in the country [45]. The sampling sites were located in five divisions lying between latitudes 7 to 10°N and 11 to 15°E and covered an area of 164,000 km2 [46]. A strong climatic gradient runs through the wet high Guinea savannah in the Adamawa up to the dry Sudano-Sahelian zone in the Far North region. The rainy season in the Guinea savannah zone is from April to October, whereas in the Sudano-Sahelian zone it is from June to September. Annual rainfall ranges from 1400 to 1700 mm in the Guinean savannah and 800–1400 mm in the Sudano-Sahelian zone (Fig. 5).

Fig. 5
figure5

Map of the study area. Geographic map showing five Agro-Ecological Zones of Cameroon (based on information from the Institute of Agricultural Research for Development, IRAD, 2009). The cattle sampling areas (red stars) were located in the climate zones Guinea wet savannah and Sudano-Sahelian dry savannah. (map depicted in Fig. 5 is from our own)

Experimental design and animal selection

A cross-sectional survey was carried out between April 2014 and June 2015. For each herd visited, about 10% of the animals were sampled using a systematic random method described by Dohoo et al. [47]. In the Faro and Mayo Tsanaga divisions only the indigenous taurine cattle breeds Doayo (Namchi) and Kapsiki, respectively, were examined and sampled. From each animal, physical examinations were made and the following variables recorded: breed, sex, body condition score (BCS) using the method described by Pullan for White Fulani [32], on a scale from 0 to 5 (0–2: poor condition, 3–4: good condition and 5 very good condition or fat), and age by dentition categorized as young (< 2.5 years), mature (> 2.5–5 years) and older (> 5 years). In many farms only very few males were present in the herds causing the random selection to be applied on the animals found in the herd without balancing for sex proportions.

Assessment of packed cell volume (PCV) and trypanosome detection

Approximatively 5 mL of blood were collected from the jugular vein of each animal, using a vacutainer tube containing potassium ethylenediaminetetraacetic acid (EDTA) anticoagulant (VACUETTE® K3 EDTA). The samples were stored in a cooler box until processing within 6 h after collection either at a stationary or mobile laboratory in the field. Plasma was separated from blood by centrifugation at 3000 rpm for 15 min. Then the buffy coat was carefully collected and stored at 4 °C for subsequent DNA extraction. To determine the PCV, blood was introduced into capillary tubes (approx. 70 μL), and after sealing one end of the capillary tube with cristoseal (Sigma Aldrich, Germany) it was centrifuged at 12,000 rpm for 5 min using a microhaematocrit centrifuge (Hawksley, UK). The PCV was measured with a haematocrit reader (Hawksley Limited, UK). Animals that had a PCV value equal or less than 25% were considered anaemic. Subsequently, the capillary tube was cut with a diamond cutter 0.5 mm below the buffy coat to transfer the layer of white blood cells containing accumulated haemoparasites [16, 17] on to a clean microscope slide. After applying a coverslip over the buffy coat, approximately 200 fields of the preparation were examined for the presence of motile trypanosomes with a compound light microscope using 400x magnification [15]. The trypanosome species were classified according to previously described morphological criteria [14].

Genomic DNA extraction, purification, PCR amplification, sequencing of ITS-1 and gGAPDH

Genomic DNA from buffy coat was extracted using the Wizard Genomic DNA Purification Kit (Promega, Germany) according to the manufacturer’s instructions, and then stored at.

− 20 °C. Generic primers were used in a nested PCR targeting kinetoplastid ITS-1 as described previously [19, 24]. Briefly, the first reaction (25 μL final volume) contained 2 μM of each outer primers (Table 6), 0.2 mM dNTP mix, 0.5 U Dream Taq DNA polymerase (Thermo Scientific, Dreieich, Germany), 1× Dream Taq buffer, and 1 μL of extracted DNA. Nuclease-free water and genomic DNA of T. brucei, T. congolense or T. grayi were used as negative and positive controls, respectively. PCR amplification was carried out as follows: initial denaturation step at 95 °C for 60 s, followed by 30 amplification cycles at 94 °C for 60 s, at 52 °C for 60 s, at 72 °C for 30 s, and final extension at 72 °C for 5 min. Thereafter, the second PCR reaction was carried out with 1 μL of first PCR product diluted 80-fold as template under the same cycling conditions as described above, except for an annealing temperature of 54 °C, and using the inner primer pairs (Table 6). 20 μL of the resulting PCR product was loaded onto a 2% TBE agarose gel stained with 0.5 μg/mL of DNA Stain G (SERVA, Heidelberg, Germany). Positive PCR amplicons of variable fragment sizes representing different trypanosome species (Table 4, Additional file 1: Figure S1) were randomly selected for Sanger sequencing. For these samples, the second reaction was carried out in a total volume of 50 μL with 2 μL of 80-fold diluted first PCR product.

Table 6 Generic Primers used for PCR amplification

An approximately 900 bp region of the gGAPDH gene was amplified by nested PCR and sequenced using the primers described by Hamilton et al. [23]. Nested PCR was carried out using 2x Red Mastermix (Genaxxon Bioscience, Ulm, Germany) to generate PCR products for direct sequencing. Briefly, the first PCR reaction with a final volume of 25 μL contained 1x mastermix, 0.5 μM of outer primers (Table 6), and 2 μL of genomic DNA template under the following conditions: initial denaturation at 95 °C for 3 min, 30 cycles at 95 °C for 1 min, annealing at 55 °C for 30 s, elongation at 72 °C for 1 min, followed by a final elongation step at 72 °C for 10 min. The first PCR products were diluted 80-fold and 2 μL transferred to the second PCR reaction with the inner primers (Table 6, Additional file 2: Figure S2) under the same conditions as the first reaction. Amplified products were subjected to electrophoresis on 2% agarose gels. The selected positive PCR products were sent for sequencing (Macrogen, Netherlands).

A subset of positive amplicons was excised from the gel and purified using GeneJet Gel Extraction Kit (Thermo Scientific, Dreieich, Germany) according to the manufacturer’s instructions. DNA concentrations were determined by photometry on a Nanodrop 1000 (Thermo Scientific, Dreieich, Germany) before submitting them to a commercial sequencing provider (Macrogen).

Statistics and phylogenetic analysis

The results from the parasitological and molecular approaches were compared by Chi-Square tests to assess the association between prevalence, breed, BCS, sampled area and age group. Fisher’s Exact Test was done to compare mean PCV values. Since only 269 (23%) samples were collected from male, no sex-differentiating analysis was performed.

Differences were tested for significance at p < 0.05 using the statistical software program SPSS v.25.0.0 (IBM, USA). Obtained sequences were analyzed using Geneious (Biomatters, Auckland, New Zealand) and aligned to sequences retrieved from data bases searches (GenBank, NCBI, https://blast.ncbi.nlm.nih.gov/genbank/), and TriTrypDBv.6.0 (http://tritrypdb.org) using nucleotide BLAST.

To investigate the genetic diversity of trypanosomes present in the study area, and to analyze their phylogenetic relationship in order to detect subpopulations of trypanosomes restricted to respective study areas, gGAPDH was used as a marker locus. Phylogenetic trees were aligned and constructed by MEGA7 software [48], and the evolutionary history was inferred using the Maximum Likelihood method (ML) based on the Kimura 2-parameter model [49]. Confidence in branching relationships was assessed using bootstrap re-sampling over 1000 replicates. The final construct nucleotide length used in this analysis was 535 bp.

Availability of data and materials

All data generated and analyzed during this study are included in this published article and its supplementary information files or available from the corresponding author on reasonable request. The sequences generated during the present study are available in the NCBI Genbank repository under the accession numbers MK674001-MK674048, MK656901-MK656904.

Abbreviations

AEZ:

Agro Ecological Zones

EDTA:

Ethylenediaminetetraacetic acid

gGAPDH:

glycosomal glyceraldehyde 3-phosphate dehydrogenase gene

ITS-1:

Internal transcribed spacer 1 region of the ribosomal RNA gene locus

PCV:

Pack cell volume

SPP:

Subspecies

References

  1. 1.

    Hoare CA. Morphological and taxonomic studies on mammalian trypanosomes. X. Revision of the systematics. J Protozool. 1964;11(2):200–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Vickerman K, Tetley L, Hendry KA, Turner CM. Biology of african trypanosomes in the tsetse fly. Biol Cell. 1988;64(2):109–19.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  3. 3.

    Schuster JP, Schaub GA. Trypanosoma cruzi: skin-penetration kinetics of vector-derived metacyclic trypomastigotes. Int J Parasitol. 2000;30(14):1475–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  4. 4.

    Bose R, Friedhoff KT, Olbrich S. Transmission of Megatrypanum trypanosomes to Cervus dama by Tabanidae. J Protozool. 1987;34(1):110–3.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  5. 5.

    Latif AA, Bakheit MA, Mohamed AE, Zweygarth E. High infection rates of the tick Hyalomma anatolicum anatolicum with Trypanosoma theileri. Onderstepoort J Vet Res. 2004;71(4):251–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  6. 6.

    Meyer A, Holt HR, Selby R, Guitian J. Past and ongoing tsetse and animal trypanosomiasis control operations in five african countries: a systematic review. PLoS Negl Trop Dis. 2016;10(12):e0005247.

    PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Mamoudou A, Zoli A, Hamadama H, Bourdanne AS, Geerts S, Clausen PH, Zessin KH, Kyule M, van den Bossche P. Seasonal distribution and abundance of tsetse flies (Glossina spp.) in the Faro and Deo division of the Adamaoua plateau in Cameroon. Med Vet Entomol. 2008;22(1):32–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  8. 8.

    Tanenbe C, Gambo H, Musongong AG, Boris O, Achukwi MD. Prévalence de la trypanosomose bovine dans les départements du Faro et Déo, et de la Vina au Cameroun - bilan de vingt années de lutte contre les glossines. Rev Elev Med Vet Pays Trop. 2010;63(3–4):63–9.

    Article  Google Scholar 

  9. 9.

    Mamoudou A, Njanloga A, Hayatou A, Suh PF, Achukwi MD. Animal trypanosomosis in clinically healthy cattle of North Cameroon: epidemiological implications. Parasit Vectors. 2016;9:206.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Mamoudou A, Zoli A, Mbahin N, Tanenbe C, Bourdanne, Clausen PH, Marcotty T, Van den Bossche P, Geerts S. Prevalence and incidence of bovine trypanosomosis on the Adamaoua plateau in Cameroon 10 years after the tsetse eradication campaign. Vet Parasitol. 2006;142(1–2):16–22.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  11. 11.

    Afewerk Y, Clausen PH, Abebe G, Tilahun G, Mehlitz D. Multiple-drug resistant Trypanosoma congolense populations in village cattle of Metekel district, north-West Ethiopia. Acta Trop. 2000;76(3):231–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  12. 12.

    Mulugeta W, Wilkes J, Mulatu W, Majiwa PAO, Masake R, Peregrine AS. Long-term occurrence of Trypanosoma congolense resistant to diminazene, isometamidium and homidium in cattle at Ghibe, Ethiopia. Acta Trop. 1997;64(3–4):205–17.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  13. 13.

    Woo PT. The Haematocrit centrifuge technique for the diagnosis of african trypanosomiasis. Acta Trop. 1970;27(4):384–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Murray M, Murray PK, McIntyre WIM. An improved parasitological technique for the diagnosis of African trypanosomiasis. Trans R Soc Trop Med Hyg. 1977;71(4):325–6.

    CAS  Article  Google Scholar 

  15. 15.

    Paris J, Murray M, McOdimba F. A comparative evaluation of the parasitological techniques currently available for the diagnosis of African trypanosomiasis in cattle. Acta Trop. 1982;39(4):307–16.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Picozzi K, Tilley A, Fevre EM, Coleman PG, Magona JW, Odiit M, Eisler MC, Welburn SC. The diagnosis of trypanosome infections: applications of novel technology for reducing disease risk. Afr J Biotechnol. 2002;1(2):39–45.

    CAS  Article  Google Scholar 

  17. 17.

    Desquesnes M, McLaughlin G, Zoungrana A, Davila AMR. Detection and identification of Trypanosoma of African livestock through a single PCR based on internal transcribed spacer 1 of rDNA. Int J Parasitol. 2001;31(5–6):610–4.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  18. 18.

    Njiru ZK, Constantine CC, Guya S, Crowther J, Kiragu JM, Thompson RC, Davila AM. The use of ITS1 rDNA PCR in detecting pathogenic African trypanosomes. Parasitol Res. 2005;95(3):186–92.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  19. 19.

    Adams ER, Malele II, Msangi AR, Gibson WC. Trypanosome identification in wild tsetse populations in Tanzania using generic primers to amplify the ribosomal RNA ITS-1 region. Acta Trop. 2006;100(1–2):103–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  20. 20.

    Cox A, Tilley A, McOdimba F, Fyfe J, Eisler M, Hide G, Welburn S. A PCR based assay for detection and differentiation of African trypanosome species in blood. Exp Parasitol. 2005;111(1):24–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Fothergill-Gilmore LA, Michels PAM. Evolution of glycolysis. Prog Biophys Mol Biol. 1993;59(2):105–235.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Hannaert V, Opperdoes FR, Michels PA. Comparison and evolutionary analysis of the glycosomal glyceraldehyde-3-phosphate dehydrogenase from different Kinetoplastida. J Mol Evol. 1998;47(6):728–38.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  23. 23.

    Hamilton PB, Stevens JR, Gaunt MW, Gidley J, Gibson WC. Trypanosomes are monophyletic: evidence from genes for glyceraldehyde phosphate dehydrogenase and small subunit ribosomal RNA. Int J Parasitol. 2004;34(12):1393–404.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  24. 24.

    SCH N, Weber JS, Ngo Bum E, Gbem TT, Kelm S, Achukwi MD. Molecular screening of tsetse flies and cattle reveal different Trypanosoma species including T. grayi and T. theileri in northern Cameroon. Parasit Vectors. 2017;10(1):631.

    Article  CAS  Google Scholar 

  25. 25.

    Fréchou H. L’élevage et le commerce du bétail dans le nord du Cameroun. In: ORSTOM Cdl. Cahiers d’outre-mer, editor. Série Sciences Humaines. Bordeaux: Ateliers de l'Union Francaise d’Impression; 1966. p. 319–20.

    Google Scholar 

  26. 26.

    Achukwi MD, Tanya VN, Hill EW, Bradley DG, Meghen C, Sauveroche B, Banser JT, Ndoki JN. Susceptibility of the Namchi and Kapsiki cattle of Cameroon to trypanosome infection. Trop Anim Health Prod. 1997;29(4):219–26.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  27. 27.

    FAO. The management of global animal genetic resources. Rome: Anim Prod Health; 1992.

    Google Scholar 

  28. 28.

    FAO. Agro ecological Zoning Guidelines. Rome: Soil Resources MaCS; 1996.

    Google Scholar 

  29. 29.

    Garcia HA, Rodrigues AC, Martinkovic F, Minervino AH, Campaner M, Nunes VL, Paiva F, Hamilton PB, Teixeira MM. Multilocus phylogeographical analysis of Trypanosoma (Megatrypanum) genotypes from sympatric cattle and water buffalo populations supports evolutionary host constraint and close phylogenetic relationships with genotypes found in other ruminants. Int J Parasitol. 2011;41(13–14):1385–96.

    PubMed  Article  PubMed Central  Google Scholar 

  30. 30.

    Adams ER, Hamilton PB, Rodrigues AC, Malele II, Delespaux V, Teixeira MM, Gibson W. New Trypanosoma (Duttonella) vivax genotypes from tsetse flies in East Africa. Parasitology. 2010;137(4):641–50.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  31. 31.

    Takeet MI, Fagbemi BO, De Donato M, Yakubu A, Rodulfo HE, Peters SO, Wheto M, Imumorin IG. Molecular survey of pathogenic trypanosomes in naturally infected Nigerian cattle. Res Vet Sci. 2013;94(3):555–61.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  32. 32.

    Pullan NB. Condition score of white Fulani cattle. Trop Anim Health Prod. 1978;10:118–20.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  33. 33.

    Gonzales JL, Jones TW, Picozzi K, Cuellar HR. Evaluation of a polymerase chain reaction assay for the diagnosis of bovine trypanosomiasis and epidemiological surveillance in Bolivia. Kinetoplastid Biol Dis. 2003;2(1):8.

    PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Matton P, Van Melckebeke H. Bovine borreliosis: comparison of simple methods for detection of the spirochaete in the blood. Trop Anim Health Prod. 1990;22:147–52.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  35. 35.

    Tibayrenc R, Gruvel J. La campagne de lutte contre les glossines dans le bassin du lac tchad. II. Contrôle de l’assainissement glossinaire: Critique technique et financière de l’ensemblée de la compagne, conclusions générales. Rev Elev Med Vet Pays Trop. 1977;30(1):31–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  36. 36.

    Lendzele SS, Abdoulmoumini M, Lydie AYG. Spatial repartition of tabanids (Diptera: Tabanidae) in different ecological zones of North Cameroon. Biodivers Int J. 2017;1(2):00010.

    Google Scholar 

  37. 37.

    Mamoudou A, Marceline M, Pierre F, Lendzele S, Oumarou F, Grabed R, Kingsley M, Mbunkah A: Tabanids (Diptera: Tabanidae) fauna composition in different sites and biotopes of Far-North, Cameroon. JOBAN. 2016;6(3):146–154 [cited 11Jul.2019]. Available from: http://www.ikprress.org/index.php/JOBAN/article/view/1248

  38. 38.

    Desquesnes M, Dia ML. Trypanosoma vivax: mechanical transmission in cattle by one of the most common African tabanids, Atylotus agrestis. Exp Parasitol. 2003;103(1–2):35–43.

    PubMed  Article  PubMed Central  Google Scholar 

  39. 39.

    Desquesnes M, Dia ML. Mechanical transmission of Trypanosoma congolense in cattle by the african tabanid Atylotus agrestis. Exp Parasitol. 2003;105(3–4):226–31.

    PubMed  Article  PubMed Central  Google Scholar 

  40. 40.

    Taioe MO, Motloang MY, Namangala B, Chota A, Molefe NI, Musinguzi SP, Suganuma K, Hayes P, Tsilo TJ, Chainey J, et al. Characterization of tabanid flies (Diptera: Tabanidae) in South Africa and Zambia and detection of protozoan parasites they are harbouring. Parasitology. 2017;144(9):1162–78.

    PubMed  Article  PubMed Central  Google Scholar 

  41. 41.

    Achukwi MD, Gillingwater J, Njan Nloga AM, Simo G. Lack of evidence for sufficiently isolated populations of Glossina morsitans submorsitans on the Adamawa plateau of Cameroon following geometric morphometric analysis. Adv Entomol. 2013;01:1–7.

    Article  Google Scholar 

  42. 42.

    Mpouam SE, Achukwi MD, Feussom-Kameni JM, Bengaly Z, Ouedraogo GA. Serological and parasitological prevalence of bovine trypanosomosis in small holder farms of the Vina division, Adamawa region of Cameroon. J Parasitol Vector Biol. 2011;3(4):44–51.

    Google Scholar 

  43. 43.

    Radostits OM, Gay CC, Hinchcliff KW, Constable PD. Veterinary Medecine: a text book of diseases of cattle, horses, sheep, pigs and Goats10 th edition: Saunders Ltd: Elservier Health Science; 2007.

  44. 44.

    Achukwi MD, Musongong GA. Trypanosomosis in the Doayo/Namchi (Bos taurus) and zebu white Fulani (Bos indicus) cattle in Faro division, North Cameroon. J App Biosc. 2009;15(1):807–14.

    Google Scholar 

  45. 45.

    Minepia. Schéma Directeur pour le développement des filières de l’élevage au Cameroun. Yaounde: Ministère de l’Elevage, des Pêches et des Industries Animales; 2003.

    Google Scholar 

  46. 46.

    The population and housing census of Cameroon [http://cameroon.opendataforafrica.org/rfdefze/census-data]. Accessed on 11 Mar 2018.

  47. 47.

    Dohoo I, Wayne M, Henrik S. Veterinary epidemiologic Reseach 2nd edition. Charlottetown Prince Edward Island Canada: AVC Inc; 2009.

    Google Scholar 

  48. 48.

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

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  49. 49.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgments

We would like to thank Dr. Abba Ferdinand from the School of Veterinary Science, University of Ngaoundéré, late Dr. Almeck K. Aboubakar Dandjouma, Dr. Madi Palou Aboubakar, Dr. Bayang Nicolas of the veterinary research laboratory, Institute of Agricultural Research for Development Wakwa Centre for logistical and assistance during the fieldwork. The research staff personnel of the Programme Onchocercoses field station in Ngaoundéré: David Ekale, Jeremie Yembo, and late George Tamenai for their support. The research group of Glycobiochemistry, University of Bremen, in particular Ibrahim Mahamat, Dr. Frank Dietz and Petra Berger for training and support of the molecular analysis.

Funding

Data collection for this study was financed by the Otto Bayer Foundation (F-2013BS522), International Foundation for Science (IFS, B/5864–1). The German Research Foundation (DFG, grant no. RE 1536/2–1 and Ke428/10–1, Ke428/13–1) and the joint RiSC program of the State Ministry of Science, Research and Arts Baden Württemberg and the University of Tübingen (PSP-no. 4041002616) funded the molecular and bioinformatics analysis. The Funding bodies had no role in study design, data collection and analysis, interpretation of results or writing of the manuscript.

Author information

Affiliations

Authors

Contributions

Conceptualisation: AP, BA, AR, AE, MDA. Laboratory analysis: AP, JSW, SCHN, MDA. Investigation: AP, BA, MDA, TKM, MA, DN. Project administration: AP, AE, AR. Resources: AP, BA, MDA, SK, AR, AE, Supervision: MDA, A R, SK. Writing – review & editing: AP, BA, MDA, AR, AE, SCHN, JSW, SK, DN, TKM, MA. All authors read and approved the final manuscript.

Authors’ information

Department of Biological Science, Faculty of Science, University of Ngaoundéré, Cameroon

Archile Paguem, Babette Abanda, Sen Claudine Henriette Ngomtcho, Dieudonné Ndjonka

Department of Parasitology, School of Veterinary Medicine and Sciences, University of Ngaoundéré, Cameroon

Mamoudou Abdoulmoumini

Ministry of public health, Yaoundé, Cameroon

Sen Claudine Henriette Ngomtcho

Veterinary Research Laboratory, Institute of Agricultural Research for Development, Ngaoundéré, Cameroon

Manchang Tanyi Kingsley

TOZARD Research Laboratory, Bambili-Tubah, Bamenda, Cameroon

Mbunkah Daniel Achukwi

Department of Comparative Zoology, Institute for Evolution and Ecology, University of Tübingen, Tübingen, Germany.

Archile Paguem, Babette Abanda, Albert Eisenbarth, Alfons Renz

Institute of Novel and Emerging Infectious Diseases, Friedrich Löffler Institut, Federal Research Institute for Animal Health, Greifswald – Insel Riems, Germany

Albert Eisenbarth

Centre for Biomolecular Interactions Bremen (CBIB), Faculty of Biology and Chemistry, University of Bremen, 28334 Bremen, Germany

Judith Sophie Weber, Sørge Kelm

Corresponding author

Correspondence to Mbunkah Daniel Achukwi.

Ethics declarations

Ethics approval

Permission for the study and ethical approval were obtained from the Ethics committee of the Institute of Agricultural Research for Development (IRAD) in Cameroon, which is the country’s government institution for animal health and livestock husbandry improvement. Furthermore, verbal consent was given by the cattle owners and herdsmen and approved by the same ethics committee. All cattle owners and herdsmen participating in the study were contact by telephone and the purpose of the study was explained. When the oral consent was obtained from participant. The veterinarians were sent for biological material collection.

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

Additional file 1: Figure S1.

PCR amplicons of Trypanosoma species in northern Cameroon. ITS amplicon sizes of different Trypanosoma species in the range of 200 to 650 bp. The amplicons were resolved on a 2% TBE agarose gel. The first lane shows the marker (M), the second lane shows the ITS-1 fragment for T. vivax (Tv) at 200 bp and a faint band at around 180 bp. The third and fourth lanes show the presence of two species, T. theileri and T. grayi (Tth/Tg) at 380 bp, the fifth line T. brucei spp. (Tb) at 400 bp and the sixth lane the presence of T. congolense forest type (Tcf) at 640 bp. C-1: Water control of 1st reaction, C-2: Water control of 2nd reaction.

Additional file 2: Figure S2.

gGAPDH amplicons of different Trypanosoma species gave one band size of 900 bp. The first lane shows the marker (M), the second (852) and third (853) lanes are positives, the fourth (884), fifth (895), sixth (898) and seventh (849) lanes are negative samples. The eighth lane is the amplicon of T. grayi genomic DNA used as a positive control and C is double distilled water as a negative control.

Additional file 3: Table S1.

Effect of cattle breed on animal body condition. Table S2. Effect of study area, cattle breed and age group on packed cell volume. Table S3. Packed cell volume of animals infected with T. theileri clade IA, IB and IIB. Table S4. GPS Coordinates of study areas and sampled herds.

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

Paguem, A., Abanda, B., Ndjonka, D. et al. Widespread co-endemicity of Trypanosoma species infecting cattle in the Sudano-Sahelian and Guinea Savannah zones of Cameroon. BMC Vet Res 15, 344 (2019). https://doi.org/10.1186/s12917-019-2111-6

Download citation

Keywords

  • Trypanosomosis
  • ITS-1
  • gGAPDH
  • T. grayi
  • T. Theileri
  • Co-endemic trypanosomes
  • Cameroon