Emergence of behavioural avoidance strategies of malaria vectors in areas of high LLIN coverage in Tanzania.

Despite significant reductions in malaria transmission across Africa since 2000, progress is stalling. This has been attributed to the development of insecticide resistance and behavioural adaptations in malaria vectors. Whilst insecticide resistance has been widely investigated, there is poorer understanding of the emergence, dynamics and impact of mosquito behavioural adaptations. We conducted a longitudinal investigation of malaria vector host choice over 3 years and resting behaviour over 4 years following a mass long-lasting insecticidal nets (LLINs) distribution in Tanzania. By pairing observations of mosquito ecology with environmental monitoring, we quantified longitudinal shifts in host-choice and resting behaviour that are consistent with adaptation to evade LLINs. The density of An. funestus s.l., declined significantly through time. In tandem, An. arabiensis and An. funestus s.l. exhibited an increased rate of outdoor relative to indoor resting; with An. arabiensis reducing the proportion of blood meals taken from humans in favour of cattle. By accounting for environmental variation, this study detected clear evidence of intra-specific shifts in mosquito behaviour that could be obscured in shorter-term or temporally-coarse surveys. This highlights the importance of mosquito behavioural adaptations to vector control, and the value of longer-term behavioural studies.

Investigating associations between biting time in the malaria vector Anopheles arabiensis Patton and single nucleotide polymorphisms in circadian clock genes: support for sub-structure among An. arabiensis in the Kilombero valley of Tanzania.

BACKGROUND: There is growing evidence that the widespread use of Long-Lasting Insecticidal Nets (LLINs) is prompting malaria vectors to shift their biting towards times and places where people are not protected, such as earlier in the evening and/or outdoors. It is uncertain whether these behavioural shifts are due to phenotypic plasticity and/or ecological changes within vector communities that favour more exophilic species, or involve genetic factors within vector species to limit their contact with LLINs. Possibly variation in the time and location of mosquito biting has a genetic basis, but as yet this phenomenon has received little investigation. Here we used a candidate gene approach to investigate whether polymorphisms in selected circadian clock genes could explain variation in the time and location of feeding (indoors versus outside) within a natural population of the major African malaria vector Anopheles arabiensis. METHODS: Host-seeking An. arabiensis were collected from two villages (Lupiro and Sagamaganga) in Tanzania by Human Landing Catch (HLC) technique. Mosquitoes were classified into phenotypes of "early" (7 pm-10 pm) or "late" biting (4 am -7 am), and host-seeking indoors or outdoors. In these samples we genotyped 34 coding SNPs in 8 clock genes (PER, TIM, CLK, CYC, PDP1, VRI, CRY1, and CRY2), and tested for associations between these SNPs and biting phenotypes. SNPs in 8 mitochondrial genes (ATP6, ATP8, COX1, COX2, COX3, ND3, ND5 and CYTB) were also genotyped to test population subdivision within An. arabiensis. RESULTS: The candidate clock genes exhibited polymorphism within An. arabiensis, but it was unrelated to variation in the timing and location of their biting activity. However, there was evidence of strong genetic structure within An. arabiensis populations in association with the TIM, which was unrelated to geographic distance. Substructure within An. arabiensis was also detected using mitochondrial markers. CONCLUSIONS: The variable timing and location of biting in An. arabiensis could not be linked to candidate clock genes that are known to influence behaviour in other Diptera. This finding does not rule out the possibility of a genetic basis to biting behaviour in this malaria vector, but suggests these are complex phenotypes that require more intensive ecological, neuronal and genomic analyses to understand.

A set of broadly applicable microsatellite markers for analyzing the structure of Culex pipiens (Diptera: Culicidae) populations.

Microsatellite markers were isolated and developed from Culex pipiens quinquefasciatus Say (Diptera: Culicidae) sampled in Johannesburg, South Africa, to identify those that are broadly useful for analyzing Cx. pipiens complex populations between continents. Suitable loci should be 1) inherited in a codominant Mendelian manner, 2) polymorphic, 3) selectively neutral, 4) randomly associated, 5) without null alleles, and 6) applicable across broad regions and between diverse biotypes. Loci in Cx. p. quinquefasciatus from Johannesburg ranged from two to 17 alleles per locus and expected heterozygosities (H(e)) were 0.02-0.87. Loci in Cx. p. pipiens L. from Johannesburg had five to 19 alleles per locus and H(e) values ranging from 0.57 to 0.93, whereas those from George, South Africa, had five to 17 alleles per locus and H(e) values ranging from 0.54 to 0.88. Loci in North American mosquitoes were more variable. Cx. p. quinquefasciatus from South Carolina had five to 19 alleles per locus and H(e) values ranging from 0.64 to 0.90, whereas Cx. p. pipiens from Massachusetts had six to 28 alleles per locus and with H(e) values ranging from 0.65 to 0.94. All loci were associated randomly. Overall, four of nine of these new loci satisfied all six criteria for broad utility for analyzing the genetic structure of Cx. pipiens populations.

