Multigene Phylogenetics Reveals Temporal Diversification of Major African Malaria Vectors

The major vectors of malaria in sub-Saharan Africa belong to subgenus Cellia. Yet, phylogenetic relationships and temporal diversification among African mosquito species have not been unambiguously determined. Knowledge about vector evolutionary history is crucial for correct interpretation of genetic changes identified through comparative genomics analyses. In this study, we estimated a molecular phylogeny using 49 gene sequences for the African malaria vectors An. gambiae, An. funestus, An. nili, the Asian malaria mosquito An. stephensi, and the outgroup species Culex quinquefasciatus and Aedes aegypti. To infer the phylogeny, we identified orthologous sequences uniformly distributed approximately every 5 Mb in the five chromosomal arms. The sequences were aligned and the phylogenetic trees were inferred using maximum likelihood and neighbor-joining methods. Bayesian molecular dating using a relaxed log normal model was used to infer divergence times. Trees from individual genes agreed with each other, placing An. nili as a basal clade that diversified from the studied malaria mosquito species 47.6 million years ago (mya). Other African malaria vectors originated more recently, and independently acquired traits related to vectorial capacity. The lineage leading to An. gambiae diverged 30.4 mya, while the African vector An. funestus and the Asian vector An. stephensi were the most closely related sister taxa that split 20.8 mya. These results were supported by consistently high bootstrap values in concatenated phylogenetic trees generated individually for each chromosomal arm. Genome-wide multigene phylogenetic analysis is a useful approach for discerning historic relationships among malaria vectors, providing a framework for the correct interpretation of genomic changes across species, and comprehending the evolutionary origins of this ubiquitous and deadly insect-borne disease.     

Mosaic genome architecture of the Anopheles gambiae species complex

Background

Attempts over the last three decades to reconstruct the phylogenetic history of the Anopheles gambiae species complex have been important for developing better strategies to control malaria transmission.

Methodology

We used fingerprint genotyping data from 414 field-collected female mosquitoes at 42 microsatellite loci to infer the evolutionary relationships of four species in the A. gambiae complex, the two major malaria vectors A. gambiae sensu stricto (A. gambiae s.s.) and A. arabiensis, as well as two minor vectors, A. merus and A. melas.

Principal Findings

We identify six taxonomic units, including a clear separation of West and East Africa A. gambiae s.s. S molecular forms. We show that the phylogenetic relationships vary widely between different genomic regions, thus demonstrating the mosaic nature of the genome of these species. The two major malaria vectors are closely related and closer to A. merus than to A. melas at the genome-wide level, which is also true if only autosomes are considered. However, within the Xag inversion region of the X chromosome, the M and two S molecular forms are most similar to A. merus. Near the X centromere, outside the Xag region, the two S forms are highly dissimilar to the other taxa. Furthermore, our data suggest that the centromeric region of chromosome 3 is a strong discriminator between the major and minor malaria vectors.

Conclusions

Although further studies are needed to elucidate the basis of the phylogenetic variation among the different regions of the genome, the preponderance of sympatric admixtures among taxa strongly favor introgression of different genomic regions between species, rather than lineage sorting of ancestral polymorphism, as a possible mechanism.     

Chromosome Inversions,Genomic Differentiation and Speciation in the African Malaria Mosquito Anopheles gambiae

The African malaria vector, Anopheles gambiae, is characterized by multiple polymorphic chromosomal inversions and has become widely studied as a system for exploring models of speciation. Near complete reproductive isolation between different inversion types, known as chromosomal forms, has led to the suggestion that A. gambiae is in early stages of speciation, with divergence evolving in the face of considerable gene flow. We compared the standard chromosomal arrangement (Savanna form) with genomes homozygous for j, b, c, and u inversions (Bamako form) in order to identify regions of genomic divergence with respect to inversion polymorphism. We found levels of divergence between the two sub-taxa within some of these inversions (2Rj and 2Rb), but at a level lower than expected and confined near the inversion breakpoints, consistent with a gene flux model. Unexpectedly, we found that the majority of diverged regions were located on the X chromosome, which contained half of all significantly diverged regions, with much of this divergence located within exons. This is surprising given that the Bamako and Savanna chromosomal forms are both within the S molecular form that is defined by a locus near centromere of X chromosome. Two X-linked genes (a heat shock protein and P450 encoding genes) involved in reproductive isolation between the M and S molecular forms of A. gambiae were also significantly diverged between the two chromosomal forms. These results suggest that genes mediating reproductive isolation are likely located on the X chromosome, as is thought to be the case for the M and S molecular forms. We conclude that genes located on the sex chromosome may be the major force driving speciation between these chromosomal forms of A. gambiae.     

