Modeling the role of environmental variables on the population dynamics of the malaria vector Anopheles gambiae sensu stricto

ABSTRACT: BACKGROUND: The impact of weather and climate on malaria transmission has attracted considerable attention in recent years, yet uncertainties around future disease trends under climate change remain. Mathematical models provide powerful tools for addressing such questions and understanding the implications for interventions and eradication strategies, but these require realistic modeling of the vector population dynamics and its response to environmental variables. METHODS: Published and unpublished field and experimental data are used to develop new formulations for modeling the relationships between key aspects of vector ecology and environmental variables. These relationships are integrated within a validated deterministic model of Anopheles gambiae s.s. population dynamics to provide a valuable tool for understanding vector response to biotic and abiotic variables. RESULTS: A novel, parsimonious framework for assessing the effects of rainfall, cloudiness, wind speed, desiccation, temperature, relative humidity and density-dependence on vector abundance is developed, allowing ease of construction, analysis, and integration into malaria transmission models. Model validation shows good agreement with longitudinal vector abundance data from Tanzania, suggesting that recent malaria reductions in certain areas of Africa could be due to changing environmental conditions affecting vector populations. CONCLUSIONS: Mathematical models provide a powerful, explanatory means of understanding the role of environmental variables on mosquito populations and hence for predicting future malaria transmission under global change. The framework developed provides a valuable advance in this respect, but also highlights key research gaps that need to be resolved if we are to better understand future malaria risk in vulnerable communities.     

Effect of temperature on the development of the aquatic stages of Anopheles gambiae sensu stricto (Diptera: Culicidae)

Global warming may affect the future pattern of many arthropod-borne diseases, yet the relationship between temperature and development has been poorly described for many key vectors. Here the development of the aquatic stages of Africa's principal malaria vector, Anopheles gambiae s.s. Giles, is described at different temperatures. Development time from egg to adult was measured under laboratory conditions at constant temperatures between 10 and 40 degrees C. Rate of development from one immature stage to the next increased at higher temperatures to a peak around 28 degrees C and then declined. Adult development rate was greatest between 28 and 32 degrees C, although adult emergence was highest between 22 and 26 degrees C. No adults emerged below 18 degrees C or above 34 degrees C. Non-linear models were used to describe the relationship between developmental rate and temperature, which could be used for developing process-based models of malaria transmission. The utility of these findings is demonstrated by showing that a map where the climate is suitable for the development of aquatic stages of A. gambiae s.s. corresponded closely with the best map of malaria risk currently available for Africa.     

Differential behaviour of Anopheles gambiae sensu stricto (Diptera: Culicidae) to human and cow odours in the laboratory

A dual port olfactometer was used to study the response of Anopheles gambiae Giles sensu stricto to odours of human and animal origin. Human odour consisted of human skin emanations collected on a nylon stocking, which was worn for 24 h. This was tested alone or together with 4.5% carbon dioxide, the concentration in human and cattle breath. Cattle odours consisted of cow skin emanations and/or carbon dioxide. Cow skin emanations were collected by tying a nylon stocking ('cow sock') around the hind leg of a cow for 12 h. Anopheles gambiae s.s. was consistently highly attracted by human odour, which is consistent with the high degree of anthropophily in this mosquito. Anopheles gambiae s.s. was not attracted by human or cattle equivalent volumes of carbon dioxide and this gas did not enhance the effect of human skin residues. Furthermore, A. gambiae s.s. showed a high degree of aversion to cow odour. When human odour and cow odour were tested together in the same port, mosquitoes were still highly attracted, indicating that whilst cattle odour may deter A. gambiae s.s., these mosquitoes can detect human odour in the presence of cattle odour. It was concluded that carbon dioxide plays a minor role in the host seeking behaviour of A. gambiae s.s., whilst host specific cues such as human skin residues play a major role and very effectively demonstrated anthropophilic behaviour in the laboratory.     

Mapping distributions of chromosomal forms of Anopheles gambiae in West Africa using climate data

