Evidence of Transmission of Plasmodium vivax 210 and Plasmodium vivax 247 by Anopheles gambiae and An. coluzzii, Major Malaria Vectors in Benin/West Africa

Simple Summary For any country to achieve an effective malaria control program, the information on causative parasite diversity and dynamics in the region is a key factor. This study demonstrates that in Benin, the occurrence of malaria cases is not solely caused by Plasmodium falciparum (common malaria parasite) and involves other non-falciparum species. The aim of this work is to assess the prevalence of various Plasmodium species in locally collected mosquito cohorts, either through human landing catches or pyrethrum spray catches. Thus, the comparative data on the prevalence of circumsporozoite protein (CSP) antibodies of P. falciparum and P. vivax in Anopheles gambiae s.l., the major insect-vector responsible for transmitting the malaria parasite to humans in Benin, were reported. Further, to delineate the possible contribution of various Plasmodium species in malaria infection in collected mosquitoes, a molecular species-specific polymerase chain reaction analysis was performed. The important finding of this study reveals the variation in the distribution of Plasmodium species prevalence in Benin. The information incurred through this study will be very helpful for malaria control stakeholders as well as planning and management agencies designing effective diagnostics and preventive measures and treatments to curb the lethal impact of malaria in Benin. Abstract Current diagnostic and surveillance systems in Benin are not designed to accurately identify or report non-Plasmodium falciparum (Pf) human malaria infections. This study aims to assess and compare the prevalence of circumsporozoite protein (CSP) antibodies of Pf and P. vivax (Pv) in Anopheles gambiae s.l. in Benin. For that, mosquito collections were performed through human landing catches (HLC) and pyrethrum spray catches (PSC). The collected mosquitoes were morphologically identified, and Pf, Pv 210, and Pv 247 CSP antibodies were sought in An. gambiae s.l. through the ELISA and polymerase chain reaction (PCR) techniques. Of the 32,773 collected mosquitoes, 20.9% were An. gambiae s.l., 3.9% An. funestus gr., and 0.6% An. nili gr. In An. gambiae s.l., the sporozoite rate was 2.6% (95% CI: 2.1–3.1) for Pf, against 0.30% (95% CI: 0.1–0.5) and 0.2% (95% CI: 0.1–0.4), respectively, for Pv 210 and Pv 247. P. falciparum sporozoite positive mosquitoes were mostly An. gambiae (64.35%), followed by An. coluzzii (34.78%) and An. arabiensis (0.86%). At the opposite, for the Pv 210 sporozoite-positive mosquitoes, An. coluzzii and An. gambiae accounted for 76.92% and 23.08%, respectively. Overall, the present study shows that P. falciparum is not the only Plasmodium species involved in malaria cases in Benin.

In communes located in the southern and central part of the country, the hot and humid climate is of the sub-equatorial type with four seasons, including a long (April to July) and short (October to November) rainy season, as well as a long (December to March) and short (August to September) dry season. In the northern communes, the climate is Sudano-Guinean, with a single rainy season (June to November). The rest of the year (December to May) is the dry season. The average annual rainfall is about 1340 mm in the south and 700 mm in the north [20], while the temperature varies little with an annual average between 26 • C and 30 • C [21].

Mosquito Collection
In all sites, the collection of adult mosquitoes was performed using both Human Landing Catches (HLCs) and morning Pyrethrum Spray Catches (PSCs). The two methods were used to increase the likelihood to collect a great number of Plasmodium sporozoitepositive mosquitoes.
HLCs were performed between 9:00 p.m. and 5:00 a.m. using eight people, with four (two positioned indoors and two outdoors) used between 9:00 p.m. and 1:00 a.m. and replaced with four others from 1:00 a.m. to 5:00 a.m. Thus, each seated collector used a flashlight as well as a sucking tube to collect the mosquitoes that landed on his/her naked lower legs. Collectors were rotated over collection of hours to limit bias related to individual variations in the attractiveness of Anopheles vectors. They were previously trained to collect the mosquitoes.
PSCs were carried out in 10 bedrooms randomly selected in each commune. The houses that were selected have ceilings and do not have eaves. They also have windows and doors that are kept closed before spraying, which prevents mosquitoes from escaping. This method, which was started at 7:00 a.m., enabled the collection of mosquitoes that entered the bedrooms overnight. Thus, the Rambo aerosol (Gongoni Company Limited, Kano, Nigeria) combining Beta-cypermethrin, tetramethrin, and permethrin was sprayed in the interior of rooms where white sheets were previously laid on the floor. This allowed the collection of indoor-resting mosquitoes.
Adult mosquitoes were then morphologically identified using a binocular loupe, as well as the morphological identification key of Gillies & Coetzee [22]. Anopheles were stored individually in eppendorf tubes containing silica gel and cotton and kept at −20 • C until processing.