DNA probes for species identification of mosquitoes in the Anopheles gambiae complex

Abstract. Identification of species within the Anopheles gambiae Giles species complex is essential for the correct evaluation of malaria vector ecology studies and control programmes. The development of DNA probes to distinguish species of the An.gambiae complex is described. Genomic libraries were prepared for four members of the An.gambiae complex. These were screened using radiolabeled DNA from different species of An. gambiae sensu lato and a number of clones selected on the basis of their species specificity. These clones could be divided into two groups, each containing homologous sequences. Sequences homologous to group 1 inserts are highly reiterated in the genomes of Anopheles arabiensis Patton and Anopheles merus Dönitz, present in low copy number in Anopheles melas Theobald, but were not detected in Anopheles gambiae sensu stricto. Studies on the organization of this sequence in the genome of An.arabiensis show that homologous sequences are male specific and interspersed within the chromatin. Sequences homologous to group 2 inserts are highly repeated in the genomes of An.merus and An.melas, but present in low copy number in An.gambiae s.s. and An.arabiensis. Group 2 homologous sequences are not sex-specific in the species tested and appear to be tandemly repeated. When used as hybridization probes, these sequences provide a sensitive means for the identification of species within the Anopheles gambiae complex.     

Localization of candidate regions maintaining a common polymorphic inversion (2La) in Anopheles gambiae

Chromosomal inversion polymorphisms are thought to play a role in adaptive divergence, but the genes conferring adaptive benefits remain elusive. Here we study 2La, a common polymorphic inversion in the African malaria vector Anopheles gambiae. The frequency of 2La varies clinally and seasonally in a pattern suggesting response to selection for aridity tolerance. By hybridizing genomic DNA from individual mosquitoes to oligonucleotide microarrays, we obtained a complete map of differentiation across the A. gambiae genome. Comparing mosquitoes homozygous for the 2La gene arrangement or its alternative (2L+^a), divergence was highest at loci within the rearranged region. In the 22 Mb included within alternative arrangements, two ~1.5 Mb regions near but not adjacent to the breakpoints were identified as being significantly diverged, a conclusion validated by targeted sequencing. The persistent association of both regions with the 2La arrangement is highly unlikely given known recombination rates across the inversion in 2La heterozygotes, thus implicating selection on genes underlying these regions as factors responsible for the maintenance of 2La. Polymorphism and divergence data are consistent with a model in which the inversion is maintained by migration-selection balance between multiple alleles inside these regions, but further experiments will be needed to fully distinguish between the epistasis (coadaptation) and local adaptation models for the maintenance of 2La.     

Genome-wide QTL mapping of saltwater tolerance in sibling species of Anopheles (malaria vector) mosquitoes

Although freshwater (FW) is the ancestral habitat for larval mosquitoes, multiple species independently evolved the ability to survive in saltwater (SW). Here, we use quantitative trait locus (QTL) mapping to investigate the genetic architecture of osmoregulation in Anopheles mosquitoes, vectors of human malaria. We analyzed 1134 backcross progeny from a cross between the obligate FW species An. coluzzii, and its closely related euryhaline sibling species An. merus. Tests of 2387 markers with Bayesian interval mapping and machine learning (random forests) yielded six genomic regions associated with SW tolerance. Overlap in QTL regions from both approaches enhances confidence in QTL identification. Evidence exists for synergistic as well as disruptive epistasis among loci. Intriguingly, one QTL region containing ion transporters spans the 2Rop chromosomal inversion that distinguishes these species. Rather than a simple trait controlled by one or a few loci, our data are most consistent with a complex, polygenic mode of inheritance.     

