Introduction {#s1}
============

Malaria and lymphatic filariasis (LF) are the two most important vector borne parasitic diseases worldwide [@pone.0028347-WHO1], [@pone.0028347-WHO2]. In Africa these diseases are both primarily transmitted by *Anopheles* species. Nigeria has the largest burden of malaria and lymphatic filariasis in Africa -- yet very little is known about the distribution of *Anopheles* mosquitoes that act as vectors for both diseases, how the species interact, overlap or differ across the country [@pone.0028347-WHO1], [@pone.0028347-WHO2]. Knowledge of the geographical distribution of the different species, their ecological parameters, role in transmission, and susceptibility to insecticide-based interventions is critical if malaria and LF are to be controlled and eliminated in the next decade [@pone.0028347-WHO3].

The World Health Organization\'s (WHO) Position Statement on Integrated Vector Management to control malaria and lymphatic filariasis promotes integrated vector management (IVM) to improve the cost effectiveness of vector-control operations, and to strengthen the capacity of programmes, partnerships and intersectoral collaboration in their efforts to control vector-borne diseases [@pone.0028347-WHO3]. There is overlapping geographical distribution of malaria and LF in large areas of Africa, and where *Anopheles* mosquitoes transmit both the malarial and lymphatic filariasis parasites, scaling up vector-control methods such of insecticide-treated mosquito nets (ITNs) and implementing indoor residual spraying (IRS) for malaria control can effectively reduce transmission of LF [@pone.0028347-WHO3]--[@pone.0028347-Pluess1].

Malaria is caused by infection with a protozoan parasite of the genus *Plasmodium* and is endemic in 106 countries and responsible for about 225 million clinical cases and 781,000 deaths annually [@pone.0028347-WHO1]. The short-term goal of WHO\'s Global Malaria Programme is to reduce the burden of malaria until it is no longer a public-health problem, while the long-term goal is to reduce the global incidence to zero by progressively eliminating the disease in endemic countries [@pone.0028347-WHO3]. The two main components of this programme are vector control and appropriate case-management through diagnosis and treatment [@pone.0028347-WHO3], [@pone.0028347-ANVR1]. Malaria vector control involves a two pronged approach: (1) use of ITNs and/or long lasting insecticide nets (LLINs); and (2) IRS with insecticides [@pone.0028347-WHO3], [@pone.0028347-ANVR1]--[@pone.0028347-ANVR2].

International funding for malaria control has risen steeply in the past decade and has led to rapid scale-up of ITNs in Africa through support from various donors, including the Global Fund, World Bank, UNITAID, UNICEF, DfID, USAID and Canadian Red Cross, as well as funds from governments and other development agencies [@pone.0028347-Katz1], [@pone.0028347-Nigeria1]. While this increasing trend in vector control is encouraging, there is still a need to reach more households given that in recent years, an estimated 42% of households in Africa owned at least one ITN, and only 10% received IRS [@pone.0028347-WHO1], [@pone.0028347-WHO3]. By the end of 2009 a total of 19,300,000 ITNs and/or LLINs had been distributed in Nigeria and 330,000 people were protected by IRS but this was still significantly below the WHO target for 2010. Only 15% of the Nigerian population (100% at risk) owned an ITN [@pone.0028347-WHO1]. Also critical is the need for monitoring and evaluation systems to assess the impact and efficacy of these interventions. Current vector control methods involving ITNs and IRS primarily use pyrethroid insecticides although bendiocarb and DDT (Dichloro-Diphenyl-Trichloroethane) are used in some areas [@pone.0028347-Ranson1]. This widespread use of a single class of insecticide could give rise to the development of insecticide resistance in the mosquito vectors and lead to a major public health problem given the limited availability of alternative insecticides [@pone.0028347-Ranson1]--[@pone.0028347-IRAC1]. In Nigeria the entire population (154,728,895) is at risk of malaria and in 2009 there were 4,295,686 confirmed cases, 658,732 inpatient malaria cases and 7,522 malaria attributed deaths [@pone.0028347-WHO1]. The National Malaria Control Programme (NMCP) of Nigeria has various intervention policies and strategies. These include the following: distribution of ITNs or LLINs; IRS with insecticides; intermittent preventive treatment (IPT) during pregnancy; and malaria case management.

Lymphatic filariasis is caused by parasitic worms that are transmitted by mosquitoes. About 90% of infections are caused by *Wuchereria bancrofti*, while most of the remainder are caused by *Brugia malayi*. In Africa, the major vectors of *W. bancrofti* are mosquitoes of the genera *Anopheles* and *Culex* [@pone.0028347-WHO2]. The Global Program to Eliminate Lymphatic Filariasis (GPELF) was launched in 2000 with the aim of (1) interrupting transmission; and (2) reducing morbidity and preventing disability [@pone.0028347-WHO2], [@pone.0028347-WHO3]. Interrupting transmission between mosquitoes and humans is possible through mass drug administration (MDA), using once-yearly treatment with a single dose of albendazole plus either ivermectin or diethylcarbamazine (DEC) for 4--6 years [@pone.0028347-WHO2], [@pone.0028347-WHO3].

The GPELF has scaled up rapidly and by the end of 2009, 52 out of 81 endemic countries were implementing MDA, and 2.7 billion treatments had been delivered to 695 million people [@pone.0028347-WHO3]. Although significant progress has been made in the use of MDA for LF control, the role of vector control is an increasingly important issue for meeting the challenges of eliminating the disease [@pone.0028347-WHO3], [@pone.0028347-Bockarie1]. Vector control is recommended as a possible strategy for controlling LF in some countries where LF and *Loa loa* are co-endemic and where the burden of LF is heaviest e.g, Nigeria, Democratic Republic of the Congo, India, Indonesia and Bangladesh [@pone.0028347-WHO3]. Reducing mosquito populations through the use of ITNs and other insecticide treated materials, as well as IRS, may help accelerate or sustain the interruption of transmission. In the year 2000, the WHO African Region launched the Programme for Elimination of Lymphatic Filariasis, and started MDA in four countries (Ghana, Nigeria, Togo and the United Republic of Tanzania). This programme has since expanded, and in 2009 a total of 19 African countries had started MDA implementation [@pone.0028347-WHO2]. However MDA coverage in Nigeria is still less than 100% of the geographical area, which is of concern given that the country has the highest burden of the disease in Africa [@pone.0028347-WHO2].

Several maps of malaria and LF vector spatial distributions in Africa have been produced, however, most are at continental or sub-regional scale with limited specific data available for further use by national programmes and researchers aiming to better understand the epidemiology of the diseases [@pone.0028347-Lindsay1]--[@pone.0028347-MARAARMA1]. National vector spatial distribution maps and small databases have been developed for Ghana [@pone.0028347-deSouza1], Kenya [@pone.0028347-Okara1], and Mali [@pone.0028347-Sogoba1], [@pone.0028347-Sogoba2]. However, to the best of our knowledge, no country has specifically developed a national database for all *Anopheles* vectors which includes extensive historical data and specific information on locations, methodologies, insecticide resistance and parasite prevalence available in the public domain.

In Nigeria, recent studies have identified mosquitoes of the *An. gambiae* (principally *An. gambiae* s.s. and *An. arabiensis*) and *Anopheles funestus* complexes as the main vectors of malaria and LF [@pone.0028347-Sinka1], [@pone.0028347-Lenhart1]. *An. melas* is found in the coastal areas and is involved in malaria transmission [@pone.0028347-Awolola1]--[@pone.0028347-Oyewole1]. Most entomological studies have focused on small area/district based collections except for a few studies [@pone.0028347-Okwa1], [@pone.0028347-Onyabe1]--[@pone.0028347-Awolola2]. The most extensive data available are by Service [@pone.0028347-Service1] who described the distribution of 29 distinct *Anopheles* species across the country in the 1960s, and more recently by Awolola et al. [@pone.0028347-Awolola3] who developed a Malaria Entomological Profile for Nigeria.

The aim of this study was to compile a national database on the *Anopheles* vectors of malaria and LF, including related information on the location, time/period of study, species abundance, sampling and collection methods, morphological and molecular species identification methods, insecticide resistance status, including evidence of the kdr allele, and *P. falciparum* sporozoite rate and *W. bancrofti* microfilaria (mf) prevalence. This database will provide an essential baseline on *Anopheles* vectors and will be an important resource for malaria and LF vector control and elimination programmes in Nigeria.

Methods {#s2}
=======

Study area {#s2a}
----------

Nigeria is a federal constitutional republic comprising thirty-six states and its Federal Capital Territory, Abuja. The states are grouped into six geopolitical zones (northwest, northeast, north central, southwest, southsouth and southeast) [@pone.0028347-FMEN1], [@pone.0028347-FMH1]. The main latitude and longitude of Nigeria is 10° North and 8° East respectively [@pone.0028347-Iloeje1]. Nigeria is approximately 923,768 sq. km and it is located in West Africa and shares borders with Benin in the west, Chad and Cameroon in the east, and Niger in the north [@pone.0028347-FMEN1]. Its coast in the south lies on the Gulf of Guinea on the Atlantic Ocean. There are two main seasons: the wet and the dry season. Most of the rainfall in Nigeria occurs between June and September. Nigeria is broadly grouped into two zones: forests and savanna [@pone.0028347-FMH1], [@pone.0028347-NGA1]. The following vegetation types are recognized in the country: the mangrove and fresh water swamps, the rain forest, the Guinea savanna, the Sudan savanna and the Sahel in a south-north transect. Between the rain forest and the Guinea savanna is a modified vegetation transition consisting of light deciduous forest and derived savanna [@pone.0028347-FMEN1].

Data collection/collation {#s2b}
-------------------------

A systematic collation of primary empirical occurrence data for *Anopheles* mosquitoes in Nigeria, in published articles, was carried out to develop a comprehensive geo-referenced database of the distribution of *Anopheles* mosquitoes. The search was conducted using electronic searches in online bibliographic archives i.e. PubMed, SCOPUS and The Walter Reed Biosystematics Unit Culicidae Systematic Literature Database. Search terms, and combinations thereof, included Nigeria, *Anopheles*, mosquito, vectors, malaria, lymphatic filariasis. All articles with information on *Anopheles* were included. Many articles were identified through this method and many of the references were obtained from the references listed within articles, and then from the references within those articles and so on. Articles that could not be obtained online were sourced from three main libraries: College of Medicine Library, University of Ibadan; Liverpool School of Tropical Medicine Library, and the National Institute of Health Library. The references of articles obtained were also searched for additional sources of information.

For each article the following information was recorded: the locality, latitude and longitude, time/period of study, year project was initiated, species, number of specimens recorded, collection method (animal baits, human baits, indoor/outdoor resting collections, bednet trap, exit traps, human landing catches, and pyrethrum spray catches), stage of collection (adult or larval), morphological identification method (cytogenetic analysis of polytene chromosome, cross mating, morphology), molecular identification method, insecticide resistance status, insecticide tested, presence or absence of the kdr allele, *P. falciparum* sporozoite rate, *W. bancrofti* microfilaria prevalence and the reference. Only articles containing information related to the aim of this study were included in the final analysis.

The locations (i.e. mosquito collection sites) were geo-referenced using the latitude and longitude coordinates obtained by cross-checking the names with data from the GEOnet Names Server [@pone.0028347-Iloeje1], Directory of Cities and Towns in the World [@pone.0028347-Directory1] databases. Some site locations were also obtained from the research articles, as provided by the authors, while others were obtained from the WHO Malaria Entomological Profile for Nigeria [@pone.0028347-Awolola3] and other internet sources. Degree/minutes/seconds were converted into decimal degrees. It is acknowledged that there are limitations in using geographical coordinates obtained retrospectively; however, we have listed the main source of coordinates in the absence of such data. All the relevant information was entered into an Excel spreadsheet and data analysis was performed using SPSS (Version 15 for Windows, SPSS Inc., Chicago, IL). All data were mapped using the geographical information systems software ArcGIS 9.2 (ESRI, Redlands, CA). The overall species distributions and levels of insecticide resistance of the main *Anopheles* species were mapped and compared between two time periods i.e. before and after 2000. The total number of studies carried out was also compared between the six geopolitical zones and between time periods. All other data were tabulated or graphed to highlight differences in sampling, collection and identification methodologies, information related to insecticides and parasite prevalence. The mean, minimum and maximum sporozoite and microfilariae prevalence rate was calculated using SPSS software and was calculated according to species.

Results {#s3}
=======

Mosquito species and distribution maps {#s3a}
--------------------------------------

In total, 110 publications reporting on a total of 484,747 *Anopheles* mosquitoes from 632 spatially unique descriptions at 142 geo-referenced locations across Nigeria were identified between 1900 and 2010 (full database available in [Table S1](#pone.0028347.s001){ref-type="supplementary-material"}). Overall, the vector species most often reported included *An. gambiae* complex (65.2%), *An. funestus* complex (17.3%), *An. gambiae* s.s. (6.5%) *An. arabiensis* (5. 0%) and *An. funestus* s.s. (2.5%). Other species (4.5%) included *An. coustani*, *An. hancocki*, *An. leesoni*, *An. nili*, *An. melas*, *An. moucheti*, *An. rivulorum* and *An. wellcomei*.

For the time period 1900--1999, a total of 420 species-specific data points (n = 422,137 *Anopheles* mosquitoes) were recorded across 66 geo-referenced locations from 63 references, while for the time period 2000--2010, a total of 212 data points (n = 62,610 *Anopheles* mosquito) were recorded across 82 geo-referenced locations from 48 references. A higher proportion of data is recorded for the *An. gambiae* and *An. funestus* complexes in 1900--1999, while more data is recorded for *An. gambiae* s.s. and *An. funestus* s.s. in 2000--2010.

The distribution of each species and locations across the two different time periods is summarised in [Table 1](#pone-0028347-t001){ref-type="table"} and [Figure 1a--f](#pone-0028347-g001){ref-type="fig"}. Of the total 632 data points recorded, 50% (316/632) and 30.5% (193/632) were from the northwest zone and southwest zone respectively ([Figure 2](#pone-0028347-g002){ref-type="fig"}). The proportion of studies that were carried out in the north central, south south, southeast and northeast zones of the country was 10.1%, 4.3%, 3.5% and 1.6% respectively. The majority of the studies in the northwest zone were carried out before, rather than after the year 2000: 72.6% (305/420) versus 5.2% (11/212). However, the reverse was the case in the southwest zone, where the majority of the studies was carried out after, rather than before the year 2000: 71.7% (152/212) versus 9.8% (41/420).

![Distribution of *Anopheles* species a. *An. gambiae* s.l. b. *An. gambiae* s.s. c. *An. arabiensis* d. *An. funestus* complex e. *An. funestus* s.s. f. Other *Anopheles*.](pone.0028347.g001){#pone-0028347-g001}

![Number of data points across six geopolitical zones by time period.](pone.0028347.g002){#pone-0028347-g002}

10.1371/journal.pone.0028347.t001

###### The number and proportion of *Anopheles* species found in studies between 1900 and 2010.

![](pone.0028347.t001){#pone-0028347-t001-1}

                                 1900--1999    2000--2010     Total                           
  ----------------------------- ------------- ------------ ------------ ------- ------------- -------
  *An. gambiae* s.l (181)          302677         71.7        13472      21.5      316149      65.2
  *An. gambiae* s.s. (156)          8546          2.0         22760      36.4       31306       6.5
  *An. arabiensis* (122)            19529         4.6          4634       7.4       24163       5.0
  *An. funestus* complex (95)       79998         19.0         4064       6.5       84062      17.3
  *An. funestus* s.s. (21)          2382          0.6          4946       1.7       7328        2.5
  Other species (57)                9005          2.1         12734      20.3       21739       4.5
  **Total (632)**                **422,137**     100.0      **62,610**   100.0   **484,747**   100.0

The maps in [Figure 1a--f](#pone-0028347-g001){ref-type="fig"} show the spatial and temporal differences of each main species, and highlight the focus in the southern regions of country in the past decade. [Figure 1d--e](#pone-0028347-g001){ref-type="fig"} also shows that very few studies have been carried out on the *An. funestus* complex and *An. funestus* s.s. The number of data points for all species by each year and time period are summarised in [Table S2](#pone.0028347.s002){ref-type="supplementary-material"} and [S3](#pone.0028347.s003){ref-type="supplementary-material"}.

The total number of *An. gambiae* s.s. molecular forms M and S form collected throughout the study period was 4784 and 5224 ([Table 2](#pone-0028347-t002){ref-type="table"}) across 59 and 61 geo-referenced locations, respectively. For the time period 1900--1999, a total of 170 and 244 species specific data points of *An. gambiae* M form and S form were collected ([Table S4](#pone.0028347.s004){ref-type="supplementary-material"}), compared with the time period 2000--2010, when a total of 4614 and 4980 data points of *An. gambiae* M form and S form were collected, respectively ([Table S5](#pone.0028347.s005){ref-type="supplementary-material"}). [Figure 3 a--c](#pone-0028347-g003){ref-type="fig"} shows the distribution of the molecular forms and the percentage contribution of each molecular form.

![a. Distribution of *An. gambiae* s.s. molecular forms. b. Percent of M form c. Percent of S form.](pone.0028347.g003){#pone-0028347-g003}

10.1371/journal.pone.0028347.t002

###### The number and proportion of *An. gambiae* s.s. molecular forms found in studies between 1900 and 2010.

![](pone.0028347.t002){#pone-0028347-t002-2}

                     1900--1999   2000--2010    Total                    
  ----------------- ------------ ------------ ---------- ------- ------- -------
  M form (61)           170          41.1        4614     48.1    4784    47.8
  S form (59)           244          58.9        4980     51.9    5224    52.2
  **Total (120)**     **414**       100.0      **9594**   100.0   10008   100.0

Mosquito sampling and collection methods {#s3b}
----------------------------------------

Overall, the most used method of collecting *Anopheles* mosquitoes was via the adult stage. Chronologically (based on the year each project was initiated) adult collection represented 87.5% (7/8), 100% (3/3), 67.1% (47/70), 96.1% (221/230), 33.9% (41/121) and 36.2% (54/149) of collections made between 1900--1920, 1921--1940, 1941--1960, 1961--1980, 1981--2000 and 2001--2010 respectively. Larval collections represented 12.5% (1/8), 0% (0/3), 4.3% (3/70), 0.9% (2/230), 48.8% (59/121) and 34.2% (51/149) for the same time periods.

1900--1999 {#s3c}
----------

The main methods of collecting *Anopheles* mosquitoes included adult collections (72.9%; 306/420) and larval collections (15.7%; 65/420), while a small proportion used larval and adult collections combined (8.6%; 36/420) or did not specify the specific methods (2.9%; 12/420) ([Table 3](#pone-0028347-t003){ref-type="table"}). In general, the frequency of methods used was similar between the main species ([Figure 4a](#pone-0028347-g004){ref-type="fig"}).

![Summary of *Anopheles* species collection methods by time period a 1900--1999 b. 2000--2010.](pone.0028347.g004){#pone-0028347-g004}

10.1371/journal.pone.0028347.t003

###### Summary of *Anopheles* species collection methods in studies between 1900 and 1999.

![](pone.0028347.t003){#pone-0028347-t003-3}

  Species                       Adult collection   Larval collection   Adult+larval collection   Unspecified    Total
  ---------------------------- ------------------ ------------------- ------------------------- ------------- ---------
  *An. gambiae* s.l                   135                  8                     10                   5          158
  *An. gambiae* s.s.                   37                 28                     11                   2          78
  *An. arabiensis*                     28                 30                     15                   2          75
  *An. funestus* complex               77                  0                      0                   3          80
  *An. funestus* s.s.                  5                   0                      0                   0           5
  *An. leesoni*                        1                   0                      0                   0           1
  *An. melas*                          1                   0                      0                   0           1
  *An. moucheti*                       2                   0                      0                   0           2
  *An. nili*                           17                  0                      0                   0          17
  *An. wellcomei*                      3                   0                      0                   0           3
  **Total**(n = data points)        **304**             **65**                 **36**              **12**      **420**

Adult collections were made primarily from pyrethrum spray catches (PSC) (17.6%), human landing catches (HLC) (12.9%), indoor resting (IR) collections (8.6%), exit trap collections (ETC) (6.9%) and animal baits (1.7%). Some collections were made using a combination of methods e.g. net catch/human bait, pyrethrum spray catch/exit trap, human landing catch/pyrethrum spray catch.

Larval collections were made from potential breeding sites, which included gutters, abandoned road sides, standing waters, vehicle tracks, tires, shallow wells, ponds, swamps, drains, rivers, small streams, irrigation ditches, hoof prints, domestic containers and empty cans ([Table S1](#pone.0028347.s001){ref-type="supplementary-material"}). Some collections were made using a combination of larval and adult collection methods, e.g. pyrethrum spray catches/larval, indoor resting collections/larval, human landing catches/exit trap and indoor resting/human landing catch/pyrethrum spray catches/larval collections.

2000--2010 {#s3d}
----------

The main methods of collecting *Anopheles* mosquitoes included adult collections (45.3%; 96/212), larval collections (32.5%; 69/212), and larval and adult collections combined (21.2%; 45/212), while a small proportion used larval and pupal collections (0.5%; 1/212), or did not specify the methods (0.5%; 1/212), as shown in [Table 4](#pone-0028347-t004){ref-type="table"}. Compared with data collection method in the 1900--1999 time period, overall there was a decrease in adult collections and an increase in the number of larval collections ([Figure 4b](#pone-0028347-g004){ref-type="fig"}).

10.1371/journal.pone.0028347.t004

###### Summary of *Anopheles* species collection methods in studies between 2000 and 2010.

![](pone.0028347.t004){#pone-0028347-t004-4}

  Species                       Adult collection   Larval collection   Larval+pupal collections   Adult+larval collection   Unspecified    Total
  ---------------------------- ------------------ ------------------- -------------------------- ------------------------- ------------- ---------
  *An. gambiae* s.l                    15                  6                      1                          1                   0          23
  *An. gambiae* s.s.                   25                 30                      0                         22                   1          78
  *An. arabiensis*                     11                 29                      0                          7                   0          47
  *An. funestus* complex               11                  3                      0                          1                   0          15
  *An. funestus* s.s.                  10                  0                      0                          6                   0          16
  *An. leesoni*                        1                   0                      0                          3                   0           4
  *An. melas*                          3                   0                      0                          2                   0           5
  *An. moucheti*                       9                   0                      0                          3                   0          12
  *An. nili*                           6                   0                      0                          0                   0           6
  *An. rivulorum*                      3                   0                      0                          0                   0           3
  *An. coustani*                       1                   1                      0                          0                   0           2
  *An. hancocki*                       1                   0                      0                          0                   0           1
  **Total**(n = data points)         **96**             **69**                  **1**                     **45**               **1**      **212**

Adult collections were made primarily from human landing catches (18.4%), indoor resting collections (9.0%), CDC light traps (5.2%) and pyrethrum spray catches (5.2%). Some methods were also used in combination: human landing collections/pyrethrum spray catches/indoor resting collections, indoor resting/human landing catch.

Similar to the 1900--1999 time period, larval collections were made from a variety of potential breeding sites ([Table S1](#pone.0028347.s001){ref-type="supplementary-material"}), and some studies used a combination of adult and larval methods, e.g. human landing catch/larval/pyrethrum spray catches, larval/indoor resting, human landing catch/indoor resting/larval, larval/outdoor baited scarecrow, and larval/pupal collections.

Species identification methods {#s3e}
------------------------------

Chronologically (based on the year each project was initiated), the method of species identification was not specified in studies initiated between 1900--1960 (n = 81). For 1961--1980 period, cytogenetic analysis of polytene chromosome represented 28.7% (66/230) of identification method employed. Morphological identification and use of PCR was the most used method of identification in studies carried out between 1981--2000 and 2001--2010 representing 60.3% (73/281) and 38.3% (57/149) respectively.

1900--1999 {#s3f}
----------

The methods of species identification were based on cytogenetic analysis of polytene chromosome (16.0%; 67/420), morphology and PCR (14.0%; 59/420), PCR (6.0%; 25/420), cross mating techniques (1.9%; 8/420) and morphology alone (0.7%; 3/420) ([Figure 5a](#pone-0028347-g005){ref-type="fig"}). The method of species identification was not specified in studies initiated between 1900 and 1963 representing 61.4% (258/420) of sites. The earliest study that identified species based on cytogenetic analysis of chromosomes was in a study initiated in 1969 [@pone.0028347-Service2]; cross mating techniques was first used in a study initiated in 1963 [@pone.0028347-Service3], morphology was first specified in a study initiated in 1965 [@pone.0028347-Pant1] while use of PCR was first mentioned in a study initiated in 1997 [@pone.0028347-Onyabe2], [@pone.0028347-Onyabe3].

![Summary of *Anopheles* species identification methods by time period a 1900--1999 b. 2000--2010.](pone.0028347.g005){#pone-0028347-g005}

The morphological keys used included those described in [@pone.0028347-DeMeillon1]--[@pone.0028347-Gillies2]. The molecular methods - polymerase chain reaction (PCR) assays for the vector species complex and molecular forms - were those of [@pone.0028347-Scott1] for *An. gambiae* complex; [@pone.0028347-Favia1] for distinguishing the molecular forms of *An. gambiae* s.s. and [@pone.0028347-Koekemoer1] for *An. funestus* complex.

For *An. gambiae* s.l. and the *An. funestus* complex, the species identification method was not specified for 93.0% (147/158) and 95.0% (76/80) of sites, respectively. The PCR method was not used at any site for the *An. gambiae* complex or other species, and only once for the *An. funestus* complex. In contrast, for *An. gambiae* s.s. and *An. arabiensis* the PCR method of [@pone.0028347-Scott1] was used for species identification and the method of [@pone.0028347-Favia1] for *An. gambiae* s.s molecular form identification. For the *An. gambiae* s.s. data, a total of 10 (n = 78) collection sites had information on both the M form and the S form.

2000--2010 {#s3g}
----------

The methods of species identification were based on morphology and PCR (47.6%; 101/212), morphology alone (21.7%; 46/212), and PCR alone (9.9%; 21/212) ([Figure 5b](#pone-0028347-g005){ref-type="fig"}). The method of species identification was not specified in 20.8% (44/212) of sites. The morphological keys used included [@pone.0028347-Gillies1]--[@pone.0028347-Gillies2], [@pone.0028347-DeMeillon2]--[@pone.0028347-Brunhes1]. The molecular PCR methods for the identification of vector species complex and molecular forms were those of [@pone.0028347-Scott1]--[@pone.0028347-Koekemoer1], [@pone.0028347-Fanello1]--[@pone.0028347-Kengne1] ([Table S1](#pone.0028347.s001){ref-type="supplementary-material"}).

For *An. gambiae* s.l. and *An. funestus* complex, morphological methods were used at 69.6% (16/23) and 73.3% (11/15) of sites, and not specified for 26.1 (6/23) and 20% (3/15) of sites respectively ([Figure 5b](#pone-0028347-g005){ref-type="fig"}). For *An. gambiae* s.s., *An. arabiensis* and *An. funestus* s.s., a combination of PCR and morphological identification was primarily used, while for other species morphological methods were primarily used (*An. coustani*, *An. leesoni*, *An. melas*, *An. moucheti*, *An. nili*, *An. hancocki* and *An. rivulorum*).

Insecticide resistance {#s3h}
----------------------

Insecticide resistance to pyrethroids and organochlorines was recorded in *An. gambiae* s.l., *An. gambiae* s.s., *An. arabiensis*, *An. funestus* complex and *An. funestus* s.s. across 45 geo-referenced locations ([Table 5](#pone-0028347-t005){ref-type="table"}). The specific locations where resistance was found for each species, is available in [Table S6](#pone.0028347.s006){ref-type="supplementary-material"} and highlighted in [Figures 6a--g](#pone-0028347-g006){ref-type="fig"}. The species most tested for insecticide resistance were *An. gambiae* s.s. 47.1% (33/70), *An. arabiensis* 25.7% (18/70) and *An. gambiae* s.l. 18.6% (13/70).

![Distribution of insecticide resistance a. *An. gambiae* s.l. b. *An. gambiae* s.s. c. *An. arabiensis* d. *An. funestus* complex e. *An. funestus* s.s.](pone.0028347.g006){#pone-0028347-g006}

10.1371/journal.pone.0028347.t005

###### Insecticides resistance recorded in *Anopheles* species.

![](pone.0028347.t005){#pone-0028347-t005-5}

  Insecticide                         *An. gambiae* s.l   *An. gambiae* s.s.   *An. arabiensis*   *An. funestus* complex   *An. funestus* s.s.   Total
  ---------------------------------- ------------------- -------------------- ------------------ ------------------------ --------------------- --------
  BHC                                         3                   0                   0                     1                       0              4
  DDT                                         1                   12                  6                     2                       1              22
  Deltamethrin, Permethrin and DDT            0                   4                   0                     0                       0              4
  Dieldrin                                    4                   0                   0                     1                       1              6
  Dieldrin and BHC                            3                   0                   0                     0                       0              3
  Dieldrin, Lindane, DDT                      1                   0                   0                     0                       0              1
  Permethrin                                  1                   12                  10                    0                       0              23
  Permethrin and DDT                          0                   5                   2                     0                       0              7
  **Total** (n = data points)              **13**               **33**              **18**                **3**                   **2**          **70**

Data were recorded for individual insecticides as well as a combination of insecticides. These were predominantly permethrin (32.9%; 23/70), DDT (31.4%; 22/70), and permethrin and DDT (10.0%; 70), as well as dieldrin (8.6%; 6/70), deltamethrin and permethrin and DDT (5.7%; 4/70), lindane (BHC) (5.7%; 4/70), dieldrin and lindane (BHC) (4.3%; 3/70), and dieldrin and lindane (BHC) and DDT (1.4%; 1/70) ([Table 5](#pone-0028347-t005){ref-type="table"}). For *An. gambiae* s.s., resistance to DDT and permethrin was mainly reported, and the kdr allele was tested 26 times and found to be present in ten locations predominantly in the southwestern region of the country ([Figure 7](#pone-0028347-g007){ref-type="fig"}). For *An. arabiensis*, resistance to permethrin, and for *An. gambiae* s.l resistance to lindane (BHC), and dieldrin were mainly reported ([Table S1](#pone.0028347.s001){ref-type="supplementary-material"}).

![Distribution of kdr alleles.](pone.0028347.g007){#pone-0028347-g007}

*Plasmodium falciparum* sporozoite rate {#s3i}
---------------------------------------

In total, *P. falciparum* sporozoite rates were recorded at 106 data points and geo-referenced locations from 20 references, and most of these were carried out in the last decade. Overall, the sporozoite rate varied by *Anopheles* species ([Table 6](#pone-0028347-t006){ref-type="table"}) and location, with most data recorded in the southwestern region of the country ([Figure 8a](#pone-0028347-g008){ref-type="fig"}). Approximately 35% (30/86) of data points were recorded in two of the studies [@pone.0028347-Okwa2], [@pone.0028347-Okwa1].

![Distribution of a. sporozoite rates and b. LF microfilariae.](pone.0028347.g008){#pone-0028347-g008}

10.1371/journal.pone.0028347.t006

###### Mean sporozoite rate recorded for the various *Anopheles* species (data point = 106).

![](pone.0028347.t006){#pone-0028347-t006-6}

  Species(n = data points)       Mean sporozoite rate (%)   Min sporozoite rate (%)   Max sporozoite rate (%)
  ----------------------------- -------------------------- ------------------------- -------------------------
  *An. gambiae* s.l (32)                   2.81                       0.4                       10
  *An. gambiae* s.s (20)                   19.4                        0                        91
  *An. arabiensis* (11)                    2.34                        0                        6.1
  *An. funestus* complex (10)              1.99                        0                       5.35
  *An. funestus* s.s (13)                 11.79                        0                        50
  *An. leesoni* (1)                         0                          0                         0
  *An. melas* (5)                          2.57                       0.4                       6.5
  *An. moucheti* (7)                       5.22                        0                       21.9
  *An. nili* (2)                            0                          0                         0
  *An. rivulorum* (3)                       0                          0                         0
  *An hancocki* (1)                         0                          0                         0
  *An. wellcomei* (1)                       0                          0                         0

The mean sporozoite rate in *An. gambiae* s.l. was 4.3% and ranged from 0.4% (Sokoto [@pone.0028347-Dodge1]) to 10% (Lagos [@pone.0028347-MuirheadThomson1]). For *An. gambiae* s.s. the mean rate was 19.4%, and ranged from 0 (Bama [@pone.0028347-Noutcha1]) to 91% (Alimosho [@pone.0028347-Okwa2]). For *An. arabiensis* the mean rate was 2.3% and ranged from 0 (Amuwo-Odofin [@pone.0028347-Okwa2] and Lemu surburb of Lagos [@pone.0028347-Awolola1]) to 6.1% (Kaduna area [@pone.0028347-Service2]). For the *An. funestus* complex, the mean rate was 2.1% and ranged from 0 (Sokoto [@pone.0028347-Dodge1]) to 5.35% (Kaduna area [@pone.0028347-Hanney1]). For *An. funestus* s.s. the mean rate was 11.8% and ranged from 0 (Agege) to 50% (Amuwo-Odofin) [@pone.0028347-Okwa2]. For *An. melas* was 2.57% and ranged from 0.4 (Bonny [@pone.0028347-Okwa1]) to 6.5% (Moba [@pone.0028347-Oyewole1]). For *An. moucheti* the mean rate was 5.2%, and ranged from 0 (Remo north and Odogbolu LGA [@pone.0028347-Oyewole2]) to 21.9% (Mushin [@pone.0028347-Okwa2]). For other species, no sporozoites were recorded for *An. leesoni*, *An. nili*, *An. rivulorum*, *An. hancocki* and *An. wellcomei* ([Table 6](#pone-0028347-t006){ref-type="table"}).

Lymphatic filariasis {#s3j}
--------------------

The prevalence of LF microfilariae (mf) in *Anopheles* species was recorded for 21 data points, across ten geo-referenced locations and eight references ([Table 7](#pone-0028347-t007){ref-type="table"} and [Figure 8b](#pone-0028347-g008){ref-type="fig"}). The species found to be infected with mf included *An. gambiae* s.l., *An. gambiae* s.s., *An. arabiensis* and *An. funestus* complex. The mean mf prevalence in *An. gambiae* s.l was 5.64% and ranged from 0 (Seri [@pone.0028347-Richards1] and Ibadan [@pone.0028347-Ogunba1]) to 21.7% (Igwun Basin [@pone.0028347-Udonsi1]). For *An. gambiae* s.s. and *An. arabiensis* mf was found only in one site and was 9.2% (Jos Plateau) and 11.1% (Jos Plateau) respectively [@pone.0028347-Lenhart1]. For *An. funestus* complex the mean mf prevalence was 3.15% and ranged from 0% (Ibadan [@pone.0028347-Ogunba1]) to 6.6% (Amassoma [@pone.0028347-Agi1]). The data from two articles [@pone.0028347-Richards1], [@pone.0028347-Richards2] had to be excluded from the analysis because the general data on mf prevalence were combined from various sites and/or species. The infection rate of mosquitoes containing *W. bancrofti* was reported to be 1.0% for *An. funestus* complex and 3.55% for *An. gambiae* s.l but the particular site was not specified [@pone.0028347-Richards1]. Also [@pone.0028347-Agi1] gave the prevalence of mf to be 4.9% and 4.7% in the untreated and ivermectin treated village respectively, but the percentage contribution of each species was not given. These two papers both presented results from studies carried out in Plateau and Nassarawa states ([Table S1](#pone.0028347.s001){ref-type="supplementary-material"}).

10.1371/journal.pone.0028347.t007

###### Prevalence of *W. bancrofti* microfilariae in various *Anopheles* species.

![](pone.0028347.t007){#pone-0028347-t007-7}

  Species (data points)         Mean mf prevalence (%)   Min mf prevalence (%)   Max mf prevalence (%)
  ---------------------------- ------------------------ ----------------------- -----------------------
  *An. gambiae* s.l (10)                 5.64                      0                     21.7
  *An. gambiae* s.s (1)                  9.2                      9.2                     9.2
  *An. arabiensis* (1)                   11.1                    11.1                    11.1
  *An. funestus* complex (9)             3.15                      0                      6.6

Discussion {#s4}
==========

This collation of data from 110 publications has produced one of the most comprehensive geo-referenced databases on *Anopheles* vectors of malaria and filariasis available for any country in the world. Importantly, it has been carried out for Nigeria, which has the largest burden of both diseases in Africa, and is a priority country for international donors (the Global Fund, World Bank, UNITAID, UNICEF, DFID, USAID and Canadian Red Cross) supporting large-scale control programmes [@pone.0028347-Nigeria1]. The data span more than 100 years and show changes in the geographical focus of research sites, field and laboratory methodologies, and they highlight critical gaps in our knowledge on *Anopheles* species distributions, transmission dynamics, species and parasite interactions, level and impact of insecticide resistance.

The results reported by the various articles show that the most abundant vectors are the *An. gambiae* complex, *An. funestus* complex, *An. gambiae* s.s., *An. arabiensis* and *An funestus* s.s., with these vectors existing in sympatry in most locations. The earliest study found was the expedition by the Liverpool School of Tropical Medicine [@pone.0028347-Annett1]. In the earlier reports *An. gambiae* s.l. was referred to as "*A. costalis*" [@pone.0028347-Johnson1], [@pone.0028347-Taylor1]. In the early days (and early reports) the term "*An. gambiae*" encompassed the species complex [@pone.0028347-Service3], [@pone.0028347-Dodge1], [@pone.0028347-Hanney1], [@pone.0028347-Archibald1]--[@pone.0028347-Fontaine1]. In subsequent reports *An. gambiae s.s.* was referred to as "species A" [@pone.0028347-Molineaux1], [@pone.0028347-Molineaux2] and *An. arabiensis* was referred to as "species B" [@pone.0028347-Molineaux1], [@pone.0028347-Molineaux2].

The total number of mosquitoes collected is under represented as the numbers of mosquitoes collected were not available in some studies. As the species abundance was not specified in some studies, the totals for each species was carried out independently and did not interfere with any of the calculations reported here. In general, there is more information provided on the *An. gambiae* s.s. and *An. funestus* s.s. collected in studies initiated in the 2000s than in those from earlier years. Also, more information was available on *Anopheles gambiae* s.l., which represented over 60% of the overall species composition compared to about 19.8% for *An. funestus* (17.5% of the studies were on *An. funestus* complex while 2.3% were on *An. funestus* s.s.). Clearly, research on the composition and distribution of the *An. funestus* group is important, as it is an important vector contributing to the transmission of both LF and malaria in Nigeria [@pone.0028347-Okwa1], [@pone.0028347-Richards1], [@pone.0028347-Molineaux2], [@pone.0028347-Awolola4]. One main reason for the lack of research on *An. funestus* may be that this vector is refractory to colonization, and very difficult to find in the field, and its larvae are very difficult to find at low densities due to their tendency to stay submerged for long periods [@pone.0028347-Gillies1]. However, new techniques for egg laying and colonising are being developed, which will facilitate research opportunities in the future [@pone.0028347-Morgan1].

Of particular interest is how the molecular forms of *An. gambiae* s.s. varied across the country. In most studies the *An. gambiae* M and S molecular forms were found in sympatry, although the ratio varied from location to location [@pone.0028347-Onyabe3], [@pone.0028347-Awolola2], [@pone.0028347-Kristan1]--[@pone.0028347-Oduola1]. At the 120 sites where M and/or S forms were reported, pure populations of the S form were recorded at 7 sites, all located in the northern region (Guinea and Sudan savanna ecological zones) [@pone.0028347-Onyabe3], [@pone.0028347-Awolola2], [@pone.0028347-Lehmann1], [@pone.0028347-Lehmann2]. Pure populations of the M form were reported at 17 sites, 14 of them located in the southern region (mangrove, forest, and transitional ecological zones) [@pone.0028347-Oyewole1], [@pone.0028347-Onyabe3], [@pone.0028347-Oduola1], [@pone.0028347-Djouaka1], [@pone.0028347-Djouaka2]. The other 3 sites with pure populations of the M form were from the north central (Guinea savanna) part of the country [@pone.0028347-Oduola1]. Geographical and ecological differences in the chromosomal and molecular forms of *An. gambiae* s.s. and their potential role in disease transmission have been examined in regions of West Africa [@pone.0028347-Bayoh1], [@pone.0028347-deSouza1], [@pone.0028347-KellyHope2] with the *An. gambiae* S form broadly associated with malaria distributions, and the *An. gambiae* M form with LF distributions. However, the role of these two sibling species in the transmission of malaria and/or LF and how they interact is largely unknown.

The series of mosquito vector distribution maps produced in this study highlight that studies in the 20th Century were carried out predominately in the northern region of the country, while studies carried out more recently have focussed mainly in the southern region. The reason for this shift may be related to limited infrastructure and an overall lack of trained staff, vector ecologists and medical entomologists, which are common problems across Africa [@pone.0028347-ANVR1], [@pone.0028347-ANVR2], [@pone.0028347-RBM1], [@pone.0028347-Greenwood1]. It may also be due to the limited accessibility of remote locations, as the majority of studies were carried out in proximity to cities and towns in two main zones, northwest and southwest, leaving large gaps in the eastern regions of the country. Furthermore, very few studies were carried out in the same locations in different years or seasons: this lack of systematic sampling/collections restricts analysis and comparisons on multiple fronts, for instance of changes in species composition and abundance over time. The data collated here can serve as a baseline for future work that includes such same-site collections and allows for comparisons within and between populations. Currently, this is particularly important given the widespread vector control interventions taking place across the country [@pone.0028347-Nigeria1], [@pone.0028347-FMH1].

There was great variability in the reporting of the mosquito sampling and collection methods used over time. For example, between 1900 and 1999, PSC, HLC, IR and ETC were the methods most frequently used either singly or in combination for adult mosquito collections, whereas in studies carried out between 2000 and 2010, the CDC light trap was mainly used; there was no record of ETC use during this period. Surprisingly, the expensive and labour intensive method of the HLC [@pone.0028347-Sikulu1] was used consistently throughout the study period from 1900 to 2010, which is probably due to a lack of a suitable substitute for this important metric. However, recent tent traps have been tested and calibrated to the HLC and their use as an alternative tool is being trialled in different locations across Africa [@pone.0028347-Sikulu1], [@pone.0028347-Govella1]. The variability in the use of different sampling and collection methods calls for simpler and more standardized methods; the value of this has been previously emphasized [@pone.0028347-Service8]--[@pone.0028347-KellyHope3].

Similarly, there were changes in species identification methods over time. Earlier studies used more cross mating techniques [@pone.0028347-Service3], [@pone.0028347-Ramsdale1], morphological [@pone.0028347-Onyabe1], [@pone.0028347-Pant1], [@pone.0028347-Anyanwu1], [@pone.0028347-Michel1] and cytogenetic methods[@pone.0028347-Service2], [@pone.0028347-Molineaux2], [@pone.0028347-Shidrawi1]--[@pone.0028347-Pombi1] for species identification, however, by the 2000s cross mating techniques and cytogenetic methods were nonexistent. Most of the missing data were from earlier studies when there were few options for species identification, For example, most of the missing data were from articles published between 1900 and 1963, when cross mating and morphological methods were predominately used. However it was interesting to note that taxonomic keys [@pone.0028347-DeMeillon2]--[@pone.0028347-Chandler1] were used for species identification in a study carried out between 2002 and 2003 [@pone.0028347-Okogun1]. In recent decades the emphasis has shifted to molecular techniques including PCR assays that target specific regions of repeat gene families, such as the ribosomal RNA (rRNA) gene family. A diagnostic assay for the identification of the *An. gambiae* complex based on IGS and ITS sequence differences was developed and applied routinely [@pone.0028347-Scott1]. Similar techniques were also developed for *An. funestus* s.l. [@pone.0028347-Koekemoer1], and *An. moucheti* s.l. [@pone.0028347-Kengne1], [@pone.0028347-Kengne2], along with more advancements allowing the simultaneous identification the *An. gambiae* complex species and *An gambiae* s.s. M and S molecular forms [@pone.0028347-Scott1], [@pone.0028347-Favia1], [@pone.0028347-Fanello1].

In contrast, there were no major changes in the methods used for the detection of insecticide resistance in Nigeria. In most of the studies, bioassays and molecular assays to detect resistance alleles (specifically the kdr mutation) were used. Bioassays are the best indicators of the presence of resistance in a field population and are widely used across Africa [@pone.0028347-ANVR2], [@pone.0028347-KellyHope1], [@pone.0028347-WHO4], [@pone.0028347-Penilla1]. However, molecular and biochemical techniques are important to verify bioassay results of resistance in wild populations and gain an understanding of the underlying mechanisms of resistance [@pone.0028347-KellyHope1], [@pone.0028347-IRAC1], [@pone.0028347-Brooke1]. In Nigeria only one study used bioassay and microarray analysis [@pone.0028347-Djouaka2], and another used bioassays, molecular assays to detect resistance alleles, and biochemical analysis and microarray analysis to characterize pyrethroid resistance mechanism in colony breed *An. gambiae* s.s. [@pone.0028347-Awolola8]. The reason for the heavy reliance on bioassays and molecular assays for resistance monitoring is probably because the biochemical assays and microarray analysis require qualified personnel and specialized and costly equipment [@pone.0028347-Ranson1].

Insecticide resistance in the *An. gambiae* complex, *An. gambiae* s.s., *An.arabiensis*, *An. funestus* complex and *An. funestus* s.s. to various insecticides (DDT, permethrin, dieldrin, deltamethrin and lindane (BHC)) was widespread across the southwestern region of the country. Resistance was most frequently reported in *An. gambiae* s.s. to permethrin and DDT, which is of major concern given that this vector was found to have among the highest malaria sporozoite and LF microfilaria rates in Nigeria, permethrin is currently the only class of insecticide used for ITNs, and permethrin and DDT are widely used for IRS [@pone.0028347-Ranson1]. Vector control can be hampered by the occurrence of insecticide resistance and as such monitoring of insecticide resistance regularly across a wide geographical area is critical [@pone.0028347-ANVR2], [@pone.0028347-KellyHope1], [@pone.0028347-IRAC1], [@pone.0028347-KellyHope4]--[@pone.0028347-IRAC3]. The findings also highlight the importance of a consensus standardization of the number and location of surveillance sites and the frequency with which resistance monitoring should occur [@pone.0028347-KellyHope1]. To date very little has been done in the northwestern region of Nigeria, especially, in the area where The Garki Project was carried out in the 1960--70s [@pone.0028347-Molineaux2]. This is disappointing given that it was one of the largest and most comprehensive control programmes to be carried out in Africa at the time, and vector and resistance monitoring may have provided some invaluable insights into the long term effects of large scale vector control programmes.

In general, there were few studies that focussed on malaria and LF parasites. The sporozoite rates recorded emphasize the role of *An. gambiae* s.s. and *An. funestus* s.s. as very efficient malaria vectors. In areas of the southern part of the country when the two species were studied in sympatry, *An. gambiae* s.s. had higher sporozoite rates than *An. funestus* s.s. At Bungudu-Gusau in the northern part of the country, however, *An. funestus* s.s. had a sporozoite rate of 2.3% compared to 0.46% recorded in *An. gambiae* s.s. [@pone.0028347-Okwa1]. The sporozoite rates recorded for *An. arabiensis* were also lower than those for *An. gambiae* s.s. in areas where they were studied at the same time [@pone.0028347-Okwa1], [@pone.0028347-Oyewole2], [@pone.0028347-Oyewole3]. These findings call for a widespread systematic sampling across the country. For LF, the results show that *An. gambiae* s.l., *An. gambiae* s.s., *An. arabiensis* and *An. funestus* complex are important vectors. Of the eight studies [@pone.0028347-Lenhart1], [@pone.0028347-Richards1]--[@pone.0028347-Agi1], [@pone.0028347-Annett1], [@pone.0028347-Taylor1], [@pone.0028347-Anosike1] carried out to determine the prevalence of LF only one study [@pone.0028347-Taylor1] reported co-infection of *Plasmodium* parasites and LF in the mosquitoes that were trapped. Annett et al [@pone.0028347-Annett1] confirmed the role played by *Anopheles* mosquitoes as malaria vectors but it was not stated whether co-infection of malaria parasites and LF was studied. The few records found make it difficult for comparisons to be made between studies, locations and species. This becomes alarming as Nigeria bears the greatest potential burden of LF in Africa, with 80 million people (19% of the total population) at risk [@pone.0028347-Manguin1].

As this study has restated that LF and malaria are transmitted by the same vectors in Nigeria, both diseases can be jointly controlled since the two diseases share a large proportion of their target population, and the national programmes have similar goals and strategies [@pone.0028347-WHO3], [@pone.0028347-Njenga1]. The result from all the articles used in the database [@pone.0028347-Lenhart1]--[@pone.0028347-Awolola1],[@pone.0028347-Awolola3],[@pone.0028347-Service3],[@pone.0028347-Pant1],[@pone.0028347-Kengne1]--[@pone.0028347-Awolola4],[@pone.0028347-Kristan1]--[@pone.0028347-Djouaka2],[@pone.0028347-Ramsdale1]--[@pone.0028347-Kengne2],[@pone.0028347-Oyewole3]--[@pone.0028347-Anosike1],[@pone.0028347-WHO5]--[@pone.0028347-Manga1] emphasizes the need for a detailed understanding of the distribution, species composition, behaviour and insecticide susceptibility levels of local vectors in order to successfully control the diseases. It reveals the areas where a dearth of data exists and emphasizes the importance of collecting data systematically, so that the impact of the interventions can be measured. These considerations are important in the light of the goal of the WHO\'s Global Malaria Programme and the WHO\'s Global Programme to Eliminate Lymphatic Filariasis to eliminate the diseases as a public health problem [@pone.0028347-WHO3]. The findings in this research shows that better linkages and partnerships between entomologists, parasitologists, national control programmes, international donors and other stakeholders in the country is needed, in order to carry out meaningful research and control on the vectors and their disease transmission across the country.

Supporting Information {#s5}
======================

###### 

Nigeria *Anopheles* vector database.

(XLSX)

###### 

Click here for additional data file.

###### 

The number of data points by each species and each year from 1900--1999.

(XLSX)

###### 

Click here for additional data file.

###### 

The number of data points by each species and each year from 2000--2010.

(XLSX)

###### 

Click here for additional data file.

###### 

Number of specimen of *Anopeles gambiae* s.s. M form and S form from 1900--1999.

(XLSX)

###### 

Click here for additional data file.

###### 

Number of specimen of *Anopeles gambiae* s.s. M form and S form from 2000--2010.

(XLSX)

###### 

Click here for additional data file.

###### 

Localities and species recording resistance and the frequency.

(XLSX)

###### 

Click here for additional data file.

**Competing Interests:**The authors have declared that no competing interests exist.

**Funding:**The authors have no support or funding to report.

[^1]: Conceived and designed the experiments: PNO FEM LK-H. Performed the experiments: PNO LK-H. Analyzed the data: PNO FEM OGA MJB LK-H. Contributed reagents/materials/analysis tools: LK-H. Wrote the paper: PNO LK-H. Gave conceptual advice and provided editorial feedback on final manuscript: FEM OGA MJB.
