# Tabanids as possible pathogen vectors in Senegal (West Africa)

Mohamed Lamine Keita, Hacene Medkour, Masse Sambou, Handi Dahmana, Oleg Mediannikov

► To cite this version:

Mohamed Lamine Keita, Hacene Medkour, Masse Sambou, Handi Dahmana, Oleg Mediannikov. Tabanids as possible pathogen vectors in Senegal (West Africa). Parasites & Vectors, 2020, 13 (1), 500 [15 p.]. 10.1186/s13071-020-04375-w. hal-03149703

HAL Id: hal-03149703

https://amu.hal.science/hal-03149703v1

Submitted on 7 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Tabanids as possible pathogen vectors in Senegal (West Africa)

Mohamed Lamine Keita<fnref n="1" />,²†, Hacène Medkour<fnref n="1" />,²†, Masse Sambou³,⁴, Handi Dahmana<fnref n="1" />,² and Oleg Mediannikov<fnref n="1" />,²*

## Abstract

**Background:** Species of the Tabanidae are potent vectors of human and animal diseases, but they have not been thoroughly investigated to date. In Senegal (West Africa), little information is available on these dipterans. Our objective in this study was to investigate Senegalese tabanids and their diversity by using molecular and proteomics approaches, as well as their associated pathogens.

**Methods:** A total of 171 female tabanids were collected, including 143 from Casamance and 28 from Niokolo-Koba. The samples were identified morphologically by PCR sequencing and by MALDI-TOF MS, and PCR analysis was employed for pathogen detection and blood-meal characterization.

**Results:** The morphological identification revealed four species concordantly with the molecular identification: *Atylo-tus fuscipes* (79.5%), *Tabanus guineensis* (16.4%), *Chrysops distinctipennis* (3.5%) and *Tabanus taeniola* (0.6%) (not identi-fied by PCR). The molecular investigation of pathogens revealed the presence of *Trypanosoma theileri* (6.6%), *Leishma-nia donovani* (6.6%), *Setaria digitata* (1.5%), *Rickettsia* spp. (5.1%) and *Anaplasmataceae* bacteria (0.7%) in *A. fuscipes*. *Tabanus guineensis* was positive for *L. donovani* (35.7%), *S. digitata* (3.6%) and *Anaplasmataceae* (17.8%). *Leishmania donovani* has been detected in 50% of *C. distinctipennis* specimens and the only *T. taeniola* specimen. No Piroplasmida, *Mansonella* spp. or *Coxeilla burnetii* DNA was detected. In addition to humans (96.43%), *Chlorocebus sabeus*, a non-human primate, has been identified as a host of (3.57%) analysed tabanids. MALDI-TOF MS enabled us to correctly identify all tabanid species that had good quality spectra and to create a database for future identification.

**Conclusions:** Tabanids in Senegal could be vectors of several pathogens threatening animal and public health. To fully characterize these dipterans, it is therefore necessary that researchers in entomology and infectiology employ molecular characterization and mass spectrometric techniques such as MALDI-TOF MS to analyse these dipterans in Senegal and West Africa.

**Keywords:** Tabanids, *Leishmania donovani*, *Trypanosoma* spp., PCR, MALDI-TOF MS, Senegal

## Background

The Tabanidae is an important dipteran family belonging to the infraorder of Tabanomorpha and to the superfamily Tabanoidae [1, 2]. There are 4400 species

grouped in 114 genera, distributed in four subfamilies: Chrysopsinae; Pangoninae; Tabaninae; and Scepsidinae [2–4]. In Africa, only the first three subfamilies are of economic and medical interest due to their impact on livestock and their role as potential vectors of diseases [1, 2, 5]. Tabanidae have a global distribution, including 410 species in South Africa, 335 in North America, 244 in India, 120 in Malaysia and nearly a hundred in Thailand [1, 3–6]. Male insects feed on nectar and do not present medical interest. In contrast, female tabanids of

<footnote>*Correspondence: olegusss1@gmail.com</footnote>

<footnote>†Mohamed Lamine Keita and Hacène Medkour contributed equally to this work</footnote>

<footnote n="1">¹ IHU Méditerranée Infection – Microbes, Evolution, Phylogeny and Infection (MEΦI), Marseille, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
most species are haematophagous, being able to obtain large stores of blood from their hosts by telmophagy [1, 2]. These arthropods are considered to be primarily mechanical vectors of various human and animal pathogens, including bacteria, viruses, helminths, and protozoans [2, 6, 7]. Foil et al. [8] identified 25 different diseases transmitted by these dipterans. Krinsky [9] listed 189 species of the Tabanidae with approximately 30 infectious agents. Despite this role, the Tabanidae remain among the least-studied dipteran species [2]. To date, our research on PubMed using “Tabanidae” as a keyword yielded only 382 references, including 200 publications in 1976 by Krinsky [8] and 46 for African Tabanidae and 182 publications over 43 years.

To prevent damage from the Tabanidae and to elucidate their vectorial roles in the transmission of infectious diseases, correctly identifying their species and associated pathogens has always been necessary [3, 5]. This identification is based on both morphological and molecular tools [7, 10]. Nonetheless, the infectious status and hosts associated with these species can only be determined by molecular biology [11, 12]. Morphological identification has the advantage of being fast and economical but requires expertise from entomologists. This method may present such difficulties as alteration of specimens, variations in developmental stages and morphological similarities [1, 3].

The molecular identification of species proposed by Hebert et al. [3] is based on the amplification of 658 bp of the mitochondrial gene coding for cytochrome c oxidase subunit 1 (cox1). This method, called “barcoding”, enables the discrimination of a wide range of species and offers the advantage of integrating taxonomic and genomic data, regardless of genital or developmental status [5, 13]. In addition, one of the fundamental tasks of barcoding is to associate genomic sequences to the species names [14]. Furthermore, this method makes it possible to reveal the genetic diversity and to solve the ambiguities at the species level with precision, constancy and celerity [5]. Also, the molecular characterization of pathogens was largely based on the amplification of conserved ribosomal gene subunits [15]. However, there are some limitations for molecular identification, such as the scarcity of sequences on GenBank DNA database, particularly with regard to the Tabanidae.

The MALDI-TOF/MS technique is another more recent arthropod identification approach but is developing rapidly in terms of cost, time, sensitivity, specificity and reproducibility [16]. Unlike molecular analysis, this technique relies on the production and analysis of protein spectra specific to each species, regardless of the stage of development [17]. Abdomen, thorax or limbs

may be used for identification [17]. The difficulty with MALDI-TOF/MS is the lack of databases and reference spectra [16, 18]. Additionally, the use of certain parts of arthropods, such as the abdomen, containing exogenous proteins of various hosts influences the quality of spectra [17]. One of the objectives of the present study was to initiate a MALDI-TOF/MS database for the Tabanidae that has not been available to date.

In the Republic of Senegal, in West Africa, livestock farming is one of the major activities, particularly in rural areas. Vector-borne diseases, such as trypanosomiasis and filariasis, are endemic [19–21]. To date, 43 species of the Tabanidae belonging to the genera *Atylotus*, *Chrysops*, *Ancala*, *Hybomitra*, *Haematopota*, *Philoliche* and *Tabanus* have been reported [22]. A series of entomological works conducted by a generation of researchers, including Oldroyd [23, 24], Ovazza and Valade [25] and Raymond et al. [22], made it possible to perform a morphological description and study the taxonomy of this notably diverse fauna. The recurrent presence of HIV, tuberculosis, malaria or dengue fever since the 1980s has diverted researchers from investigating tabanids in Senegal. Except for some studies, such as Touré's in 1968 on cattle trypanosomiasis in Casamance [26], knowledge on pathogens carried by tabanids in Senegal is non-existent and marked by a scarcity of bibliographic resources.