The mosquito Anopheles gambiae Giles sensu stricto (Diptera: Culicidae), the principal vector of malaria in West Africa, comprises several chromosomal forms (e.g. Bissau, Forest, Mopti, Savanna) associated with climatic zones. Here we show how climate data can be used to map the geographical distribution of these chromosomal forms. The climate at 144 sites surveyed for mosquitoes in West Africa between 1971 and 92 was determined using computerized climate surfaces. Forest and Bissau forms occurred at relatively wet sites: median annual precipitation 1325 mm and 1438 mm, respectively, interquartile ranges (IQR) 1144-1858 mm and 1052-1825 mm), whilst the Mopti form was found at dry sites (annual 938 mm, IQR 713-1047 mm) and the Savanna form at sites intermediate between the wet and dry forms (annual 1067 mm, IQR 916-1279). Logistic regression analyses of the climate variables were carried out on a stratified random sample of half the sites. The resulting models correctly classified over 80% of the sites for presence or absence of each chromosomal form. When these models were tested against excluded sites they were also correct at over 80% of sites. The combined data produced models that were correct at over 86% of sites. Mean annual precipitation, evapotranspiration, minimum temperature and maximum temperature were the most important climate variables correlated with the distribution of these forms of An. gambiae. We used the logistic models to map the distribution of each chromosomal form within the reported range for An. gambiae s.s. in West Africa employing a geographical information system. Our maps indicate that each chromosomal form favours particular climate envelopes in well-defined ecoclimatic zones, although these forms are sympatric at the edges of their ranges. This study demonstrates that climate can be used to map the distribution of chromosomal forms of insects across large areas.     

Islands and Stepping-Stones: Comparative Population Structure of Anopheles gambiae sensu stricto and Anopheles arabiensis in Tanzania and Implications for the Spread of Insecticide Resistance

Population genetic structures of the two major malaria vectors Anopheles gambiae s.s. and An. arabiensis, differ markedly across Sub-Saharan Africa, which could reflect differences in historical demographies or in contemporary gene flow. Elucidation of the degree and cause of population structure is important for predicting the spread of genetic traits such as insecticide resistance genes or artificially engineered genes. Here the population genetics of An. gambiae s.s. and An. arabiensis in the central, eastern and island regions of Tanzania were compared. Microsatellite markers were screened in 33 collections of female An. gambiae s.l., originating from 22 geographical locations, four of which were sampled in two or three years between 2008 and 2010. An. gambiae were sampled from six sites, An. arabiensis from 14 sites, and both species from two sites, with an additional colonised insectary sample of each species. Frequencies of the knock-down resistance (kdr) alleles 1014S and 1014F were also determined. An. gambiae exhibited relatively high genetic differentiation (average pairwise F_ST = 0.131), significant even between nearby samples, but without clear geographical patterning. In contrast, An. arabiensis exhibited limited differentiation (average F_ST = 0.015), but strong isolation-by-distance (Mantel test r = 0.46, p = 0.0008). Most time-series samples of An. arabiensis were homogeneous, suggesting general temporal stability of the genetic structure. An. gambiae populations from Dar es Salaam and Bagamoyo were found to have high frequencies of kdr 1014S (around 70%), with almost 50% homozygote but was at much lower frequency on Unguja Island, with no. An. gambiae population genetic differentiation was consistent with an island model of genetic structuring with highly restricted gene flow, contrary to An. arabiensis which was consistent with a stepping-stone model of extensive, but geographically-restricted gene flow.     

Molecular divergence of the soil yeasts Williopsis sensu stricto

Fifty-three strains of saturn-spored yeasts were analyzed by means of restriction analysis of the amplified fragment of rDNA which comprised the 5.8S rRNA gene and the internal transcribed spacers ITS1 and ITS2. The use of endonucleases HaeIII and MspI enabled clear differentiation of yeast species Williopsis mucosa, W. salicorniae, Zygowilliopsis californica, Komagataea pratensis, and the Williopsis sensu stricto complex. Minisatellite primer M13 was proposed for the differentiation between twin species of Williopsis sensu stricto, which have identical restriction profiles. PCR with primer M13 enabled reidentification of a number of collection strains, species identification of saturn-spored isolates from the Far East, and detection of three strains affiliated to novel taxa. The latter have unique PCR profiles and differ in the nucleotide sequences of ITS1 and ITS2 fragments of rDNA. Possible variations in the results obtained with different molecular methods are discussed.     

Molecular identification of a myosuppressin receptor from the malaria mosquito Anopheles gambiae

