Study sites
The study was carried out in three active foci situated in the south and centre of the Republic of Congo.
- Mpouya (2°36'57''S, 16°12'43''E): this foci is located in the Plateau administrative division about 300 km north of Brazzaville. The foci is compose of several villages bordering the Congo river. Its vegetation is essentially made of savannah with sparse vegetation with few high trees. The population mainly practice fishing and farming for living.
- Ngabé (3°12'52''S, 16°10'1''E) is situated about 200 km north of Brazzaville in the Pool region. This foci is also compose of several villages bordering river Congo. The vegetation is made of savannah and forest gallery along the edges of rivers. The population mainly practice farming, hunting and fishing for living.
- Loudima (4°6'45''S, 13°3'30''E) is located in the Bouenza division about 300km south of Brazzaville. The vegetation is dominated by grassfield. Main activities practiced by the population is farming.
In each foci tsetse flies sampling was conducted in four to five villages.
(Figure 1)
Tsetse flies sampling
Entomological surveys were conducted in November 2019, November 2020 (rainy season) and June 2021 (dry season). During each survey, pyramid traps[10]were set up in five villages in Ngabe, four villages in Mpouya and Loudima for at least three days in different places (close to the river or swamps, shaded areas and in the village). A total of 251 traps were set, including 106 in November 2019, 94 in November 2020 and 51 in june 2021. For each trap, the geographical coordinates were recorded using a GARMIN GPS (eTrex® 10). The minimum distance between traps was 100m. The traps were visited twice daily (10 A.M. and 4 P.M.). Once the tsetse flies were captured, they were identified morphologically to determine the sex and species[11] and then preserved in microtubes containing alcohol. In the field, these microtubes were kept in a refrigerated cooler (at +4°C) and in the laboratory at -20°C until DNA extraction.
DNA extraction
DNA was extracted from whole tsetse fly by the cethyl trimethyl ammonium bromide (CTAB) method[12]. In the laboratory, the alcohol was evaporated at room temperature and 600 µl of CTAB buffer (5% CTAB; 1 M Tris, pH 8; 0.5 M EDTA, pH 8; 5 M NaCl) was added, the tsetse fly was ground with a pestle and then incubated in a water bath at 60 °C for 30 min. Then, 600 µl of chloroform/isoamyl alcohol (24/1) was added to each tube. This mixture was slowly homogenised for 15 min and then centrifuged at 13 000 rpm for 15 min. The upper aqueous phase was removed and transferred to another 1.5 mL microtube. The DNA was precipitated by the addition of 600 µl isopropanol. After gentle homogenisation of each microtube for 5 min and incubation overnight at -20 °C, each microtube was centrifuged at 13,000 rpm for 15 min. The DNA pellet was then washed twice with cold 70% ethanol and dried overnight at room temperature. The resulting DNA pellet was resuspended in 50 µL of sterile nuclease-free water and stored at -20 °C or used directly for PCR.
Molecular identification of tsetse fly species
The principle of this PCR is based on the amplification of the ITS1 region ("internal transcribed spacer 1") of the tsetse genome, using the primer pair "Diag forward" (5'-TGGACTTCGGATTAAGTACAACA-3') and "Diag reverse" (5'-TCATTATGCGCTATTAAGGTAAGC-3')[13]. For this identification, only one PCR cycle was performed; the cycle was performed in a final volume of 25 µl containing 1× PCR buffer (10 mM Tris-HCl (pH 9.0), 50 mM KCl), 3 mM MgCl2, 1 µl (15 pmol) of each primer, 0.5 μl (200 mM) of dNTPs, 0.1 μl (one unit) of Taq DNA polymerase, 2.5 µl of DNA and 17.4 μl of sterile water. The amplification programme began with a denaturation step at 95 °C for 3 min and 30 s, followed by 30 cycles of amplification; each cycle contained a denaturation step at 95 °C for 30 s, a hybridisation step at 56 °C for 1 min and an extension step at 72 °C for 1 min, followed by a final extension step at 72 °C for 5 min.
The amplification products were run on a 2% agarose gel to enable species identification.
Molecular identification of different trypanosome species
The identification of trypanosomes was performed by amplification of the internal transcribed spacer 1 (ITS1) of the ribosomal DNA of different trypanosome species, as described by Ravel et al [14]. For this identification, two PCR cycles were performed; the first cycle was performed in a final volume of 25 µl containing 1× PCR buffer (10 mM Tris-HCl (pH 9. 0), 50 mM KCl), 2 mM MgCl2, 1 µl (10 pmol) of each primer (5′-CAA ATT GCC CAA TGT CG-3′ and 5′-GCT GCG TTC TTC AAC GAA-3′), 0. 5 µl (200 mM) of dNTPs, 1 µl (one unit) of Taq DNA polymerase (5 U/µl), 5 μl of DNA and 14 µl of nuclease-free water. The amplification programme began with a denaturation step at 94 °C for 3 min and 30 s, followed by 30 cycles of amplification; each cycle contained a denaturation step at 94 °C for 30 s, a hybridisation step at 58 °C for 1 min and an extension step at 72 °C for 1 min, followed by a final extension step at 72 °C for 5 min.
The amplified products from the first round of PCR were diluted 10-fold and 3 µl of each dilution was used as template for the second round of PCR. The second round of PCR was performed with two different primers (5′-CCT GCA GCT GGA TCA T-3′ and 5′-ATC GCG ACA CGT TGT G-3′). The amplification program was identical to that of the first round of PCR. After nested PCR, the amplicons were separated by electrophoresis on a 2% agarose gel which was then stained with midori green advance and visualised under UV light.
Trypanosomes species were identified on the basis of the length of ITS1 fragments. For example, T. congolense strains generate DNA fragments of about 650 bp (630 bp for T. congolense forest and 610 bp for T. congolense savannah) whereas fragments of about 150 bp and 400 bp, respectively, are expected for T. vivax and all trypanosomes belonging to the Trypanozoon subgenus (T. brucei (s.l.), T. evansi and T. equiperdum).
Identification of Trypanosoma congolense forest and Trypanosoma congolense savannah
After amplification of the ITS1 sequences, all samples that had a DNA fragment between 600-650 bp, corresponding to the expected size of T. congolense, were subjected to a further PCR where specific primers were used to separate T. congolense forest to T. congolense savannah. These specific identifications were performed as described by Simo et al. [15]using primers TCF1 (5′-GGA CAC ACG CCA GAA GGT ACT T-3′) and TCF2 (5′-GTT CTC TCG CAC CAA ATC CAA C-3′) for T. congolense forest "type"[16], and TCS1 (5′-CGA GCG AGA ACG GGC AC-3′) and TCS2 (5′-GGG ACA AAC AAA TCC CGC-3′) for T. congolense savannah "type"[17]. PCR reactions were performed in a final volume of 25 μl containing 1× PCR buffer (10 mM Tris-HCl (pH 9.0), 50 mM KCl), 3 mM MgCl2, 1 µl (15 pmol) of each primer, 0.5 μl (200 mM) of dNTPs, 1 μl (one unit) of Taq DNA polymerase, 3 µl of DNA and 16 μl of sterile water. The amplification program consisted of a denaturation step at 94°C for 3 min 30 s, followed by 40 cycles of amplification including a denaturation step at 94°C for 30 s, a hybridisation step at 60°C for 1 min and an elongation step at 72°C for 1 min, followed by a final elongation step at 72°C for 5 min.
The amplified products were separated by electrophoresis on a 2% agarose gel containing midori green advance. The DNA bands were visualised under ultraviolet (UV) light and a picture of the gel was taken and stored.
Data analysis
The apparent density of tsetse flies was calculated using the following formula ADT = C/(T*D) where C represents the number of flies collected, T the number of traps used and D the number of days of trapping. The Anova was used to compare the ADT between foci. Statistical analysis were conducted using the software MedCalc[18]. The infection rate of tsetse flies between sites was compared using the pearson chi-square test. The threshold for significance was P < 0.05.