Moreover, despite the existence of vector-borne diseases, such as filariasis and trypanosomiasis, to the best of our knowledge, no reports have described the molecular and proteomics identification of tabanid flies in Senegal. The aim of this study was to identify samples of the Tabanidae from Senegal and characterize their associated pathogens and hosts.

## Methods

### Study sites

Sample collection was carried out at two sites in southern Senegal: Lower Casamance in December 2014 and Niokolo-Koba National Park in August 2016.

Lower Casamance is located in southern Senegal and is isolated from the rest of the country by Gambia. Ziguinchor, the capital of the region, borders the Casamance River. The Lower Casamance covers an area of 7339 km² and is composed of three departments: Ziguinchor (regional capital), Bignona and Oussoye. The catches took place in several localities of the department of Oussoye. The climate is tropical of type subguineen, characterized by a rainfall of approximately 1500 mm/year, as well as a high hygrometry related to both the marine influence and a low thermal amplitude.

Niokolo-Koba National Park (PNNK) is a park located 650 km east of Dakar in the region of Tambacounda (southeastern Senegal) near the Guinean border. The
Gambia River runs through the park, as well as its two tributaries: the Koulountou River and the Niokolo-Koba River, which gave the park its name. The relief is quite flat with an altitude of 16–311 m. Mount Assirik is the highest point. Large plains, sometimes marshy in the rainy season, separate small hills of barely 200 m. The terrain becomes more rugged in the extreme southeast on the foothills of the Fouta-Djalon mountain range where the Niger, Senegal and Gambia rivers originate. These foothills form a natural border between Senegal and Guinea. The Simenti pond, the site of our collection, is the only permanent pond in the park and attracts savannah animals that drink very salty water and roll in the mud.

## Methods and period of collection

Lower Casamance tabanids were caught using tsetse monocone traps. Ten traps were placed at 6:00-12:00 h the morning between pastures and houses, close to water points. Those of Niokolo-Koba were collected manually around camp lights during the evening. Samples were pooled by location and day of collection and then stored in 70% ethanol at room temperature until morphological identification was performed.

## Morphological identification

All samples were subjected to double morphological identification; first, at the genus level in the collection sites, then at the species level in the University Hospital Institute-Méditerranée Infection insectarium, Marseille, France, using the Zeiss Axio Zoom V16 microscope (ZEISS, Oberkochen, Germany) and according to the different characteristics described by Oldroyd [23, 24]. Pictures were taken on three specimens of each species (Fig. 1).

## Dissection and sample preparation

The samples were rinsed with RNase-free water to eliminate ethanol and dissected longitudinally into two equal parts. One half was used for the molecular characterization of pathogens, species and hosts and the other half was stored in 70% ethanol for further analysis. Legs were used for analysis by MALDI-TOF MS; we preferred the legs to avoid interference during identification due to non-specific exogenous proteins from tabanids. These proteins are found in some parts of arthropods and generally affect the quality of MALDI-TOF MS spectra [27].

## DNA extraction

DNA extraction was performed using the QIAamp DNA Mini Kit® (Qiagen, Courtaboeuf, France) on BIOROBOT EZ1 (Qiagen) according to the manufacturer's protocol.

## Molecular identification of Senegalese tabanids and their hosts

To confirm the morphological identification of the tabanids, a standard PCR targeting the *cox1* gene using LCO/HCO primers developed by Folmer et al. (Table 1) was performed. The reaction mixture was 25 µl per reaction, including 5 µl of tabanid DNA template and 20 µl of mix (12.5 µl of Amplitaq Gold Master Mix; Eurogentec, Liège, Belgium; 6 µl of RNase-DNase-free water and 0.75 µl of each primer). No positive controls were used. After visualization by electrophoresis on a 2% agarose gel, the PCR product was purified using NucleoFast 96 PCR plates (Macherey Nagel EURL, Hoerdt, France) according to the manufacturer's instructions. The ampli-cons were sequenced using the Big Dye Terminator Cycle Sequencing Kit (Perkin Elmer Applied Biosystems, Foster City, CA, USA) with an ABI automated sequencer (Applied Biosystems). The obtained electropherograms were assembled and edited using ChromasPro software (ChromasPro 1.7, Technelysium Pty Ltd., Tewantin, Australia) and checked against the GenBank database by NCBI BLASTn (https://blast.ncbi.nlm.nih.gov/Blast.cgi).

An identical process using primers targeting the vertebrate-specific 16S rRNA gene was performed to identify hosts (Table 1). The *Agama agama* DNA extract from the IHU collection represented the positive control. The molecular identification involved 6 specimens of *Atylotus fuscipes*, 6 of *Tabanus guineensis*, 6 of *Chrysops distingutipennis* and the only specimen of *T. taeniola*.

## Pathogen study

Real-time quantitative PCR (qPCR) was used in the initial screening for pathogens. Reactions were performed in a CFX96 Real-Time system (Bio-Rad Laboratories, Foster City, CA, USA) according to the manufacturer's instructions. The reaction mixture consisted of 5 µl of DNA and 15 µl of mix, including 10 µl Eurogentec Master Mix (Eurogentec, Liège, Belgium), 0.5 µl of each primer and uracile-DNA glycosylase (UDG), and 0.5 µl of the FAM-labelled probe. qPCR amplification was performed using the following thermal profile: incubation at 50 °C for 2 min for UDG action (eliminating PCR amplicon contaminants), followed by an activation step at 95 °C for 3 min, followed by 40 cycles of denaturation at 95 °C for 15 s and annealing-extension at 60 °C for 30 s.

The 171 specimens of the Tabanidae were screened using qPCR systems targeting the following genes: 5.8S rRNA for *Piroplasmida*; IS1111 for *Coxiella burnetii*; 28S rRNA for Filaroida; 5.8S rRNA for *Trypanosoma* spp.; kinetoplast DNA (kDNA) for *Leishmania donovani* complex (*L. donovani* and *L. infantum*), ITS1 for *Mansonella* spp.; citrate synthase (*gltA*) specific for *Rickettsia* spp.; and 23S rRNA for *Anaplasmataceae* (Table 1). Then, a
standard PCR amplifying larger amplicons followed by sequencing was carried out for the qPCR-positive samples (Table 1). PCR assays and sequencing were performed as described above. In addition, semi-nested PCR was applied to optimize the detection capacity for Filaroida and *Rickettsiae*. The primers and probes used in this study are listed in Table 1.

## Phylogenetic analysis

