Urban agriculture is a widely practiced farming activity involving more than 800 million people, worldwide . Growing population and rural migration to the cities of African countries have resulted in the practice of urban and peri-urban agriculture to provide food security, to overcome poverty and create employment. Although most of these urban inhabitants are engaged in subsistence gardening, more than 200 million practice market-oriented farming on undeveloped urban spaces .
The same trend is observed in Ghana where the population of Accra metropolis has grown from 1.6 million people in 2000 to about 2.27 million in 2015 . While it is critical meeting urban challenges such as food security and unemployment, these practices also present a serious public health problem . The economic and social value of urban and peri-urban agriculture is hindered by a number of factors including the proliferation of mosquito breeding sites   with potentially higher malaria epidemiological risk . Malaria transmission in urban settings such as in Accra is influenced by several factors including housing, availability of health care, knowledge on vector control strategies and vector susceptibility to insecticides.
Intensive use of insecticides for both public health and agricultural pest control has led to reports of global spread of insecticide resistance. The trend is threatening the efficacy of malaria vector control tools and even worsened by the use of the same insecticide or class of insecticides     . Ghana is not exempted from this phenomenon particularly where pyrethroid resistance has been reported in many areas in the country    . Recently, there is a scaling up of Long Lasting Insecticidal Net (LLIN) distribution and Indoor Residual Spraying (IRS) coverage in the country . Resistance to pyrethroids by malaria vectors is known to be driven by selection pressure resulting from the increased use of LLINs and IRS   and has also been attributed to the contamination of mosquito breeding habitats during the application of agricultural pesticides. Agricultural environments can create mosquito breeding habitats, especially when irrigated, thereby increasing vector density   . Anopheles gambiae s.s. and Anopheles coluzzii constitute the main malaria vectors in Ghana and particularly in the Southern part of the country where Accra is located. Anopheles gambiae s.s., (formerly named A. gambiae s.s. molecular form) prefers rain-dependent pools, sunny and temporally puddles while A. coluzzii (formerly A. gambiae s.s. M molecular form) exploits preferentially immature sites that exist across seasons and also linked to human activities such as rice cultivation, urbanization and irrigation .
Malaria incidence in urban areas has generally been low compared to rural areas. However, data collected throughout the decade indicate that urban malaria remains a major public health problem      . It is therefore important to understand the factors that promote urban malaria to help formulate vector control strategies and help manage insecticide resistance in urban areas. The aim of this study was to determine the distribution of malaria vectors and the contribution of urban agriculture to the development of insecticide resistance in some selected areas within Accra Metropolis.
2. Material and Methods
2.1. Study Sites
The study was conducted in Accra, the capital of Ghana (5˚32'59.99''N and 0˚12'60.00''E). Accra is the country’s biggest metropolis, most diverse and most cosmopolitan with an estimated population of 2.27 million of people as in 2015. Accra is in the coastal savanna ecological zone, characterized by two rainfall peaks, the first occurring from April to July and the second from September to October. The region receives an average rainfall of 740 mm to 890 mm per year. The lowest mean monthly temperature (about 26˚C) is recorded during August and the highest (about 30˚C) between March and April. The mean relative humidity throughout the year ranges between 65% and 75%.
For the study, most of the urban and peri-urban vegetable growing areas were mapped using a 3D global positioning system (GPS) to determine the specific geographical coordinates or geo-location of the sites. Two vegetable growing sites were selected among them and paired with two non-vegetable growing areas close to each vegetable growing site and served as controls. The sites included Korle Bu (5˚32'20.124''N; 0˚14'10.115''E) and Opeibea (5˚35'54.54''N; 0˚10'53.609''E) as vegetable growing sites and Mataheko (5˚34'0.4''N; 0˚15'13.8''W), and Achimota (5˚36'59.99''N; 0˚13'60.00''E) as non-vegetable growing sites (Figure 1). A questionnaire survey on insecticides and herbicides classes and type was conducted to determine the extent of usage of insecticides.
Figure 1. Map of the study sites including the A. gambiae s.l. species composition and the frequencies of the target site mutations within each population.
2.2. Mosquito Sampling and WHO Susceptibility Test
Larvae and pupae were collected from the study sites between January and April 2017. The mosquito larvae were reared to adults in the insectary of Vestergaard-NMIMR Vector Labs (VNVL) at the Noguchi Memorial Institute for Medical Research (NMIMR) using standard rearing methods .
Susceptibility tests of WHO were conducted according to WHO standard protocols (WHO, 2016). Mosquitoes were assayed using WHO discriminating dosages of five insecticides: 0.05% deltamethrin, 0.75% permethrin, 4% DDT, 5% malathion and 0.1% bendiocarb. The insecticides were selected to match with the four classes of insecticides. In addition, synergist assays with 4% Piperonyl butoxide (PBO) impregnated papers were conducted to determine the probable presence of metabolic mechanisms such as P450 enzymes. Twenty to twenty-five non-blood fed female A. gambiae s.l., aged 3 - 5 days were exposed for one hour to the different insecticides before and holding period of 24 hours for mortality. For the synergist assay, mosquitoes were pre-exposed to PBO for one hour before exposure to the insecticide for an additional hour. The number of mosquitoes knocked down was recorded at 60 minutes and mortality recorded after 24 hours post exposure for both tests. Tests with silicone and olive oil impregnated papers were run in parallel and served as controls. After the susceptibility tests, all mosquitoes (including controls) were kept at −20˚C for further identification of A. gambiae species complex and characterization of the voltage-gated sodium channel 1014F (Vgsc-1014F) and acetylcholinesterase G119S (ace-1) mutations.
2.3. DNA Extraction and Species Characterization
Forty mosquito samples kept at −20 ˚C after WHO susceptibility test were randomly selected per site for molecular characterization. The whole DNA of every single mosquito was extracted using a modified CTAB method described by Collins et al. . Thereafter, the A. gambiae species identification was done using the protocol of Fanello et al.  and the primers UN [GTGTGCCGCTTCCTCGATGT], AG [CTGGTTTGGTCGGCACGTTT], AR [AAGTGTCCTTCTCCATCCTA] and AM [GTGACCAACCCACTCCCTTGA] and Santolamazza et al.  with primers F6.1a [5’-TCGCCTTAGACCTTGCGTTA-3’] and R6.1b [5’-CGCTTCAAGAATTCGAGATAC-3’]. The Vgsc-1014F and ace-1 allele of each mosquito were characterized following the protocol of Martinez-Torres et al.  using the primers AGD1 [5’-ATAGATTCCCCGACCATG-3’]; AGD2 [5’-AGACAAGGATGATGAACC-3’], AGD3 [5’-AATTTGCATTACTTACGACA-3’] and AGD4 [5’-CTGTAGTGATAGGAAATTTA-3’], and Weill et al.  with primers EX3AGdir [5’-GATCGTGGACACCGTGTTCG-3’] and EX3AGrev [5’-AGGATGGCCCGCTGGAACAG-3’] respectively.
Furthermore, 50 non-exposed mosquitoes per site were kept at −85˚C and the enzyme activities of each population were characterized following the methods described in MR4  , and compared with the susceptible strain of A. gambiae Kisumu. These enzymes included oxidases (P450s), esterases (α and β), acetylcholine esterase (Ache), proteins and Glutathione S-Transferases (GSTs).
2.4. Data Analysis
The resistance status of each A. gambiae mosquito population was determined following the WHO criteria, where mortality ≤ 90% indicates that the colony is resistant to the insecticide tested, >90% and <97%, the resistance is suspected and ≥98%, the colony is susceptible (WHO, 2016).
The Vgsc-1014F and ace-1 frequencies were calculated using Hardy Weinberg formula and compared among populations using the z-test of proportion.
The enzyme activities were estimated as the amount of product per minute per milligram of protein. The mean values were plotted and compared using the non-parametric test of Kruskal Wallis of the Graph Pad Prism 5.0 software.
3.1. Insecticide Usage and A. gambiae s.l. Resistance Status
The results of the insecticide usage survey are described in Appendix. Almost all the classes of insecticides were found within the pesticides and herbicides been used for controlling pests in the vegetable growing areas. However, the organophosphates were the most common insecticides found followed by the pyrethroids.
The WHO susceptibility test results are described in Figure 2. The resistance to deltamethrin was significantly higher in Opeibea and Korle-Bu (2.6% and 5.7% respectively), the vegetable growing sites than the non-vegetable growing
Figure 2. Mortality of A. gambiae s.l. of the different sites against various insecticides.
sites Achimota and Mataheko (36.2% and 38.9% respectively) (p < 0.05). A similar trend was observed with permethrin though the mortalities were relatively lower than deltamethrin in all the sites. The mortalities against permethrin ranged between 2.5% and 13.1% for both vegetable and non-vegetable sites. When A. gambiae s.l. was pre-exposed to PBO, a significant enhancement of susceptibility was observed with both insecticides and in all the sites. However, the synergistic action of PBO was significantly higher for deltamethrin than permethrin in all the sites. The mortalities against PBO + deltamethrin in the vegetable growing sites were 56.5% and 61.5% for Opeibea and Korle-Bu respectively, while PBO + permethrin showed mortalities of 46.8% and 14% respectively in the same sites. Similarly, the non-vegetable sites Achimota and Mataheko showed 76.5% and 93.2% for PBO + deltamethrin respectively compared to 44.1% and 66.2% for PBO + permethrin respectively. All the populations of A. gambiae s.l. from all the sites were highly resistant to DDT with mortalities ranging between 1.2% at Korle-Bu and 7% in Achimota. The other insecticide classes including the carbamate (bendiocarb) and organophosphate (malathion) recorded high resistance in all the sites. Though the mortality against bendiocarb was low, malathion in contrast recorded a relatively higher mortality of 64.6% at Korle-Bu.
Overall, no significant difference was observed between the percentage mortalities of the mosquitoes exposed to insecticides across sites (p = 0.493).
3.2. Molecular Characterization of A. gambiae s.l. Populations of the Sites
The results of the species identification and resistance allele mutations are shown in Figure 1. The majority (95%) of A. gambiae s.s. occurred in Opeibea compared to 97.5% of A. coluzzii in Korle-Bu. The Vgsc-1014F frequency of A. gambiae in Opeibea and Korle-Bu were similar (0.91 vs 0.88). In contrast, the ace-1 mutation was significantly higher (p = 0.001) within the population of Opeibea (0.79) than Korle-Bu (0.45).
For the non-vegetable growing areas, the populations were living in sympatry with A. gambiae s.s. and A. coluzzii constituting 62.5% and 37.5% respectively in Achimota. In Mataheko, A. coluzzii and A. gambiae s.s. formed 75% and 12.5% of the species respectively. Hybrids of both species (12.5%) were found in Mataheko. The frequency of the Vgsc-1014F is also similar in both non-vegetable growing sites with 0.88 and 0.81 observed in Achimota and Mataheko respectively. The ace-1 frequencies were 0.63 for Achimota and 0.42 for Mataheko.
3.3. Biochemical Analysis
Figure 3. Comparison of the different enzyme activities of A. gambiae s.l. from the study sites.
Table 1. Enzyme activities of A. gambiae s.l. populations of the different sites compared to the susceptible A. gambiae s.s. Kisumu.
Degree of significance of the p value is 5%. (p value in bold are significantly lower than the reference A. gambiae Kisumu).
A. gambiae s.l. from Opeibea recorded the highest level of enzyme activity compared to the other sites. P450s, α-esterase and acetylcholinesterase were significantly higher in Opeibea (p < 0.05) than both non-vegetable growing sites and Korle-Bu populations. However, the expression of β-esterase activity was relatively higher mostly in the A. gambiae s.l. population of the non-vegetable growing areas.
Farming activities in West African urban areas (capitals and big cities) have an important impact on the density of mosquito populations   . Nowadays, urban and peri-urban vegetable growing activities contribute significantly to the gross domestic product (GDP) of countries and serve as sources of income for citizens who are engaged in these activities. Nonetheless, these urban vegetable-growing activities involving indiscriminate use of insecticides has led to the creation of favorable breeding sites for Anopheles mosquitoes and the development of insecticide resistance. These are major factors that contribute to urban malaria transmission in cities such in Accra .
In this study, high insecticide resistance of A. gambiae s.l. to all the classes of insecticides was observed in the vegetable growing areas and particularly to pyrethroids and organochlorine (Figure 2). The resistance to deltamethrin and permethrin was more than ten-fold greater in the vegetable growing areas (Opeibea and Korle-Bu) than thenon-vegetable growing areas (Achimota and Mataheko).
This observation confirms work done by Yadouleton et al. in Benin and Kabula et al. in Ghana   showing high insecticide resistance of A. gambiae s.l. within urban vegetable growing areas. Additionally, the demand of freshly produced vegetables by the urban market has increased the development of vegetable growing spaces and the use of insecticides and pesticides. These factors might have contributed to the development of insecticide resistance of A. gambiae s.l. mosquitoes that is observed. However, the synergistic effect using Piperonyl butoxide (PBO) and deltamethrin was very effective in enhancing susceptibility of these resistant mosquitoes, thereby increasing the mortality of A. gambiae s.l. collected from both vegetable and non-vegetable growing areas. Nevertheless, there were differences in the percentage increment of susceptibility by PBO in both collection sites. In the vegetable growing areas, the mortalities were increased by five and twenty folds using PBO + deltamethrin in Korle-Bu and Opeibea respectively while it doubled in the non-vegetable growing areas (Mataheko and Achimota).
The same trend was observed with PBO + permethrin where the percentage of susceptibility enhancement in A. gambiae s.l. from Opeibea, a vegetable growing area was significantly higher compared to the non-vegetable growing sites. The results also showed low level of enhancement of susceptibility when mosquitoes from Korle-Bu were pre-exposed to PBO before permethrin compared to Opeibea. PBO + permethrin showed three folds enhancement of susceptibility in A. gambiae s.l. from Korle-Bu, whilst that of Opeibea was about twenty folds (Figure 2). The difference observed in the ability of PBO to increase the mortality of the two populations of A. gambiae s.l. might be due to difference in mosquito species composition and the non-involvement of metabolic resistance mechanisms. Anopheles coluzzii was abundant at Korle-Bu while A. gambiae s.s. dominated the population of Opeibea (Figure 1). Furthermore, the latter species was observed to exhibit higher P450 enzyme activity than A. coluzzii from this study (Figure 3, Table 1) and this observation also confirms the fact that A. gambiae s.s. was more resistant than A. coluzzii  . Piperonyl butoxide (PBO) is known as a synergist which acts as an inhibitor of the P450 enzyme activities .
The level and probable involvement of P450s in the resistance status of the mosquitoes from all the study sites have been first described using the synergistic effect of PBO  and was further confirmed by the characterization of the level of the different enzyme activities including P450s, esterases, Glutathione S-transferase (GST) and acetylcholine esterase (Ache). The high level of P450s activities observed mainly in A. gambiae s.l. from Opeibea confirmed the involvement of oxidases in the insecticide resistance level of the mosquitoes (Figure 3 and Table 1).
Relatively high resistance to other insecticides including bendiocarb and malathion was observed in both vegetable and non-vegetable growing areas. The slightly high level of resistance observed in the non-vegetable growing areas is an indicator of the spread of insecticide resistance which may be sustained by the household use of vector control products such as mosquito coils and aerosols as well as insecticide treated nets. The study also showed that the two major malaria vectors A. gambiae s.s. and A. coluzzii occurred in sympatry though in significantly different proportions within the various sites (Figure 1). The two species were specially found in larger proportion in the vegetable growing areas with A. gambiae s.s. representing more than 95% of the total population in Opeibea while 98% of the anopheles from Korle-Bu were A. coluzzii. The difference in the species composition could be due to the geographical position of the two vegetable growing sites. Korle-Bu is located at a lower altitude than Opeibea and the vegetation is relatively dense while Opeibea is a sunny and opened area. The difference in the ecological niche might have shaped their adaptation since A. gambiaes.s.is more adapted to sunny areas than A. coluzzii . In addition, the species distribution in the non-vegetable growing areas seemed to have been driven by the proximity of the population to the vegetable growing sites. A higher proportion of A. gambiae s.s. was observed at Achimota which was close to Opeibea, while Mataheko showed high percentage of A. coluzzii due to its proximity to Korle-Bu. However, this observation need to be further investigated and confirmed using a genomic tree to characterize the different population structures and their association .
The voltage-gated sodium channel (Vgsc-1014F) which is known to confer resistance to pyrethroids and DDT was highly observed in all the sites including the non-vegetable growing areas (Figure 1). This is a call for concern since this observation is indicative of widespread of Vgsc mutation within the A. gambiae s.l. populations in Accra. Over the years, various studies have shown gradual increase of the frequencies of Vgsc mutation in Anopheles gambiae s.l. from Accra and the whole country      and probably sustained by the intense use of pyrethroids in most of the vector control tools. It is therefore not surprising that high frequencies of resistance alleles were observed in the non-vegetable growing sites purported to be residential areas.
The ace-1 mutation, known to confer resistance to carbamates and organophosphates showed a relatively higher frequency in vegetable growing areas than the non-vegetable growing sites (p < 0.05). This could be due to the massive use of carbamate and organophosphate-based pesticides to control agricultural pests as confirmed by the results of the pesticide use survey in the study areas (Appendix). It is therefore important to properly manage urban vegetable growing activities and use of pesticides to prevent further spread of insecticide resistance and eventual failure of vector control strategies using these insecticides. The use of insecticides to control pest in agriculture obviously contributes to the spread of insecticide resistance of malaria vectors  , which was also observed within urban and peri-urban vegetable growing activities in Accra. Therefore, the monitoring of insecticide resistance and associated mutations in vector populations are relevant for early detection and prompt decision making by appropriate authorities.
The distribution of A. gambiae s.l. indicated that both A. gambiae s.s. and A. coluzzii occurred in sympatry in both study areas. The Vgsc-1014F is highly and widely spread within the A. gambiae s.l. populations in Accra, irrespective of whether the site is a vegetable growing area or not. The ace-1 was also found in both settings but in a significantly higher proportion within the vegetable growing areas. Agricultural activities using pesticides and household vector control influence insecticide resistance in malaria vector populations. The findings of this study will be useful in understanding the contribution of urban agriculture to the increase in urban malaria transmission and help formulate policies for future vector control in urban areas.
Limitation: While the study was able to emphasize the impact of the use of agricultural pesticides on the resistant status of the malaria vectors within Accra, it was in contrast difficult to estimate the correlation between that trend and household use of protective measures using long lasting insecticidal nets (LLINs) and other treated materials.
This work was financially supported by Vestergaard and the Department of Parasitology of NMIMR. We are also grateful to Mr. Duncan K. Athinya for the map of the study site.
JC, KKF and SKD designed, implemented the study, analyzed, interpreted data and drafted the manuscript. MCE, RP, JJ, DO, GA, DAB, CA, AI, SA, and SG carried out mosquito collections and all laboratory experiments. SDK, KKF, MPH and HPJ revised the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Table A1. Class and active ingredient of agrochemicals used by the farmers in both vegetable growing areas.
 Drechsel, P., Graefe, S., Sonou, M. and Cofie, O.O. (2006) Informal Irrigation in Urban West Africa: An Overview. IWMI Research Report 102, International Water Management Institute, Colombo, 40 p.
 Yadouleton, A.W., Asidi, A., Djouaka, R.F., Braima, J., Agossou, C.D. and Akogbeto, M.C. (2009) Development of Vegetable Farming: A Cause of the Emergence of Insecticide Resistance in Populations of Anopheles gambiae in Urban Areas of Benin. Malaria Journal, 8, 103.
 Afrane, Y.A., Klinkenberg, E., Drechsel, P., Owusu-Daaku, K., Garms, R. and Kruppa, T. (2004) Does Irrigated Urban Agriculture Influence the Transmission of Malaria in the City of Kumasi, Ghana? Acta Tropica, 89, 125-134.
 Afrane, Y.A., Lawson, B.W., Brenya, R., Kruppa, T. and Yan, G. (2012) The Ecology of Mosquitoes in an Irrigated Vegetable Farm in Kumasi, Ghana: Abundance, Productivity and Survivorship. Parasites & Vectors, 5, 233.
 Donnelly, M.J., McCall, P.J., Lengeler, C., Bates, I., D’Alessandro, U., Barnish, G., et al. (2005) Malaria and Urbanization in Sub-Saharan AFRICA. Malaria Journal, 4, 12.
 Czeher, C., Labbo, R., Arzika, I. and Duchemin, J.B. (2008) Evidence of Increasing Leu-Phe Knockdown Resistance Mutation In Anopheles gambiae from Niger Following a Nationwide Long-Lasting Insecticide-Treated Nets Implementation. Malaria Journal, 7, 189.
 Diabate, A., Baldet, T., Chandre, F., Akogbeto, M., Guiguemde, T.R., Darriet, F., et al. (2002) The Role of Agricultural Use of Insecticides in Resistance to Pyrethroids in Anopheles gambiae s.l. in Burkina Faso. The American Journal of Tropical Medicine and Hygiene, 67, 617-622.
 N’Guessan, R., Corbel, V., Akogbeto, M. and Rowland, M. (2007) Reduced Efficacy of Insecticide-Treated Nets and Indoor Residual Spraying for Malaria Control in Pyrethroid Resistance Area, Benin. Emerging Infectious Diseases, 13, 199-206.
 Adasi, K. and Hemingway, J. (2008) Susceptibility to Three Pyrethroids and Detection of Knockdown Resistance Mutation in Ghanaian Anopheles gambiae sensu stricto. Journal of Vector Ecology, 33, 255-262.
 Chabi, J., Baidoo, P.K., Datsomor, A.K., Okyere, D., Ablorde, A., Iddrisu, A., et al. (2016) Insecticide Susceptibility of Natural Populations of Anopheles coluzzii and Anopheles gambiae (sensu stricto) from Okyereko Irrigation Site, Ghana, West Africa. Parasites & Vectors, 9, 182.
 Muller, P., Donnelly, M.J. and Ranson, H. (2007) Transcription Profiling of a Recently Colonised Pyrethroid Resistant Anopheles gambiae Strain from Ghana. BMC Genomics, 8, 36.
 Yawson, A.E., McCall, P.J., Wilson, M.D. and Donnelly, M.J. (2004) Species Abundance and Insecticide Resistance of Anopheles gambiae in Selected Areas of Ghana and Burkina Faso. Medical and Veterinary Entomology, 18, 372-377.
 Matowo, J., Kitau, J., Kaaya, R., Kavishe, R., Wright, A., Kisinza, W., et al. (2015) Trends in the Selection of Insecticide Resistance in Anopheles gambiae s.l. Mosquitoes in Northwest Tanzania during a Community Randomized Trial of Longlasting Insecticidal Nets and Indoor Residual Spraying. Medical and Veterinary Entomology, 29, 51-59.
 Keiser, J., Utzinger, J., Caldas de Castro, M., Smith, T.A., Tanner, M. and Singer, B.H. (2004) Urbanization in Sub-Saharan Africa and Implication for Malaria Control. The American Journal of Tropical Medicine and Hygiene, 71, 118-127.
 Kone, A.B., Tiembre, I., Cisse, G., Diallo, A., Tanner, M. and N’Goran, K.E. (2015) The Impact of Urbanization on Malaria Infection Rate and Parasite Density in Children in the Municipality of Yopougon, Abidjan (Cote d’Ivoire). Medecine Et Sante Tropicales, 25, 69-74.
 Robert, V., Macintyre, K., Keating, J., Trape, J.F., Duchemin, J.B., Warren, M. and Beier, J.C. (2003) Malaria Transmission in Urban Sub-Saharan Africa. The American Journal of Tropical Medicine and Hygiene, 68, 169-176
 Wang, S.J., Lengeler, C., Smith, T.A., Vounatsou, P., Cisse, G., Diallo, D.A., et al. (2005) Rapid Urban Malaria Appraisal (RUMA) in Sub-Saharan Africa. Malaria Journal, 4, 40.
 Siri, J.G., Lindblade, K.A., Rosen, D.H., Onyango, B., Vulule, J., Slutsker, L. and Wilson, M.L. (2008) Quantitative Urban Classification for Malaria Epidemiology in Sub-Saharan Africa. Malaria Journal, 7, 34.
 Collins, F.H., Mendez, M.A., Rasmussen, M.O., Mehaffey, P.C., Besansky, N.J. and Finnerty, V. (1987) A Ribosomal RNA Gene Probe Differentiates Member Species of the Anopheles gambiae Complex. The American Journal of Tropical Medicine and Hygiene, 37, 37-41.
 Fanello, C., Santolamazza, F. and della Torre, A. (2002) Simultaneous Identification of Species and Molecular Forms of the Anopheles gambiae Complex by PCR-RFLP. Medical and Veterinary Entomology, 16, 461-464.
 Santolamazza, F., Mancini, E., Simard, F., Qi, Y., Tu, Z. and della Torre, A. (2008) Insertion Polymorphisms of SINE200 Retrotransposons within Speciation Islands of Anopheles gambiae Molecular Forms. Malaria Journal, 7, 163.
 Martinez-Torres, D., Chandre, F., Williamson, M.S., Darriet, F., Berge, J.B., Devonshire, A.L., et al. (1998) Molecular Characterization of Pyrethroid Knockdown Resistance (kdr) in the Major Malaria Vector Anopheles gambiae s.s. Insect Molecular Biology, 7, 179-184.
 Weill, M., Malcolm, C., Chandre, F., Mogensen, K., Berthomieu, A., Marquine, M. and Raymond, M. (2004) The Unique Mutation in ace-1 Giving High Insecticide Resistance Is Easily Detectable in Mosquito Vectors. Insect Molecular Biology, 13, 1-7.
 Dossou-yovo, J., Doannio, J., Riviere, F. and Duval, J. (1994) Rice Cultivation and Malaria Transmission in Bouake City (Cote d’Ivoire). Acta Tropica, 57, 91-94.
 Audibert, M., Brun, J., Mathonnat, J. and Henry, M. (2009) Malaria and Agricultural Production: Are There Bidirectional Effects? The Case of Coffee and Cocoa in Côte D’Ivoire. In: De Boeck Supérieur, Ed., Revue D’économie du Développement, Vol. 17, Cairn. Info, France, 107-126.
 Klinkenberg, E., McCall, P., Wilson, M.D., Amerasinghe, F.P. and Donnelly, M.J. (2008) Impact of Urban Agriculture on Malaria Vectors in Accra, Ghana. Malaria Journal, 7, 151.
 Kabula, B.I., Atta, P.K., Wilson, M.D. and Boakye, D.A. (2011) Characterization of Anopheles gambiae s.l. and Insecticide Resistance Profile Relative to Physicochemical Properties of Breeding Habitats within Accra Metropolis, Ghana. Tanzania Journal of Health Research, 13, 3.
 della Torre, A., Tu, Z. and Petrarca, V. (2005) On the Distribution and Genetic Differentiation of Anopheles gambiae s.s. Molecular Forms. Insect Biochemistry and Molecular Biology, 35, 755-769.
 Simard, F., Ayala, D., Kamdem, G.C., Pombi, M., Etouna, J., Ose, K., et al. (2009) Ecological Niche Partitioning between Anopheles gambiae Molecular Forms in Cameroon: The Ecological Side of Speciation. BMC Ecology, 9, 17.
 Young, S.J., Gunning, R.V. and Moores, G.D. (2006) Effect of Pretreatment with Piperonyl Butoxide on Pyrethroid Efficacy against Insecticide-Resistant Helicoverpa armigera (Lepidoptera: Noctuidae) and Bemisia tabaci (Sternorrhyncha: Aleyrodidae). Pest Management Science, 62, 114-119.
 Chouaibou, M., Zivanovic, G.B., Knox, T.B., Jamet, H.P. and Bonfoh, B. (2014) Synergist Bioassays: A Simple Method for Initial Metabolic Resistance Investigation of Field Anopheles gambiae s.l. Populations. Acta Tropica, 130, 108-111.
 Dadzie, S.K., Chabi, J., Asafu-Adjaye, A., Owusu-Akrofi, O., Baffoe-Wilmot, A., Malm, K., et al. (2017) Evaluation of Piperonyl Butoxide in Enhancing the Efficacy of Pyrethroid Insecticides against Resistant Anopheles gambiae s.l. in Ghana. Malaria Journal, 16, 342.