Evidence for subdivision within the M molecular form of Anopheles gambiae.

The principal vector of malaria in sub-Saharan Africa, Anopheles gambiae is subdivided into two molecular forms M and S. Additionally, several chromosomal forms, characterized by the presence of various inversion polymorphisms, have been described. The molecular forms M and S each contain several chromosomal forms, including the Savanna, Mopti and Forest forms. The M and S molecular forms are now considered to be the reproductive units within A. gambiae and it has recently been argued that a low recombination rate in the centromeric region of the X chromosome has facilitated isolation between these forms. The status of the chromosomal forms remains unclear however. Therefore, we studied genetic differentiation between Savanna S, Forest S, Forest M and Mopti M populations using microsatellites. Genetic differentiation between Savanna S and Forest S populations is very low (F(ST) = 0.0053 +/- 0.0049), even across large distances. In comparison, the Mopti M and Forest M populations show a relatively high degree of genetic differentiation (F(ST) = 0.0406 +/- 0.0054) indicating that the M molecular form may not be a single entity, but could be subdivided into at least two distinct chromosomal forms. Previously it was proposed that inversions have played a role in the origin of species within the A. gambiae complex. We argue that a possible subdivision within the M molecular form could be understood through this process, with the acquisition of inversions leading to the expansion of the M molecular form into new habitat, dividing it into two distinct chromosomal forms.

An unusual distribution of the kdr gene among populations of Anopheles gambiae on the island of Bioko, Equatorial Guinea.

In West Africa, Anopheles gambiae exists in discrete subpopulations known as the M and S molecular forms. Although these forms occur in sympatry, pyrethroid knock-down resistance (kdr) is strongly associated with the S molecular form. On the island of Bioko, Equatorial Guinea we found high frequencies of the kdr mutation in M form individuals (55.8%) and a complete absence of kdr in the S form. We also report the absence of the kdr allele in M and S specimens from the harbour town of Tiko in Cameroon, representing the nearest continental population to Bioko. The kdr allele had previously been reported as absent in populations of An. gambiae on Bioko. Contrary to earlier reports, sequencing of intron-1 of this sodium channel gene revealed no fixed differences between M form resistant and susceptible individuals. The mutation may have recently arisen independently in the M form on Bioko due to recent and intensive pyrethroid application.

Effect of seminal fluids in mating between M and S forms of Anopheles gambiae.

Previous studies have shown that sympatric populations of M and S molecular forms of Anopheles gambiae sensu stricto exhibit strong assortative mating. In the few documented cases of cross-mating between M and S forms, females that mated with amale of the alternative form were often also mated with a male of their own form. A potential explanation for the association between cross-mating and double mating could be that male accessory gland or sperm proteins that are responsible for inducing refractoriness to further mating by females have diverged between the M and S forms. This mechanism of postmating reproductive isolation would have important implications for our understanding of the speciation processes in the An. gambiae complex. We tested for this mechanism, by comparing the likelihood of mating, feeding, and laying eggs, as well as the fertility of females presented with males of their own form or the alternate form in the laboratory. We also compared the likelihood of remating in cross-mated and assortatively-mated females, and we analyzed their progeny to unravel patterns of sperm precedence. We found that cross-mated females differed from assortatively-mated females only in terms of egg-hatching rate and larval survival but that these effects could be attributed to hybrid vigor rather than differential response to seminal products. Cross-mating between forms was not associated with remating behavior. These results indicate that the sex proteins responsible for inhibiting further insemination and triggering the gonotrophic cycle in females have not diverged between these M and S populations. We discuss alternative explanations for the patterns of cross-mating and multiple mating observed in the field.

Genetic variation in the sand fly salivary protein, SP-15, a potential vaccine candidate against Leishmania major.

SP-15 is a sandfly salivary protein that provides strong protection against cutaneous leishmaniasis, caused by Leishmania major, and has been proposed as a potential vaccine against this disease. To investigate possible antigenic variation in this protein, we examined genetic polymorphism of SP-15 in 100 Phlebotomus papatasi sandflies, from a natural population from Sudan and four laboratory colonies from Egypt, Jordan, Israel and Saudi Arabia. We found that although many variants of SP-15 may be found in nature, differences among them are minimal (mean+/-SD pairwise differences=1.69+/-0.83% for forty nucleotide sequences and 3.06+/-1.13% for thirty amino acid sequence variants). Analysis of proportions of synonymous and non-synonymous substitutions indicated that SP-15 is not under diversifying selection. Our results suggest that a vaccine based on SP-15 protein should result in a uniform immune response.

