Malaria transmission dynamics surrounding the first nationwide long-lasting insecticidal net distribution in Papua New Guinea
© Reimer et al. 2016
Received: 24 July 2015
Accepted: 22 December 2015
Published: 12 January 2016
The major malaria vectors of Papua New Guinea exhibit heterogeneities in distribution, biting behaviour and malaria infection levels. Long-lasting, insecticide-treated nets (LLINs), distributed as part of the National Malaria Control Programme, are the primary intervention targeting malaria transmission. This study evaluated the impact of LLINs on anopheline density, species composition, feeding behaviour, and malaria transmission.
Mosquitoes were collected by human landing catch in 11 villages from East Sepik Province and Madang Province. Mosquitoes were collected for 3 years (1 year before distribution and 2 years after), and assayed to determine mosquito species and Plasmodium spp. infection prevalence. The influence of weather conditions and the presence of people and animals on biting density was determined. Determinants of biting density and sporozoite prevalence were analysed by generalized estimating equations (GEE).
Mosquito biting rates and entomological inoculation rates decreased significantly after the distribution. Plasmodium falciparum and P. vivax sporozoite prevalence decreased in year 2, but increased in year 3, suggesting the likelihood of resurgence in transmission if low biting rates are not maintained. An earlier shift in the median biting time of Anopheles punctulatus and An. farauti s.s. was observed. However, this was not accompanied by an increase in the proportion of infective bites occurring before 2200 hours. A change in species composition was observed, which resulted in dominance of An. punctulatus in Dreikikir region, but a decrease in An. punctulatus in the Madang region. When controlling for village and study year, An. farauti s.s., An. koliensis and An. punctulatus were equally likely to carry P. vivax sporozoites. However, An. punctulatus was significantly more likely than An. farauti s.s. (OR 0.14; p = 0.007) or An. koliensis (OR 0.27; p < 0.001) to carry P. falciparum sporozoites.
LLINs had a significant impact on malaria transmission, despite exophagic and crepuscular feeding behaviours of dominant vectors. Changes in species composition and feeding behaviour were observed, but their epidemiological significance will depend on their durability over time.
Malaria transmission in Papua New Guinea (PNG) is highly variable across environmentally diverse zones, ranging from intense perennial transmission in the northern coastal lowlands to seasonal moderate transmission in the southern coast and unstable transmission at higher altitudes . A recent survey reported weighted malaria parasite prevalence of 12 % nationwide, with substantial heterogeneity, ranging from 0 to 49.7 % . Four malaria species are endemic to PNG (Plasmodium falciparum, P. ivax, P. ovale, and P. malariae) with the majority of infections caused by P. falciparum and P. vivax [3, 4]. The major malaria and filariasis vectors in PNG are members of the Anopheles punctulatus group. This group comprises 13 species, each exhibiting different degrees of exophily and anthropophily and different habitat preferences [5–8]. The five major malaria vectors in this group, due to their widespread distribution and high abundance, include An. punctulatus, An. farauti s.s., An. koliensis, An. hinesorum, and An. farauti 4 . In addition to the Anopheles punctulatus group, An. bancroftii, An. longirostris, An. karwari, and An. subpictus have been incriminated as malaria vectors.
Prior to the development of more sensitive molecular diagnostics in the 1990s [9–11], identification was restricted to the Punctulatus clade (now known to include An. punctulatus, An. farauti 4 and An. sp. nr punctulatus), Farauti clade (now known to include An. farauti s.s., An. hinesorum, An. torresiensis, An. farauti 5, An. farauti 6, and An. farauti 7) and An. koliensis. Morphological identification, based on proboscis scale patterns and the presence of a sector spot on the costal wing vein, was often unreliable in distinguishing An. koliensis due to variable scaling patterns, even within isofemale lines [12–15].
Although restricted by morphological identification, early studies highlight heterogeneities in habitat preference, seasonality, human blood index, and transmission potential among members of this complex. Differences were observed in the spatial distribution of species among and within villages, larval habitats and vegetation. An. farauti s.s. is primarily found in coastal villages with a high tolerance for breeding in brackish water. An. farauti s.s. and An. hinesorum are more commonly found in natural breeding sites, such as ground pools, while An. punctulatus is commonly found in areas disturbed by human activity [14, 16]. An. punctulatus is abundant in the hills and An. koliensis has a patchy distribution in lowland inland areas , and both exploit different larval habitats . Abundance is correlated with recent rainfall, with An. koliensis showing greater temporal variability than An. punctulatus and An. farauti s.l. Peak outdoor biting times vary with the majority of An. farauti s.l. biting in the early evening and An. koliensis and An. punctulatus biting in the late night and early hours of the morning .
Pilot projects in PNG in the 1950s demonstrated the likelihood that DDT could successfully control malaria given the proper resources. It was not until the 1970s, through the support of United Nations Development Programme, that the DDT indoor residual spraying campaign was scaled up to cover over 50 % of the population . In Madang, DDT spraying was ineffective against An. farauti s.l., but very effective against An. punctulatus and moderately effective against An. koliensis . In the Solomon Islands, DDT residual spraying impacted members of the complex differently [20, 21] with a stronger impact on An. punctulatus and An. koliensis than An. farauti s.s.
Long-lasting, insecticide-treated nets (LLINs) can be a powerful tool in reducing malaria-associated morbidity, particularly when high community coverage is achieved . In PNG, where the vector population is 100 % susceptible to pyrethroids  LLINs remain effective for up to 5 years of in-home use , and are an attractive choice for malaria control. However, the tendency of the vector population to spend the majority of time outdoors  may render them less susceptible to indoor vector control due to low biological coverage. Furthermore, home-based interventions may have an unequal impact in areas with multiple vector species exhibiting a range of exophily and anthropophily.
The PNG Malaria Control Programme, supported by the Global Fund to Fight AIDS, Tuberculosis and Malaria, launched a countrywide, free, LLIN distribution between 2005 and 2009. PermaNet 2.0 LLINs, treated with 55 mg/cu m deltamethrin and manufactured by Vestergaard Frandsen, were distributed at a target ratio of one net per 2.5 household members. Independent household surveys reported 68.7 % of households surveyed in the Momase region (Morobe, Madang and Sepik Provinces included in this study) owned at least one LLIN (with 47 % sleeping underneath a LLIN the night before), while 95 % of households had at least one net of any type (with 74 % sleeping underneath the night before) . In addition, the malaria control programme introduced rapid diagnostic tests and artemisinin-based combination therapy in 2012 .
The purpose of this study was to determine the impact of the LLIN distribution on mosquito abundance, species composition, peak biting times, and malaria transmission in three geographical regions of PNG with multiple, sympatric vector species.
Mosquitoes were collected by the outdoor human landing catch method for 1 year prior to and 1 year following the nationwide LLIN distribution. Four representative villages (Matukar, Dimer, Nanaha, Yauatong) were chosen for further surveillance throughout year 3. Collectors were trained to aspirate host-seeking mosquitoes before biting in order to minimize exposure to malaria. Pairs of collectors worked in teams, with one individual collecting all mosquitoes landing on exposed legs from 1800 to 2400 hours and the second individual collecting from 24.00 to 06.00. The collectors rotated from the early shift to the late shift on subsequent nights. Mosquitoes were stored in cups according to the date, location and hour of collection. Following morphological identification [27, 28], mosquitoes were stored dry on silica gel and returned to the laboratory for mosquito species confirmation and enzyme-linked immunosorbent assay to detect Plasmodium spp. circumsporozoite protein. In the Madang villages, additional data were recorded hourly by each collector, including: the number of animals and number of additional people present in the hamlet, as well as a qualitative record of wind, rainfall and cloud cover (none, light, moderate or heavy).
In both inland and coastal regions of Madang, the first year of mosquito collections began in August 2008 until LLINs were distributed in July and August 2009. Year 2 collections began after LLIN distribution through September 2010, and year 3 collections concluded in November 2011. In the Dreikikir region, the first year of collections began in September 2008 until LLINs were distributed in late August 2009. Year 2 collections continued until July 2010 and year 3 collections concluded in July 2011. Mosquitoes were collected monthly in every village for the first 2 years. In the third year, collections occurred only in select villages with bimonthly collections in Dimer and Matukar and quarterly collections in Nanaha and Yauatong. Details of collection effort per year are shown in Additional file 1: Table S1.
Mosquitoes that were morphologically identified as members of the An. punctulatus group were confirmed to species by PCR-restriction fragment length polymorphism of the ITS2 region (4) using either an individual leg or extracted DNA (QIAGEN, Maryland, USA). Lysates from whole mosquitoes were screened for P. falciparum, P. vivax 210 and P. vivax 247 circumsporozoite proteins by enzyme-linked immunosorbent assay .
Rainfall data were collected by the PNG National Weather Service from Madang airport, 37–52 km from the Madang study villages. Although variations are to be expected between inland and coastal villages, regional rainfall records demonstrate typical seasonality and a lack of aberrant rainfall during the collection period.
This study was approved by the institutional review boards at University Hospitals Case Medical Center in Cleveland, Ohio, USA and the Institute of Medical Research and Medical Research Advisory Council of PNG.
Generalized estimating equations (GEE) were used to identify the determinants of mosquito abundance and infection prevalence. When the dependent variable was binary (sporozoite prevalence), a binomial distribution with a log link function was used. For count data (mosquito abundance), a Poisson distribution with log link function was used. In both cases, an exchangeable working correlation matrix was assumed. For sporozoite prevalence, separate models were constructed with either P. falciparum or P. vivax positivity as the dependent variable. Mosquito species (only An. punctulatus complex members were included) and year of the study (years 1, 2, 3) were covariates and village was the subject variable. For mosquito abundance, LLINs (presence or absence), rain in mm the previous month, the number of people and animals present outdoors within the compound, and the qualitative variables of rain, wind and cloud cover were included as covariates. Subject variables were village and collector, and within-subject variables were the date and hour of collection.
Species composition before and after the distribution are presented for the four intensively surveyed villages. Due to a significant association between P. falciparum infection and An. punctulatus, the proportion of this species of the total catch between year 1 and years 2 and 3 was compared by Fisher’s exact test.
Median biting times were calculated based on the entire catch of a given species per village and year. Not all collected species could be PCR-confirmed, therefore only villages that showed 95 % concordance between morphological and molecular identification in a given year are presented, and the data include all morphologically identified species. Mann–Whitney U tests were performed to determine if the medians between years were the same.
Mean man-biting rates were calculated for each village and compared between years 1 and 2 or 1 and 3 (t test with Bonferroni correction for multiple comparisons). Entomological inoculation rates, a measure of the number of infective bites per person per year, were calculated based on the average total number of bites per person per year and the sporozoite prevalence in mosquitoes.
The effect of environmental variables on hourly catch of anopheline mosquitoes
95 % CI
Mosquito populations have been extensively studied in Madang and East Sepik Provinces, where P. falciparum and P. vivax are both highly endemic. In Dreikikir, pre-LLIN sporozoite rates and man biting rates were similar to those observed in the late 1950s . In Madang Province, sporozoite rates are similar to previous studies [13, 31] but observed man biting rates were lower than historically reported [31, 32].
Overall, man biting rates dropped significantly in the year following LLIN distribution and remained low in representative villages through the third year. Even though members of the An. punctulatus group do not exhibit preferences for endophagy , LLINs may still have a large communal impact if LLIN coverage and usage is high . LLIN usage in study sites was higher than regionally reported . A survey 3 months following the distribution in Dreikikir revealed that 79 % had slept under a LLIN the previous night . A separate survey conducted in Madang study sites 3 years after the distribution revealed that 79 % had slept under a LLIN in Mirap and 54 % had slept under a LLIN in Wasab (JBK, unpublished). A separate study has shown that individual LLIN use and high community LLIN coverage were independently and strongly associated with reduced odds of malaria infection in PNG .
Although biting rates decreased, a trend towards earlier biting was observed in An. punctulatus and An. farauti s.s. Changes in biting behaviour were previously reported in An. farauti and An. koliensis, in a nearby village, following distribution of insecticide-treated nets in 1985 . In PNG, this shift could be indicative of behavioural flexibility or the impact of LLINs on population age structure; it has been observed that nulliparous females bite earlier than parous females . The former scenario has been suggested in an An. farauti s.l. population in the Solomon Islands following a DDT spraying campaign . Additional behavioural and genetic studies could elucidate whether selection is occurring in populations in this study. In the latter scenario, early biting patterns would not be accompanied by increased exposure even if biting densities returned to pre-control levels. No increase in the proportion of infective bites occurring before 22.00 in this study was observed, although sample sizes were limited. Continued surveillance can determine if the observed shifts in biting times are epidemiologically significant, and whether they will limit the protective effect of LLINs.
Differences in species composition were observed after LLIN distribution. However, there was not a consistent trend across villages towards species dominance. For example, in the Madang villages, An. farauti s.s. increased in proportion over An. punctulatus while in Dreikikir, An. punctulatus increased in proportion over An. koliensis. It is likely that some species will be more resilient to this intervention than others, and changes in species composition will depend on behaviour and physiology of the entire population. The observed pattern may be attributed to a greater impact of LLINs on the later biting species in the community. The observed increase in the proportion of An. punctulatus could be epidemiologically significant given the close association between this species and P. falciparum transmission. Further work is needed to determine whether these changes in species composition are durable, and whether they are associated with malaria prevalence.
While LLINs were very effective in reducing man biting rate and entomological inoculation rate in the three regions, sporozoite prevalence did not remain low. The observed increase in sporozoite prevalence between years 2 and 3 may be attributable to a low sample size, a challenge in post-intervention evaluations. However, it is clear that residual malaria transmission is still occurring, and that any rebound in mosquito density will be accompanied by an increase in transmission intensity. In both Madang regions, P. falciparum sporozoite rate was higher in year 3 than in year 1. This result is surprising because vector control is expected to have a greater impact on P. falciparum than P. vivax, due in part to the longer extrinsic incubation period for P. falciparum and the role of dormant hypnozoites in P. vivax transmission. These higher infection rates in year 3 may be influenced by low mosquito numbers and further work to determine long-term trends is needed.
This study has demonstrated a strong community impact of LLINs on exophagic and early biting vector populations. Success may be attributable to high LLIN ownership as well as a portion of the vector population that continues host seeking during the time when people are indoors. Individual malaria exposure will be affected by house construction, personal bed net use and other human behavioural patterns. Although LLINs had a clear impact on vector populations, the impact on sporozoite prevalence did not extend to year 3. Given the plasticity of mosquito behaviour and natural fluctuations in vector densities, continued surveillance is needed to determine changes in the effectiveness of LLINs.
PAZ, LJR, IM, MH, JWK, and PS conceived and designed the study; GK, JBK, LJR, and EKT supervised the field activities and conducted the research. LJR and EKT analysed the data; LJR wrote and revised the manuscript. All authors read and approved the final manuscript.
Funding provided by Fogarty International Center Ecology of Infectious Diseases TW007872 (PAZ), National Institutes of Health K01 TW008778 (LJR), NHMRC Senior Research Fellowship 1043345 (IM). We are grateful to data manager Thomas Adiguma, to laboratory assistants Absalom Mai and Adela Keven, and to field officers Lemen Kilepak, Siub Yabu, Yule E’ele, Muker Sakur, George Satokan, Daphne Sepe and Melinda Susapu.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Müller I, Bockarie M, Alpers M, Smith T. The epidemiology of malaria in Papua New Guinea. Trends Parasitol. 2003;19:253–9.PubMedView ArticleGoogle Scholar
- Hetzel MW, Morris H, Tarongka N, Barnadas C, Pulford J, Makita L, Siba PM, Mueller I. Prevalence of malaria across Papua New Guinea after initial roll-out of insecticide-treated mosquito nets. Trop Med Int Heal. 2015;20:1745–55. doi:10.1111/tmi.12616.View ArticleGoogle Scholar
- Mehlotra RK, Gray LR, Blood-Zikursh MJ, Kloos Z, Henry-Halldin CN, Tisch DJ, Thomsen E, Reimer L, Kastens W, Baea M, Baea K, Baisor M, Tarongka N, Kazura JW, Zimmerman PA. Short report: molecular-based assay for simultaneous detection of four Plasmodium spp. and Wuchereria bancrofti infections. Am J Trop Med Hyg. 2010;82:1030–3.PubMedPubMed CentralView ArticleGoogle Scholar
- Mehlotra RK, Kasehagen LJ, Baisor M, Lorry K, Kazura JW, Bockarie MJ, Zimmerman PA. Malaria infections are randomly distributed in diverse holoendemic areas of Papua New Guinea. Am J Trop Med Hyg. 2002;67:555–62.PubMedPubMed CentralGoogle Scholar
- Cooper RD, Waterson DGE, Frances SP, Beebe NW, Pluess B, Sweeney AW. Malaria vectors of Papua New Guinea. Int J Parasitol. 2009;39:1495–501.PubMedView ArticleGoogle Scholar
- Cooper RD, Waterson DG, Frances SP, Beebe NW, Sweeney AW. The anopheline fauna of Papua New Guinea. J Am Mosq Control Assoc. 2006;22:213–21.PubMedView ArticleGoogle Scholar
- Beebe NW, Russell T, Burkot TR, Cooper RD. Anopheles punctulatus group: evolution, distribution, and control. Annu Rev Entomol. 2015;60:335–50.PubMedView ArticleGoogle Scholar
- Burkot TR, Russell TL, Reimer LJ, Bugoro H, Beebe NW, Cooper RD, et al. Barrier screens: a method to sample blood-fed and host-seeking exophilic mosquitoes. Malar J. 2013;12:49.PubMedPubMed CentralView ArticleGoogle Scholar
- Cooper L, Cooper RD, Burkot TR. The Anopheles punctulatus complex: DNA probes for identifying the Australian species using isotopic, chromogenic, and chemiluminescence detection systems. Exp Parasitol. 1991;73:27–35.PubMedView ArticleGoogle Scholar
- Foley DH, Paru R, Dagoro H, Bryan JH. Allozyme analysis reveals six species within the Anopheles punctulatus complex of mosquitoes in Papua New Guinea. Med Vet Entomol. 1993;7:37.PubMedView ArticleGoogle Scholar
- Beebe NW, Saul A. Discrimination of all members of the Anopheles punctulatus complex by polymerase chain reaction—restriction fragment length polymorphism analysis. Am J Trop Med Hyg. 1995;53:478–81.PubMedGoogle Scholar
- Rozeboom LE, Knight KL. The punctulatus complex of Anopheles (diptera: culicidae). J Parasitol. 1946;32:95–131.PubMedView ArticleGoogle Scholar
- Benet A, Mai A, Bockarie F, Lagog M, Zimmerman P, Alpers MP, et al. Polymerase chain reaction diagnosis and the changing pattern of vector ecology and malaria transmission dynamics in Papua New Guinea. Am J Trop Med Hyg. 2004;71:277–84.PubMedGoogle Scholar
- Cooper RD, Waterson DGE, Frances SP, Beebe NW, Sweeney AW. Speciation and distribution of the members of the Anopheles punctulatus (diptera: culicidae) group in Papua New Guinea. J Med Entomol. 2002;39:16–27.PubMedView ArticleGoogle Scholar
- Woodhill AR. Observations on the morphology and biology of the subspecies of Anopheles punctulatus Donitz. Proc Linn Soc NSW. 1946;70:276–87.Google Scholar
- Hii JL, Smith T, Mai A, Mellor S, Lewis D, Alexander N, et al. Spatial and temporal variation in abundance of Anopheles (diptera:culicidae) in a malaria endemic area in Papua New Guinea. J Med Entomol. 1997;34:193–205.PubMedView ArticleGoogle Scholar
- Charlwood JD, Graves PM, Alpers MP. The ecology of the Anopheles punctulatus group of mosquitoes from Papua New Guinea: a review of recent work. P N G Med J. 1986;29:19–26.PubMedGoogle Scholar
- Charlwood JD, Bryan JH. A mark-recapture experiment with the filariasis vector Anopheles punctulatus in Papua New Guinea. Ann Trop Med Parasitol. 1987;81:429–36.PubMedGoogle Scholar
- Parkinson AD. Malaria in Papua New Guinea 1973. PNG Med J. 1974;17:8–16.Google Scholar
- Taylor B. Observations on malaria vectors of the Anopheles punctulatus complex in the British Solomon Islands Protectorate. J Med Entomol. 1975;11:677–87.PubMedView ArticleGoogle Scholar
- Bugoro H, Iro’ofa C, Mackenzie DO, Apairamo A, Hevalao W, Corcoran S, et al. Changes in vector species composition and current vector biology and behaviour will favour malaria elimination in Santa Isabel Province, Solomon Islands. Malar J. 2011;10:287.PubMedPubMed CentralView ArticleGoogle Scholar
- Lengeler C. Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004;(2):CD000363. doi:10.1002/14651858.CD000363.pub2.
- Keven JB, Henry-Halldin CN, Thomsen EK, Mueller I, Siba PM, Zimmerman PA, et al. Pyrethroid susceptibility in natural populations of the Anopheles punctulatus group (diptera: culicidae) in Papua New Guinea. Am J Trop Med Hyg. 2010;83:1259–61.PubMedPubMed CentralView ArticleGoogle Scholar
- Katusele M, Gideon G, Thomsen EK, Siba PM, Hetzel MM, Reimer LJ. Long-lasting insecticidal nets remain efficacious after five years of use in Papua New Guinea. PNG Med J. 2014;57:86–93.Google Scholar
- Hetzel MW, Gideon G, Lote N, Makita L, Siba PM, Mueller I. Ownership and usage of mosquito nets after four years of large-scale free distribution in Papua New Guinea. Malar J. 2012;11:192.PubMedPubMed CentralView ArticleGoogle Scholar
- Hetzel MW, Pulford J, Maraga S, Barnadas C, Reimer LJ, Tavul L, et al. Evaluation of the global fund-supported National malaria control program in Papua New Guinea, 2009-2014. PNG Med J. 2014;57:2009–14.Google Scholar
- Belkin JN. The mosquitoes of the South Pacific (diptera: culicidae), vol. 1. Berkeley: University of California Press; 1962.Google Scholar
- Lee DJ, Hicks MM, Griffiths M, Debenham ML, Bryan JH, Russell RC, Geary M, Marks EN. The culicidae of the Australasian Region, vol. 5. Canberra: Australian Government Publishing Service; 1987.Google Scholar
- Wirtz RA, Burkot TR, Graves PM, Andre RG. Field evaluation of enzyme-linked immunosorbent assays for Plasmodium falciparum and Plasmodium vivax sporozoites in mosquitoes (diptera: culicidae) from Papua New Guinea. J Med Entomol. 1987;24:433–7.PubMedView ArticleGoogle Scholar
- Peters W, Standfast HA. Studies on the epidemiology of malaria in New Guinea. Part II. Holoendemic malaria—the entomological picture. Trans R Soc Trop Med Hyg. 1960;54:249–54.PubMedView ArticleGoogle Scholar
- Burkot TR, Graves PM, Paru R, Wirtz RA, Heywood PF. Human malaria transmission studies in the Anopheles punctulatus complex in Papua New Guinea: sporozoite rates, inoculation rates, and sporozoite densities. Am J Trop Med Hyg. 1988;39:135–44.PubMedGoogle Scholar
- Afifi SE, Spencer M, Hudson PB, Tavil NW. Biting prevalence and malaria transmission patterns in the Anopheles punctulatus complex (diptera: culicidae) in Papua New Guinea. Aust J Exp Biol Med Sci. 1980;58:1–17.PubMedView ArticleGoogle Scholar
- Govella NJ, Okumu FO, Killeen GF. Insecticide-treated nets can reduce malaria transmission by mosquitoes which feed outdoors. Am J Trop Med Hyg. 2010;82:415–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Reimer LJ, Thomsen EK, Tisch DJ, Henry-Halldin CN, Zimmerman PA, Baea ME, Dagoro H, Susapu M, Hetzel MW, Bockarie MJ, Michael E, Siba PM, Kazura JW. Insecticidal bed nets and filariasis transmission in Papua New Guinea. N Engl J Med. 2013;369:745–53.PubMedView ArticleGoogle Scholar
- Charlwood JD, Graves PM. The effect of permethrin-impregnated bednets on a population of Anopheles farauti in coastal Papua New Guinea. Med Vet Entomol. 1987;1:319–27.PubMedView ArticleGoogle Scholar
- Bockarie MJ, Alexander N, Bockarie F, Ibam E, Barnish G, Alpers M. The late biting habit of parous Anopheles mosquitoes and pre-bedtime exposure of humans to infective female mosquitoes. Trans R Soc Trop Med Hyg. 1996;90:23–5.PubMedView ArticleGoogle Scholar