Onchocerca volvulus transmission in the Mbam valley of Cameroon following 16 years of annual community-directed treatment with ivermectin, and the description of a new cytotype of Simulium squamosum
Parasites & Vectors volume 14, Article number: 563 (2021)
The onchocerciasis focus surrounding the lower Mbam and Sanaga rivers, where Onchocerca volvulus is transmitted by Simulium damnosum s.l. (Diptera: Simuliidae), was historically the largest in the southern regions of Cameroon. Annual community-directed treatment with ivermectin (CDTI) has been taking place since 2000, but recent studies have shown that new infections are occurring in children. We aimed to investigate blackfly biting and O. volvulus transmission rates along the lower Mbam river 16 years after the formal onset of annual CDTI.
Black flies were collected for three consecutive days each month between July 2016 and June 2017 at two riverside villages and two inland sites situated 4.9 km and 7.9 km from the riverside. Specimens collected at each site were dissected on one of the three collection days each month to estimate parity rates and O. volvulus infection rates, while the remaining samples were preserved for pool screening.
In total, 93,573 S. damnosum s.l. black flies were recorded biting humans and 9281 were dissected. Annual biting rates of up to 606,370 were estimated at the riverside, decreasing to 20,540 at 7.9 km, while, based on dissections, annual transmission potentials of up to 4488 were estimated at the riverside, decreasing to 102 and 0 at 4.9 km and 7.9 km, respectively. However, pool screening showed evidence of infection in black flies at the furthest distance from the river. Results of both methods demonstrated the percentage of infective flies to be relatively low (0.10–0.36%), but above the WHO threshold for interruption of transmission. In addition, a small number of larvae collected during the dry season revealed the presence of Simulium squamosum E. This is the first time S. squamosum E has been found east of Lake Volta in Ghana, but our material was chromosomally distinctive, and we call it S. squamosum E2.
Relatively low O. volvulus infection rates appear to be offset by extremely high densities of biting black flies which are sustaining transmission along the banks of the lower Mbam river.
Onchocerciasis was hyperendemic in five main foci in Cameroon before community-directed treatment with ivermectin (CDTI) was introduced to control the disease in 1997 [1, 2]. The foci extended across savannahs in the northern regions, through forest–savannah transition zones in the centre, and into dense humid forests further south [1, 3]. Fast-flowing rivers running through these areas provide suitable breeding habitats for Simulium damnosum (Diptera: Simuliidae) complex black flies, the only vectors of Onchocerca volvulus in the country [4, 5].
Eight cytoforms (both cytospecies and cytotypes) of the S. damnosum complex have been described from Cameroon , and their geographical distribution has been well mapped [4, 6]. Simulium damnosum sensu stricto (s.s.) and Simulium sirbanum are cytospecies common in savannah habitats , while Simulium yahense and Simulium squamosum are associated with forest and transitional zones [4, 6], although the latter may spread north into savannah rivers during the rainy season . Simulium squamosum is divided into five cytotypes (chromosomally distinct populations of unconfirmed taxonomic status) designated A–E, of which cytotypes A–D are present in Cameroon, while cytotype E occurs west of Lake Volta in Ghana [4, 7, 8]. Simulium squamosum A is the typical form described by Vajime and Dunbar . It is found throughout most of Cameroon except for the Sanaga river, which is the only known breeding locality of S. squamosum B . Cytotypes C and D are known from the areas around Mount Cameroon where A and C also appear to interbreed [4, 7]. Another member of the complex, Simulium mengense, has a scattered distribution throughout the north where it occupies similar habitats to S. damnosum s.s. and S. sirbanum . It is also present in rivers around Mount Cameroon and in the Centre Region, where it often breeds sympatrically with S. squamosum [10,11,12]. All S. damnosum sensu lato (s.l.) cytospecies present in Cameroon are known or suspected vectors of O. volvulus [13,14,15].
The onchocerciasis focus surrounding the lower Mbam and Sanaga rivers was historically the largest in the southern regions of Cameroon . Disease prevalence was particularly high in villages along the lower Mbam, where infection was associated with severe ocular pathologies and high rates of epilepsy [2, 11, 16]. The latter condition was found to be related to high community microfilarial loads following a study by Boussinesq et al. in 1991/1992 . Furthermore, two cohort studies in these river valleys have shown a temporal relationship with a dose effect between microfilarial load in children and the development of epilepsy [17, 18]. Prior to the implementation of CDTI, an investigation of blackfly biting rates and O. volvulus transmission was conducted by Barbazan et al. . This took place over 12 months from April 1993 until March 1994 and included two transects perpendicular to the Mbam river and close to the town of Bafia. The study showed that blackfly biting rates were high throughout the year, but O. volvulus transmission was seasonal, occurring mainly between February and May and peaking in February and March. Breeding populations of S. damnosum s.l. were restricted to the main river, where 90% of larvae were found to be S. squamosum s.s. (presumably cytotype A) and 10% were S. mengense .
There is no history of vector control along the Mbam river, where onchocerciasis is controlled solely through annual CDTI [3, 19]. The first large-scale ivermectin treatments commenced in 1994 in selected villages as part of a clinical trial to evaluate the macrofilaricidal potential of the drug [11, 20]. CDTI initially commenced in 1997 but encountered problems due to severe adverse events associated with Loa loa co-infection . Treatment coverage during the first years was consequently low, and the programme was relaunched in 2000 [2, 3]. Epidemiological studies conducted in 2011 showed that the CDTI project, Centre 1, which includes Bafia Health District, had a mean microfilarial prevalence of 52.3% based on a survey of 12 villages, and was progressing more slowly than expected towards elimination [3, 21]. After corrective measures had been made to the CDTI project, Kamga et al.  conducted a follow-up survey in 2015 which revealed a mean microfilarial prevalence of 41.6% in four villages. While this was still higher than expected, community microfilarial loads had decreased dramatically from pre-control levels, but microfilaridermia and the presence of nodules in children aged less than 10 years implied that transmission was ongoing .
Onchocerciasis is still a country-wide public health problem in Cameroon, and at the onset of this study there was evidence of ongoing O. volvulus transmission despite long-term public health interventions [3, 22]. Between July 2016 and June 2017, we investigated blackfly biting rates and seasonal patterns of O. volvulus transmission near villages surrounding the lower Mbam river, almost 23 years after the pre-control survey by Barbazan et al.  and 16 years after the formal onset of annual CDTI [3, 11].
The study was conducted in villages close to the perennially flowing Mbam river near Bafia (4.73939°N, 11.22073°E) in the Centre Region of Cameroon (Fig. 1). The Mbam originates in the northern savannah region and flows through the transitional forest–savannah mosaic surrounding Bafia before joining the Sanaga river as its main tributary . As it passes near Bafia, the river is characterised by a series of rapids that provide ideal sites for blackfly breeding . The climate is equatorial and the area receives 1700–1850 mm in annual rainfall, occurring mainly in two peaks (Fig. 2) . The lesser peak occurs between March and June which precedes a brief dry period in July and is followed by heavier rainfall between August and October. River discharge steadily increases across the rainy seasons before it abruptly decreases in November at the start of the long dry season (Fig. 2) . blackfly biting occurs throughout the year but is most intense between February and May, and September and October. Transmission of O. volvulus mainly occurs between February and May . The estimated population of Bafia Health District was 226,073 in 2014 . Many residents engage in subsistence farming, fishing, and sand mining along the Mbam river , while nomadic pastoralists (Bororo herdsmen) migrate to the area at varying times between November and May each year (M Ronse pers. comm.). Small numbers of cattle are also kept locally .
blackfly collection sites were established at or near four villages along the lower Mbam river (Table 1). The sites along the transect extending from Nyamongo I to Ondouano were chosen to enable comparison with data collected by Barbazan et al. , while the site near Bayomen was chosen due to historically high rates of epilepsy in the area . Nyamongo I and Bayomen sites were both at the riverside. The former was next to a ferry crossing approximately 1 km downstream from major rapids, while the latter was used by local fishermen and was known for its high densities of biting black flies. Collection sites at Egona II and Ondouano were located 4.9 km and 7.9 km from the riverside, respectively.
Collection and preservation of black flies
blackfly larvae and pupae were collected from rocks and trailing vegetation in large rapids 1 km upstream from the Nyamongo I ferry crossing and 300 m downstream from the collection site at Bayomen in January 2017 (dry season), when the river was most accessible (Fig. 1). Samples were placed in three changes of Carnoy’s fixative (≈3:1 ethanol: glacial acetic acid) to preserve S. damnosum complex larval polytene chromosomes for cytotaxonomic study. All pupae and non-damnosum larvae were transferred to absolute ethanol in the laboratory, where they were stored at −20 °C until needed.
Adult black flies were collected simultaneously at the four sites for three consecutive days each month, beginning in July 2016 and ending in June 2017. At each site, catches were made between 07:00 and 18:00 daily by two individuals working alternate hours who were recruited from the nearest village and trained in standard human landing collection methods . Due to anticipated high landing rates, black flies were collected using mouth aspirators which were labelled and changed hourly. Adult black flies were either dissected in the field to determine parity rates and to detect Onchocerca spp. larvae or were preserved in absolute ethanol for subsequent pool screening. Preserved specimens were maintained at ambient temperatures in the field until stored at 4 °C in the laboratory. Landing rates were interpreted as being representative of exposure to biting and are, therefore, referred to as biting rates.
Identification of S. damnosum complex
Larvae of the S. damnosum complex were morphologically identified by the presence of dorsal abdominal tubercles and scales on the prothoracic proleg . Late-instar larvae were used for cytotaxonomy. Prior to chromosome preparation, the head and thorax of each larva was separated from the remainder of the specimen and stored individually in absolute ethanol for morphological identification. Salivary glands were then dissected from the abdominal cavity of each respective specimen and were prepared following a Feulgen staining method outlined by Adler et al. . Larvae were identified with reference to the cytotaxonomic key in Post et al. , and chromosome maps in Vajime and Dunbar , Boakye , and Mustapha et al. . Inversion nomenclature follows Post et al. . Adult black flies were identified by their enlarged fore-tarsi bearing crests of dark hair and the presence of white bands on the hind basi-tarsi .
Dissection of adult black flies
Field stations were established at Bayomen and Nyamongo I (Fig. 1) to estimate S. damnosum s.l. parity rates and Onchocerca spp. infection rates at each site. To achieve this, a subsample of adult black flies (up to ≈30 per hour, and ≈330 per day) were dissected from each site on one of the three collection days each month. Black flies were exposed to chloroform vapour inside a closed box until anaesthetised. Each specimen was then placed in a drop of saline solution on a microscope slide and dissected to determine whether it was parous or nulliparous using standard methods . Nulliparous flies were discarded while parous flies were further dissected and examined for the presence of O. volvulus larvae [35, 36]. If parasites were found, the number and developmental stages (L1–L3, and L3H) were recorded. L3 larvae were defined as infective stages found anywhere within the blackfly, while L3H were those present in the head. Dissections were carried out by experienced technicians using Wild M5 stereo microscopes (Wild Heerbrugg, Switzerland) illuminated with 6v bulbs.
The objective of pool screening was to confirm the identity of the Onchocerca species present in S. damnosum s.l. and to provide a robust estimate of infection rates using a larger sample of black flies. Pools of heads and bodies (maximum pool size n = 200) of ethanol-preserved black flies were prepared separately, by site and month of collection, following methods outlined by Hendy et al. . DNA was extracted and samples were checked for PCR inhibitors as described previously . Samples positive for PCR inhibitors were serially diluted (1:10 steps) until no inhibition remained. Dilutions of up to 1:100 and 1:10,000 were necessary to remove inhibition in pools of bodies and heads, respectively. Diluted samples were then used in a triplex real-time PCR that differentiates O. volvulus from Onchocerca ochengi, a morphologically similar bovine parasite also transmitted by S. damnosum s.l., based on differences in their respective ND5 genes (GenBank: AY462885.1 and FM206483.1) . Samples were not screened for Onchocerca ramachandrini, a parasite of warthogs (Phacochoerus africanus) also transmitted by S. damnosum s.l. , since the study area was beyond the limit of its host range .
Triplicate reactions were carried out in 20 μl total volumes containing 2 μl template DNA or water as a negative control, 1× HotStar Taq Buffer (Qiagen, N.V.), 4.5 mM MgCl2, 200 µM dNTPs (50 µM each, Fisher Scientific GmbH, Schwerte, Germany), 0.5 units HotStar Taq (Qiagen), and primers (Microsynth AG, Balgach, Switzerland) and hybridization probes (biomers.net GmbH, Ulm, Germany) at the concentrations listed in Table 2 . The reactions were performed in a Rotor Gene Q cycler (Qiagen, Hilden, Germany) with the following cycling conditions: Taq polymerase activation at 95 °C for 15 min and then 45 cycles at 95 °C for 10 s and 61 °C for 30 s, with fluorescence acquisition on the Fam, Hex, and Cy5 channels at the end of each cycle. Plasmids containing the respective sequences were used as PCR-positive controls in every run . Due to the need for such strong dilutions (1:10,000), the samples were run in triplicate up to four times. Samples were accepted as true positives for O. volvulus or O. ochengi infection if at least two of the three runs were positive, or if one run was O. volvulus ND5-positive alongside a positive 16S result (pan-Onchocerca probe). Bodies positive for either filarial species were considered to contain microfilariae or developing larvae, while heads positive for either species were considered to contain L3 larvae.
The Wilson method was used to calculate confidence intervals for all proportions (biting rates, parity rates, and L1–L2 and L3H infection rates), as recommended by Agresti and Coull . The coverage probabilities were close to the nominal confidence levels. All the confidence intervals used later were calculated at the 95% significance level. Biting rates were analysed using a generalized negative binomial regression to take into account overdispersion and to test possible differences across collection sites and seasons. A logistic regression model was used to test the effect of seasonality on parity and infection rates between the different collection sites. The above analyses were performed in R v4.1.0 . blackfly collection and dissection data were used to estimate monthly and annual S. damnosum s.l. biting rates and the monthly (MTP) and annual (ATP) transmission potentials at each site following methods described by Walsh et al. . Poolscreen v2.0  was used to estimate O. volvulus infection rates (maximum likelihood) with 95% confidence intervals (Bayes credibility estimate) at each site. Infection rates were only estimated in pools of heads since almost all pools of bodies were positive and the sensitivity of analysis is not known when pool positivity exceeds 25% .
Identification of S. damnosum complex
Morphological and cytotaxonomic identification of S. damnosum complex larvae revealed the presence of S. squamosum s.l. and S. mengense (Figs. 3, 4), both of which have previously been reported from the area . However, S. squamosum did not conform to descriptions of the cytotypes A, B, C, or D, known from Cameroon [4, 5, 7]. All specimens analysed from Nyamongo I (n = 39) and Bayomen (n = 10) contained homozygous inversions 1S-1 and 1L-3, which are generally fixed in the S. squamosum subcomplex (Fig. 4a) . They also possessed a new 1L inversion from sections p34 to p39 which was homozygous and fixed within the population, and which we designate 1L-57. In addition, 16/17 male specimens from Nyamongo I and 3/3 male specimens from Bayomen possessed a sex-linked band dimorphism (3C-Sp) and a heterozygous inversion (3L/82) near the centromere of chromosome 3 (Fig. 4b). This was absent in all 29 females collected from both sites. The involvement of 3C in sex determination among S. squamosum is diagnostic for S. squamosum E (= type III) [32, 33]. Since only a small number of specimens were examined from a narrow geographic range, we propose the name S. squamosum E2 for this variant possessing inversion 1L-57.
An additional five larval specimens collected at Bayomen were identified as S. mengense based on the 1CER expanded centromere of chromosome 1 (Fig. 4c) and tufts of hair-like scales on the anterior dorsum of the larval thorax (Fig. 3b) [10, 15, 44]. Adult S. mengense have hairs on the subcostal wing vein [10, 44], but no adult flies examined from collections made between July 2016 and January 2017 (Bayomen n = 254, Nyamongo I n = 74, and Ondouano n = 77) possessed this characteristic.
Adult blackfly collections and biting rates
A total of 93,573 adult female S. damnosum s.l. were recorded biting humans across the four sites during the sampling period (Additional file 1: Tables S1–S4, Additional file 2: Dataset). Of the specimens collected, 9281 were dissected for parity and evidence of Onchocerca species infection. The remainder were preserved in ethanol for pool screening. There was some discrepancy between the number of black flies recorded on field collection forms and the actual numbers dissected and preserved for pool screening. These manual count errors were mainly caused by high blackfly densities but do not alter patterns of biting and transmission. Biting rates were highest at the two riverside sites, but a generalized negative binomial regression showed that monthly catches were significantly higher at Bayomen than Nyamongo I (P < 0.001) (Additional file 1: Table S5). The mean daily biting rate at Bayomen remained consistently high throughout the year and although it decreased in July, December, and February, it remained above 700 black flies/person/day (Fig. 5). Annual biting rates (ABR) of 606,370 and 233,167 were estimated at Bayomen and Nyamongo I, respectively. There were two clear peaks in biting at Nyamongo I. The first coincided with the months of highest average rainfall in September and October, before rates decreased abruptly in November at the onset of the long dry season. A second increase occurred in December, and biting rates remained high until March, before decreasing in April at the onset of the new rainy season. Most of the biting at Nyamongo I, therefore, occurred during the long dry season. Similar seasonal patterns were observed at Egona II (ABR 89,849) and Ondouano (ABR 20,540), although blackfly activity decreased with increasing distance from the river (Fig. 5). Peaks in biting were recorded at both sites in October and March, but daily biting rates at Ondouano otherwise remained relatively low (< 85). Monthly biting rates at Egona II and Ondouano were both significantly lower than at Nyamongo I (P < 0.001) (Additional file 1: Table S5).
The Wilson method showed that the overall percentage of parous flies was higher at Bayomen (36.3%, 95% CI 34.7–37.8%) than at Nyamongo I (20.7%, 95% CI 19.3–22%) (P < 0.001) (Additional file 1: Tables S1–S2). Parity rates were lower at Egona II (11.1%, 95% CI 9.6–12.6%) (P < 0.001) and Ondouano (9.4%, 95% CI 0.7–12.3%) (P < 0.001) than at Nyamongo I, but there was no difference between Egona II and Ondouano (P = 0.29) (Additional file 1: Tables S3–S4). Parity rates were higher in the March–June rainy season than in the August–October rainy season based on data from both riverside sites combined (P < 0.001) (Fig. 6a). The odds ratio for being bitten by a parous fly was 1.8 times higher in the March–June rainy season than in the August–October rainy season. Parity rates were < 10% at Egona II and < 5% at Ondouano when biting rates peaked at these sites in October 2016 and March 2017.
Onchocerca species transmission intensity as assessed by dissection
Parasites that were morphologically similar to O. volvulus were mainly transmitted between January and May, overlapping the dry and rainy seasons during this period. The ATP was higher and peaks in transmission occurred slightly later at Bayomen (ATP 4488) than at Nyamongo I (2360) (Fig. 6b; Additional file 1: Table S6). The mean number of L3H per infective fly was higher at Bayomen (6.4) than Nyamongo I (2.4), and the highest estimated MTP (2406) at Bayomen was calculated based on just two infective flies carrying 17 L3H parasites. When data were combined for the two riverside sites, there was no effect of rainy season on the proportion of flies carrying only L1 and L2 stage parasites (P = 0.528), suggesting that they were present throughout the year (Fig. 6a). However, the proportion of flies carrying L3H parasites was higher in the March–June rainy season (P = 0.043), corresponding with higher parity rates and suggesting that parasites were developing to transmissible forms during this period (Fig. 6a). At Egona II, the ATP (102) was substantially lower than at the riverside sites (Fig. 6b; Additional file 1: Table S6). Infective L3H parasites were only found at Egona II in two flies collected in April 2017, at the beginning of the rainy season and shortly after biting rates had peaked (Additional file 1: Table S3). No parasites of any stage were found in the 489 flies dissected at Ondouano during the study (Additional file 1: Table S4). When considering all L3 larvae, regardless of their location in the blackfly, the ATP at Bayomen increased to 5479, but all other sites remained unchanged.
Onchocerca spp. transmission intensity as assessed by pool screening
Onchocerca volvulus was the predominant species developing to infective stages in S. damnosum s.l. (106/434 head pools), and while O. ochengi was present in the area (7/417 body pools), only 1/434 head pools were positive for infection. Pool screening revealed estimated infection rates in blackfly heads at Bayomen (0.10%), Nyamongo I (0.33%), and Egona II (0.13%), similar to the percentage of infective (L3H) flies encountered by dissection during the study at these sites (0.13%, 0.36%, and 0.13%, respectively) (Fig. 7). However, 4/15 pools containing 1453 black flies were positive at Ondouano, 7.9 km from the riverside, despite parasites being absent from the 489 flies dissected at this site.
Our study confirmed ongoing transmission of O. volvulus by black flies along the lower Mbam river in 2016/2017. Breeding site surveys were limited to a single period in the dry season and revealed the presence of S. squamosum s.l. and S. mengense. The involvement of chromosomal inversion 3C in sex determination and the presence of a new inversion in 1L indicates that the former is a variant of S. squamosum E, and not S. squamosum A, which was previously reported in the Mbam river basin, upstream from our study sites . The E cytotype is currently known from western Côte d’Ivoire, Ghana, Guinea, Liberia, and Sierra Leone [5, 33], while previous reports of its presence in Benin, Central African Republic, and Togo were made in error [8, 33]. Our collection of S. squamosum E2 from the Mbam river is, therefore, the first record of this cytotype (or a variant of it) east of Lake Volta, more than 1200 km away.
While there is little reason to suggest that S. squamosum A was absent from the area, it is possible that the composition of cytoforms breeding along the lower Mbam has changed over time. Reductions in rainfall and river discharge have occurred since the early 1970s, and dams were built upstream from Bafia on the Noun and Mbam rivers in 1974 and 1987, respectively . Environmental changes such as these alter the physical properties of river water, including rates of discharge, which could affect the blackfly species present [24, 45]. However, these events pre-date collections by Barbazan et al.  in 1993/1994. Another possible explanation for the different cytotypes encountered between the studies is that their composition changes seasonally. Traoré-Lamizana and Lemasson  have already shown that S. squamosum s.l. spreads north from the forest–savannah transition zones in Cameroon and into savannah rivers during the rainy season, but this was before the different S. squamosum cytotypes had been described. Spatiotemporal studies of the S. squamosum cytotypes present in Cameroon are needed to document the extent and dynamics of their distributions. Additionally, the breeding range and vector competence of S. squamosum E2 needs to be investigated in view of the severe clinical presentation of onchocerciasis in the Mbam valley.
We encountered high densities of biting black flies throughout the year at riverside sites in agreement with findings recently reported by Abong et al.  from nearby Biatsota village. Biting rates were higher at Bayomen than at Nyamongo I, possibly due to the proximity of collection sites to the nearest known breeding sites (< 300 m and ≈1 km, respectively). At Nyamongo I, seasonal biting patterns were similar to those recorded in 1993/1994 although ABRs were considerably higher in the current study (ABRs of 233,167 and 98,208, respectively) . This was partly due to the extended biting peak that occurred during the long dry season compared with the same period in 1993/1994 . The dry season peak at Nyamongo I occurred at a time (January–March) when the average river discharge was relatively low and stable (< 500 m3/s) (Fig. 2) . These conditions minimise the impact of fluctuating water levels on blackfly breeding sites and may expose previously submerged oviposition substrates, leading to population increases [24, 47, 48]. Similar biting patterns were observed at Egona II and Ondouano, although biting peaked at these sites in October and March. The estimated ABR of 89,849 at Egona II was almost as high as the 98,208 estimated by Barbazan et al.  at the riverside in 1993/1994. However, the ABR declined markedly at Ondouano (20,540), 7.9 km from the riverside, whereas it was still 43,790 at a similar distance (7.2 km) in the previous study. In addition, Barbazan et al. had two more collection sites along the same transect at 13.5 km and 23 km, both of which had ABRs > 43,000 . Considering that they reported no breeding along tributaries, it may be that black flies dispersed greater distances at the time of the earlier study. Recent population growth in villages close to the river has also increased the availability of potential blood hosts which might decrease the need for dispersal . Alternatively, decreased dispersal may indicate the presence of a different S. damnosum cytotype with different behavioural characteristics.
Differences in the proximity of our riverside collection sites to the nearest breeding sites might also explain why parity rates were higher at Bayomen than at Nyamongo I. Duke  showed that the number of parous flies can decrease rapidly with increasing distance to breeding sites, and our results would appear to agree. It is, therefore, worth noting that the Bayomen collection site was approximately 1.5 km from Bayomen village, where blackfly parity rates might have been lower. Differences in the presence and relative abundance of non-human blood hosts might also affect parity rates between populations. It is well known that black flies feed on a variety of vertebrates including cattle [51, 52], and it has been shown that O. ochengi appears in black flies at times coinciding with cattle migrations in northern Cameroon . Our blackfly pools that tested positive for O. ochengi were collected at Bayomen between October 2016 and April 2017, coinciding with the transhumance period in Cameroon . The number of migratory cattle arriving in the area is unclear (M. Ronse pers. comm.), however, a 2015 report by the Programme National de Développement Participatif  revealed there to be fewer than 100 cows kept in the Kon-Yambetta canton (population 18,268) which includes Bayomen, and which may explain why only small numbers of O. ochengi were detected.
We documented a gradual increase in blackfly parity rates throughout the sampling period despite our dissection data being relatively limited (Fig. 6a). Combined parity rates at the two riverside sites were higher in the March–June rainy season than in the August–October rainy season, which agrees with previous results . The parity rates recorded at our riverside collection sites were similar to those recently reported from Biatsota (18.9–41.4%)  but were lower than the 78% average reported from the Mbam/Sanaga basin in 1993/1994 . Barbazan et al. also reported that there was little difference in parity between black flies on the shoreline and at sites further inland . Contrastingly, we found significantly lower parity rates at Egona II and Ondouano than at the riverside. At our inland sites, parity rates decreased to < 10% during the months of peak biting indicating relatively little exposure to parous flies among these communities.
The combined results of blackfly catches, dissections, and pool screening showed that O. volvulus transmission occurred mostly at the riverside, although the estimated percentage of infective flies was lower than recently documented at Biatsota . Based on dissections, transmission peaked slightly earlier at Nyamongo I than at Bayomen (Fig. 6b), but generally occurred between January and May, overlapping the dry and rainy seasons during this period. This is in broad agreement with Barbazan et al.  who documented peaks in O. volvulus transmission in February and March along their Ngoro transect, extending from Nyamongo I. The ATP of 3113 reported at the Ngoro riverside site was similar to our finding at Nyamongo I (ATP 2360). However, the higher biting rates that we used to calculate this ATP indicate that O. volvulus infection rates in black flies have decreased since the earlier study.
We found that peaks in O. volvulus infection coincided with peak biting rates at Nyamongo I, whereas at Bayomen, biting rates were continuously high throughout the year. Overall, 90% of all L3 larvae were found in the heads of dissected S. damnosum s.l., which is higher than in many other studies . The high densities of biting black flies and a higher mean number of L3H per infective fly resulted in a higher estimated ATP at Bayomen than at Nyamongo I. Nevertheless, the collection site at Nyamongo I was more likely to be a point of significant human–vector contact and parasite transmission. It was adjacent to a busy but slow-moving ferry crossing where people often gathered for lengthy periods of time, while the site at Bayomen was 1.5 km from the village and accessed mainly by farmers and fishermen. Building a bridge at the ferry crossing might alleviate some of the congestion and reduce human–vector contact. Parasite transmission declined rapidly with increasing distance from the riverside, indicating that front line communities and others working on or near the river are at highest risk of exposure to infective bites. However, pool screening detected the presence of transmissible O. volvulus parasites at each of our collection sites, up to 7.9 km from the riverside. Importantly, the upper bound of the confidence interval calculated at each site was above 0.05% (equivalent to 1 in 2000 flies infected) and consequently above the WHO threshold for interruption of O. volvulus transmission .
At the time of our study, CDTI was taking place around July each year [2, 3]. Ivermectin is most effective at reducing the uptake of microfilariae by black flies during the first 6 months after treatment . If O. volvulus transmission along the Mbam river was constant throughout the year, this might explain the lower transmission rates encountered between July and December based on the results of dissections. However, similar patterns of transmission were reported in the area in 1993/1994, and at additional sites along the Sanaga river  where S. squamosum B breeds and bites perennially . These entomological studies were conducted at a time pre-dating ivermectin mass treatment in the Mbam valley [2, 58]. Even the earliest mass treatments (1994–1997) only covered approximately 10% of the population, and high (> 65%) therapeutic coverage was not achieved until several years after CDTI commenced [2, 3]. Accordingly, ivermectin would not have affected the results of the earlier study  and these patterns of transmission probably reflect true seasonal cycles. Similar seasonality has been reported for perennially breeding S. squamosum in Togo , where peaks in O. volvulus transmission coincided with high dry-season biting rates and high parity rates despite black flies being present throughout the year. Favourable weather conditions are also thought to increase the probability of blackfly survival and human–vector contact resulting in dry season peaks in transmission involving the forest cytoform, S. yahense [59, 60]. Consideration should, therefore, be given to seasonal patterns of O. volvulus transmission along the Mbam river when planning CDTI activities. Whereas communities normally decide on the period of ivermectin distribution, optimizing the timing of CDTI to coincide with the period prior to peak transmission may enhance the programme outcome at no additional cost . Alternatively, annual moxidectin treatment might overcome the effects of seasonal transmission since it is effective at maintaining low skin microfilarial levels for more than 12 months after administration [61, 62]. Complementary treatment strategies such as seasonal larviciding during times of peak O. volvulus transmission  or slash and clear-based vector control methods  should also be considered in view of the high blackfly densities encountered along the Mbam.
Our study was mainly limited by problems encountered during pool screening which prevented longitudinal analyses of infection rates or estimates of transmission potentials based on our larger sample of black flies. High dilutions were needed (> 1:1000) to remove the effect of PCR inhibitors from DNA samples of specimens collected between February and June 2017, which suggests a problem with the DNA extraction. The dilution of samples to remove inhibitors also dilutes the target DNA and may have pushed the template below our detection limit. Smaller pool sizes might have reduced this source of error and will be considered in future studies. As a result, many of our inferences relied upon data obtained from a comparatively small number of dissected black flies, and our study does not take into consideration year-to-year variation in biting rates and parasite transmission. Nevertheless, our findings agree with previous investigations [11, 46] and show that O. volvulus transmission is sustained by high blackfly biting rates along the lower Mbam river where onchocerciasis continues to be a public health concern [3, 21, 58].
After 16 years of annual CDTI in the lower Mbam valley, O. volvulus infection rates in black flies were above the 0.05% threshold for interruption of transmission. Although the percentage of infective black flies was relatively low, extremely high riverside biting rates were likely a major contributing factor to ongoing transmission in an area where a new variant of S. squamosum E is described.
Availability of data and materials
All data generated or analysed during this study are included in this published article [and its Additional files].
Macé JM, Boussinesq M, Ngoumou P, Enyegue Oye J, Koeranga A, Godin C. Country-wide rapid epidemiological mapping of onchocerciasis (REMO) in Cameroon. Ann Trop Med Parasitol. 1997;91(4):379–91.
Pion SD, Clement MC, Boussinesq M. Impact of four years of large-scale ivermectin treatment with low therapeutic coverage on the transmission of Onchocerca volvulus in the Mbam valley focus, central Cameroon. Trans R Soc Trop Med Hyg. 2004;98(9):520–8. https://doi.org/10.1016/j.trstmh.2003.11.010.
Kamga GR, Dissak-Delon FN, Nana-Djeunga HC, Biholong BD, Mbigha-Ghogomu S, Souopgui J, et al. Still mesoendemic onchocerciasis in two Cameroonian community-directed treatment with ivermectin projects despite more than 15 years of mass treatment. Parasites Vectors. 2016;9(1):581. https://doi.org/10.1186/s13071-016-1868-8.
Traoré-Lamizana M, Somiari S, Mafuyai HB, Vajime CG, Post RJ. Sex chromosome variation and cytotaxonomy of the onchocerciasis vector Simulium squamosum in Cameroon and Nigeria. Med Vet Entomol. 2001;15(2):219–23.
Adler PH, Crosskey RW. World blackflies (Diptera: Simuliidae): a comprehensive revision of the taxonomic and geographical inventory . 2014. http://biomia.sites.clemson.edu/pdfs/blackflyinventory.pdf.
Traoré-Lamizana M, Lemasson JJ. Participation à une étude de faisabilité d'une campagne de lutte contre l’onchocercose dans la région du bassin du Logone. Répartition des espèces du complexe Simulium damnosum dans la zone camerounaise du projet. Cah ORSTOM Sér Ent Méd Parasit. 1987;25:171–86.
Mustapha M, Post RJ, Enyong P, Lines J. A new cytotype of Simulium squamosum from south-west Cameroon. Med Vet Entomol. 2004;18(3):296–300. https://doi.org/10.1111/j.0269-283X.2004.00501.x.
Post RJ, Cheke RA, Boakye DA, Wilson MD, Osei-Atweneboana MY, Tetteh-Kumah A, et al. Stability and change in the distribution of cytospecies of the Simulium damnosum complex (Diptera: Simuliidae) in southern Ghana from 1971 to 2011. Parasites Vectors. 2013;6:205. https://doi.org/10.1186/1756-3305-6-205.
Vajime CG, Dunbar RW. Chromosomal identification of eight species of the subgenus Edwardsellum near and including Simulium (Edwardsellum) damnosum Theobald (Diptera: Simuliidae). Tropenmed Parasitol. 1975;26(1):111–38.
Mustapha M, Post RJ, Krüger A. The cytotaxonomy and morphotaxonomy of Simulium mengense (Diptera: Simuliidae). Ann Trop Med Parasitol. 2004;98(5):509–23. https://doi.org/10.1179/000349803225003523.
Barbazan P, Escaffre H, Mbentengam R, Boussinesq M. Entomologic study on the transmission of onchocerciasis in a forest-savanna transition area of Cameroon. Bull Soc Pathol Exot. 1998;91(2):178–82.
Vajime CG, Dunbar RW. The chromosomal identification of Simulium (Edwardsellum) mengense new species (Diptera: Simuliidae). Parassitologia. 1977;19(1–2):95–102.
Renz A, Barthelmess C, Eisenbeiß W. Vectorial capacity of Simulium damnosum s.l. populations in Cameroon. Trop Med Parasitol. 1987;38:344–5.
Demanou M, Enyong P, Pion SD, Basanez MG, Boussinesq M. Experimental studies on the transmission of Onchocerca volvulus by its vector in the Sanaga valley (Cameroon): Simulium squamosum B. Intake of microfilariae and their migration to the haemocoel of the vector. Ann Trop Med Parasitol. 2003;97(4):381–402. https://doi.org/10.1179/000349803235002254.
Post RJ, Onyenwe E, Somiari SAE, Mafuyai HB, Crainey JL, Ubachukwu PO. A guide to the Simulium damnosum complex (Diptera: Simuliidae) in Nigeria, with a cytotaxonomic key for the identification of the sibling species. Ann Trop Med Parasitol. 2011;105(4):277–97. https://doi.org/10.1179/136485911X12987676649700.
Boussinesq M, Pion SD, Demanga-Ngangue KJ. Relationship between onchocerciasis and epilepsy: a matched case–control study in the Mbam Valley, Republic of Cameroon. Trans R Soc Trop Med Hyg. 2002;96(5):537–41.
Chesnais CB, Nana-Djeunga HC, Njamnshi AK, Lenou-Nanga CG, Boullé C, Bissek AZ, et al. The temporal relationship between onchocerciasis and epilepsy: a population-based cohort study. Lancet Infect Dis. 2018;18(11):1278–86. https://doi.org/10.1016/s1473-3099(18)30425-0.
Chesnais CB, Bizet C, Campillo JT, Njamnshi WY, Bopda J, Nwane P, et al. A second population-based cohort study in Cameroon confirms the temporal relationship between onchocerciasis and epilepsy. Open Forum Infect Dis. 2020;7(6):206. https://doi.org/10.1093/ofid/ofaa206.
Nana-Djeunga HC, Bourguinat C, Pion SD, Bopda J, Kengne-Ouafo JA, Njiokou F, et al. Reproductive status of Onchocerca volvulus after ivermectin treatment in an ivermectin-naive and a frequently treated population from Cameroon. PLoS Negl Trop Dis. 2014;8(4):e2824. https://doi.org/10.1371/journal.pntd.0002824.
Gardon J, Boussinesq M, Kamgno J, Gardon-Wendel N, Demanga N, Duke BO. Effects of standard and high doses of ivermectin on adult worms of Onchocerca volvulus: a randomised controlled trial. Lancet. 2002;360(9328):203–10. https://doi.org/10.1016/s0140-6736(02)09456-4.
Tekle AH, Zouré HGM, Noma M, Boussinesq M, Coffeng LE, Stolk WA, et al. Progress towards onchocerciasis elimination in the participating countries of the African Programme for Onchocerciasis Control: epidemiological evaluation results. Infect Dis Poverty. 2016;5(1):66. https://doi.org/10.1186/s40249-016-0160-7.
World Health Organization: African Programme for Onchocerciasis Control (APOC): Cameroon. http://www.who.int/apoc/countries/cmr/en/ (2015). Accessed 25 Sep 2015.
Bird MI. A seasonal cycle in the carbon isotope composition of organic carbon in the Sanaga River, Cameroon. Limnol Oceanogr. 1998;43(1):143–6.
Ndam Ngoupayou JR, Dzana JG, Kpoumie A, Ghogomu RT, Fouepe Takounjou A, Braun JJ, et al. Present-day sediment dynamics of the Sanaga catchment (Cameroon): from the total suspended sediment (TSS) to erosion balance. Hydrolog Sci J. 2016;61(6):1080–93. https://doi.org/10.1080/02626667.2014.968572.
Programme National de Développement Participatif. Plan Communal de Développement de la commune de Kon Yambetta - Région du Centre. PNDP. 2015. https://www.pndp.org/documents/15_PCD_KON_YAMBETTA1.pdf. Accessed 10 Oct 2021.
KNMI Climate Explorer: Time series monthly Bafia GHCN v2 precipitation. https://climexp.knmi.nl (2017). Accessed 21 Aug 2017.
National Center for Atmospheric Research: Augmented Monthly Flow Rates of World Rivers (except former Soviet Union). https://rda.ucar.edu/datasets/ds552.1/ (2017). Accessed 21 Aug 2017.
Programme National de Développement Participatif. Plan Communal de Développement de Ngoro. PNDP. 2013. https://www.pndp.org/documents/13_PCD_NGORO1.pdf. Accessed 10 Oct 2021.
Walsh JF, Davies JB, Le Berre R, Garms R. Standardization of criteria for assessing the effect of Simulium control in onchocerciasis control programmes. Trans R Soc Trop Med Hyg. 1978;72(6):675–6.
Crosskey RW. The natural history of blackflies. UK: Wiley; 1990.
Adler PH, Currie DC, Wood DM. The Black Flies (Simuliidae) of North America. New York: Cornell University Press; 2004.
Boakye DA. A pictorial guide to the chromosomal identification of members of the Simulium damnosum Theobald complex in West Africa with particular reference to the Onchocerciasis Control Programme Area. Trop Med Parasitol. 1993;44(3):223–44.
Post RJ, Mustapha M, Krüger A. Taxonomy and inventory of the cytospecies and cytotypes of the Simulium damnosum complex (Diptera: Simuliidae) in relation to onchocerciasis. Trop Med Int Health. 2007;12(11):1342–53. https://doi.org/10.1111/j.1365-3156.2007.01921.x.
Freeman P, de Meillon B. Simuliidae of the Ethiopian Region. London: British Museum (Natural History); 1953.
World Health Organization. Training module for entomology technicians in charge of entomological activities in onchocerciasis control. Ouagadougou: Onchocerciasis Control Programme in West Africa; 2002.
Davies JB, Crosskey RW: Simulium-vectors of onchocerciasis. In: Vector Control Series. Geneva: World Health Organization; 1991:115.
Hendy A, Krüger A, Pfarr K, De Witte J, Kibweja A, Mwingira U, et al. The blackfly vectors and transmission of Onchocerca volvulus in Mahenge, south eastern Tanzania. Acta Trop. 2018;181:50–9. https://doi.org/10.1016/j.actatropica.2018.01.009.
Colebunders R, Mandro M, Mokili JL, Mucinya G, Mambandu G, Pfarr K, et al. Risk factors for epilepsy in Bas-Uele Province, Democratic Republic of the Congo: a case–control study. Int J Infect Dis. 2016;49:1–8. https://doi.org/10.1016/j.ijid.2016.05.018.
Verocai GG, Hassan HK, Lakwo T, Habomugisha P, Katabarwa MN, Begumisa S, et al. Molecular identification of Onchocerca spp. larvae in Simulium damnosum sensu lato collected in northern Uganda. Am J Trop Med Hyg. 2017;97(6):1843–5. https://doi.org/10.4269/ajtmh.17-0525.
Jori F, Bastos AD. Role of wild suids in the epidemiology of African swine fever. EcoHealth. 2009;6(2):296–310. https://doi.org/10.1007/s10393-009-0248-7.
Agresti A, Coull BA. Approximate is better than “exact” for interval estimation of binomial proportions. Am Stat. 1998;52(2):119–26. https://doi.org/10.1080/00031305.1998.10480550.
R Core Team: R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2021.
Barker JT. Statistical estimators of infection potential based on PCR pool screening with unequal pool sizes [Ph.D. dissertation]. Birmingham, AL: University of Alabama at Birmingham; 2000.
Elsen P. Définition d’un groupe particulier d’espèces d’Afrique Centrale au sein du complexe Simulium damnosum et création de trois groupes dans le sous-genre Edwardsellum (Diptera, Simuliidae). Rev Zool Afric. 1983;97:633–9.
Grunewald J. The hydro-chemical and physical conditions of the environment of the immature stages of some species of the Simulium (Edwardsellum) damnosum complex (Diptera). Tropenmed Parasitol. 1976;27(4):438–54.
Abong RA, Amambo GN, Hamid AA, Enow BA, Beng AA, Nietcho FN, et al. The Mbam drainage system and onchocerciasis transmission post ivermectin mass drug administration (MDA) campaign, Cameroon. PLoS Negl Trop Dis. 2021;15(1):e0008926-e. https://doi.org/10.1371/journal.pntd.0008926.
Poff NL, Ward JV. Drift responses of benthic invertebrates to experimental streamflow variation in a hydrologically stable stream. Can J Fish Aquat Sci. 1991;48(10):1926–36. https://doi.org/10.1139/f91-229.
Car M. The influence of water-level fluctuation on the drift of Simulium chutteri Lewis, 1965 (Diptera, Nematocera) in the Orange River, South Africa. Onderstepoort J Vet Res. 1983;50(3):173–7.
Boullé C, Njamnshi AK, Dema F, Mengnjo MK, Siewe Fodjo JN, Bissek A-CZ, et al. Impact of 19 years of mass drug administration with ivermectin on epilepsy burden in a hyperendemic onchocerciasis area in Cameroon. Parasites Vectors. 2019;12(1):114. https://doi.org/10.1186/s13071-019-3345-7.
Duke BO. The differential dispersal of nulliparous and parous Simulium damnosum. Tropenmed Parasitol. 1975;26(1):88–97.
Lamberton PHL, Cheke RA, Walker M, Winskill P, Crainey JL, Boakye DA, et al. Onchocerciasis transmission in Ghana: the human blood index of sibling species of the Simulium damnosum complex. Parasites Vectors. 2016;9(1):432. https://doi.org/10.1186/s13071-016-1703-2.
Adler PH, Cheke RA, Post RJ. Evolution, epidemiology, and population genetics of black flies (Diptera: Simuliidae). Infect, Genet Evol. 2010;10(7):846–65. https://doi.org/10.1016/j.meegid.2010.07.003.
Wahl G, Enyong P, Ngosso A, Schibel JM, Moyou R, Tubbesing H, et al. Onchocerca ochengi: epidemiological evidence of cross-protection against Onchocerca volvulus in man. Parasitology. 1998;116(4):349–62.
Motta P, Porphyre T, Hamman SM, Morgan KL, Ngwa VN, Tanya VN, et al. Cattle transhumance and agropastoral nomadic herding practices in Central Cameroon. BMC Vet Res. 2018;14(1):214. https://doi.org/10.1186/s12917-018-1515-z.
Cheke RA, Garms R. Indices of onchocerciasis transmission by different members of the Simulium damnosum complex conflict with the paradigm of forest and savanna parasite strains. Acta Trop. 2013;125(1):43–52. https://doi.org/10.1016/j.actatropica.2012.09.002.
World Health Organization. Guidelines for stopping mass drug administration and verifying elimination of human onchocerciasis: criteria and procedures. Geneva: World Health Organization; 2016.
Cupp EW, Bernardo MJ, Kiszewski AE, Collins RC, Taylor HR, Aziz MA, et al. The effects of ivermectin on transmission of Onchocerca volvulus. Science. 1986;231(4739):740–2.
Dissak-Delon FN, Kamga GR, Humblet PC, Robert A, Souopgui J, Kamgno J, et al. Adherence to ivermectin is more associated with perceptions of community directed treatment with ivermectin organization than with onchocerciasis beliefs. PLoS Negl Trop Dis. 2017;11(8):e0005849. https://doi.org/10.1371/journal.pntd.0005849.
Cheke RA, Sowah SA, Avissey HS, Fiasorgbor GK, Garms R. Seasonal variation in onchocerciasis transmission by Simulium squamosum at perennial breeding sites in Togo. Trans R Soc Trop Med Hyg. 1992;86(1):67–71.
Barbiero VK, Trpis M. Transmission of onchocerciasis by local black flies on the Firestone Rubber Plantation, Harbel, Liberia. Am J Trop Med Hyg. 1984;33(4):586–94.
Boussinesq M, Fobi G, Kuesel AC. Alternative treatment strategies to accelerate the elimination of onchocerciasis. Int Health. 2018;10(Suppl_1):i40–8. https://doi.org/10.1093/inthealth/ihx054.
Turner HC, Walker M, Attah SK, Opoku NO, Awadzi K, Kuesel AC, et al. The potential impact of moxidectin on onchocerciasis elimination in Africa: an economic evaluation based on the Phase II clinical trial data. Parasites Vectors. 2015;8(1):167. https://doi.org/10.1186/s13071-015-0779-4.
Siewe Fodjo JN, Vieiri MK, Ngarka L, Njamnshi WY, Nfor LN, Mengnjo M, et al. 'Slash and clear’ vector control for onchocerciasis elimination and epilepsy prevention: a protocol of a cluster randomised trial in Cameroonian villages. BMJ Open. 2021:0:e050341. https://doi.org/10.1136/bmjopen-2021-050341.
We thank Peter Adler and Andreas Krüger for advice regarding the cytotaxonomy; Emilia Agbor and Akem Mbi for blackfly dissections; Oben Bruno for technical assistance in the field; Julia Irani, Jospeh Nelson Siewe Fodjo, and Maya Ronse for input when preparing and conducting the work and while preparing the manuscript. We also thank the residents of the Mbam valley who supported the work throughout.
Work was funded by the Department of Economy, Science, and Innovation (EWI) of the Flemish government through a Structural Research Funding (SOFI) grant awarded via the Institute of Tropical Medicine, Antwerp; the European Research Council (Grant number 671055); and the German Center for Infection Research (Grant number TI 03.907).
Ethics approval and consent to participate
blackfly collections involving human participants were subject to review and approval by the Institutional Review Board at the Institute of Tropical Medicine, Antwerp, Belgium (1041/15) and the Comité National d’Ethique de la Recherche pour la Santé Humaine (CNERSH), Cameroon (2016/03/677/L/CNERSH/SP). Collectors received training and were not considered to be at a higher risk of exposure than others living in local communities. Participants were adults over the age of 18 years who were receiving ivermectin treatment in accordance with the national onchocerciasis control programme and provided written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Summary of blackfly collection and dissection data from catches made at Bayomen between July 2016 and June 2017. Table S2. Summary of blackfly collection and dissection data from catches made at Nyamongo I between July 2016 and June 2017. Table S3. Summary of blackfly collection and dissection data from catches made at Egona II between July 2016 and June 2017. Table S4. Summary of blackfly collection and dissection data from catches made at Ondouano between July 2016 and June 2017. Table S5. Estimated monthly (MBR) and annual (ABR) biting rates at the four collection sites. Table S6. Estimated monthly (MTP) and annual (ATP) transmission potentials at the four collection sites.
Dataset. File contains pool screen, dissection, and biting rate data.
About this article
Cite this article
Hendy, A., Krit, M., Pfarr, K. et al. Onchocerca volvulus transmission in the Mbam valley of Cameroon following 16 years of annual community-directed treatment with ivermectin, and the description of a new cytotype of Simulium squamosum. Parasites Vectors 14, 563 (2021). https://doi.org/10.1186/s13071-021-05072-y
- Onchocerca volvulus
- Simulium damnosum
- Simulium squamosum