Chromosome arm‐specific patterns of polymorphism associated with chromosomal inversions in the major African malaria vector,Anopheles funestus

Chromosomal inversions facilitate local adaptation of beneficial mutations and modulate genetic polymorphism, but the extent of their effects within the genome is still insufficiently understood. The genome of Anopheles funestus, a malaria mosquito endemic to sub‐Saharan Africa, contains an impressive number of paracentric polymorphic inversions, which are unevenly distributed among chromosomes and provide an excellent framework for investigating the genomic impacts of chromosomal rearrangements. Here, we present results of a fine‐scale analysis of genetic variation within the genome of two weakly differentiated populations of Anopheles funestus inhabiting contrasting moisture conditions in Cameroon. Using population genomic analyses, we found that genetic divergence between the two populations is centred on regions of the genome corresponding to three inversions, which are characterized by high values of F_ST, absolute sequence divergence and fixed differences. Importantly, in contrast to the 2L chromosome arm, which is collinear, nucleotide diversity is significantly reduced along the entire length of three autosome arms bearing multiple overlapping chromosomal rearrangements. These findings support the idea that interactions between reduced recombination and natural selection within inversions contribute to sculpt nucleotide polymorphism across chromosomes in An. funestus.     

Dieldrin resistance in the malaria vector Anopheles gambiae in Ghana

Anopheles gambiae Giles s.s. (Diptera: Culicidae) is one of the principal vectors of malaria in the Ashanti region of central Ghana. High levels of resistance to dieldrin were recorded in a wild-caught sample from Obuasi (south of Kumasi) as well as a laboratory colony established using material from the wild population. Cytogenetic analysis of wild-caught and laboratory samples revealed chromosomal polymorphism for inversions 2La and 2Rb. Although inversion 2La has previously been shown to be associated with dieldrin resistance in certain other laboratory strains originating from West Africa, there was no obvious association between inversion karyotype assortment and the resistance phenotype in the Obuasi population. In addition, polymerase chain reaction analysis indicated the presence of the alanine296 to glycine mutation in the GABA (gamma amino-butyric acid) receptor (which has been mapped to a chromosomal position within inversion 2La). This mutation has previously been shown to be associated with dieldrin resistance in the same An. gambiae laboratory strains of West African origin. Our data show only a weak association between the dieldrin resistance phenotype and the presence of this mutation, suggesting that another dieldrin resistance mechanism is operational in the Obuasi population. Biochemical and synergist exposure assays suggest a metabolic component, probably mediated by monooxygenase P450 enzymes. We conclude that dieldrin resistance in the An. gambiae population of the Obuasi region occurs at a high level - most likely in the absence of selection - and that control of the resistance phenotype is polyfactorial and must include components other than mutations in the GABA receptor locus.     

Dynamics of gene introgression in the African malaria vector Anopheles gambiae

Anopheles gambiae is a major malaria vector in Africa and a popular model species for a variety of ecological, evolutionary, and genetic studies on vector control. Genetic manipulation of mosquito vectorial capacity is a promising new weapon for the control of malaria. However, the release of exotic transgenic mosquitoes will bring in novel alleles in addition to the parasite-inhibiting genes, which may have unknown effects on the local population. Therefore, it is necessary to develop methodologies that can be used to evaluate the spread rate of introduced genes in A. gambiae. In this study, the effects and dynamics of genetic introgression between two geographically distinct A. gambiae populations from western Kenya (Mbita) and eastern Tanzania (Ifakara) were investigated with amplified fragment length polymorphisms (AFLPs) and microsatellite markers. Microsatellites and polymorphic cDNA markers revealed a large genetic differentiation between the two populations (average F(ST) = 0.093, P < 0.001). When the two strains were crossed in random mating between the two populations, significant differences in the rate of genetic introgression were found in the mixed populations. Allele frequencies of 18 AFLP markers (64.3%) for Mbita and of 26 markers (92.9%) for Ifakara varied significantly from F5 to F20. This study provides basic information on how a mosquito release program would alter the genetic makeup of natural populations, which is critical for pilot field testing and ecological risk evaluation of transgenic mosquitoes.     

Distribution of genetic diversity in relation to chromosomal inversions in the malaria mosquito Anopheles gambiae