Shifts in malaria vector species composition and transmission dynamics along the Kenyan coast over the past 20 years

Background

Over the past 20 years, numerous studies have investigated the ecology and behaviour of malaria vectors and Plasmodium falciparum malaria transmission on the coast of Kenya. Substantial progress has been made to control vector populations and reduce high malaria prevalence and severe disease. The goal of this paper was to examine trends over the past 20 years in Anopheles species composition, density, blood-feeding behaviour, and P. falciparum sporozoite transmission along the coast of Kenya.

Methods

Using data collected from 1990 to 2010, vector density, species composition, blood-feeding patterns, and malaria transmission intensity was examined along the Kenyan coast. Mosquitoes were identified to species, based on morphological characteristics and DNA extracted from Anopheles gambiae for amplification. Using negative binomial generalized estimating equations, mosquito abundance over the period were modelled while adjusting for season. A multiple logistic regression model was used to analyse the sporozoite rates.

Results

Results show that in some areas along the Kenyan coast, Anopheles arabiensis and Anopheles merus have replaced An. gambiae sensu stricto (s.s.) and Anopheles funestus as the major mosquito species. Further, there has been a shift from human to animal feeding for both An. gambiae sensu lato (s.l.) (99% to 16%) and An. funestus (100% to 3%), and P. falciparum sporozoite rates have significantly declined over the last 20 years, with the lowest sporozoite rates being observed in 2007 (0.19%) and 2008 (0.34%). There has been, on average, a significant reduction in the abundance of An. gambiae s.l. over the years (IRR?=?0.94, 95% CI 0.90–0.98), with the density standing at low levels of an average 0.006 mosquitoes/house in the year 2010.

Conclusion

Reductions in the densities of the major malaria vectors and a shift from human to animal feeding have contributed to the decreased burden of malaria along the Kenyan coast. Vector species composition remains heterogeneous but in many areas An. arabiensis has replaced An. gambiae as the major malaria vector. This has important implications for malaria epidemiology and control given that this vector predominately rests and feeds on humans outdoors. Strategies for vector control need to continue focusing on tools for protecting residents inside houses but additionally employ outdoor control tools because these are essential for further reducing the levels of malaria transmission.     

Transcriptomic differences between euryhaline and stenohaline malaria vector sibling species in response to salinity stress

Evolution of osmoregulatory systems is a key factor in the transition of species between fresh‐ and saltwater habitats. Anopheles coluzzii and Anopheles merus are stenohaline and euryhaline malaria vector mosquitoes belonging to a larger group of sibling species, the Anopheles gambiae complex, which radiated in Africa within the last 2 million years. Comparative ecological genomics of these vector species can provide insight into the mechanisms that permitted the rapid radiation of this species complex into habitats of contrasting salinity. Here, we use RNA‐Seq to investigate gene expression differences between An. coluzzii and An. merus after briefly exposing both young and old larval instars of each species to either saltwater (SW) or freshwater (FW). Our study aims to identify candidate genes and pathways responsible for the greater SW tolerance of An. merus. Our results are congruent with the ability of gene induction to mediate salinity tolerance, with both species showing increasing amounts of differential gene expression between SW and FW as salt concentrations increase. Besides ion transporters such as AgAE2 that may serve as effectors for osmoregulation, we also find mitogen‐activated protein kinases that may serve in a phosphorylation signalling pathway responding to salinity, and report potential cross‐talk between the mosquito immune response and osmoregulation. This study provides a key step towards applying the growing molecular knowledge of these malaria vectors to improve understanding of their ecological tolerances and habitat occupancy.     

Prominent intraspecific genetic divergence within Anopheles gambiae sibling species triggered by habitat discontinuities across a riverine landscape