The insect myosuppressins (X1DVX2HX3FLRFamide) are neuropeptides that generally block insect muscle activities. We have used the genomic sequence information from the malaria mosquito Anopheles gambiae Genome Project to clone a G protein-coupled receptor that was closely related to the two previously cloned and characterized myosuppressin receptors from Drosophila [Proc. Natl. Acad. Sci. USA 100 (2003) 9808]. The mosquito receptor cDNA was expressed in Chinese hamster ovary cells and was found to be activated by low concentrations of Anopheles myosuppressin (TDVDHVFLRFamide; EC50, 1.6 x 10(-8)M). The receptor was not activated by a library of 35 other insect neuropeptides and monoamines, including neuropeptides that resembled myosuppressin in their C-terminal moiety, such as PDRNFLRFamide (Anopheles FMRFamide-3), other Anopheles FMRFamide peptides, or neuropeptide F-like peptides, showing that the receptor was quite selective for myosuppressin. These results also showed that the myosuppressin receptor needs a much larger portion than the C-terminal FLRFamide sequence for its activation. The insect myosuppressins are often grouped together with the insect FMRFamides under the name FaRPs (FMRFamide-related peptides). However, this is not justified anymore, because the insect myosuppressin receptor/ligand couple is both functionally and evolutionarily fully unrelated to the insect FMRFamide receptor/ligand couple. To our knowledge, this is the first report on the molecular identification of a mosquito neuropeptide receptor.     

Variation in cuticular hydrocarbons among strains of the Anopheles gambiae sensu stricto by analysis of cuticular hydrocarbons using gas liquid chromatography of larvae

Cuticular hydrocarbons of larvae of individual strains of the Anopheles gambiae sensu stricto were investigated using gas liquid chromatography. Biomedical discriminant analysis involving multivariate statistics suggests that there was clear hydrocarbon difference between the Gambian(G3), the Nigerian (16CSS and, its malathion resistant substrain, REFMA) and the Tanzanian (KWA) strains. The high degree of segregation (95%) in hydrocarbons among the four strains investigated indicates that further analysis is needed to enable understanding of hydrocarbon variation in samples of An. gambiae especially from areas where these populations co-exist.     

Genetic identification of biological species in theSaccharomyces sensu stricto complex

Studies on taxonomic and evolutionary genetics of theSaccharomyces sensu stricto complex are considered in light of the biological species concept. Genetic variability of some physiological properties traditionally used in yeast taxonomy is discussed. Genetic hybridization analysis and molecular karyotyping revealed six biological species in theSaccharomyces sensu stricto complex. DNA-DNA reassociation data are concordant with the data obtained by genetic analysis. A new system for naming the cultivatedSaccharomyces yeast (groups of cultivars) is proposed.This paper is dedicated to Danish scientists Ö Winge and V Jensen in recognition of their contributions to zymology.     

Molecular analysis of alpha-galactosidase MEL genes from Saccharomyces sensu stricto

To infer the molecular evolution of yeast Saccharomyces sensu stricto from analysis of the alpha-galactosidase MEL gene family, two new genes were cloned and sequenced from S. bayanus var. bayanus and S. pastorianus. Nucleotide sequence homology of the MEL genes of S. bayanus var. bayanus (MELb), S. pastorianus (MELpt), S. bayanus var. uvarum (MELu), and S. carlsbergensis (MELx) was rather high (94.1-99.3%), comparable with interspecific homology (94.8-100%) of S. cerevisiae MEL1-MEL11. Homology of the MEL genes of sibling species S. cerevisiae (MEL1), S. bayanus (MELb), S. paradoxus (MELp), and S. mikatae (MELj) was 76.2-81.7%, suggesting certain species specificity. On this evidence, the alpha-galactosidase gene of hybrid yeast S. pastorianus (S. carlsbergensis) was assumed to originate from S. bayanus rather than from S. cerevisiae.     

The "far-west" of Anopheles gambiae molecular forms