Blood Meal Origin in An. gambiae s.l. Collected by PSC
The origin of blood meal in females An. gambiae s.l. collected through PSC in the different communes were determined using the ELISA test, the protocol of which was described by Beier et al. [23]. Antibodies used for this testing were from human (Sigma, St. Louis, MO, USA, 11 mg/mL), cow (Sigma, 11 mg/mL), sheep (Sigma, 4 mg/mL), and pig (Sigma, 11 mg/mL). The blood meal origin was in Anopheles to ensure that the detected Plasmodium species were from humans living in the commune. Abdomens of blood-fed mosquitoes were ground in 250 µL of Phosphate Buffer Saline (PBS). For each well, 50 µL of each blood eluate were introduced (one mosquito per column) and incubated for 3 h at room temperature. The sensitized plates were washed twice with PBS-Tween 20. For each well and each row, 50 µL of labeled IgG solution were distributed. After a 1 h incubation at room temperature, the plates were emptied and washed with PBS/Tween 20. Then, 100 µL of the peroxidase enzyme substrate were added to each well and incubated in the dark for 30 min. A blue color was obtained, and the reaction was blocked with 50 µL of a 4 N sulfuric acid solution. Results were read at 450-620 nm with the BioTek Elx808 ELISA plate reader (BioTek, Winsooki, VT, USA ). The negative (PBS 1×) and positive controls, in particular the sera of the different hosts sought (human, cow, goat, pig), were placed on each line.

Circumsporozoite
Protein ELISA for the Detection of P. falciparum, P. vivax 210 and P. vivax 247 In this study, the detection of P. falciparum and P. vivax 210 and 247 was done using the heads and thorax of these sampled mosquitoes according to the method described by Burkot et al. [24]. A sampling of the collected Anopheles was done by commune in the laboratory before the tests were performed. The heads-thoraxes of the sampled mosquitoes were individually crushed, after a 1 h incubation in 50 µL of blocking buffer (BB) containing IGEPAL CA-630 (Sigma, Darmstadt, Germany). During the grinding, two times 100 µL of BB was supplemented to each ground material. Then, the 96 wells of a plate were first sensitized with circumsporozoite (CS) protein capture monoclonal antibodies (0.5 mg/mL, BEI Resources). For each well of the plate, 50 µL of this solution were introduced and left overnight on the bench. The following day, the wells were emptied without being washed and received 200 µL of BB for 1 h at room temperature. These wells were again emptied without being washed, and each received 50 µL of mosquito ground material.
After a 2 h incubation at room temperature, the wells were emptied and washed twice with PBS/Tween 20. Each well then received 50 µL of peroxidase-conjugated monoclonal antibodies (0.5 mg/mL, BEI Resources) for 1 h on the bench. Then, 7.5 µL of the substrate solution was mixed with 5 mL of BB. After 1 h of incubation at room temperature, the wells were emptied and washed 4 times with PBS/Tween 20. Then, 100 µL of the peroxidase substrate were placed in each well. The reaction occurs in the dark for 30 min. At the end, the positive wells were colored in blue. This reaction was blocked with 4 N sulfuric acid. The positive wells were colored in yellow. Plates were read using the BioTek Elx808 ELISA plate reader at 405 nm and 650 nm (BioTek, Winsooki, VT, USA). Laboratory-susceptible strain served as a negative control. P. falciparum, P. vivax 210, and P. vivax 247 at 2 pg/µL (BEI Resources) were used as positive controls. The testing was simultaneously conducted on 3 plates to detect the presence of the three Plasmodium species (P. falciparum, P. vivax 210, and P. vivax 247).