The Anopheles gambiae complex of mosquitoes includes malaria vectors at different stages of speciation, whose study enables a better understanding of how adaptation to divergent environmental conditions leads to evolution of reproductive isolation. We investigated the population genetic structure of closely related sympatric taxa that have recently been proposed as separate species (An. coluzzii and An. gambiae), sampled from diverse habitats along the Gambia river in West Africa. We characterized putatively neutral microsatellite loci as well as chromosomal inversion polymorphisms known to be associated with ecological adaptation. The results revealed strong ecologically associated population subdivisions within both species. Microsatellite loci on chromosome‐3L revealed clear differentiation between coastal and inland populations, which in An. coluzzii is reinforced by a unusual inversion polymorphism pattern, supporting the hypothesis of genetic divergence driven by adaptation to the coastal habitat. A strong reduction of gene flow was observed between An. gambiae populations west and east of an extensively rice‐cultivated region apparently colonized exclusively by An. coluzzii. Notably, this ‘intraspecific’ differentiation is higher than that observed between the two species and involves also the centromeric region of chromosome‐X which has previously been considered a marker of speciation within this complex, possibly suggesting that the two populations may be at an advanced stage of differentiation triggered by human‐made habitat fragmentation. These results confirm ongoing ecological speciation within these most important Afro‐tropical malaria vectors and raise new questions on the possible effect of this process in malaria transmission.     

Arm-specific dynamics of chromosome evolution in malaria mosquitoes

Background  

The malaria mosquito species of subgenus Cellia have rich inversion polymorphisms that correlate with environmental variables. Polymorphic inversions tend to cluster on the chromosomal arms 2R and 2L but not on X, 3R and 3L in Anopheles gambiae and homologous arms in other species. However, it is unknown whether polymorphic inversions on homologous chromosomal arms of distantly related species from subgenus Cellia nonrandomly share similar sets of genes. It is also unclear if the evolutionary breakage of inversion-poor chromosomal arms is under constraints.     

Spatially explicit analyses of anopheline mosquitoes indoor resting density: implications for malaria control

Background

The question of sampling and spatial aggregation of malaria vectors is central to vector control efforts and estimates of transmission. Spatial patterns of anopheline populations are complex because mosquitoes'' habitats and behaviors are strongly heterogeneous. Analyses of spatially referenced counts provide a powerful approach to delineate complex distribution patterns, and contributions of these methods in the study and control of malaria vectors must be carefully evaluated.

Methodology/Principal Findings

We used correlograms, directional variograms, Local Indicators of Spatial Association (LISA) and the Spatial Analysis by Distance IndicEs (SADIE) to examine spatial patterns of Indoor Resting Densities (IRD) in two dominant malaria vectors sampled with a 5×5 km grid over a 2500 km² area in the forest domain of Cameroon. SADIE analyses revealed that the distribution of Anopheles gambiae was different from regular or random, whereas there was no evidence of spatial pattern in Anopheles funestus (Ia = 1.644, Pa<0.05 and Ia = 1.464, Pa>0.05, respectively). Correlograms and variograms showed significant spatial autocorrelations at small distance lags, and indicated the presence of large clusters of similar values of abundance in An. gambiae while An. funestus was characterized by smaller clusters. The examination of spatial patterns at a finer spatial scale with SADIE and LISA identified several patches of higher than average IRD (hot spots) and clusters of lower than average IRD (cold spots) for the two species. Significant changes occurred in the overall spatial pattern, spatial trends and clusters when IRDs were aggregated at the house level rather than the locality level. All spatial analyses unveiled scale-dependent patterns that could not be identified by traditional aggregation indices.

Conclusions/Significance

Our study illustrates the importance of spatial analyses in unraveling the complex spatial patterns of malaria vectors, and highlights the potential contributions of these methods in malaria control.     

Endogenous siRNAs and piRNAs derived from transposable elements and genes in the malaria vector mosquito Anopheles gambiae

Background

The siRNA and piRNA pathways have been shown in insects to be essential for regulation of gene expression and defence against exogenous and endogenous genetic elements (viruses and transposable elements). The vast majority of endogenous small RNAs produced by the siRNA and piRNA pathways originate from repetitive or transposable elements (TE). In D. melanogaster, TE-derived endogenous siRNAs and piRNAs are involved in genome surveillance and maintenance of genome integrity. In the medically relevant malaria mosquito Anopheles gambiae TEs constitute 12-16% of the genome size. Genetic variations induced by TE activities are known to shape the genome landscape and to alter the fitness in An. gambiae.