The epidemiology of malaria in Africa is complicated by the fact that its principal vector, the mosquito Anopheles gambiae, constitutes a complex of six sibling species. Each species is characterized by a unique array of paracentric inversions, as deduced by karyotypic analysis. In addition, most of the species carry a number of polymorphic inversions. In order to develop an understanding of the evolutionary histories of different parts of the genome, we compared the genetic variation of areas inside and outside inversions in two distinct inversion karyotypes of A. gambiae. Thirty-five cDNA clones were mapped on the five arms of the A. gambiae chromosomes with divisional probes. Sixteen of these clones, localized both inside and outside inversions of chromosome 2, were used as probes in order to determine the nucleotide diversity of different parts of the genome in the two inversion karyotypes. We observed that the sequence diversity inside the inversion is more than threefold lower than in areas outside the inversion and that the degree of divergence increases gradually at loci at increasing distance from the inversion. To interpret the data we present a selectionist and a stochastic model, both of which point to a relatively recent origin of the studied inversion and may suggest differences between the evolutionary history of inversions in Anopheles and Drosophila species.Correspondence to: K.D. Mathiopoulos     

Coalescent patterns for chromosomal inversions in divergent populations

Chromosomal inversions allow genetic divergence of locally adapted populations by reducing recombination between chromosomes with different arrangements. Divergence between populations (or hybridization between species) is expected to leave signatures in the neutral genetic diversity of the inverted region. Quantitative expectations for these patterns, however, have not been obtained. Here, we develop coalescent models of neutral sites linked to an inversion polymorphism in two locally adapted populations. We consider two scenarios of local adaptation: selection on the inversion breakpoints and selection on alleles inside the inversion. We find that ancient inversion polymorphisms cause genetic diversity to depart dramatically from neutral expectations. Other situations, however, lead to patterns that may be difficult to detect; important determinants are the age of the inversion and the rate of gene flux between arrangements. We also study inversions under genetic drift, finding that they produce patterns similar to locally adapted inversions of intermediate age. Our results are consistent with empirical observations, and provide the foundation for quantitative analyses of the roles that inversions have played in speciation.     

Indirect evidence that agricultural pesticides select for insecticide resistance in the malaria vector Anopheles gambiae

We investigated the possible relationship between the agricultural use of insecticides and the emergence of insecticide resistance. Bioassays were conducted using simulated mosquito larval habitats and well known Anopheles gambiae strains. Soil samples were collected from vegetable production areas in Benin, including one site with insecticide use, one site where insecticides had not been used for two months, and a third where insecticides had not been used. Pupation and emergence rates were very low in pyrethroid‐susceptible strains when exposed to soil that had been recently exposed to insecticides. Pupation and emergence rates in strains with the kdr mutation alone or both the kdr and Ace‐1 mutations were much higher. Overall, strains with the kdr mutation survived at higher rates compared to that without kdr mutation. Although this study is observational, we provide indirect evidence indicating that soils from agricultural areas contain insecticide residues that can play a role in the emergence of insecticide resistance in Anopheles. This aspect should be taken into account to better utilize the insecticide in the context of integrated pest management programs.     

Adaptation to aridity in the malaria mosquito Anopheles gambiae: chromosomal inversion polymorphism and body size influence resistance to desiccation

Chromosomal inversions are thought to confer a selective advantage in alternative habitats by protecting co-adapted alleles from recombination. The frequencies of two inversions (2La and 2Rb) of the afro-tropical malaria mosquito Anopheles gambiae change gradually along geographical clines, increasing in frequency with degree of aridity. Such clines can result from gene flow and local selection acting upon alternative karyotypes along the cline, suggesting that these inversions may be associated with tolerance to xeric conditions. Since water loss represents a major challenge in xeric habitats, it can be supposed that genes inside these inversions are involved in water homeostasis. To test this hypothesis, we compared the desiccation resistance of alternative karyotypes from a colonised 2Rb/2La polymorphic population of A. gambiae from Cameroon. The strain included only the molecular form S, one of the genetic units marking incipient speciation in this taxon. Day-old mosquitoes of both sexes were assayed individually for time to death in a dry environment and the karyotype of each was determined post-mortem using molecular diagnostic assays for each inversion. In agreement with expectations based on their eco-geographical distribution, we found that 2La homokaryotypes survived significantly longer (1.3 hours) than the other karyotypes. However, there was weak support for the effect of 2Rb on desiccation resistance. Larger mosquitoes survived longer than smaller ones. Median survival of females was greater than males, but the effect of sex on desiccation resistance was weakly supported, indicating that differential survival was correlated to differences between sexes in average size. We found weak evidence for a heterotic effect of 2La karyotype on size in females. These results support the notion that genes located inside the 2La inversion are involved in water balance, contributing towards local adaptation of A. gambiae to xeric habitats, beyond the adaptive value conferred by a larger body size.     

Viral paratransgenesis in the malaria vector Anopheles gambiae