The main Afrotropical malaria vector, Anopheles gambiae sensu stricto, is undergoing a process of sympatric ecological diversification leading to at least two incipient species (the M and S molecular forms) showing heterogeneous levels of divergence across the genome. The physically unlinked centromeric regions on all three chromosomes of these closely related taxa contain fixed nucleotide differences which have been found in nearly complete linkage disequilibrium in geographic areas of no or low M-S hybridization. Assays diagnostic for SNP and structural differences between M and S forms in the three centromeric regions were applied in samples from the western extreme of their range of sympatry, the only area where high frequencies of putative M/S hybrids have been reported. The results reveal a level of admixture not observed in the rest of the range. In particular, we found: i) heterozygous genotypes at each marker, although at frequencies lower than expected under panmixia; ii) virtually all possible genotypic combinations between markers on different chromosomes, although genetic association was nevertheless detected; iii) discordant M and S genotypes at two X-linked markers near the centromere, suggestive of introgression and inter-locus recombination. These results could be indicative either of a secondary contact zone between M and S, or of the maintenance of ancestral polymorphisms. This issue and the perspectives opened by these results in the study of the M and S incipient speciation process are discussed.     

Mitochondria--tool for taxonomic identification of yeasts from Saccharomyces sensu stricto complex

Mitochondrial genomes of Saccharomyces and close relatives previously used for transplacement of mitochondria to S. cerevisiae were examined. The origins of replication in mitochondrial DNA, the presence of nuclear and mitochondrial polymorphic loci and the ability to produce mitochondrial respiration-deficient mutants were used to reclassify some collection yeasts and to assign others into four separate subgroups. The first included isolates identical to Saccharomyces cerevisiae (S. italicus, S. oviformis, S. chevalieri and S. capensis) which possess 5 or more replication origins. The second group consists of S paradoxus (var douglasii) mitochondrial genome with the equal number of ori sequences but incompatible mitochondria. The third group represents Saccharomyces sensu stricto petite-positive species (S. carlsbergensis, S. heterogenicus, S. uvarum, S. willianus) with 1-2 origins of replication significantly different from S. cerevisiae. In addition, the locus between tRNA(fMet) and tRNA(Pro) is about one-half of the 1400 bp members of S. cerevisiae complex. The last group includes isolates that do not belong to Saccharomyces sensu stricto group as they are petite-negative and devoid of any S. cerevisiae-like replication origins.     

Molecular and metabolic characterization of cold-tolerant alpine soil Pseudomonas sensu stricto

Alpine soils undergo dramatic temporal changes in their microclimatic properties, suggesting that the bacteria there encounter uncommon shifting selection gradients. Pseudomonads constitute important members of the alpine soil community. In order to characterize the alpine Pseudomonas community and to assess the impact of shifting selection on this community, we examined the ability of cold-tolerant Pseudomonas isolates to grow on a variety of carbon sources, and we determined their phylogenetic relationships based on 16S ribosomal DNA sequencing. We found a high prevalence of Pseudomonas in our soil samples, and isolates from these soils exhibited extensive metabolic diversity. In addition, our data revealed that many of our isolates form a unique cold-adapted clade, representatives of which are also found in the Swedish tundra and Antarctica. Our data also show a lack of concordance between the metabolic properties and 16S phylogeny, indicating that the metabolic diversity of these organisms cannot be predicted by phylogeny.     

Characterization ofYersinia enterocolitica sensu stricto

The speciesYersinia enterocolitica is definedsensu stricto on the bases of biochemical and other phenotypic characteristics. Biochemically,Y. enterocolitica contains five major biotypes: 1 through 4 of Niléhn and of Wauters, and the trehalose-negative, metabolically inactive, socalled hare strains in biotype 5 of Niléhn and of Wauters, and biochemically atypical strains, including urease-negative, Simmons' citrate-positive, and lactose-and raffinose-positive strains.Y. enterocolitica sensu stricto was distinguishable from the newly described speciesYersinia kristensenii by sucrose and Voges-Proskauer reactions (negative inY. kristensenii). These species were previously separated by DNA relatedness.Y. enterocolitica was also separable biochemically and by DNA relatedness from the two newly proposed rhamnose-positive species,Yersinia intermedia andYersinia frederiksenii. Strain 161(=CIP 80-27=ATCC 9610) is proposed as the neotype forY. enterocolitica.     

Mapping the ranges and relative abundance of the two principal African malaria vectors,Anopheles gambiae sensu stricto and An. arabiensis,using climate data.