Results

Here, using bioinformatics approaches we analyzed the small RNA data sets from 6 libraries formally reported in a previous study and examined the expression of the mixed germline/somatic siRNAs and piRNAs produced in adult An. gambiae females. We characterized a large population of TE-derived endogenous siRNAs and piRNAs, which constitutes 56-60% of the total siRNA and piRNA reads in the analysed libraries. Moreover, we identified a number of protein coding genes producing gene-specific siRNAs and piRNAs that were generally expressed at much lower levels than the TE-associated small RNAs. Detailed sequence analysis revealed that An. gambiae piRNAs were produced by both “ping-pong” dependent (TE-associated piRNAs) and independent mechanisms (genic piRNAs). Similarly to D. melanogaster, more than 90% of the detected piRNAs were produced from TE-associated clusters in An. gambiae. We also found that biotic stress as blood feeding and infection with Plasmodium parasite, the etiological agent of malaria, modulated the expression levels of the endogenous siRNAs and piRNAs in An. gambiae.

Conclusions

We identified a large and diverse set of the endogenously derived siRNAs and piRNAs that share common and distinct aspects of small RNA expression across insect species, and inferred their impact on TE and gene activity in An. gambiae. The TE-specific small RNAs produced by both the siRNA and piRNA pathways represent an important aspect of genome stability and genetic variation, which might have a strong impact on the evolution of the genome and vector competence in the malaria mosquitoes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1436-1) contains supplementary material, which is available to authorized users.     

Optimized Pan-species and Speciation Duplex Real-time PCR Assays for Plasmodium Parasites Detection in Malaria Vectors

Background

An accurate method for detecting malaria parasites in the mosquito’s vector remains an essential component in the vector control. The Enzyme linked immunosorbent assay specific for circumsporozoite protein (ELISA-CSP) is the gold standard method for the detection of malaria parasites in the vector even if it presents some limitations. Here, we optimized multiplex real-time PCR assays to accurately detect minor populations in mixed infection with multiple Plasmodium species in the African malaria vectors Anopheles gambiae and Anopheles funestus.

Methods

Complementary TaqMan-based real-time PCR assays that detect Plasmodium species using specific primers and probes were first evaluated on artificial mixtures of different targets inserted in plasmid constructs. The assays were further validated in comparison with the ELISA-CSP on 200 field caught Anopheles gambiae and Anopheles funestus mosquitoes collected in two localities in southern Benin.

Results

The validation of the duplex real-time PCR assays on the plasmid mixtures demonstrated robust specificity and sensitivity for detecting distinct targets. Using a panel of mosquito specimen, the real-time PCR showed a relatively high sensitivity (88.6%) and specificity (98%), compared to ELISA-CSP as the referent standard. The agreement between both methods was “excellent” (κ = 0.8, P<0.05). The relative quantification of Plasmodium DNA between the two Anopheles species analyzed showed no significant difference (P = 0, 2). All infected mosquito samples contained Plasmodium falciparum DNA and mixed infections with P. malariae and/or P. ovale were observed in 18.6% and 13.6% of An. gambiae and An. funestus respectively. Plasmodium vivax was found in none of the mosquito samples analyzed.

Conclusion

This study presents an optimized method for detecting the four Plasmodium species in the African malaria vectors. The study highlights substantial discordance with traditional ELISA-CSP pointing out the utility of employing an accurate molecular diagnostic tool for detecting malaria parasites in field mosquito populations.     

Environmental Factors Associated with the Distribution of Anopheles gambiae s.s in Ghana; an Important Vector of Lymphatic Filariasis and Malaria