Paratransgenesis, the genetic manipulation of insect symbiotic microorganisms, is being considered as a potential method to control vector-borne diseases such as malaria. The feasibility of paratransgenic malaria control has been hampered by the lack of candidate symbiotic microorganisms for the major vector Anopheles gambiae. In other systems, densonucleosis viruses (DNVs) are attractive agents for viral paratransgenesis because they infect important vector insects, can be genetically manipulated and are transmitted to subsequent generations. However, An. gambiae has been shown to be refractory to DNV dissemination. We discovered, cloned and characterized the first known DNV (AgDNV) capable of infection and dissemination in An. gambiae. We developed a flexible AgDNV-based expression vector to express any gene of interest in An. gambiae using a two-plasmid helper-transducer system. To demonstrate proof-of-concept of the viral paratransgenesis strategy, we used this system to transduce expression of an exogenous gene (enhanced green fluorescent protein; EGFP) in An. gambiae mosquitoes. Wild-type and EGFP-transducing AgDNV virions were highly infectious to An. gambiae larvae, disseminated to and expressed EGFP in epidemiologically relevant adult tissues such as midgut, fat body and ovaries and were transmitted to subsequent mosquito generations. These proof-of-principle data suggest that AgDNV could be used as part of a paratransgenic malaria control strategy by transduction of anti-Plasmodium peptides or insect-specific toxins in Anopheles mosquitoes. AgDNV will also be extremely valuable as an effective and easy-to-use laboratory tool for transient gene expression or RNAi in An. gambiae.     

Chromosome inversions and ecological plasticity in the main African malaria mosquitoes

Chromosome inversions have fascinated the scientific community, mainly because of their role in the rapid adaption of different taxa to changing environments. However, the ecological traits linked to chromosome inversions have been poorly studied. Here, we investigated the roles played by 23 chromosome inversions in the adaptation of the four major African malaria mosquitoes to local environments in Africa. We studied their distribution patterns by using spatially explicit modeling and characterized the ecogeographical determinants of each inversion range. We then performed hierarchical clustering and constrained ordination analyses to assess the spatial and ecological similarities among inversions. Our results show that most inversions are environmentally structured, suggesting that they are actively involved in processes of local adaptation. Some inversions exhibited similar geographical patterns and ecological requirements among the four mosquito species, providing evidence for parallel evolution. Conversely, common inversion polymorphisms between sibling species displayed divergent ecological patterns, suggesting that they might have a different adaptive role in each species. These results are in agreement with the finding that chromosomal inversions play a role in Anopheles ecotypic adaptation. This study establishes a strong ecological basis for future genome‐based analyses to elucidate the genetic mechanisms of local adaptation in these four mosquitoes.     

Polymorphic chromosomal inversions in Anopheles moucheti,a major malaria vector in Central Africa

Anopheles moucheti Evans (Diptera: Culicidae) is a major vector of malaria in forested areas of Central Africa. However, few genetic tools are available for this species. The present study represents the first attempt to characterize chromosomes in An. moucheti females collected in Cameroon. Ovarian nurse cells contained polytene chromosomes, which were suitable for standard cytogenetic applications. The presence of three polymorphic chromosomal inversions in An. moucheti was revealed. Two of these inversions were located on the 2R chromosome arm. The homology between the 2R chromosome arms of An. moucheti and Anopheles gambiae Giles was established by fluorescent in situ hybridization of six An. gambiae genic sequences. Mapping of the probes on chromosomes of An. moucheti detected substantial gene order reshuffling between the two species. The presence of polytene chromosomes and polymorphic inversions in An. moucheti provides a new basis for further population genetic, taxonomic and ecological studies of this neglected malaria vector.     

Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes

Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients.     

Identifying targets of selection in mosaic genomes with machine learning: applications in Anopheles gambiae for detecting sites within locally adapted chromosomal inversions

Chromosomal inversions are important structural changes that may facilitate divergent selection when they capture co‐adaptive loci in the face of gene flow. However, identifying selection targets within inversions can be challenging. The high degrees of differentiation between heterokaryotypes, as well as the differences in demographic histories of collinear regions compared with inverted ones, reduce the power of traditional outlier analyses for detecting selected loci. Here, we develop a new approach that uses discriminant functions informed from inversion‐specific expectations to classify loci that are under selection (or drift). Analysis of RAD sequencing data we collected in a classic dipteran species with polymorphic inversion clines—Anopheles gambiae, a malaria vector species from sub‐Saharan Africa—demonstrates the benefits of the approach compared with traditional outlier analyses. We focus specifically on two polymorphic inversions, the 2La and 2Rb arrangements that predominate in dry habitats and the 2L+^a and 2R+^b arrangements in wet habitats, which contrast with the minimal geographic structure of SNPs from collinear regions. With our approach, we identify two strongly selected regions within 2La associated with dry habitat. Moreover, we also show that the prevalence of selection is greater in the arrangement 2L+^a that is associated with wet habitat (unlike presumed importance of selective divergence associated with the shift of the mosquitoes into dry habitats). We discuss the implications of these results with respect to studies of rapid adaptation in these malaria vectors, and in particular, the insights our newly developed approach offers for identifying not only potential targets of selection, but also the population that has undergone adaptive change.     