Members of the Anopheles gambiae complex are major malaria vectors in Africa. We tested the hypothesis that the range and relative abundance of the two major vectors in the complex, An. gambiae sensu stricto and An. arabiensis, could be defined by climate. Climate was characterized at mosquito survey sites by extracting data for each location from climate surfaces using a Geographical Information System. Annual precipitation, together with annual and wet season temperature, defined the ranges of both vectors and were used to map suitable climate zones. Using data from West Africa, we found that where the species were sympatric, An. gambiae s.s. predominated in saturated environments, and An. arabiensis was more common in sites subject to desiccation (r2 = 0.875, p < 0.001). We used the nonlinear equation that best described this relationship to map habitat suitability across Africa. This simple model predicted accurately the relative abundance of both vectors in Tanzania (rs = 0.745, p = 0.002), where species composition is highly variable. The combined maps of species'' range and relative abundance showed very good agreement with published maps. This technique represents a new approach to mapping the distribution of malaria vectors over large areas and may facilitate species-specific vector control activities.     

Molecular analysis of photic inhibition of blood-feeding in Anopheles gambiae

Background

Anopheles gambiae mosquitoes exhibit an endophilic, nocturnal blood feeding behavior. Despite the importance of light as a regulator of malaria transmission, our knowledge on the molecular interactions between environmental cues, the circadian oscillators and the host seeking and feeding systems of the Anopheles mosquitoes is limited.

Results

In the present study, we show that the blood feeding behavior of mosquitoes is under circadian control and can be modulated by light pulses, both in a clock dependent and in an independent manner. Short light pulses (~2–5 min) in the dark phase can inhibit the blood-feeding propensity of mosquitoes momentarily in a clock independent manner, while longer durations of light stimulation (~1–2 h) can induce a phase advance in blood-feeding propensity in a clock dependent manner. The temporary feeding inhibition after short light pulses may reflect a masking effect of light, an unknown mechanism which is known to superimpose on the true circadian regulation. Nonetheless, the shorter light pulses resulted in the differential regulation of a variety of genes including those implicated in the circadian control, suggesting that light induced masking effects also involve clock components. Light pulses (both short and long) also regulated genes implicated in feeding as well as different physiological processes like metabolism, transport, immunity and protease digestions. RNAi-mediated gene silencing assays of the light pulse regulated circadian factors timeless, cryptochrome and three takeout homologues significantly up-regulated the mosquito's blood-feeding propensity. In contrast, gene silencing of light pulse regulated olfactory factors down-regulated the mosquito's propensity to feed on blood.

Conclusion

Our study show that the mosquito's feeding behavior is under circadian control. Long and short light pulses can induce inhibition of blood-feeding through circadian and unknown mechanisms, respectively, that involve the chemosensory system.     

Genetic identification of naturalSaccharomyces sensu stricto yeasts from Finland,Holland and Slovakia

Genetic and karyotypic studies of naturalSaccharomyces sensu stricto yeasts from Finland, Holland and Slovakia revealed three wild sibling-species:Saccharomyces cerevisiae, Saccaromyces bayanus andSaccharomyces paradoxus.     

Gene flow-dependent genomic divergence between Anopheles gambiae M and S forms

Anopheles gambiae sensu stricto exists as two often-sympatric races termed the M and S molecular forms, characterized by fixed differences at an X-linked marker. Extreme divergence between M and S forms at pericentromeric "genomic islands" suggested that selection on variants therein could be driving interform divergence in the presence of ongoing gene flow, but recent work has detected much more widespread genomic differentiation. Whether such genomic islands are important in reproductive isolation or represent ancestral differentiation preserved by low recombination is currently unclear. A critical test of these competing hypotheses could be provided by comparing genomic divergence when rates of recent introgression vary. We genotyped 871 single nucleotide polymorphisms (SNPs) in A. gambiae sensu stricto from locations of M and S sympatry and allopatry, encompassing the full range of observed hybridization rates (0-25%). M and S forms were readily partitioned based on genomewide SNP variation in spite of evidence for ongoing introgression that qualitatively reflects hybridization rates. Yet both the level and the heterogeneity of genomic divergence varied markedly in line with levels of introgression. A few genomic regions of differentiation between M and S were common to each sampling location, the most pronounced being two centromere-proximal speciation islands identified previously but with at least one additional region outside of areas expected to exhibit reduced recombination. Our results demonstrate that extreme divergence at genomic islands does not simply represent segregating ancestral polymorphism in regions of low recombination and can be resilient to substantial gene flow. This highlights the potential for islands comprising a relatively small fraction of the genome to play an important role in early-stage speciation when reproductive isolation is limited.