Anopheles gambiae s.s mosquitoes are important vectors of lymphatic filariasis (LF) and malaria in Ghana. To better understand their ecological aspects and influence on disease transmission, we examined the spatial distribution of the An. gambiae (M and S) molecular forms and associated environmental factors, and determined their relationship with disease prevalence. Published and current data available on the An. gambiae species in Ghana were collected in a database for analysis, and the study sites were georeferenced and mapped. Using the An. gambiae s.s sites, environmental data were derived from climate, vegetation and remote-sensed satellite sources, and disease prevalence data from existing LF and malaria maps in the literature. The data showed that An. gambiae M and S forms were sympatric in most locations. However, the S form predominated in the central region, while the M form predominated in the northern and coastal savanna regions. Bivariate and multiple regression analyses identified temperature as a key factor distinguishing their distributions. An. gambiae M was significantly correlated with LF, and 2.5 to 3 times more prevalent in the high LF zone than low to medium zones. There were no significant associations between high prevalence An. gambiae s.s locations and malaria. The distribution of the An. gambiae M and S forms and the diseases they transmit in Ghana appear to be distinct, driven by different environmental factors. This study provides useful baseline information for disease control, and future work on the An. gambiae s.s in Ghana.     

The Population Genomics of Trans-Specific Inversion Polymorphisms in Anopheles gambiae