The obtained sequences for the Tabanidae and their associated pathogens were assembled and edited as mentioned above. Then, the sequences were compared to sequences available on GenBank (https://blast.ncbi.nlm.nih.gov/Blast.cgi). Molecular phylogenetic and evolutionary analyses were conducted in MEGA7 (https://www.megasoftware.net/) using the maximum likelihood (ML) or neighbour-joining (NJ) methods with 100 replicates.

## MALDI-TOF MS analysis

Tabanids were individually rinsed with distilled water, and the 6 legs were removed and homogenized in an Eppendorf tube containing a pinch of glass powder (Sigma-Aldrich, Lyon, France) and 40 µl of a 70% (v/v) formic acid mixture and 50% (v/v) acetonitrile (Fluka, Buchs, Switzerland) using the TissueLyser device (Qiagen, Hilden, Germany) and following configuration parameters as previously described [18].

The homogenized tabanid legs were centrifuged at 2000× g for 30 s, and 1 µl of the supernatant from each sample was carefully dropped onto the MALDI-TOF target plate. Each spot was then recovered with 1 µl of matrix solution composed of saturated α-cyano-4-hydroxycinnamic acid (Sigma-Aldrich), 50% acetonitrile

(v/v), 2.5% trifluoroacetic acid (v/v) (Sigma-Aldrich, Dorset, UK) and HPLC-grade water. The target plate, after several minutes of drying at room temperature, was introduced into the Microflex LT MALDI-TOF Mass Spectrometer device (Bruker Daltonics, Bremen, Germany) for analysis. The loading of the MS target plate, the matrix quality and the performance of the MALDI-TOF were performed as previously described [18]. The obtained MALDI-TOF MS spectra were analysed using Flex 3.3 analysis software and ClinPro-tools 2.2. Only good-quality spectra were selected according to their intensity and reproducibility for subsequent analyses. The spectra of two specimens of *T. guineensis*, *C. distinctipennis* and *A. fuscipes* were introduced into the database. These species have been previously identified by morphology and confirmed by molecular biology (except for *T. taeniola*, which could not be identified molecularly). The spectra of the remaining specimens were used for the blind test to finally assess the ability of MALDI-TOF MS to identify them.

## Results

### Collection and morphological characterization of tabanids

A total of 171 samples of the Tabanidae were captured. Due to the quality of the bait, all the specimens were females, including 143 in Casamance and 28 in Niokolo-Koba. The morphological identification revealed three species in Casamance, including *A. fuscipes* Ricardo, 1908 for 136 samples (95%), *C. distinctipennis* Austen, 1906 for 6 samples (4%) and *T. taeniola* Palisto de Beavios 1807 for 1 sample (0.69%). Niokolo-Koba samples consisted only of *T. guineensis* Wiedemann, 1824 (*Tabanus subangustus* Ricardo, 1908), 28 samples (100%) [28]. In the present study, the abundant species was *A. fuscipes* (79.5%).

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Morphotypes of the collected tabanid species. **a** *Tabanus taeniola*. **b** *Chrysops distinctipennis*. **c** *Tabanus guineensis*. **d** *Atylotus fuscipes*</figcaption></figure>
<table><caption>**Table 1** PCR systems used in this study and their sources</caption><thead><tr><th>Targeted organism</th><th>Target gene</th><th>Name</th><th>Primers (5'-3')</th><th>Tm* (°C)</th><th>Reference</th></tr></thead><tbody><tr><td rowspan="2">Tabanidae</td><td rowspan="2">cox1</td><td>LCO1490</td><td>GGTCAACAAATCATAAAGATATTGG</td><td rowspan="2">50</td><td rowspan="2">[49]</td></tr><tr><td>HCO2198</td><td>TAAACTTCAGGGTGACCAAAAAAATCA</td></tr><tr><td rowspan="2">Hosts</td><td rowspan="2">16S</td><td>16SA-2290</td><td>CGCCTGTTTACCAAAAACAT</td><td rowspan="2">50</td><td rowspan="2">[50]</td></tr><tr><td>16SB-2860</td><td>CCGGTCTGAACTCAGATC ACGT</td></tr><tr><td rowspan="2">Leishmania donovani complex</td><td rowspan="2">kDNA</td><td>RV1</td><td>CTTTTCTGGTCCTCCGGGTAGG</td><td rowspan="2">60</td><td rowspan="2">[51]</td></tr><tr><td>RV2</td><td>CCACCCGGCCCTATTTTACACCAA</td></tr><tr><td rowspan="2">Leishmania Leish. spp.</td><td rowspan="2">ITS1</td><td>Probe. Leish*</td><td>FAM-TTTTCGCAGAACGCCCCTACCCGC-TAMRA</td><td rowspan="2">55</td><td rowspan="2">[52]</td></tr><tr><td>rDNA-10 F</td><td>CAATACAGGTGATCGGACAGG</td></tr><tr><td rowspan="4">Piroplasmida</td><td rowspan="4">5.8S</td><td>rDNA-14 R</td><td>CACGGGGATGACACAATAGAG</td><td rowspan="4">60</td><td rowspan="4">[53]</td></tr><tr><td>5.8S-F5</td><td>AYYKTYAGCGRTGGATGTC</td></tr><tr><td>5.8S-R</td><td>TCGCAGRAGTCTKCAAGTC</td></tr><tr><td>5.8S-S*</td><td>FAM-TTYGCTGCGTCCTTCATCGTTGT-MGB</td></tr><tr><td rowspan="4">Trypanosoma Leish. spp.</td><td rowspan="4">5.8S</td><td>S. 5.8 S Tryp sp.</td><td>FAM-GTTGAAGAACGCAGCAAAGIGCGAT-TAMRA</td><td rowspan="4">60</td><td rowspan="4">[37]</td></tr><tr><td>F. 5.8 S Tryp sp.</td><td>CAACGTGTCGCGATGGATGA</td></tr><tr><td>R. 5.8 S Tryp sp.</td><td>ATTCTGCAATTGATACCACTTATC</td></tr><tr><td>28S</td><td>F2 28S</td><td>ACCAAGGAGTCAAACAGACG</td></tr><tr><td rowspan="4">Kinetoplastida Leish. spp.</td><td rowspan="4">28S</td><td>R1 28S</td><td>GACGCCACATATCCCTAAG</td><td rowspan="4">58</td><td rowspan="4">[54]</td></tr><tr><td>ITS1-CF</td><td>CCGGAAGTTCACCGATATTG</td></tr><tr><td>ITS1-BR</td><td>TTGCTGCGTTCTTCAACGAA</td></tr><tr><td>F2 28S</td><td>ACCAAGGAGTCAAACAGACG</td></tr><tr><td rowspan="4">Coxiella burnetii</td><td rowspan="4">IS1111</td><td>R1 28S</td><td>GACGCCACATATCCCTAAG</td><td rowspan="4">58</td><td rowspan="4">[55]</td></tr><tr><td>F 720</td><td>GTTAAAGGGTTCGTAGTTGAA</td></tr><tr><td>R1495</td><td>GACTACAATGGTCTCTAATCA</td></tr><tr><td>F</td><td>CAAGAAACGTATCGCTGTGGC</td></tr><tr><td rowspan="4">Anaplasma Leish. spp.</td><td rowspan="4">23S</td><td>R</td><td>CACAGAGCCACCGTATGAATC</td><td rowspan="4">60</td><td rowspan="4">[56]</td></tr><tr><td>S</td><td>FAM-CCGAGTTCGAAACAATGAGGGCTG-TAMRA</td></tr><tr><td>TtAna-F</td><td>TGACAGCGTACCTTTTGCAT</td><td rowspan="2">60</td></tr><tr><td>TtAna-R</td><td>GTAACAGGTTCGGTCCTCCA</td></tr><tr><td rowspan="2">Borrelia Leish. spp.</td><td rowspan="2">16S</td><td>TtAna-S*</td><td>FAM-CTTGGTTTCGGGTCTAATCC-TAMRA</td><td rowspan="2">60</td><td rowspan="2">[58]</td></tr><tr><td>Bor_16S_3F</td><td>AGCCTTTAAAGCTTCGCTTGTAG</td></tr><tr><td rowspan="4">Rickettsia Leish. spp.</td><td rowspan="4">Citrate synthase (gltA)</td><td>Bor_16S_3R</td><td>GCCTCCCGTAGGAGTCTGG</td><td rowspan="4">60</td><td rowspan="4">[46]</td></tr><tr><td>Bor_16S_3P</td><td>FAM-CCGGCCTGAGAGGGTGAACGG-TAMRA</td></tr><tr><td>RKNDO3_F</td><td>GTGAATGAAAGATTACACTATTTAT</td></tr><tr><td>RKNDO3_R</td><td>GTATCTTAGCAATCATTCTAATAGC</td></tr><tr><td rowspan="4">gltA (First PCR)</td><td rowspan="4">Citrate synthase (gltA)</td><td>RKNDO3_S</td><td>FAM-CTATTATGCTTGCGGCTGTCGGTTC-TAMRA</td><td rowspan="4">54</td><td rowspan="4">[59]</td></tr><tr><td>CS2D</td><td>ATGACCAATGAAAATAATAAT</td></tr><tr><td>CSEndR</td><td>CTTATACTCTCTATGTACA</td></tr><tr><td>409D</td><td>CCTATGGCTATTATGCTTGC</td></tr><tr><td rowspan="2">gltA (nested PCR)</td><td rowspan="2">Citrate synthase (gltA)</td><td>1258R</td><td>ATTGCAAAAAAGTACAGTGAACA</td><td rowspan="2">54</td><td rowspan="2">[59]</td></tr><tr><td>Manso_ITS_F</td><td>CCTGCGGAAGGATCATTAAC</td></tr><tr><td rowspan="3">Mansonella Leish. spp.</td><td rowspan="3">ITS</td><td>Manso_ITS_R</td><td>ATCGACGGTTTAGGCGATAA</td><td rowspan="3">60</td><td rowspan="3">[20]</td></tr><tr><td>Manso_ITS_P</td><td>FAM- CGGTGATATTCGTTGGTGTCT-TAMRA</td></tr><tr><td>Fwd.18S.631</td><td>TCGTCATTGCTGCGGTTAAA</td></tr><tr><td rowspan="2">Filaria Leish. spp.</td><td rowspan="2">18S</td><td>Rwd.18S.1825r</td><td>GGTTCAAGCCACTGCGATTAA</td><td rowspan="2">54</td><td rowspan="2">[60]</td></tr><tr><td></td><td></td></tr></tbody></table>

**Molecular characterization and phylogeny of the Tabanidae**

BLAST analysis of the *cox1* sequences revealed that *T*.

*guineensis* specimens from Niokola-Koba exhibited approximately 93% identity with *Tabanus subangustus* described in East Africa (Kenya) (GenBank: KX946537)
[28]. *Atylotus fuscipes* had 99% identity with *A. nigromaculatus* (GenBank: KX946515) and *A. diurnus* (GenBank: KX946508) from East Africa. The identity of *C. distingutipennis* was up to 98% with that described in Uganda (GenBank: KX946519). We were unable to confirm the molecular identity of *T. taeniola* due to the poor quality of the sequences. The phylogenetic analysis using the neighbour-joining method revealed the existence of three clades: Clade 1 (*Tabanus* spp.); Clade 2 (*Atylotus* spp.); and Clade 3 (*Chrysops* spp.). The *cox1* sequence of *Acalles guadarramaensis* (GenBank: GU988096) was used as the out-group (Fig. 2). *cox1* sequences of tabanids in the present study were deposited in the GenBank database under the accession numbers MW013783-MW013793.

## Molecular characterization and phylogeny of hosts

Molecular and phylogenetic analysis of mitochondrial 16S rRNA genes amplified from the blood meal of tabanids using vertebrate-specific primers identified *Homo sapiens* as the host of *A. fuscipes*, *C. distincipennis* and *T. taeniola*. The hosts of *T. guineensis* were *C. sabeus* and *H. sapiens*. The results obtained showed that 96.43% of all tabanids analysed had their *H. sapiens* meals compared to 3.57% for *C. sabeus*.

## Molecular detection, identification and phylogenetic analysis of pathogens

### Detection of Trypanosoma spp.

Among the analysed samples, 9/136 (6.6%) *A. fuscipes* were positive for *Trypanosoma* spp. by qPCR for the 5.8S gene (Table 2). Amplifications were obtained for the 28S and 18S genes (550 bp and 600 bp, respectively) using Kinetoplastida-specific primers; and for the ITS1 region (250–710 bp) using *Trypanosoma* spp.-specific primers. No positive amplification was found in the other tabanid species (Table 2). Sequencing revealed that 8 of these samples presented 97% identity for the 28S gene with *T. minasense* (GenBank: AB362411), > 99% identity with *T. theileri* (GenBank: KR024688; LC385952) for the 18S gene and 92–93% identity with *Trypanosoma* sp. (GenBank: AB569248) and *T. theileri* (GenBank: KX569248) for the ITS1, respectively. One sample (TBC 05) exhibited divergent similarity for the ITS1 sequence with the other trypanosomes in the GenBank database, but sequencing for the other genes (18S and 28S) was not possible. The sequences of the specimen were of remarkable specificity. The phylogenetic analysis based on ITS1 sequences (Fig. 3) showed the formation of two clusters, one authentic *T. theileri* cluster and another cluster atypical of all Gen-Bank kinetoplastid sequences. *Homo sapiens* has been

identified as the host of all positive samples for *Trypanosoma* spp. Sequences were deposited in the GenBank database under the accession numbers MN244142-MN244147 (ITS1); MN244154-MN244159 (28S) and MN244160-MN244163 (18S).

### Detection of Leishmania spp.

Molecular detection of *Leishmania* showed that 13.45% (23/171) of tabanids, consisting of 9/136 (6.61%) samples of *A. fuscipes*, 3/6 (50%) of *C. distingutipennis*, 1/1 of *T. taeniola* and 10/28 (35.7%) of *T. guineensis*, were positive by qPCR. Only 4/9 positive samples of *A. fuscipes* and 2/10 positive samples *T. guineensis* were positive by standard PCR for the ITS1 region (480 bp) and the ampli-cons were successfully sequenced (Table 2). The obtained sequences were deposited to the GenBank database under the accession numbers MN244150-MN244153. BLAST and phylogenetic analysis confirmed the presence of *L. donovani* (GenBank: KU680955) (Fig. 4). *Homo sapiens* was identified as host for the tabanids positive to *Leishmania* spp.

### Filaroid nematode detection

The search for filaria yielded, by qPCR for 28S and standard PCR for 18S, 2 positive samples out of 136 (1.5%) *A. fuscipes* and 1/28 (3.5%) *T. guineensis* (Table 2). Phylogenetic analysis confirmed the detection of *S. digitata* (GenBank: DQ094175) (Fig. 5). These sequences were deposited to the GenBank database under the accession numbers MN244148-MN244149. The host identified for all *Setaria digitata*-positive samples was *H. sapiens*.

### Rickettsia spp. detection

Among all the tabanids tested, only *A. fuscipes* was positive for *Rickettsiae*, with 7/136 (5.1%) samples by qPCR. Four of them were amplified by nested PCR and successfully sequenced (Table 2). Two obtained sequences were nearly identical (99.33%) to *R. slovaca* isolate WB07 (GenBank: MK624992). One sequence had 99% identity with *R. africana* (GenBank: KX819298), which was detected on *Hyalomma marginatum* (camel tick) from Egypt [29]. Another sequence showed 100% similarity with *R. montanensis* str. OSU 85-930 (GenBank: CP003340) (Fig. 6).

### Detection of Anaplasmataceae

Using a 23S-based qPCR, we detected DNA of *Anaplasmataceae* in 1/136 (0.73%) *A. fuscipes*, 1/6 (16.66%) *C. distincipennis* and 5/28 (17.85%) *T. guineensis*. The only specimen of *T. taeniola* tested negative (Table 2). We were unable to identify the positives by 23S-based standard PCR.
<figure><img src="image_4.png" /></figure>

## MALDI-TOF/MS analysis

A total of 22 specimens, including 14 *T. guineensis*, 6 *C. distinctipennis*, one *T. taeniola* and one *A. fuscipes*, were submitted for MALDI-TOF MS analysis. We obtained good quality spectra for 20 specimens and poor quality for 1 one *C. distinctipennis* and one *T. taeniola*. Intraspecies reproducibility and species specificity by analysing MS spectra obtained from different species using Flex analysis software are shown in Fig. 7. The blind test showed that the 16 specimens with good quality spectra (12 *T. guineensis* and 4 *C. distinctipennis*) were correctly identified, with scores ranging from 1.800 to 2.519.

## Discussion

In this period, marked by the increasing success of genomics, the molecular and proteomic identification of Senegalese tabanids as vectors of pathogens is a reliable approach. Indeed, this approach makes it possible to further characterize these important arthropods, to better define their roles in the transmission of pathogens and to improve the tools intended for their elimination.

All of the specimens captured in this study were females maybe due to the quality of the trap used in Casamance. In Niokolo-Koba, capture of females only was a coincidence. Ovazza & Valade [25], observing the diurnal
<table><caption>**Table 2** Summary of pathogen research results</caption><thead><tr><th rowspan="2">Pathogen</th><th colspan="3">Atylotus fuscipes (n = 136)</th><th colspan="3">Chrysops distinctipennis (n = 6)</th><th colspan="3">Tabanus taeniola (n = 1)</th><th colspan="3">Tabanus guineensis (n = 28)</th><th rowspan="2">Total (n = 171)</th></tr><tr><th>qPCR n (%)</th><th>PCR n (%)</th><th>BLAST n (%)</th><th>qPCR n (%)</th><th>PCR n (%)</th><th>BLAST n (%)</th><th>qPCR n (%)</th><th>PCR n (%)</th><th>BLAST n (%)</th><th>qPCR n (%)</th><th>PCR n (%)</th><th>BLAST n (%)</th></tr></thead><tbody><tr><td>Trypanosoma spp.</td><td>9 (6.6)</td><td>9 (6.6)</td><td>T. theileri</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>9 (5.26)</td></tr><tr><td>Leishmania spp.</td><td>9 (6.6)</td><td>4 (2.9)</td><td>L. donovani</td><td>3 (50.0)</td><td>0 (0)</td><td>-</td><td>1 (100)</td><td>0 (0)</td><td>-</td><td>10 (35.7)</td><td>2 (7.1)</td><td>L. donovani</td><td>23 (13.4)</td></tr><tr><td>Filaria</td><td>2 (1.47)</td><td>2 (1.5)</td><td>S. digitata</td><td>0 (0)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1 (3.5)</td><td>1 (3.57)</td><td>S. digitata</td><td>3 (1.7)</td></tr><tr><td>Rickettsia spp.</td><td>7 (5.1)</td><td>4 (2.7)</td><td>R. slovaca; R. africana; R. mon-tanensis</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>7 (4.1)</td></tr><tr><td>Anaplasmataceae</td><td>1 (0.7)</td><td>0 (0)</td><td>-</td><td>1 (16.7)</td><td>0 (0)</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>5 (17.8)</td><td>0 (0)</td><td>-</td><td>7 (4.1)</td></tr><tr><td>Coxiella burnetii</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td></tr><tr><td>Mansonia spp.</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td></tr><tr><td>Piroplasmida sp.</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td><td>-</td><td>-</td><td>0 (0)</td></tr></tbody></table>

Abbreviations: n, number of samples by species; (%), percentage positive; –, not realized
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Maximum likelihood phylogeny of *Trypanosoma theileri* and *Kinetoplastida* sp. detected on tabanids in the present study. The evolutionary history, based on the ITS1 region, was inferred by using the maximum likelihood method based on the Tamura 3-parameter model. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. The analysis involved 17 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 206 positions in the final dataset</figcaption></figure>

activity of certain tabanid species in Senegal, reached the same conclusions. In addition, the haematophagous nature of females with a preference to draw their meals on humans or animals, particularly domestic ones, could also explain these findings [30].

According to the results, there are more samples and greater species diversity in Casamance than in Niokolo-Koba, probably due to the lengths of the capture period (one week in Casamance and one day in Niokolo-Koba) and the capture methods. We should also consider the park of Niokolo-Koba, a wildlife reserve, with low human density. In agreement with Ovazza [30], the uniqueness of the specimen of *T. taeniola* we captured could be explained by the fact that this species is more frequently harvested by day than by night on cattle and humans.

In the present study, four species were morphologically characterized, i.e. *A. fuscipes*, *C. distingutipennis*, *T. taeniola* and *T. guineensis*. Hamon et al. in 1956 [31], Ovazza in 1967 [30] and Raymond et al. in 1980 [22] reported in separate works their existence in Senegal. The description of specimens at the species level performed in our study was consistent with that of Ricardo [4] and Crosskey [28] for *T. guineensis*, Mugasa et al. [7] for *C. distingutipennis*

and Baldacchino et al. [2] and Hamon et al. [31] and Al-Talatha et al. [1] for *A. fuscipes*.

The molecular identification revealed three species, i.e. *A. fuscipes*, *C. distinctipennis* and *T. guineensis*. To date, there are no sequences of *A. fuscipes* in the GenBank database. However, phylogenetic analysis involving *cox1* sequences of several *Atylotus* species has shown that *A. fuscipes* is distinct from *A. diurnus* and *A. nigromaculatus* described in East Africa. It exhibited 99% identity with them. To our knowledge, this is the first study where genomic sequences were obtained for *A. fuscipes*. It seems that morphological differences in the genus *Atylotus* are not supported by genetic differences.

The sequences of *C. distingutipennis* had up to 98% identity with those in GenBank (KX946519) [7]. BLAST analysis of *T. guineensis* sequences against those described by Mugasa et al. [7] in Kenya (GenBank: KX946537, KX946566) showed an identity of only 93%. Phylogenetic analysis revealed two distinct clades, confirming that they most likely represent two genetically different species. Based on arguments from the literature and phylogeny, it appears that our samples are more likely to be authentic *T. guineensis* than species captured
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Molecular phylogenetic tree showing the position of the *Leishmania donovani* complex identified in the present study according to sequences from GenBank. The evolutionary history, based on the ITS1 region, was inferred by using the maximum likelihood method based on the Tamura 3-parameter model. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. The analysis involved 19 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were approximately 300 positions in the final dataset</figcaption></figure>

in Kenya, as the holotypes and lectotypes described in 1908 by Ricardo [4] and again in 1963 by Crosskey [28] originated in West Africa (Nigeria and Guinea). In addition, according to Oldroyd [24], the distribution of this species is limited exclusively to West Africa.

Proteomic analysis by MALDI-TOF MS made it possible to determine, to the best of our knowledge, for the first time, specific and reproducible spectra of four tabanid species, *A. fuscipes*, *C. distingutipennis*, *T. guineensis* and *T. taeniola*. Differential peaks were found to be distinct by species (Fig. 7). However, not without being unusable, *T. taeniola* spectra were slightly defective. The uniqueness of the specimen collected made it impossible to repeat the tests.

Since this study is the first proteomics analysis of tabanids by MALDI-TOF MS, it will not be possible to compare its results here. However, the data provided in our study may establish a foundation for future work concerning the identification of tabanid flies by MALDI-TOF MS.

The analysis of the blood meals revealed that *H. sapiens* is the principal host of the tabanids from Casamance (*A. fuscipes*, *C. distinctipennis* and *T. taeniola*). The hosts

of the tabanids from Niokolo-Koba (*T. guineensis*) were *H. sapiens* (96.42%) and *C. sabeus* (3.57%). These two host types occupied different ecosystems, specifically the village agglomerations of Casamance for the first and Niokolo-Koba Park for the second. In the present study, humans were identified as the main source of blood meals for the tabanid females in Senegal. The host preference cannot be established conclusively, and we suggest that this finding is most likely due to sampling sites. Ovazza [30] also found similar results. The exploration of other sites can reveal more hosts. However, we may not exclude that the identification of humans as associated hosts may be influenced by some possible sample contamination during collection.

In this study, *C. sabeus*, a non-human primate, was characterized as a host of tabanid flies for the first time. Previously, in a study conducted in Congo, Gouteux et al. [28] found non-human hosts, including the hippopotamus, ruminants and lizards. These results suggest possible transmissions by tabanids of pathogens harboured by these hosts.

Pathogens detected in *A. fuscipes* belonged to the Kinetoplastida: *T. theileri*, a high percentage of the *L.*
<figure><img src="image_7.png" /><figcaption>**Fig. 5** Phylogenetic analysis showing the position of *Setaria digitata* detected in the present study. The evolutionary history based on 18S rRNA gene was inferred using the neighbour-joining method. The evolutionary distances were computed using the Tamura-Nei method and are in units of the number of base substitutions per site. The differences in the composition bias among sequences were considered in evolutionary comparisons. The analysis involved 24 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 708 positions in the final dataset</figcaption></figure>

*donovani* complex; and an unknown Kinetoplastida sp. The obtained 28S rDNA sequences of *T. theileri* exhibited 97% similarity with *T. minasense* (GenBank: AB362411). This finding is due to the absence of 28S rDNA sequences for *T. theileri* on GenBank. *Trypanosoma theileri* has been confirmed by amplifying the 18S gene and the ITS1 region. The existence of the *L. donovani* complex in tabanids has already been mentioned [32], but its involvement in the transmission cycle and the occurrence of disease has not been determined. Our study could not provide an explanation for the strong preponderance of *Leishmania* among Senegalese tabanids. This molecular detection could also be attributed to the fact that tabanids feed on infected animals/humans and cannot maintain pathogen multiplication [33]. However, visceral leishmaniosis (VL) cases reported in West Africa are rare compared to other African endemic regions, such as north and eastern Africa. In Senegal, high seroprevalence of canine leishmaniasis (> 40%) confirmed by PCR and more than 30% seropositive people have been reported [34]. Niger reported 21 cases, Côte d'Ivoire reported eight cases, and Gambia and Burkina Faso reported one case each.

Most of these cases occurred more than 10 years ago; one Togolese case was diagnosed in Lama-Kara (northern region) in 1994; and in Nigeria, suspected cases were reported between 1936 and 1947 [35], and 60 cases were reported between 2005 and 2012 [36]. A recent survey in Ivorian dogs showed a seroprevalence of 15%, and *L. infantum* DNA was present in 4% of dogs [37]. On the other hand, VL vectors are not described in the literature. *Phlebotomus orientalis* and *P. alexandri*, the confirmed vectors of *L. donovani* and *L. infantum*, respectively, have not been identified in Senegal, Côte d'Ivoire, Mali, and Burkina Faso [38–40]. Experimental models are required to clarify this ambiguity concerning the *Leishmania*-vectorial role of tabanids.

Phylogenetic analysis based on the ITS1 region of the detected trypanosome sequences showed that the species named Kinetoplastida sp. was most likely a new species of *Trypanosoma* sp. (Fig. 3). Studies had previously described *T. theileri* as pathogens in tabanids [9, 41–43]. Similarly, the presence of *T. theileri* in Casamance was noted by Touré [26] in a study describing this species. A high (25%) prevalence of *T. theileri* in the
<figure><img src="image_8.png" /><figcaption>**Fig. 6** Molecular phylogenetic analysis for *Rickettsia* spp. detected on Senegalese tabanids. The evolutionary history based on the *gltA* gene was inferred by using the maximum likelihood method based on the Tamura 3-parameter model. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. The analysis involved 17 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 441 positions in the final dataset</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 7** MALDI-TOF specific spectra for four species of Senegalese tabanids</figcaption></figure>
studied horse flies was detected in several West African countries [44]. Investigating trypanosomes in Senegal is of great interest in the prevention and control of human and especially animal trypanosomiasis in view of the size of the livestock available in this country. Other protozoans, such as piroplasmids, were not detected in the present study.

A nematode belonging to the family Setaridae was detected in *T. guineensis* and *A. fuscipes*. It is well-known that filariae, such as *Loa loa*, are transmitted by tabanids [9, 43], but *S. digitata*, mainly recognized as the aetiological agent of animal filariasis, has not been investigated in tabanids. The present study demonstrates the existence of a new genus of the Tabanidae, including *Atylotus* (*A. fuscipes*), as potential carriers of wireworms. To date, only *Chrysops* spp. were known as worm vectors in the Tabanidae. The revelation concerning the presence of filariae in Casamance and Niokolo-Koba provides important data for understanding the many filariases that affect human populations and their animals in Senegal. This research is part of the government's efforts to prevent and control locusts. In the past, only two species of *Onchocerca* spp. has been reported in the cattle of Kedougou by Vassiliades et al. [21]. Finally, in the present study, positive tabanids for the different pathogens had *H. sapiens* as the principal host.

We detected bacteria belonging to the family *Anaplasmataceae* and to the genus *Rickettsia*. We were not able to amplify and characterize *Anaplasmataceae*-positive samples, probably due to the low bacterial load reflected by the high cycle threshold (34–36) found by qPCR. Several bacterial species have been detected in tabanids [2, 45]. *Tabanus* spp. and *Chrysops* spp. have been reported as mechanical vectors of *Anaplasma marginale*, while *C. vittatus* and *Tabanus* spp. were recognized as vectors for *Ehrlichia risticii* [2]. In this study, DNA of *Anaplasmataceae* was detected in *A. fuscipes*, *C. distinctipennis* and *T. guineensis*. In addition, three *Rickettsia* genotypes have been detected in *A. fuscipes*. Unfortunately, we could not better characterize the identified rickettsial species because of the low bacterial load in the samples studied. Senegal, however, is an endemic area for *Anaplasmatacea* [45] and *Rickettsia* spp. [46] infections. In the same study area, Mediannikov et al. [47] detected a new *Rickettsia* sp. in tsetse flies. *Anaplasma marginale* and rickettsial endosymbionts have been molecularly identified in tabanids from Hungary [48]. The role of tabanids as vectors of *Anaplasmatacea* and/or *Rickettsia* spp. has not been characterized to date. In addition, no *Coxiella burnetii* DNA was detected, while *C. burnetii* has already been detected in *T. staegeri* and was not captured in the present study [2].

## Conclusions

Previously overlooked and therefore insufficiently examined, tabanids from Senegal in particular and from West Africa in general, are nevertheless potential carriers of multiple pathogens. Investigating these dipterans using innovative tools and genomic and proteomics methods is of considerable interest to animal and human public health. Several limitations to this research are notable, particularly with regard to the duration and uniqueness of some specimens characterized in this study. However, this study represents a strong contribution to the understanding and prevention of vector-borne diseases caused by tabanids in West Africa. Furthermore, this report is the first to describe the diversity among Senegalese tabanids, their hosts and associated pathogens. The molecular results and sequences of *A. fuscipes* and *T. guineensis* and the presence of *L. donovani*, *Trypanosoma* spp., *S. digitata* and *Rickettsia* spp. are notable. To thoroughly characterize these species, it is important that studies in entomology and infectiology be pursued to achieve the molecular characterization and MALDI-TOF MS analysis of tabanids from Senegal and West Africa. These efforts, if realized, may enrich the genomic and proteomic databases for tabanids. Screening for other pathogens not studied in this research is equally important. Finally, experimental models should be considered to better understand the ability and competence of Senegalese tabanids to carry or transmit the pathogens detected in this study.

## Abbreviations

MALDI-TOF MS: Matrix assisted laser desorption ionization-time of flight mass spectrometry; PCR: Polymerase chain reaction.

## Acknowledgements

We are grateful for Jean-Michel Berenger for providing technical assistance in laboratory entomological work, Cheikh Sokhna for his help in field mission organization, Maxence Aubadie-Ladrix and Alpha Kabinet Keita for their help in arthropod collection, and Adama Zan Diarra and Maureen Laroche for their help.

## Authors' contributions

MLK, HM and OM designed the study. MS and OM collected the samples. HM and MLK carried out the data analysis. MLK, HM and HD performed the laboratory work. HM, MLK and HD drafted the manuscript. OM reviewed and edited the manuscript. All authors read and approved the final manuscript.

## Funding

This study was supported by the Institut Hospitalo-Universitaire (IHU) Méditerranée Infection, the National Research Agency under the programme «Investissements d'avenir», reference ANR-10IAHU-03, the Région Provence-Alpes-Côte d'Azur and European funding FEDER PRIMI.

## Availability of data and materials

All data are included in the manuscript. The newly generated sequences were deposited in the GenBank database under the accession numbers: MW013783-MW013793 (cox1) of tabanids; MN244142-MN244147 (ITS1); MN244154-MN244159 (28S) and MN244160-MN244163 (18S) of *Trypanosma* spp., MN244150-MN244153 (ITS1) of *Leishmania* spp., MN244148-MN244149 (18S) of *Setaria digitata*.
**Ethics approval and consent to participate**

Not applicable.

**Consent for publication**

Not applicable.

**Competing interests**

The authors declare that they have no competing interests.

**Author details**

¹ IHU Méditerranée Infection – Microbes, Evolution, Phylogeny and Infection (MEΦI), Marseille, France. ² UMR Aix-Marseille University, IRD, APHM -19-21, Bd Jean Moulin, 13385 Marseille Cedex 05, France. ³ IHU Méditerranée Infection - Vecteurs - Infections Tropicales et Méditerranéennes (VITROME), Marseille, France. ⁴ Vectors-Tropical and Mediterranean Infections (VITROME), Campus International, UCAD-IRD, Dakar, Sénégal.

Received: 3 February 2020 Accepted: 23 September 2020

Published online: 01 October 2020

**References**

1. Al-Talatha HA, Yaakop S, Idris AB. Two new species and seven new records of horse fly (Diptera: Tabanidae) from Malaysia, including a description of new species and modified keys. J Med Entomol. 2018;55:112-21.

2. Baldacchino F, Desquesnes M, Mihok S, Foil LD, Duvallet G, Jittapalapong S. Tabanids: neglected subjects of research, but important vectors of disease agents! Infect Genet Evol. 2014;28:596-615.

3. Changbunjong T, Bhusri B, Sedwisai P, Weluwanarak T, Nitiyamawat E, Chareonviriyaphap T, et al. Species identification of horse flies (Diptera: Tabanidae) in Thailand using DNA barcoding. Vet Parasitol. 2018;259:35-43.

4. Ricardo G. Descriptions of thirty new species of Tabani from Africa and Madagascar. Ann Mag Nat Hist. 1908;1:268-78.

5. Banerjee D, Kumar V, Maity A, Ghosh B, Tyagi K, Singha D, et al. Identification through DNA barcoding of Tabanidae (Diptera) vectors of surra disease in India. Acta Trop. 2015;150:52-8.

6. Cywinska A, Hannan MA, Kevan PG, Roughley RE, Iranpour M, Hunter FF. Evaluation of DNA barcoding and identification of new haplomorphs in Canadian deerflies and horseflies. Med Vet Entomol. 2010;24:382-410.

7. Mugasa CM, Villinger J, Gitau J, Ndungu N, Ciosi M, Masiga D. Morphological re-description and molecular identification of tabanidae (Diptera) in East Africa. ZooKeys. 2018;2018:117-44.

8. Foil LD. Tabanids as vectors of disease agents. Parasitol Today. 1989;5:88-96.

9. Krinsky WL. Animal disease agents transmitted by horse flies and deer flies (Diptera: Tabanidae). J Med Entomol. 1976;13:225-75.

10. Wilson JJ. DNA barcodes for insects. Methods Mol Biol. 2012;858:17-46.

11. Taioe MO, Motloang MY, Namangala B, Chota A, Molefe NI, Musinguzi SP, et al. Characterization of tabanid flies (Diptera: Tabanidae) in South Africa and Zambia and detection of protozoan parasites they are harbouring. Parasitology. 2017;144:1-17.

12. Gouteux JP, Noireau F, Staak C. The host preferences to *Chrysops silacea* and *C. dimidiata* (Diptera: Tabanidae) in an endemic area of *Loa loa* in the Congo. Ann Trop Med Parasitol. 1989;83:167-72.

13. Goldstein PZ, DeSalle R. Integrating DNA barcode data and taxonomic practice: determination, discovery, and description. Bioessays. 2011;33:135-47.

14. Kress WJ, Erickson DL. DNA barcodes: methods and protocols. Methods Mol Biol. 2012;858:3-8.

15. Chakraborty C, Doss CGP, Patra BC, Bandyopadhyay S. DNA barcoding to map the microbial communities: current advances and future directions. Appl Microbiol Biotechnol. 2014;98(8):3425-36.

16. Nebbak A, El Hamzaoui B, Berenger JM, Bitam I, Raoult D, Almeras L, et al. Comparative analysis of storage conditions and homogenization methods for tick and flea species for identification by MALDI-TOF MS. Med Vet Entomol. 2017;31:438-48.

17. Nebbak A, Willcox AC, Bitam I, Raoult D, Parola P, Almeras L. Standardization of sample homogenization for mosquito identification using an innovative proteomic tool based on protein profiling. Proteomics. 2016;16:3148-60.

18. Diarra AZ, Almeras L, Laroche M, Berenger JM, Koné AK, Bocoum Z, et al. Molecular and MALDI-TOF identification of ticks and tick-associated bacteria in Mali. PLoS Negl Trop Dis. 2017;24(11):e005762.

19. Bouyer F, Bouyer J, Seck MT, Sall B, Dicko AH, Lancelot R, et al. Importance des infections à transmission vectorielle selon le système de production: trypanosomose bovine et dynamique d'innovation des éleveurs au Sénégal. Rev Sci Tech. 2015;34:199-225.

20. Bassene H, Sambou M, Fenollar F, Clarke S, Djiba S, Mourembou G, et al. High prevalence of *Mansonia perstans* filariasis in rural Senegal. Am J Trop Med Hyg. 2015;93:601-6.

21. Vassiliades G, Delbove P, Bain O. Onchocercoses bovines au Sénégal. Note préliminaire. Rev Elev Med Vet Pays Trop. 1983;36:351-3.

22. Raymond HL, Taufflieb R, Cornet M, Camicas J-L, Chateau R, Dieng PY. Annotated list of Tabanidae (Diptera) of Senegal and the Gambia. Bull Inst Fondam Afr Noire. 1980;42:812-22.

23. Oldroyd H. The horse-flies (Diptera: Tabanidae) of the Ethiopian region. London: The Trustees of the British Museum; 1952.

24. Oldroyd H. The horse-flies (Diptera: Tabanidae) of the Ethiopian region: Vol. II. *Tabanus* and related genera. https://catalogue.nla.gov.au/Recor d/1108859.

25. Ovazza M, Valade M. Contribution to the study of Tabanidae (Diptera) in French West Africa. II. On two species of *Tabanus*: *T. hamoni* n. sp. and *T. sugens* Wiedemann, 1828. Bull Soc Pathol Exot Fil. 1958;51:992-8.

26. Touré S. Description complémentaire de Trypanosoma theileri Laveran, 1902 Mention particulière des formes observées en Casamance (Sénégal). Rev Elev Med Vet Pays Trop. 1968;21:365-73.

27. Yssouf A, Socolovschi C, Flaudrops C, Ndiath MO, Sougoufara S, Dehecq JS, et al. Matrix-assisted laser desorption ionization-time of flight mass spectrometry: an emerging tool for the rapid identification of mosquito vectors. PLoS ONE. 2013;8:e72380.

28. Crosskey RW. Lectotype selection for *Tabanus guineensis* Wiedemann, 1824, and neotype designation for *Tabanus sugens* Wiedemann, 1828 (diptera: Tabanidae). Proc R Entomol Soc London Ser B, Taxon. 2009;30:79-82.

29. Abdullah HHAM, El-sMolla A, Salib FA, Ghazy AA, Allam NAT, Sanad YM, Abdel-Shafy S. Molecular characterization of *Rickettsiae* infecting camels and their ticks vectors in Egypt. Anim Husb Dairy Vet Sci. 2018;2:2-10.

30. Ovazza M. Observations sur l'activité nocturne de certaines espèces de Tabanidae (Diptera), en savane d'Afrique de l'Ouest. Cah ORSTOMSérie Entomol Méd Parasitol. 1967;5:53-61.

31. Hamon J, Morel J, Ovazza M, Rickenbach A. Tabanidae (Diptera) of French West Africa. Ann Parasitol Hum Comp. 1956;31:436-48.

32. Coelho WMD, Buzetti WAS, Bresciani KDS. Histochemical and molecular evaluation of the prevalence of *Leishmania* spp. in hematophagous insects. Parasite Epidemiol Control. 2016;1:85-9.

33. Baldacchino F, Muenworn V, Desquesnes M, Desoli F, Charoenviriyaphap T, Duvallet G. Transmission of pathogens by *Stomoxys* flies (Diptera, Muscidae): a review. Parasite. 2013;20:26.

34. Faye B, Bañuls AL, Bucheton B, Dione MM, Bassanganan O, Hide M, et al. Canine visceral leishmaniasis caused by *Leishmania infantum* in Senegal: risk of emergence in humans? Microbes Infect. 2010;12:1219-25.

35. Desjeux P, WHO. Division of Control of Tropical Diseases. Information sur l'épidémiologie des leishmanioses et la lutte contre ces maladies par pays ou territoire. Geneva: World Health Organization; 1991. https://apps.who.int/iris/handle/10665/61576.

36. Kone AK, Niaré DS, Piarroux M, Izri A, Marty P, Laurens MB, et al. Visceral leishmaniasis in West Africa: clinical characteristics, vectors, and reservoirs. J Parasitol Res. 2019;2019:9282690.

37. Medkour H, Laidoudi Y, Athias E, Bouam A, Dizoé S, Davoust B, et al. Molecular and serological detection of animal and human vector-borne pathogens in the blood of dogs from Côte d'Ivoire. Comp Immunol Microbiol Infect Dis. 2020;69:101412.

38. Berdjane-Brouk Z, Koné AK, Djimdé AA, Charrel RN, Ravel C, Delaunay P, et al. First detection of *Leishmania major* DNA in *Sergentomyia* (*Spe- laemyia*) darlingi from cutaneous leishmaniasis foci in Mali. PLoS ONE. 2012;7:e28266.
39. Sangare L, Gantier JC, Koalaga G, Deniau M, Ouari A, Guiguemdé RT. Sandflies of the South part of Ouagadougou City. Burkina Faso. *Parasite*. 2009;16:231–3.

40. Senghor MW, Faye MN, Faye B, Diarra K, Elguero E, Gaye O, et al. Ecology of phlebotomine sand flies in the rural community of Mont Rolland (Thiès region, Senegal): area of transmission of canine leishmaniasis. PLoS ONE. 2011;6:e14773.

41. Böse R, Friedhoff KT, Olbrich S, Büscher G, Domeyer I. Transmission of *Trypanosoma theileri* to cattle by Tabanidae. Parasitol Res. 1987;73:421–4.

42. Desquesnes M, Dia ML. Mechanical transmission of *Trypanosoma congolense* in cattle by the African tabanid *Atylotus agrestis*. Exp Parasitol. 2003;105:226–31.

43. Salim B, Bakheit MA, Kamau J, Nakamura I, Sugimoto C. Molecular epidemiology of camel trypanosomiasis based on ITS1 rDNA and RoTat 1.2 VSG gene in the Sudan. Parasit Vectors. 2011;4(4):31.

44. Votýpka J, Brzoňová J, Ježek J, Modrý D. Horse flies (Diptera: Tabanidae) of three West African countries: a faunistic update, barcoding analysis and trypanosome occurrence. Acta Trop. 2019;197:105069.

45. Dahmani M, Davoust B, Sambou M, Bassene H, Scandola P, Ameur T, et al. Molecular investigation and phylogeny of species of the *Anaplasma*-ceae infecting animals and ticks in Senegal. Parasit Vectors. 2019;12:495.

46. Mediannikov O, Diatta G, Fenollar F, Sokhna C, Trape JF, Raoult D. Tick-borne rickettsioses, neglected emerging diseases in rural Senegal. PLoS Negl Trop Dis. 2010;14:e821.

47. Mediannikov O, Audoly G, Diatta G, Trape JF, Raoult D. New *Rickettsia* sp. in tsetse flies from Senegal. Comp Immunol Microbiol Infect Dis. 2012;35:145–50.

48. Hornok S, Földvári G, Elek V, Naranjo V, Farkas R, de la Fuente J. Molecular identification of *Anaplasma marginale* and rickettsial endosymbionts in blood-sucking flies (Diptera: Tabanidae, Muscidae) and hard ticks (Acari: Ixodidae). Vet Parasitol. 2008;154:354–9.

49. Folmer O, Black MB, Hoeh WR, Lutz R, Vrijenhoek RC. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3:294–9.

50. Kocher TD, Thomas WK, Meyer A, Edwards SV, Paabo S, Villablanca FX, et al. Dynamics of mitochondrial DNA evolution in animals: amplification and sequencing with conserved primers. Proc Natl Acad Sci USA. 1989;86:6196–200.

51. Mary C, Faraut F, Lascombe L, Dumon H. Quantification of *Leishmania infantum* DNA by a real-time PCR assay with high sensitivity. J Clin Microbiol. 2004;42:5249–55.

52. Ogado Ceasar Odiwuor S, Ageed Saad A, De Doncker S, Maes I, Laurent T, El Safi S, et al. Universal PCR assays for the differential detection of all Old World *Leishmania* species. Eur J Clin Microbiol Infect Dis. 2011;30:209–18.

53. Dahmana H, Amanzougaghene N, Davoust B, Normand T, Carette O, Demoncheaux JP, et al. Great diversity of *Piroplasmina* in Equidae in Africa and Europe, including potential new species. Vet Parasitol Reg Stud Reports. 2019;18:100332.

54. Njiru ZK, Constantine CC, Guya S, Crowther J, Kiragu JM, Thompson RCA, et al. The use of ITS1 rDNA PCR in detecting pathogenic African trypanosomes. Parasitol Res. 2005;95:186–92.

55. Medkour H, Varloud M, Davoust B, Mediannikov O. New molecular approach for the detection of kinetoplastida parasites of medical and veterinary interest. Microorganisms. 2020;8:356.

56. Mediannikov O, Fenollar F, Socolovschi C, Diatta G, Bassene H, Molez JF, et al. *Coxiella burnetii* in humans and ticks in rural Senegal. PLoS Negl Trop Dis. 2010;4:e654.

57. Dahmani M, Loudahi A, Mediannikov O, Fenollar F, Raoult D, Davoust B. Molecular detection of *Anaplasma platys* and *Ehrlichia canis* in dogs from Kabylie. Algeria. Ticks Tick Borne Dis. 2015;6:198–203.

58. Parola P, Diatta G, Socolovschi C, Mediannikov O, Tall A, Bassene H, et al. Tick-borne relapsing fever borreliosis, rural senegal. Emerg Infect Dis. 2011;17:883–5.

59. Mediannikov OY, Sidelnikov Y, Ivanov L, Mokretsova E, Fournier PE, Tarasevich I, et al. Acute tick-borne rickettsiosis caused by *Rickettsia heilongjiangensis* in Russian Far East. Emerg Infect Dis. 2004;10:810–7.

60. Laidoudi Y, Davoust B, Varloud M, Niang EHA, Fenollar F, Mediannikov O. Development of a multiplexed qPCR-based approach for the diagnosis of *Dirofilaria immitis*, *D repens*, *Acanthocheilonema reconditum* and the others filariosis. Parasit Vectors. 2020;13:319.

## Publisher's Note

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

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_10.png" /></figure>