Genetic differentiation of Anopheles gambiae s.s. populations in Mali, West Africa, using microsatellite loci.

Anopheles gambiae sensu stricto is a principal vector of malaria through much of sub-Saharan Africa, where this disease is a major cause of morbidity and mortality in human populations. Accordingly, population sizes and gene flow in this species have received special attention, as these parameters are important in attempts to control malaria by impacting its mosquito vector. Past measures of genetic differentiation have sometimes yielded conflicting results, in some cases suggesting that gene flow is extensive over vast distances (6000 km) and is disrupted only by major geological disturbances and/or barriers. Using microsatellite DNA loci from populations in Mali, West Africa, we measured genetic differentiation over uniform habitats favorable to the species across distances ranging from 62 to 536 km. Gene flow was strongly correlated with distance (r(2) = 0.77), with no major differences among chromosomes. We conclude that in this part of Africa, at least, genetic differentiation for microsatellite DNA loci is consistent with traditional models of isolation by distance.

Phylogeography of the neotropical sand fly Lutzomyia longipalpis inferred from mitochondrial DNA sequences.

Sand flies in the Lutzomyia longipalpis species complex include the primary vector of Leishmania chagasi, the etiologic agent of visceral leishmaniasis in the Neotropics. Twelve L. longipalpis populations from South and Central America were compared using the cytochrome c oxidase I (COI) gene from the mitochondrial genome. The haplotype profiles for each population revealed that the majority of sequence variation was inter-population (98%) rather than intra-population, suggesting that sequence polymorphisms at the COI locus should provide excellent characters for the study of phylogenetic relationships among populations. Phylogenetic reconstruction using distance (neighbor-joining) and maximum parsimony analysis revealed the existence of four clades among the L. longipalpis populations studied: (1) Laran, (2) Brazilian, (3) cis-Andean and (4) trans-Andean. We suggest that these clades represent species. A biogeographical interpretation of the molecular phylogeny suggests that the process of speciation in the L. longipalpis complex began in the Pliocene, from a sub-Andean-Amazonian gene pool resulting from the Andean orogeny (formation of the East Andean Cordillera). The four clades probably diverged as a result of vicariance events that occurred throughout the late Pliocene and Pleistocene. We propose and discuss several historical scenarios, based on the biogeography and historical geology of Central and South America.

Sandfly maxadilan exacerbates infection with Leishmania major and vaccinating against it protects against L. major infection.

Bloodfeeding arthropods transmit many of the world's most serious infectious diseases. Leishmania are transmitted to their mammalian hosts when an infected sandfly probes in the skin for a bloodmeal and injects the parasite mixed with its saliva. Arthropod saliva contains molecules that affect blood flow and modulate the immune response of the host. Indeed, sandfly saliva markedly enhances the infectivity of L. major for its host. If the salivary molecule(s) responsible for this phenomenon was identified, it might be possible to vaccinate the host against this molecule and thereby protect the host against infection with Leishmania. Such an approach represents a novel means of controlling arthropod-borne disease transmission. Here, we report that a single molecule, maxadilan, in sandfly saliva can exacerbate infection with L. major to the same degree as whole saliva, and that vaccinating against maxadilan protects mice against infection with L. major.

DNA analysis of transferred sperm reveals significant levels of gene flow between molecular forms of Anopheles gambiae.

Anopheles gambiae populations in west Africa are complex, being composed of multiple, sympatric subpopulations. Recent studies have failed to reveal significant genetic differences among subpopulations, stimulating a debate regarding the levels of gene flow among them. The observed homogeneity may be the consequence of substantial contemporary gene flow or it may be that reproductive isolation is complete, but too recent for the accumulation of significant levels of genic divergence. Here, we report the results of a study estimating contemporary levels of gene flow between An. gambiae subpopulations by analysing females and transferred sperm removed from their reproductive systems. A total of 251 female and associated sperm extracts was analysed from a single site in Mali. Two molecular forms of An. gambiae, the M- and S-forms, occurred in sympatry at this site. Overall, we found very strong positive assortative mating within forms, however, we did observe significant hybridization between forms. In the M subpopulation 2/195 females (1.03%) contained sperm from S-form males and in 55 S-form females we found one female containing M-form sperm (1.82%). We also identified a mated M xS hybrid adult female. From mating frequencies, we estimate the Nem between the M- and S-form at 16.8, and from the adult hybrid frequency at 5.6. These values are consistent with our earlier estimate, based on FST for 21 microsatellite loci in which Nem = 5.8. We conclude that the general lack of genetic divergence between the M and S subpopulations of An. gambiae can be explained entirely by contemporary gene flow.