In the malaria mosquito Anopheles gambiae polymorphic chromosomal inversions may play an important role in adaptation to environmental variation. Recently, we used microarray-based divergence mapping combined with targeted resequencing to map nucleotide differentiation between alternative arrangements of the 2La inversion. Here, we applied the same technique to four different polymorphic inversions on the 2R chromosome of An. gambiae. Surprisingly, divergence was much lower between alternative arrangements for all 2R inversions when compared to the 2La inversion. For one of the rearrangements, 2Ru, we successfully mapped a very small region (∼100 kb) of elevated divergence. For the other three rearrangements, we did not identify any regions of significantly high divergence, despite ample independent evidence from natural populations of geographic clines and seasonal cycling, and stable heterotic polymorphisms in laboratory populations. If these inversions are the targets of selection as hypothesized, we suggest that divergence between rearrangements may have escaped detection due to retained ancestral polymorphism in the case of the youngest 2R rearrangements and to extensive gene flux in the older 2R inversion systems that segregate in both An. gambiae and its sibling species An. arabiensis.MORE than 70 years ago Dobzhansky and Sturtevant (1938) first discovered polymorphic inversion arrangements carried by various Drosophila pseudoobscura populations. After observing correlations between environmental conditions and inversion frequencies, Dobzhansky proposed that inversions are under strong selection due to their role in promoting local adaptation to the heterogeneous conditions a species encounters both spatially and temporally (Dobzhansky 1944, 1948; Powell 1997). More recent studies have implicated chromosomal inversions in the adaptation of a diversity of eukaryotes including humans (Coluzzi et al. 1979; Feder et al. 2003; Hoffmann et al. 2004; Stefansson et al. 2005). Long known to be common in dipteran insects, more recent HapMap data suggest that polymorphic inversions may be numerous in human populations and by extension other mammals (Bansal et al. 2007). Given their potential importance in facilitating adaptation, surprisingly little is known about the mechanism(s) or the genes responsible for maintaining inversion polymorphisms in natural populations.Gene exchange between inverted and standard arrangements, although reduced, can still occur through gene flux: the action of gene conversion and multiple crossovers in inversion heterozygotes (heterokaryotypes) (Chovnick 1973; Navarro et al. 1997; Schaeffer and Anderson 2005). Over time allelic variation unrelated to ecological adaptation should become homogenized between arrangements, while alleles which are under divergent selection pressures should remain in linkage disequilibrium with each other and with the inversion itself, leading to heightened differentiation between standard and inverted arrangements at and near the target loci. In principle, this process allows the identification of specific loci involved in adaptive divergence (Schaeffer et al. 2003; Schaeffer and Anderson 2005; Storz 2005). Consistent with this model, previous low-resolution studies of Drosophila inversions revealed heterogeneous patterns of nucleotide diversity relative to divergence, as well as the interspersion of regions of high and low genetic association potentially due to the interaction of selection and gene flux (Schaeffer et al. 2003; Kennington et al. 2006; but see Munte et al. 2005). The application of high-resolution tools flowing from completely sequenced genomes will facilitate the mapping of genes that are the targets of divergent natural selection within gene arrangements.Although Drosophila has been the favored model, the African malaria vector Anopheles gambiae sensu stricto also provides an excellent system for studying the maintenance of inversion polymorphisms, not only within a species but across speciation events of different ages in the An. gambiae sibling species complex. The nominal species An. gambiae s.s. (hereafter, An. gambiae) is synanthropic: almost exclusively biting humans, resting indoors, and exploiting anthropogenic larval habitats (Coluzzi 1999). This close association with humans, vital to making An. gambiae one of the most proficient vectors of malaria, is likely to have been facilitated by chromosomal inversions thought to confer adaptive benefits in heterogeneous climatic and ecological settings in Africa. Seven common polymorphic inversions exist on the second chromosome. Six of these are located on the right arm (2R): j, b, c, u, d, and k, while 2La is the only inversion on the left arm (Coluzzi et al. 2002). Facilitated by the sequenced reference genome (Holt et al. 2002), some of the breakpoints for these polymorphic inversions have been localized to small genomic regions (Sharakhov et al. 2006; Coulibaly et al. 2007; Sangare 2007). Most of these inversions appear to be the targets of strong selection. Five of the inversions (2La and 2Rb, -c, -d, and -u) are nonrandomly associated with degree of aridity; each cycles seasonally with rainfall, and all except 2Ru form stable geographic clines in frequency from mesic forest to xeric regions bordering the Sahara (Coluzzi et al. 1979; Toure et al. 1994, 1998; Powell et al. 1999). Inversion 2Rj is not clinal, but its distribution in Mali is consistent with adaptation to novel rockpool niches (Coluzzi et al. 1985; Manoukis et al. 2008).In the An. gambiae species complex, inversion polymorphisms can be maintained across the boundaries of emerging and even full species. An. gambiae and its sibling An. arabiensis, strictly sympatric throughout most of their extensive ranges in sub-Saharan Africa, differ by multiple fixed chromosomal rearrangements on the X but share three chromosome 2 inversions: 2La, fixed in An. arabiensis and polymorphic in An. gambiae; and 2Rb and -c, polymorphic in both species (Coluzzi et al. 1979, 2002). Moreover, these same inversions and all other common An. gambiae inversions with the exception of 2Rj are shared and polymorphic in two lineages apparently undergoing ecological speciation within An. gambiae—the assortatively mating M and S molecular forms (della Torre et al. 2002, 2005). Inversion frequencies are correlated with climatic and ecological conditions in parallel in both lineages (Costantini et al. 2009; Simard et al. 2009). Unlike the full species, the M and S incipient species are not distinguished by any fixed inversion differences. Indeed, genomewide divergence mapping between the M and S forms revealed that significant differentiation was confined to two small low-recombination regions adjacent to the centromeres of 2L and X which are distant from any inversions (Turner et al. 2005). Thus, in distinction to models of speciation invoking inversions as facilitating the persistence of hybridizing species (Noor et al. 2001; Rieseberg 2001; Ortiz-Barrientos et al. 2002; Navarro and Barton 2003), the An. gambiae data suggest that chromosome 2 inversions are not directly responsible for reproductive isolation. Instead, the same chromosome 2 inversion polymorphisms appear to confer similar ecological benefits, within and across species boundaries. A long-term research goal is to identify the mechanisms and the genes controlling these processes.Previously we conducted the first high-density genomic scan of divergence across a chromosomal inversion (2La) in An. gambiae (White et al. 2007). By hybridizing genomic DNA from S form mosquitoes homokaryotypic for alternate gene arrangements on chromosome 2L (2La or 2L+^a) to oligonucleotide microarrays we were able to measure divergence across the 22-Mb inversion at nearly 14,000 markers. Differentiation in the rearranged region was significantly higher than in collinear portions of chromosome 2L. Between breakpoints the pattern of differentiation was heterogeneous: two genomic clusters of significantly higher divergence were identified near but not adjacent to the breakpoints. Directed resequencing within the S form confirmed these results and suggested that both clusters contained genes targeted by selection. Observed levels of linkage disequilibrium between the 2La breakpoints and markers in the clusters are highly unlikely under a neutral scenario, in light of known recombination rates and plausible estimates of the age of the inversion.The present study characterizes the patterns of genetic variation in polymorphic rearrangements on the opposite (right) arm of chromosome 2: 2Rj, -b, -c, and -u. With the goal of identifying candidate genes maintaining these inversions in natural populations, we applied microarray-based divergence mapping to measure differentiation between alternative 2R arrangements. Because three of four inversions have taxonomic distributions that span incipient and/or completed speciation events, we validated the microarray findings by targeted sequencing in multiple taxa: sympatric Malian populations of An. gambiae M and S forms, and the sibling species An. arabiensis.     

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.     