DNA Extraction of Plasmodium Species and Detection of P. vivax Serotypes
To characterize the Plasmodium species (P. falciparum, P. vivax 210, and 247) in the different communes, Plasmodium DNA was extracted from the grindings of the heads-thoraxes of the mosquitoes used for the ELISA CSP testing. In each commune, a total of 50 mosquitoes, including all Plasmodium sporozoite-positive samples identified through ELISA CSP testing, plus a subset of randomly selected negative ones, were analyzed by PCR according to the protocol described in Methods in Anopheles Research [25].
Plasmodium DNA was extracted with 2% CTAB (Cetyl trimethyl ammonium bromide). The grindings of the heads-thoraxes of the mosquitoes used for the ELISA CSP testing were ground in 200 µL of 2% CTAB. After 5 min in a water bath at 65 • C, the grindings were mixed with 200 µL of chloroform and then centrifuged at 14,000 rpm for 5 min. The supernatant was gently transferred into another tube containing 200 µL of isopropanol then centrifuged at 12,000 rpm for 15 min. The DNA was purified with 200 µL of 70% ethanol. Another centrifugation occurred at 14,000 rpm for 5 min. The content of each tube was gently inverted to preserve the DNA, which will then be dried for at least 3 h on the bench. Finally, 20 µL of bi-distilled water was added to the DNA, which is left on the bench overnight or half a day.
Then, PCR amplification of the Plasmodium species was performed in a thermocycler (Eppendorf, Hamburg, Germany) following PCR conditions described in Methods in Anopheles Research [25].
The identification of Plasmodium species was performed using the primers presented in Table 1. UNR and PLF are universal primers that hybridize to a particular nucleotide sequence in all Plasmodium species. FAR forms with UNR and PLF, the pair of primers that amplify a specific nucleotide fragment in P. falciparum. The same applies for the pair of primers formed by UNR-PLF and VIR, which amplifies a specific nucleotide fragment in P. vivax. Then, the PVCSP1 and PVCSP2 primers were used to identify the P. vivax variants ( Table 1). DNA of P. falciparum, P. vivax, P. vivax 210, and P. vivax 247 (BEI Resources) was used as positive controls, while sterilized water was used as a negative control. Total of 18 µL of each amplified product for each reaction was analyzed with electrophoresis gel using a 2% agarose gel (Sigma). When successfully amplified, the gel images were recorded and visualized on transilluminator EBOX 1000 (Vilber, Marne-la-Vallée, France).

Molecular Species Identification in the Anopheles gambiae Complex
The abdomens, legs, and wings of the same mosquito individuals previously tested through PCR for detection of Plasmodium infection were molecularly speciated using the protocol of Santolamazza et al. [26]. PCR amplification of Anopheles gambiae complex species was carried out in a 25 µL reaction mixture containing 5 µL genomic DNA for An. gambiae s.l., 1 × PCR buffer 5×, 1.5 mM MgCl 2 , 0.2 mM of each dNTP, 10 pmoles of each primer (Table 1), and 5 U of Taq polymerase (Eurogentec, Liège, Belgium). The amplification was performed in a thermocycler (Eppendorf, Hamburg, Germany) as follows: an initial denaturation of the DNA material at 94 • C for 5 min, followed by 35 cycles of 94 • C for 30 s, 54 • C for 30 s, and 72 • C for 30 s, with a final extension at 72 • C for 10 min. DNA of An. gambiae, An. coluzzii, and An. arabiensis were used as positive controls, while sterilized water was used as a negative control. The amplification product was migrated on a 1.5% agarose gel (Sigma) in the presence of ethidium bromide (Sigma) used as an intercalant. This allows for the assessing of the involvement of each species of the An. gambiae complex in the transmission of the different Plasmodium species.