Microsatellite DNA polymorphism and heterozygosity among field and laboratory populations of Anopheles gambiae ss (Diptera: Culicidae).

We compared microsatellite polymorphism at nine loci located on chromosome 3 among two colonies and a field population of Anopheles gambiae sensu stricto Giles mosquitoes. Numbers of microsatellite alleles observed at each locus and mean heterozygosities were drastically reduced among laboratory colonies. Genetic analysis of the field population used in this study revealed an unprecedented frequency of rare alleles (<0.05). In contrast, colony samples revealed large numbers of alleles with frequencies >0.50. Partitioning of field data to assess the impact of rare alleles, null alleles, and sample size on estimates of mean heterozygosity revealed the plasticity of this measurement and suggests that heterozygosity may be reliably estimated from relatively small collections using microsatellites.

Distribution of genetic variation among chromosomal forms of Anopheles gambiae s.s: introgressive hybridization, adaptive inversions, or recent reproductive isolation?

A series of four papers in this issue explores the reproductive status of the five chromosomal forms of An. gambiae s.s. using molecular techniques to examine the variation among twelve different genes located throughout the An. gambiae s.s. genome. Results of these and previous studies are consistent with a hypothesis of at least partial barriers to gene flow between some chromosomal forms in the Ivory Coast and other West African countries to the north and west, but introgression between S and M types in Benin and countries to the east. Collectively, these studies indicate the need for a broader geographical sampling of An. gambiae s.s., increased research on mechanisms of prezygotic reproductive isolation and field-based studies of survival and fecundity in hybrids to test for postzygotic reproductive isolation.

Gene flow among populations of the malaria vector, Anopheles gambiae, in Mali, West Africa.

The population structure of the Anopheles gambiae complex is unusual, with several sibling species often occupying a single area and, in one of these species, An. gambiae sensu stricto, as many as three "chromosomal forms" occurring together. The chromosomal forms are thought to be intermediate between populations and species, distinguishable by patterns of chromosome gene arrangements. The extent of reproductive isolation among these forms has been debated. To better characterize this structure we measured effective population size, N(e), and migration rates, m, or their product by both direct and indirect means. Gene flow among villages within each chromosomal form was found to be large (N(e)m > 40), was intermediate between chromosomal forms (N(e)m approximately 3-30), and was low between species (N(e)m approximately 0.17-1.3). A recently developed means for distinguishing among certain of the forms using PCR indicated rates of gene flow consistent with those observed using the other genetic markers.

Sibling species in the Llutzomyia longipalpis complex differ in levels of mRNA expression for the salivary peptide, maxadilan.

Maxadilan is a small ( approximately 7 kDa) protein found in the saliva of sand fly species in the Lutzomyia longipalpis complex, vectors of the parasite causing visceral leishmaniasis, Leishmania chagasi. It is a potent vasodilator and also has immunomodulatory affects. Maxadilan recovered from different sibling species of the Lu. longipalpis complex differ in amino acid content by as much as 23%, however all variants possess equivalent vasodilatory activity. Therefore, the dramatic differences in vasodilatory activity of the saliva from different sibling species is probably due to differences in the amounts of maxadilan in their saliva. This is significant because it has been suggested that maxadilan may influence the pathogenesis of leishmanial infections. In this study we measured the amount of maxadilan messenger RNA (mRNA) per pair of salivary glands from individual sand flies by quantitative reverse transcription polymerase chain reaction (RT-PCR) using a competitive method. We report a method using the gene of interest, in this case maxadilan, amplified by the PCR from genomic DNA, as a competitor in the quantitative RT-PCR, taking advantage of differences in the size of these products due to the presence of an intron. Significant differences in amounts of maxadilan mRNA among colonies from Central and South America are described. We found a strong correlation between the amount of maxadilan mRNA detected in salivary glands of different Lu. longipalpis sibling species and previously described differences in the size of erythemas produced by the bite of these species. Therefore, variation in the amount of mRNA suggests that differences in the vasodilatory properties of saliva among the different sibling species are the result of differences in the amount of maxadilan present in the saliva and not differences in the potency of maxadilan peptide variants. The geographical distribution of species with high or low levels of maxadilan gene expression are concordant with the distribution of atypical cutaneous leishmaniasis resulting from infection with Le. chagasi, lending credence to earlier suggestions that maxadilan may be involved with visceralization of this parasite.