Molecular organization of heterochromatin in malaria mosquitoes of the Anopheles maculipennis subgroup

Although heterochromatin makes up a significant portion of the malaria mosquito genome, its organization, function, and evolution are poorly understood. Sibling species of the Anopheles maculipennis subgroup, the European malaria mosquitoes, are characterized by striking differences in the morphology of pericentric heterochromatin; however, the molecular basis for the rapid evolutionary transformation of heterochromatin is not known. This study reports an initial survey of the molecular organization of the pericentric heterochromatin in nonmodel species from the A. maculipennis subgroup. Molecular identity and chromosomal localization were established for short DNA fragments obtained by microdissection from the pericentric diffuse β-heterochromatin of A. atroparvus. Among 102 sequenced clones of the Atr2R library, twenty had sequence similarity to transposable elements (TEs) from the Anopheles gambiae and Aedes aegypti genomes. At least six protein-coding single-copy genes from A. gambiae and four single-copy genes from Drosophila melanogaster were homologous to eight clones from the library. Most of these conserved genes were heterochromatic in A. gambiae but euchromatic in D. melanogaster. The remaining 74 clones were characterized as noncoding repetitive DNA. Comparative chromosome mapping of twelve clones in the sibling species A. atroparvus and A. messeae demonstrated that the noncoding repetitive sequences and the TEs have undergone independent chromosome-specific and species-specific gains and losses in the morphologically different pericentric heterochromatic regions, in accordance with the “library model.”     

Select Small Core Structure Carbamates Exhibit High Contact Toxicity to ��Carbamate-Resistant�� Strain Malaria Mosquitoes, Anopheles gambiae (Akron)

Acetylcholinesterase (AChE) is a proven target for control of the malaria mosquito (Anopheles gambiae). Unfortunately, a single amino acid mutation (G119S) in An. gambiae AChE-1 (AgAChE) confers resistance to the AChE inhibitors currently approved by the World Health Organization for indoor residual spraying. In this report, we describe several carbamate inhibitors that potently inhibit G119S AgAChE and that are contact-toxic to carbamate-resistant An. gambiae. PCR-RFLP analysis was used to confirm that carbamate-susceptible G3 and carbamate-resistant Akron strains of An. gambiae carry wild-type (WT) and G119S AChE, respectively. G119S AgAChE was expressed and purified for the first time, and was shown to have only 3% of the turnover number (k _cat) of the WT enzyme. Twelve carbamates were then assayed for inhibition of these enzymes. High resistance ratios (>2,500-fold) were observed for carbamates bearing a benzene ring core, consistent with the carbamate-resistant phenotype of the G119S enzyme. Interestingly, resistance ratios for two oxime methylcarbamates, and for five pyrazol-4-yl methylcarbamates were found to be much lower (4- to 65-fold). The toxicities of these carbamates to live G3 and Akron strain An. gambiae were determined. As expected from the enzyme resistance ratios, carbamates bearing a benzene ring core showed low toxicity to Akron strain An. gambiae (LC₅₀>5,000 μg/mL). However, one oxime methylcarbamate (aldicarb) and five pyrazol-4-yl methylcarbamates (4a–e) showed good to excellent toxicity to the Akron strain (LC₅₀ = 32–650 μg/mL). These results suggest that appropriately functionalized “small-core” carbamates could function as a resistance-breaking anticholinesterase insecticides against the malaria mosquito.