Salivary gland-specific gene expression in the malaria vector Anopheles gambiae

Molecular studies on the tissue-specific gene expression in the salivary glands of Anopheles gambiae may provide useful tools for the development of new strategies for the control of the most efficient malaria vector in the sub-Saharan Africa. We summarize here the results of a recent investigation focused on the isolation of secreted factors and putative receptors from the salivary glands of An. gambiae. Using the Signal Sequence Trap technique we have identified the first cDNAs specifically expressed in the An. gambiae salivary glands. Among these, four are exclusively expressed in female glands and encode factors presumably involved in blood-feeding, whereas two other cDNAs seem to be expressed both in male and in female glands and are likely implicated in sugar-feeding. Homologues of genes previously identified in the yellow fever mosquito Aedes aegypti, like the apyrase and D7, as well as novel salivary gland-specific cDNAs, were identified. The isolation and characterization of promoter sequences from the corresponding genes may prove useful for the expression of anti parasitic agents in the salivary glands of transgenic mosquitoes.     

New repellent effective against African malaria mosquito Anopheles gambiae: implications for vector control

Anopheles gambiae Giles sensu stricto (Diptera: Culicidae) is a vector for Plasmodium, the causative agent of malaria. Current control strategies to reduce the impact of malaria focus on reducing the frequency of mosquito attacks on humans, thereby decreasing Plasmodium transmission. A need for new repellents effective against Anopheles mosquitoes has arisen because of changes in vector behaviour as a result of control strategies and concern over the health impacts of current repellents. The response of A. gambiae to potential repellents was investigated through an electroantennogram screen and the most promising of these candidates (1‐allyloxy‐4‐propoxybenzene, 3c {3,6}) chosen for behavioural testing. An assay to evaluate the blood‐host seeking behaviour of A. gambiae towards a simulated host protected with this repellent was then performed. The compound 3c {3,6} was shown to be an effective repellent, causing mosquitoes to reduce their contact with a simulated blood‐host and probe less at the host odour. Thus, 3c {3,6} may be an effective repellent for the control of A. gambiae.     

Host-specific cues cause differential attractiveness of Kenyan men to the African malaria vector Anopheles gambiae

Background  

Many studies have suggested that variability in the attractiveness of humans to host-seeking mosquitoes is caused by differences in the make-up of body emanations, and olfactory signals in particular. Most investigations have either been laboratory-based, utilising odour obtained from sections of the body, or have been done in the field with sampling methods that do not discriminate between visual, physical and chemical cues of the host. Accordingly, evidence for differential attractiveness based on body emanations remains sparse in spite of the far-reaching epidemiological implications of this phenomenon.     

Sex-biased expression of odorant receptors in antennae and palps of the African malaria vector Anopheles gambiae

At the heart of the odor recognition process in all animals are G-protein-coupled receptors, which are seven-transmembrane domain proteins that initiate G-protein-mediated signaling cascades when activated by their ligands. Odorant receptors (ORs) are a large, diverse family of proteins with some 80 members in the mosquito Anopheles gambiae. With the assumption that more sensilla on female antennae are tuned to human odors than on male antennae, comparison of specific OR mRNA levels in male and female antennae can provide an indication as to which receptors may be stimulated by host odors. We have used RT PCR and quantitative real-time PCR (qRT PCR) to investigate sex-biased expression levels of 80 A. gambiae ORs in male and female antennae and maxillary palps. On the basis of prevalence of expression in female antennae and on a strong female relative to male expression bias we identified a short list of ORs that are likely involved in host odor recognition by female mosquitoes.     

Cas9-mediated maternal effect and derived resistance alleles in a gene-drive strain of the African malaria vector mosquito,Anopheles gambiae

CRISPR/Cas9 technologies are important tools for the development of gene-drive systems to modify mosquito vector populations to control the transmission of pathogens that cause diseases such as malaria. However, one of the challenges for current Cas9-based drive systems is their ability to produce drive-resistant alleles resulting from insertions and deletions (indels) caused principally by nonhomologous end-joining following chromosome cleavage. Rapid increases in the frequency of such alleles may impair gene-drive dynamics. We explored the generation of indels in the germline and somatic cells in female gene-drive lineages using a series of selective crosses between a gene-drive line, AgNosCd-1, and wild-type mosquitoes. We find that potential drive-resistant mutant alleles are generated largely during embryonic development, most likely caused by deposition of the Cas9 endonuclease and guide RNAs in oocytes and resulting embryos by homozygous and hemizygous gene-drive mothers.