Breeding structure of the sand fly Lutzomyia longipalpis (Lutz & Neiva) in Brazil.

Eleven populations of Lutzomyia longipalpis (Lutz & Neiva), the sand fly vector of Leishmania chagasi, from different areas of Brazil were analyzed for genetic variation at 16 enzyme loci. In this region, the prevalence of visceral leishmaniasis (VL) caused by L. chagasi is spotty and reproductive isolation among populations of Lu. longipalpis has been reported. It is thought that morphologically similar cryptic species with varying vectorial capacity may be responsible for the discontinuous distribution of VL. The aim was to study the genetic structure of populations within this region and to identify demes that may represent sibling species. Genotypic frequencies within populations were in close compliance to Hardy-Weinberg expectations, suggesting there are no sympatric species among these 11 populations. Levels of genetic distance between pairs of populations were very low (< 0.03), consistent with local populations within a single sand fly species. When genotypic frequency data for all populations were pooled, 9 of the 13 polymorphic loci deviated from Hardy-Weinberg expectations, indicating some degree of genetic substructuring. Estimates of effective migration rates (N(e)m) among all populations were low, 2.73, suggesting that gene flow is restricted among populations, which is probably the reason for the observed genetic substructuring.

Variation in the salivary peptide, maxadilan, from species in the Lutzomyia longipalpis complex.

Maxadilan is an approximately 7kDa peptide that occurs in the saliva of the sand fly Lutzomyia longipalpis. This peptide is a potent vasodilator and may also have immunomodulatory effects related to the pathogenesis of leishmanial infections. Variation in the primary DNA and inferred amino acid sequence of maxadilan is reported. Differences were found within and among natural field populations as well as among sibling species. Extensive amino acid sequence differentiation, up to 23%, was observed among maxadilan from different populations. This is a remarkable degree of polymorphism considering the small size of this peptide. The vasodilatory activity of maxadilan was equivalent among recombinant maxadilan variants. All maxadilan variants induce interleukin-6. Predicted secondary structure and hydrophobicity plots suggest that these characteristics are conserved among variant peptides. However, profiles based on the antigenic index do differ among peptides.

Metaphase karyotypes and G-banding in sandflies of the Lutzomyia longipalpis complex.

Mitotic metaphase chromosomes (2n = 8) from brain cells of fourth instar sandfly larvae of four geographical strains of the Lutzomyia longipaplis complex were examined microscopically, with bright-field illumination, after staining by a new G-banding technique involving exposure of air-dried chromosome preparations to quinacrine and ultraviolet light. Differences of G-banding and/or position of the centromere on chromosome 4 (the smallest chromosome pair) distinguished four putative sibling species from Costa Rica, Colombia and Brazil (distinctive populations from Jacobina and Lapinha Caves). The karyotype of the population from Jacobina, Brazil, showed an apparently plesiomorphic pattern of G-banding. On the basis of their recognizably different mitotic karyotypes, cytogenetic identification of separate taxa in the L. longipalpis complex should be useful for specific female vector competence and ecology studies.

Complexities in the genetic structure of Anopheles gambiae populations in west Africa as revealed by microsatellite DNA analysis.

Chromosomal forms of Anopheles gambiae, given the informal designations Bamako, Mopti, and Savannah, have been recognized by the presence or absence of four paracentric inversions on chromosome 2. Studies of karyotype frequencies at sites where the forms occur in sympatry have led to the suggestion that these forms represent species. We conducted a study of the genetic structure of populations of An. gambiae from two villages in Mali, west Africa. Populations at each site were composed of the Bamako and Mopti forms and the sibling species, Anopheles arabiensis. Karyotypes were determined for each individual mosquito and genotypes at 21 microsatellite loci determined. A number of the microsatellites have been physically mapped to polytene chromosomes, making it possible to select loci based on their position relative to the inversions used to define forms. We found that the chromosomal forms differ at all loci on chromosome 2, but there were few differences for loci on other chromosomes. Geographic variation was small. Gene flow appears to vary among different regions within the genome, being lowest on chromosome 2, probably due to hitchhiking with the inversions. We conclude that the majority of observed genetic divergence between chromosomal forms can be explained by forces that need not involve reproductive isolation, although reproductive isolation is not ruled out. We found low levels of gene flow between the sibling species Anopheles gambiae and Anopheles arabiensis, similar to estimates based on observed frequencies of hybrid karyotypes in natural populations.