Data Analysis
Data were analyzed using R Core Team software (version 4.1.  and Graph Pad Prism software version 5.0. The graphs were made using Graph Pad software version 5.0. The exact binomial method was used to calculate the prevalence and their confidence intervals with the R Core Team software (version 4.1.3-2022).

Ethical Consideration
The present study received ethical approval from the Institutional Ethics Committee of Centre de Recherche Entomologique de Cotonou (N • 04/IECC). Mosquito collectors gave their consent prior to being involved in the study. They were trained to collect the mosquito before being bitten and were regularly subjected to medical check-ups. When they had malaria, they were taken care of and received anti-malarial medication like Artemisinin-based Combination Therapy if required. They were also vaccinated against yellow fever. Table 2 shows the species composition of mosquito populations collected through HLCs in the different surveyed regions over the study period. Overall, a great diversity was observed in mosquito fauna collected on human bait, regardless of the commune. In total, 32,773 specimens of mosquitoes belonging to 18 species were collected using HLCs. Overall, the most predominant species in the whole study area were Culex quinquefasciatus (40.48%), followed by Anopheles gambiae s.l. (20.92%), Culex gr decens (10.23%), Mansonia africana (18.27%), Anopheles funestus (3.94%), and Aedes aegypti (2.61%) ( Table 1). At the region level, Anopheles gambiae s.l. was the most abundant mosquito species in Zou (63.9%), Donga (56.9%), Atacora (50.8%), Mono (39.1%), and Couffo (35%), while Anopheles funestus was the most abundant one in Collines (60.9%). In the Ouémé, Littoral, Atlantique, and Borgou regions, Cx quinquefasciatus was the major species with frequencies of 51.7%, 78.6%, 43.8%, and 48.3%, respectively. Mansonia africana was predominant in the Plateau (59%) and Alibori (57.4%) regions (Table 2).

Blood Meal Origin in An. gambiae s.l. Collected through PSC
Of the 1392 mosquitoes collected by PSC, 514 were specimens of An. gambiae s.l., of which 311 (60.5%) were blood-fed. Out of the 311 blood-fed specimens of An. gambiae s.l. that were tested through ELISA blood meal in the study area, 95.5% took their blood meal on humans (Table 3).
Approximately 16.67% and 7.14% of the tested individuals of An. gambiae s.l. had successful blood meal intake on cow in the communes of Kouandé and Karimama, respectively. In Athiémé and Karimama, the proportion of of An. gambiae s.l that blood-fed on sheep was 4.87% and 3.57%, respectively. Blood meal intake on pigs occurred in 1.69%, 8.33%, and 1.61% of specimens of An. gambiae s.l analyzed in Aguégués, Sô-Ava, and

Blood Meal Origin in An. gambiae s.l. Collected through PSC
Of the 1392 mosquitoes collected by PSC, 514 were specimens of An. gambiae s.l., of which 311 (60.5%) were blood-fed. Out of the 311 blood-fed specimens of An. gambiae s.l. that were tested through ELISA blood meal in the study area, 95.5% took their blood meal on humans (Table 3).
Approximately 16.67% and 7.14% of the tested individuals of An. gambiae s.l. had successful blood meal intake on cow in the communes of Kouandé and Karimama, respectively. In Athiémé and Karimama, the proportion of of An. gambiae s.l that blood-fed on sheep was 4.87% and 3.57%, respectively. Blood meal intake on pigs occurred in 1.69%, 8.33%, and 1.61% of specimens of An. gambiae s.l analyzed in Aguégués, Sô-Ava, and Athiémé, respectively. In the commune of Sô-Ava, 4.17% of tested individuals of An. gambiae s.l took their blood meal on both humans and pigs (Table 3). 3.4. Prevalence of P. falciparum, P. vivax 210 and 247 CSP Antibodies in An. gambiae s.l.

Plasmodium Detection by ELISA
A total of 3300 head-thoraxes of An. gambiae s.l. out of the 6857 collected were tested through ELISA CSP for each of the three Plasmodium species, equating to a total of 9900 samples analyzed. The number of sporozoite-positive mosquitoes for each Plasmodium species at the commune level is mentioned on Figure 3. The combined data showed that 86 mosquitoes were positive for P. falciparum, 9 for P. vivax 210, and 7 for P. vivax 247. Thus, the average prevalence (sporozoite index: SI) for P. falciparum was 2.6% (CI: 2.1-3.1) compared to 0.3% (CI: 0.1-0.5) for P. vivax 210 and 0.2% (CI: 0.1-0.4) for P. vivax 247 in the whole study area.
Out of the 228 specimens of An. funestus analyzed, only two (one in Ouèssè and one in Savè) were positive for P. falciparum. No specimen of An. funestus was found positive for P. vivax 210 and P. vivax 247.

Distribution of the Plasmodium Species in the Sibling Species of the Anopheles gambiae Complex
Among the 115 specimens of Anopheles gambiae s.l found positive for P. falciparum through PCR, 64.35% (n = 74) were An. gambiae, 34.78% (n = 40) were An. coluzzii, and 0.86% (n = 1) were An. arabiensis. Similarly, the P. falciparum infectivity rate was significantly higher in An. gambiae than in An. coluzzii (p < 0.001).
At the opposite, most P. vivax 210-positive mosquitoes were An. coluzzii (76.92%). The rest were An. gambiae. The infectivity rates of P. vivax 210 were significantly higher in An. coluzzii than in An. gambiae (p  0.05).
Concerning P. vivax 247, only An. coluzzii was positive to this variant (Table 4). As observed with the CSP ELISA test, no co-infection was recorded. In Grand Popo, Aplahoué, and Banikoara, P. falciparum sporozoite-positive mosquitoes were found through PCR, which was not the case with the ELISA CSP test. The same trend was observed with P. vivax 210 in Karimama.

Distribution of the Plasmodium Species in the Sibling Species of the Anopheles gambiae Complex
Among the 115 specimens of Anopheles gambiae s.l found positive for P. falciparum through PCR, 64.35% (n = 74) were An. gambiae, 34.78% (n = 40) were An. coluzzii, and 0.86% (n = 1) were An. arabiensis. Similarly, the P. falciparum infectivity rate was significantly higher in An. gambiae than in An. coluzzii (p < 0.001).
At the opposite, most P. vivax 210-positive mosquitoes were An. coluzzii (76.92%). The rest were An. gambiae. The infectivity rates of P. vivax 210 were significantly higher in An. coluzzii than in An. gambiae (p < 0.05).
Concerning P. vivax 247, only An. coluzzii was positive to this variant (Table 4).

Discussion
In the current context of malaria elimination, the identification of Plasmodium species in humans and vectors is a prerequisite not only to understand the epidemiology of the disease, but also to implement an effective control strategy against these parasites. The present study shows the presence of antibodies of P. falciparum and P. vivax in the An. gambiae complex, the main malaria vector complex in Benin. This result corroborates the work of Poirier et al. [16], who showed a surprisingly high exposure of Beninese to P. vivax, resulting in submicroscopic infections.
The molecular species identification performed in An. gambiae s.l. revealed the presence of three sibling species, namely An. gambiae, An. coluzzii, and An. arabiensis. These species were found in sympatry but in variable proportions. Overall, An. gambiae accounted for 53.68% of all specified mosquitos, followed by An. coluzzii (45.88%) and An. Arabiensis, which was found at a very low frequency. This result is similar to the one of Aikpon et al. [33] that showed the same trend in Atacora department, the cotton zone of Benin. The opposite result, displaying An. coluzzii as the predominant species, was observed by Gnanguenon et al. [30] in Allada, Dassa, and Malanville. Although these species often have different geographical and seasonal distributions, they can also live in sympatry. In the present study, An. arabiensis was found only in Savè, a wet savanna area. This is consistent with works by Gnanguenon et al. [30], who found An. arabiensis in Dassa, a commune that is quite close to Savè. Over past years, An. arabiensis was considered a dry savanna species found in abundance between the extreme north of Benin and Parakou [34]. It is, therefore, necessary to assess if the species migrated from its preferred area to a more southern area. A possible migration could be explained by factors such as climate change, the destruction of the vector's ecological zones for agricultural and hunting purposes, the harvesting of timber, and pastoralism.
The findings of the present study revealed spatial heterogeneity in the distribution of Plasmodium species prevalent in Benin. The mean prevalence of antibodies of P. falciparum, P. vivax 210, and P. vivax 247 in An. gambiae s.l. were, respectively, 2.6%, 0.30%, and 0.2% in the study area. In Cotonou, the three Plasmodium species were found in An. coluzzii with a high prevalence for P. vivax 247 (4%) compared to P. vivax 210 (2%) and P. falciparum (2%). In the communes of Aguégués, Porto-Novo, Ouidah, Agbangnizoun, Athiémé, Karimama, and Sinendé, P. falciparum and P. vivax 210 were detected with similar frequencies across communes (p > 0.05). In Sô-Ava and Dogbo, only P. falciparum and P. vivax 247 were registered with similar frequencies (p > 0.05). In the rest of the communes, P. falciparum was the sole species found except for Pobè, where no infection to three searched Plasmodium species was observed. The null infection observed in Pobè could be due to the low number of mosquitos tested in this commune. This variation of the prevalence of the Plasmodium species across communes can be explained by the period of collection of the mosquitoes, as well as the ecological and evolutionary factors associated with mosquitoes, and that can promote divergence and restrict gene flow between vector-adapted parasite lineages [35]. Furthermore, contrary to the findings of Sandeu et al. [36], no co-infection with different species of Plasmodium cells was recorded in Anopheles mosquitoes analyzed in the communes surveyed over the present study. In addition, over the trial of Sandeu et al. [36] where the different species of Plasmodium were searched and quantified in malaria vectors, P. vivax was not found in any of the mosquito samples analyzed, contrary to what is shown in the present study. Indeed, Howes et al. [37] previously showed that P. vivax was the most prevalent malaria parasite globally but clearly "rare" in Africa. The majority of African populations do not express the Duffy blood group antigen, which is the only known receptor for the P. vivax infection. Since this discovery in the 1970s, the low clinical incidence of P. vivax in Africa led to a perception that this Plasmodium species is completely absent, and that all apparent cases are from a misdiagnosis. Thus, no public health consideration is made for this parasite in terms of diagnosis, treatment, or surveillance. As more sensitive diagnosis methods become available, P. vivax infection in Africa is increasingly being reported through various types (entomological, serological, community prevalence) of surveys, as well as clinical data on infections from local residents and international travelers.
It would have been interesting to separate the infection data of mosquitoes according to the biting locations (indoor and outdoor), which we did not do during the trial. Thus, it would be important to take into account the microclimatic conditions (ambient temperature and relative humidity) of the indoor and outdoor environment where collections were undertaken as it matters to the vector availability and density. This constitutes a limitation for the study.
In the present study, the distribution of the different Plasmodium species in the sibling species of An. gambiae s.l. was also assessed across the surveyed communes. Of the 115 samples positive for P. falciparum through PCR, An. gambiae was predominant (n = 74), followed by An. coluzzii (n = 40) and An. arabiensis (n = 1). The infectivity of An. coluzzii to P. vivax (variant 210 and 247) was significantly higher than that of An. gambiae. Indeed, most (76.92%) and all mosquitoes that were positive for P. vivax 210 and P. vivax 247, respectively, were An. coluzzii. This result contradicts findings from Goupeyou-Youmsi et al. [38] obtained in two neighboring villages of the rural commune of Andriba in Madagascar. Indeed, they found that An. arabiensis, which was the sole species detected within the An. gambiae complex, carried P. falciparum, P. vivax, and P. malariae. Contrary to the present study, only An. arabiensis was the vector of three Plasmodium species (P. falciparum, P. vivax, and P. malariae). These differences could be explained by weather conditions, the ecological environment, and/or the geographical distribution of vectors that are different between countries [38]. These results suggest that more screening of An. coluzzii has to be undertaken, as P. vivax is difficult to eliminate, and more ultra-sensitive methods are required to detect such infections in the community. This is a valid point that epidemiologists and entomologists need to worry about since the urban malaria vector, An. stephensi, which invaded many parts of sub-Saharan Africa, such as Ethiopia and Sudan [39,40], Somalia, and more recently Nigeria [41], could also be infected with P. vivax. For that, the larval habitats of the malaria vectors in Benin need to be documented as efforts for larval source management, along with adult control interventions (LLINs and targeted IRS) that can be strengthened in order to have positive results.
In the present study, PCR allows for the detection of sporozoites of both P. falciparum and P. vivax in some mosquitoes, while ELISA-CSP gave negative results. A possible explanation would be that the number of sporozoites present in the salivary glands of the mosquito is too low to be detected by the ELISA-CSP test. In addition, several studies have also mentioned the difficulty of ELISA in detecting the antigens of other Plasmodium species such as P. malariae and P. ovale, but also the problem of antigenic variation of the CSP antigen between geographical areas [42,43]. Finally, an extension of this study in all geographical areas of the country and other West African regions could help generate valuable epidemiological data on the two variants of P. vivax and allow for a more elaborated malaria pre-elimination policy. The evaluation of P. vivax variants is a prerequisite to the development of an effective vaccine. Other aspects that are worth evaluation include the geographical distribution the variants, the intensity of their transmission, the ability of the vectors to allow for their development, and the variability of the human immune response towards them.

Conclusions
The present study constitutes a first in Benin that detects P. falciparum and P. vivax in An. gambiae s.l. using both CSP ELISA and PCR techniques. The results showed that the risk of exposure of populations to P. falciparum, P. vivax 210, and P. vivax 247 varies from species to species. In addition, both variants of P. vivax are present in Benin. It would, therefore, be necessary to evaluate the impact of currently used antimalarial drugs on this Plasmodium species for a better therapeutic management of malaria in the country. These results will help the National Malaria Control Program to better plan for the therapeutic management strategy of malaria cases to move towards the pre-elimination of the disease.