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

Yearly variations of the genetic structure of Aedes aegypti (Linnaeus) (Diptera: Culicidae) in the Philippines (2017-2019)

Dengue is the fastest emerging arboviral disease in the world, imposing a substantial health and economic burden in the tropics and subtropics. The mosquito, Aedes aegypti, is the primary vector of dengue in the Philippines. We examined the genetic structure of Ae. aegypti populations collected from the Philippine major islands (Luzon, Visayas and Mindanao), each with highland (Baguio city, Cebu city mountains and Maramag, Bukidnon, respectively) and lowland sites (Quezon city; Liloan, Cebu and Cagayan de Oro [CDO] city, respectively) during the wet (2017-2018 and 2018-2019) and dry seasons (2018 and 2019). Mosquitoes (n = 1800) were reared from field-collected eggs and immatures, and were analyzed using 12 microsatellite loci. Generalized linear model analyses revealed yearly variations between highlands and lowlands in the major islands as supported by Bayesian clustering analyses on: 1) stronger selection (inbreeding coefficient, F_IS = 0.52) in 2017-2018 than in 2018-2019 (F_IS = 0.32) as influenced by rainfall, 2) the number of non-neutral loci indicating selection, and 3) differences of effective population size although at p = 0.05. Across sites except Baguio and CDO cities: 1) F_IS varied seasonally as influenced by relative humidity (RH), and 2) the number of non-neutral loci varied as influenced by RH and rainfall indicating selection. Human-mediated activities and not isolation by distance influenced genetic differentiations of mosquito populations within (F_ST = 0.04) the major islands and across sites (global F_ST = 0.16). Gene flow (Nm) and potential first generation migrants among populations were observed between lowlands and highlands within and across major islands. Our results suggest that dengue control strategies in the epidemic wet season are to be changed into whole year-round approach, and water pipelines are to be installed in rural mountains to prevent the potential breeding sites of mosquitoes.

Perfect association between spatial swarm segregation and the X-chromosome speciation island in hybridizing Anopheles coluzzii and Anopheles gambiae populations

The sibling species An. coluzzii and An. gambiae s.s. are major malaria vectors thought to be undergoing sympatric speciation with gene flow. In the absence of intrinsic post-zygotic isolation between the two taxa, speciation is thought possible through the association of assortative mating and genomic regions protected from gene flow by recombination suppression. Such genomic islands of speciation have been described in pericentromeric regions of the X, 2L and 3L chromosomes. Spatial swarm segregation plays a major role in assortative mating between sympatric populations of the two species and, given their importance for speciation, genes responsible for such pre-mating reproductive barriers are expected to be protected within divergence islands. In this study 2063 male and 266 female An. coluzzii and An. gambiae s.s. individuals from natural swarms in Burkina Faso, West Africa were sampled. These were genotyped at 16 speciation island SNPs, and characterized as non-hybrid individuals, F₁ hybrids or recombinant F_1+n backcrossed individuals. Their genotypes at each speciation island were associated with their participation in An. coluzzii and An. gambiae-like swarms. Despite extensive introgression between the two species, the X-island genotype of non-hybrid individuals (37.6%), F₁ hybrids (0.1%) and F_1+n recombinants (62.3%) of either sex perfectly associated to each swarm type. Associations between swarm type and the 3L and 2L speciation islands were weakened or broken down by introgression. The functional demonstration of a close association between spatial segregation behaviour and the X speciation island lends further support to sympatric speciation models facilitated by pericentric recombination suppression in this important species complex.

Novel and Broadly Applicable Microsatellite Markers in Identified Chromosomes of the Philippine Dengue Mosquitoes, Aedes aegypti (Diptera: Culicidae)

Dengue is the leading arboviral infection in the Philippines. Its endemicity in the country is due to the presence of its primary mosquito vector, Aedes aegypti (L.). This species has limited microsatellite markers. This study characterized microsatellite markers screened in silico from intergenic regions of the updated reference genome of Ae. aegypti from Liverpool, U.K. Criteria for good markers are: polymorphic, inherited in a Mendelian codominant manner, no null alleles, selectively neutral, randomly associated, and broadly applicable across different regions. Genotypes were scored using ABI Peak Scanner and were screened for the presence of null alleles. Hardy-Weinberg equilibrium, linkage disequilibrium, and robustness of the markers were determined by GENEPOP using Ae. aegypti samples from selected highland and lowland sites (n = 30 each) in the Philippines and outgroups (Thailand and Vietnam). Mendelian codominant inheritance was examined using F1 offspring of Ae. aegypti family (n = 30 each) derived from samples collected from Cebu city highlands and Maramag, Bukidnon. From the 63 randomly selected markers, nine were polymorphic. Two markers (Aaeg1-3D of chromosome 1 and Aaeg3-4C of chromosome 3) satisfied all criteria, hence, are good broadly useful microsatellite markers. Two other markers (Aaeg2-2E of chromosome 2 and Aaeg3-2A of chromosome 3) met all criteria but deviated from Mendelian codominant inheritance. These new markers of the Philippine Ae. aegypti with their chromosomal locations relative to the other published markers are presented, and will ultimately be useful in a variety of population genetic studies of Ae. aegypti to protect the public health.

Monitoring Insecticide Resistance and Target Site Mutations of L1014 Kdr And G119 Ace Alleles in Five Mosquito Populations in Korea

Mosquitoes are globally distributed and important vectors for the transmission of many human diseases. Mosquito control is a difficult task and the cost of preventing mosquito-borne diseases is much lower than that for curing the associated diseases. Thus, chemical control remains the most effective tool for mosquito. Due to the long-term intensive use of insecticides to control mosquito vectors, resistance to most chemical insecticides has been reported. This study aimed to investigate the relationship between insecticide resistance and target site mutation of L1014 kdr and G119 ace alleles in 5 species/species group of mosquitoes (Aedes vexans, Ae. albopictus, Anopheles spp., Culex pipiens complex, and Cx. tritaeniorhynchus) obtained from 6 collection sites. For Anopheles spp., the proportion of mosquitoes with mutated alleles in L1014 was 88.4%, homozygous resistant genotypes were observed in 46.7%, and heterozygous resistant genotypes were observed in 41.8%. For the Cx. pipiens complex and Cx. tritaeniorhynchus species, homozygous resistant genotypes were found in 25.9% and 9.8%, respectively. However, target site mutation of L1014 in the Ae. vexans nipponii and Ae. albopictus species was not observed. Anopheles spp., Cx. pipiens complex, and Cx. tritaeniorhynchus mosquitoes were resistant to deltamethrin and chlorpyriphos, whereas Ae. vexans nipponii and Ae. albopictus were clearly susceptible. We also found a correlation between the resistance phenotype and the presence of the L1014 kdr and G119 ace mutations only in the Anopheles spp. population. In this study, we suggest that insecticide resistance poses a growing threat and resistance management must be integrated into all mosquito control programs.

Chromosome and Genome Divergence between the Cryptic Eurasian Malaria Vector-Species Anopheles messeae and Anopheles daciae

Chromosomal inversions are important drivers of genome evolution. The Eurasian malaria vector Anophelesmesseae has five polymorphic inversions. A cryptic species, An. daciae, has been discriminated from An. messeae based on five fixed nucleotide substitutions in the internal transcribed spacer 2 (ITS2) of ribosomal DNA. However, the inversion polymorphism in An. daciae and the genome divergence between these species remain unexplored. In this study, we sequenced the ITS2 region and analyzed the inversion frequencies of 289 Anopheles larvae specimens collected from three locations in the Moscow region. Five individual genomes for each of the two species were sequenced. We determined that An. messeae and An. daciae differ from each other by the frequency of polymorphic inversions. Inversion X1 was fixed in An. messeae but polymorphic in An. daciae populations. The genome sequence comparison demonstrated genome-wide divergence between the species, especially pronounced on the inversion-rich X chromosome (mean Fst = 0.331). The frequency of polymorphic autosomal inversions was higher in An. messeae than in An. daciae. We conclude that the X chromosome inversions play an important role in the genomic differentiation between the species. Our study determined that An. messeae and An. daciae are closely related species with incomplete reproductive isolation.

Cadmium tolerance pathway in Anopheles gambiae senso stricto

Cadmium is one of the widely used heavy metals (HM) in commercial and industrial products and contributes to environmental contamination in an urban setting. In our previous studies, we established that An. gambiae sensu stricto, a vector of malaria, had adapted to HM pollutants in nature despite their proclivity for unpolluted aquatic habitats. We further demonstrated that heavy metal tolerance adaptation process impacts a biological cost to the fitness of the mosquito and potentially involves the induction of specific HM-responsive transcripts and proteins. Here we interrogated differential proteomic profiles of the cadmium tolerant vs. naïve strains of An. gambiae to shed light on proteomic processes that underpinned biological cost to fitness. We identified a total of 1067 larval proteins and observed significant down-regulation of proteins involved in larval immune responses, energy metabolism, antioxidant enzymes, protein synthesis, and proton transport. Our results suggest that mosquitoes can adjust their biological program through proteome changes to counter HM pollution. Since our study was done in controlled laboratory settings, we acknowledge this may not wholly represent the conditions HM polluted environments. Nevertheless, mosquitoes deploying this strategy have the potential of creating an urban enclave for breeding and thrive and become agents of sporadic malaria epidemics.

Microsatellite Markers Transferability in Culex nigripalpus (Diptera: Culicidae)

Microsatellite loci specifically designed and optimized for secondary vectors and neglected species of mosquitoes are lacking. Cross-species transferability of microsatellite loci between close taxa has been used as a reliable option, being successfully done for Aedes taeniorhynchus Wiedemann (Diptera: Culicidae), Triatoma sordida Stål (Hemiptera: Reduviidae), and Aedes fluviatilis Lutz (Diptera: Culicidae). Nonetheless, prior to the utilization of transferred microsatellites, they must be assessed for their performance and reliability on the target species. Therefore, the goal of this study was to transfer and validate SSR primers previously developed for the Culex pipiens complex to Culex (Culex) nigripalpus Theobald (Diptera: Culicidae). Culex nigripalpus were collected in seven different locations throughout the city of São Paulo using CDC traps baited with dry ice between 2011 and 2013, yielding a total of 223 specimens. Six SSR primers were successfully transferred to Cx. nigripalpus yielding satisfactory results for allele richness, ranging from 5 to 12 alleles per locus. The expected heterozygosity values were higher than the observed heterozygosity (P < 0.05) for the loci CxqCTG10, CxqGT51, and Cxpq59. Linkage disequilibrium results were statistically significant only for the comparison between the loci CxqCA115 and CxqCTG10, and the estimation for the probability of null alleles was also low, ranging from 0 to 0.052. These results substantially validate the SSRs transference to Cx. nigripalpus as well as their subsequent use in microevolution studies comprising this species.

Genetic variation associated with increased insecticide resistance in the malaria mosquito, Anopheles coluzzii

BACKGROUND

Malaria mortality rates in sub-Saharan Africa have declined significantly in recent years as a result of increased insecticide-treated bed net (ITN) usage. A major challenge to further progress is the emergence and spread of insecticide resistance alleles in the Anopheles mosquito vectors, like An. coluzzii. A non-synonymous mutation in the para voltage-gated sodium channel gene reduces pyrethroid-binding affinity, resulting in knockdown resistance (kdr). Metabolic mechanisms of insecticide resistance involving detoxification genes like cytochrome P450 genes, carboxylesterases, and glutathione S-transferases are also important. As some gene activity is tissue-specific and/or environmentally induced, gene regulatory variation may be overlooked when comparing expression from whole mosquito bodies under standard rearing conditions.

RESULTS

We detected complex insecticide resistance in a 2014 An. coluzzii colony from southern Mali using bottle bioassays. Additional bioassays involving recombinant genotypes from a cross with a relatively susceptible 1995 An. coluzzii colony from Mali confirmed the importance of kdr and associated increased permethrin resistance to the CYP9K1 locus on the X chromosome. Significant differential expression of CYP9K1 was not observed among these colonies in Malpighian tubules. However, the P450 gene CYP6Z1 was overexpressed in resistant individuals following sublethal permethrin exposure and the carboxylesterase gene COEAE5G was constitutively overexpressed.

CONCLUSIONS

The significant P450-related insecticide resistance observed in the 2014 An. coluzzii colony indicates that ITNs treated with the P450 inhibitor piperonyl butoxide (PBO) would be more effective in this region. The known insecticide resistance gene CYP6Z1 was differentially expressed exclusively in the context of sublethal permethrin exposure, highlighting the importance of tissue-specificity and environmental conditions in gene expression studies. The increased activity of the carboxylesterase COEAE5G in the resistant An. coluzzii colony suggests resistance to other insecticides like organophosphates. Additional gene expression studies involving other tissues (e.g. fat body) would provide a more comprehensive view of genes underlying metabolic insecticide resistance in An. coluzzii from Mali. Identifying genetic markers linked to these regulatory alleles is an important next step that would substantially improve insecticide resistance surveillance and population genetic studies in this important vector species.

Microgeographical structure in the major Neotropical malaria vector Anopheles darlingi using microsatellites and SNP markers

BACKGROUND

In recent decades, throughout the Amazon Basin, landscape modification contributing to profound ecological change has proceeded at an unprecedented rate. Deforestation that accompanies human activities can significantly change aspects of anopheline biology, though this may be site-specific. Such local changes in anopheline biology could have a great impact on malaria transmission. The aim of this study was to investigate population genetics of the main malaria vector in Brazil, Anopheles darlingi, from a microgeographical perspective.

METHODS

Microsatellites and ddRADseq-derived single nucleotide polymorphisms (SNPs) were used to assess levels of population genetic structuring among mosquito populations from two ecologically distinctive agricultural settlements (~60 km apart) and a population from a distant (~700 km) urban setting in the western Amazon region of Brazil.

RESULTS

Significant microgeographical population differentiation was observed among Anopheles darlingi populations via both model- and non-model-based analysis only with the SNP dataset. Microsatellites detected moderate differentiation at the greatest distances, but were unable to differentiate populations from the two agricultural settlements. Both markers showed low polymorphism levels in the most human impacted sites.

CONCLUSIONS

At a microgeographical scale, signatures of genetic heterogeneity and population divergence were evident in Anopheles darlingi, possibly related to local environmental anthropic modification. This divergence was observed only when using high coverage SNP markers.

Characterisation of Species and Diversity of Anopheles gambiae Keele Colony

Anopheles gambiae sensu stricto was recently reclassified as two species, An. coluzzii and An. gambiae s.s., in wild-caught mosquitoes, on the basis of the molecular form, denoted M or S, of a marker on the X chromosome. The An. gambiae Keele line is an outbred laboratory colony strain that was developed around 12 years ago by crosses between mosquitoes from 4 existing An. gambiae colonies. Laboratory colonies of mosquitoes often have limited genetic diversity because of small starting populations (founder effect) and subsequent fluctuations in colony size. Here we describe the characterisation of the chromosomal form(s) present in the Keele line, and investigate the diversity present in the colony using microsatellite markers on chromosome 3. We also characterise the large 2La inversion on chromosome 2. The results indicate that only the M-form of the chromosome X marker is present in the Keele colony, which was unexpected given that 3 of the 4 parent colonies were probably S-form. Levels of diversity were relatively high, as indicated by a mean number of microsatellite alleles of 6.25 across 4 microsatellites, in at least 25 mosquitoes. Both karyotypes of the inversion on chromosome 2 (2La/2L+^a) were found to be present at approximately equal proportions. The Keele colony has a mixed M- and S-form origin, and in common with the PEST strain, we propose continuing to denote it as an An. gambiae s.s. line.

Patterns of Genome-Wide Variation in Glossina fuscipes fuscipes Tsetse Flies from Uganda

The tsetse fly Glossina fuscipes fuscipes (Gff) is the insect vector of the two forms of Human African Trypanosomiasis (HAT) that exist in Uganda. Understanding Gff population dynamics, and the underlying genetics of epidemiologically relevant phenotypes is key to reducing disease transmission. Using ddRAD sequence technology, complemented with whole-genome sequencing, we developed a panel of ∼73,000 single-nucleotide polymorphisms (SNPs) distributed across the Gff genome that can be used for population genomics and to perform genome-wide-association studies. We used these markers to estimate genomic patterns of linkage disequilibrium (LD) in Gff, and used the information, in combination with outlier-locus detection tests, to identify candidate regions of the genome under selection. LD in individual populations decays to half of its maximum value (r(2) max/2) between 1359 and 2429 bp. The overall LD estimated for the species reaches r(2) max/2 at 708 bp, an order of magnitude slower than in Drosophila Using 53 infected (Trypanosoma spp.) and uninfected flies from four genetically distinct Ugandan populations adapted to different environmental conditions, we were able to identify SNPs associated with the infection status of the fly and local environmental adaptation. The extent of LD in Gff likely facilitated the detection of loci under selection, despite the small sample size. Furthermore, it is probable that LD in the regions identified is much higher than the average genomic LD due to strong selection. Our results show that even modest sample sizes can reveal significant genetic associations in this species, which has implications for future studies given the difficulties of collecting field specimens with contrasting phenotypes for association analysis.

Characteristics of asymptomatic Plasmodium spp. parasitaemia in Kwahu-Mpraeso, a malaria endemic mountainous district in Ghana, West Africa

Background

Malaria control efforts in Ghana have reduced the countrywide average malaria prevalence from 71 % in 2000 to about 51 % in 2012; however, its main focus is on symptomatic malaria. If further progress is to be made, parasite reservoirs in asymptomatic carriers need to be moved into focus. This study profiles asymptomatic Plasmodium spp. parasitaemia amongst residents of mountainous Kwahu-Mpraeso in the Eastern region of Ghana.

Methods

A cross-sectional study of 360 residents was carried out from October to December 2013. This included recording demographics, malaria testing of asymptomatic residents, and gathering of their malaria history. Assessment of malaria transmission was done with molecular identification of vectors, determination of sporozoite rate, insecticide resistance status and biting pattern. Univariate and multivariate analysis were used to establish risk determinants.

Results

In Mpraeso, in the Kwahu highland of Eastern Region, children were at higher risk of asymptomatic parasitaemia, thereby contributing to the parasite reservoir and hence sustained malaria transmission. As well, findings suggested Hb AC genotype influenced susceptibility to asymptomatic malaria with 8.03-fold increase in odds (univariate) and 11.92-fold higher odds (multivariate) than the normal Hb AA. The mosquito vector predominant in the area was Anopheles gambiae sensu stricto of the homozygous pyrethroid resistant form (RR); with biting mainly occurring indoors.

Conclusion

For an effective malaria control in this area, interventions should be formulated and implemented to target asymptomatic parasite reservoirs; especially in children and people with Hb AC. The dominant vector species An. gambiae s.s. and its feeding patterns of biting indoors should also be considered.

Colonization of malaria vectors under semi-field conditions as a strategy for maintaining genetic and phenotypic similarity with wild populations

Background

Malaria still accounts for an estimated 207 million cases and 627,000 deaths worldwide each year. One proposed approach to complement existing malaria control methods is the release of genetically-modified (GM) and/or sterile male mosquitoes. As opposed to laboratory colonization, this requires realistic semi field systems to produce males that can compete for females in nature. This study investigated whether the establishment of a colony of the vector Anopheles arabiensis under more natural semi-field conditions can maintain higher levels of genetic diversity than achieved by laboratory colonization using traditional methods.

Methods

Wild females of the African malaria vector An. arabiensis were collected from a village in southern Tanzania and used to establish new colonies under different conditions at the Ifakara Health Institute. Levels of genetic diversity and inbreeding were monitored in colonies of An. arabiensis that were simultaneously established in small cage colonies in the SFS and in a large semi-field (SFS) cage and compared with that observed in the original founder population. Phenotypic traits that determine their fitness (body size and energetic reserves) were measured at 10^th generation and compared to founder wild population.

Results

In contrast to small cage colonies, the SFS population of An. arabiensis exhibited a higher degree of similarity to the founding field population through time in several ways: (i) the SFS colony maintained a significantly higher level of genetic variation than small cage colonies, (ii) the SFS colony had a lower degree of inbreeding than small cage colonies, and (iii) the mean and range of mosquito body size in the SFS colony was closer to that of the founding wild population than that of small cage colonies. Small cage colonies had significantly lower lipids and higher glycogen abundances than SFS and wild population.

Conclusions

Colonization of An. arabiensis under semi-field conditions was associated with the retention of a higher degree of genetic diversity, reduced inbreeding and greater phenotypic similarity to the founding wild population than observed in small cage colonies. Thus, mosquitoes from such semi-field populations are expected to provide more realistic representation of mosquito ecology and physiology than those from small cage colonies.

First report of an exophilic Anopheles arabiensis population in Bissau City, Guinea-Bissau: recent introduction or sampling bias?

Background

The malaria vector Anopheles arabiensis exhibits greater behavioural and ecological plasticity than the other major vectors of the Anopheles gambiae complex, which presents challenges for major control methods. This study reports for the first time the presence of An. arabiensis in Antula, a suburb of Bissau city, the capital of Guinea Bissau, where high levels of hybridization between Anopheles coluzzii and An. gambiae have been reported. Given that previous surveys in the area, based on indoor collections, did not sample An. arabiensis, the possibility of a recently introduced exophilic population was investigated.

Methods

Larval and adult mosquito collections were carried out in Antula at the end of the rainy season of 2010. Anopheles gambiae species composition, determined by rDNA-IGS and SINE200X6.1 markers, was compared with four previously collected samples dating back to 1993. Analysis of ten microsatellites was used to estimate levels of genetic diversity, relatedness and to investigate demographic stability.

Results

Anopheles arabiensis comprised 54.0% of larvae and 25.6% of adults collected in 2010, but was absent in all previous collections, a highly unlikely observation by chance if the population was stable. This species had the lowest levels of genetic diversity, highest relatedness and, along with An. gambiae, exhibited evidence of a recent population expansion.

Conclusions

Results point to the presence of a previously undetected outdoor population of An. arabiensis in Antula, which appears to have expanded recently, highlighting the importance of complementing indoor-based mosquito collections with sampling methods targeting outdoor adults and immature stages for a more complete assessment of mosquito biodiversity. A change in temporal dynamics in the species complex composition was also detected. Coupled with previous evidence of asymmetric introgression from An. coluzzii to An. gambiae, this suggests that the study area may be subject to ecological changes with a potential impact on both the genetics of these species and on malaria transmission.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-423) contains supplementary material, which is available to authorized users.

Diversity, differentiation, and linkage disequilibrium: prospects for association mapping in the malaria vector Anopheles arabiensis

Association mapping is a widely applied method for elucidating the genetic basis of phenotypic traits. However, factors such as linkage disequilibrium and levels of genetic diversity influence the power and resolution of this approach. Moreover, the presence of population subdivision among samples can result in spurious associations if not accounted for. As such, it is useful to have a detailed understanding of these factors before conducting association mapping experiments. Here we conducted whole-genome sequencing on 24 specimens of the malaria mosquito vector, Anopheles arabiensis, to further understanding of patterns of genetic diversity, population subdivision and linkage disequilibrium in this species. We found high levels of genetic diversity within the An. arabiensis genome, with ~800,000 high-confidence, single- nucleotide polymorphisms detected. However, levels of nucleotide diversity varied significantly both within and between chromosomes. We observed lower diversity on the X chromosome, within some inversions, and near centromeres. Population structure was absent at the local scale (Kilombero Valley, Tanzania) but detected between distant populations (Cameroon vs. Tanzania) where differentiation was largely restricted to certain autosomal chromosomal inversions such as 2Rb. Overall, linkage disequilibrium within An. arabiensis decayed very rapidly (within 200 bp) across all chromosomes. However, elevated linkage disequilibrium was observed within some inversions, suggesting that recombination is reduced in those regions. The overall low levels of linkage disequilibrium suggests that association studies in this taxon will be very challenging for all but variants of large effect, and will require large sample sizes.

Spatiotemporal dynamics of gene flow and hybrid fitness between the M and S forms of the malaria mosquito, Anopheles gambiae

The M and S forms of Anopheles gambiae have been the focus of intense study by malaria researchers and evolutionary biologists interested in ecological speciation. Divergence occurs at three discrete islands in genomes that are otherwise nearly identical. An "islands of speciation" model proposes that diverged regions contain genes that are maintained by selection in the face of gene flow. An alternative "incidental island" model maintains that gene flow between M and S is effectively zero and that divergence islands are unrelated to speciation. A "divergence island SNP" assay was used to explore the spatial and temporal distributions of hybrid genotypes. Results revealed that hybrid individuals occur at frequencies ranging between 5% and 97% in every population examined. A temporal analysis revealed that assortative mating is unstable and periodically breaks down, resulting in extensive hybridization. Results suggest that hybrids suffer a fitness disadvantage, but at least some hybrid genotypes are viable. Stable introgression of the 2L speciation island occurred at one site following a hybridization event.

A new multiplex SNP genotyping assay for detecting hybridization and introgression between the M and S molecular forms of Anopheles gambiae

The M and S forms of Anopheles gambiae have been the subject of intense study, but are morphologically indistinguishable and can only be identified using molecular techniques. PCR-based assays to distinguish the two forms have been designed and applied widely. However, the application of these assays towards identifying hybrids between the two forms, and backcrossed hybrids in particular, has been problematic as the currently available diagnostic assays are based on single locus and/or are located within a multicopy gene. Here, we present an alternative genotyping method for detecting hybridization and introgression between M and S molecular forms based on a multilocus panel of single-nucleotide polymorphisms (SNPs) fixed between the M and S forms. The panel of SNPs employed is located in so-called islands of divergence leading us to describe this method as the 'Divergence Island SNP' (DIS) assay. We show this multilocus SNP genotyping approach can robustly and accurately detect F1 hybrids as well as backcrossed individuals.

A preliminary investigation of the relationship between water quality and Anopheles gambiae larval habitats in Western Cameroon

Background

Water quality and anopheline habitat have received increasing attention due to the possibility that challenges during larval life may translate into adult susceptibility to malaria parasite infection and/or insecticide resistance.

Methods

A preliminary study of Anopheles gambiae s.s. larval habitats in the north-west and south-west regions of Cameroon was conducted in order to detect associations between An. gambiae s.s. molecular form and 2La inversion distributions with basic water quality parameters. Water quality was measured by temperature, pH, conductivity, total dissolved solids (TDS) at seven sites in Cameroon and one site in Selinkenyi, Mali.

Results

Principal components and correlation analyses indicated a complex relationship between 2La polymorphism, temperature, conductivity and TDS. Cooler water sites at more inland locations yielded more S form larvae with higher 2La inversion polymorphism while warmer water sites yielded more M form larvae with rare observations of the 2La inversion.

Discussion

More detailed studies that take into account the population genetics but also multiple life stages, environmental data relative to these life stages and interactions with both humans and the malaria parasite may help us to understand more about how and why this successful mosquito is able to adapt and diverge, and how it can be successfully managed.

Geographic population structure of the African malaria vector Anopheles gambiae suggests a role for the forest-savannah biome transition as a barrier to gene flow

The primary Afrotropical malaria mosquito vector Anopheles gambiae sensu stricto has a complex population structure. In west Africa, this species is split into two molecular forms and displays local and regional variation in chromosomal arrangements and behaviors. To investigate patterns of macrogeographic population substructure, 25 An. gambiae samples from 12 African countries were genotyped at 13 microsatellite loci. This analysis detected the presence of additional population structuring, with the M-form being subdivided into distinct west, central, and southern African genetic clusters. These clusters are coincident with the central African rainforest belt and northern and southern savannah biomes, which suggests restrictions to gene flow associated with the transition between these biomes. By contrast, geographically patterned population substructure appears much weaker within the S-form.

Analysis of population genetic structure of Indian Anopheles culicifacies species A using microsatellite markers

Background

Anopheles culicifacies sensu lato is an important vector of malaria in Southeast Asia contributing to almost 70% of malaria cases in India. It exists as morphologically similar sibling species A, B, C, D and E with varied geographical distribution patterns. Vector control measures have been difficult for this important vector as the sibling species have developed varying levels of resistance to the currently used insecticides. In view of the importance of this vector, we developed and validated a set of microsatellite markers and the same were used to analyze the population genetic structure of five different geographical populations of An. culicifacies A.

Methods

Anopheles culicifacies A samples were collected from different localities across India, and genotyping was performed using eight microsatellite markers on ABI Prism 310 Genetic Analyzer. Several statistical analyses were performed to ascertain the genetic diversity that exists within and between the populations.

Results

The markers were found to be moderately polymorphic in the populations. Genetic analysis indicated significant genetic differentiation between the majority of the population pairs analyzed and was not found to be related to the geographical distances between populations.

Conclusion

This is the first and successful attempt to test the microsatellite markers developed for population genetic analysis of An. culicifacies A. Host feeding and breeding habits of species A suggest that factors other than ecological and geographical barriers were responsible for the genetic differentiation that has been observed between the populations.

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.

Chromosomal localization of actin genes in the malaria mosquito Anopheles darlingi

Physical and genetic maps have been used for chromosomal localization of genes in vectors of infectious diseases. The availability of polytene chromosomes in malaria mosquitoes provides a unique opportunity to precisely map genes of interest. We report the physical mapping of two actin genes on polytene chromosomes of the major malaria vector in the Amazon, Anopheles darlingi (Diptera: Culicidae). Clones with actin gene sequences were obtained from a cDNA library constructed from RNA isolated from adult females and males of An. darlingi. Each of the two clones was mapped to a unique site on chromosomal arm 2L in subdivisions 21A (clone pl05-A04) and 23B (clone pl17-G06). The obtained results, together with previous mapping data, provide a suitable basis for comparative genomics and for establishing chromosomal homologies among major malaria vectors.

No evidence for biased co-transmission of speciation islands in Anopheles gambiae

Genome-scale scans have revealed highly heterogeneous levels of divergence between closely related taxa in many systems. Generally, a small number of regions show high differentiation, with the rest of the genome showing no or only low levels of divergence. These patterns have been interpreted as evidence for ongoing speciation-with-gene-flow, with introgression homogenizing the whole genome except loci involved in reproductive isolation. However, as the number of selected loci increases, the probability of introgression at unselected loci decreases unless there is a transmission ratio distortion causing an over-representation of specific combinations of alleles. Here we examine the transmission of three 'speciation islands' that contain fixed differences between the M and S forms of the mosquito, Anopheles gambiae. We made reciprocal crosses between M and S parents and genotyped over 2000 F(2) individuals, developing a hierarchical likelihood model to identify specific genotypes that are under- or over-represented among the recombinant offspring. Though our overall results did not match the expected number of F(2) genotypes, we found no biased co-transmission among M or S alleles in the three islands. Our likelihood model did identify transmission ratio distortion at two of the three islands, but this distortion was small (approx. 3%) and in opposite directions for the two islands. We discuss how our results impinge on hypotheses of current gene flow between M and S and ongoing speciation-with-gene-flow in this system.

Can Wolbachia be used to control malaria?

Malaria is a mosquito-borne infectious disease caused by Plasmodium parasites transmitted by the infectious bite of Anopheles mosquitoes. Vector control of malaria has predominantly focused on targeting the adult mosquito through insecticides and bed nets. However, current vector control methods are often not sustainable for long periods so alternative methods are needed. A novel biocontrol approach for mosquito-borne diseases has recently been proposed, it uses maternally inherited endosymbiotic Wolbachia bacteria transinfected into mosquitoes in order to interfere with pathogen transmission. Transinfected Wolbachia strains in Aedes aegypti mosquitoes, the primary vector of dengue fever, directly inhibit pathogen replication, including Plasmodium gallinaceum, and also affect mosquito reproduction to allow Wolbachia to spread through mosquito populations. In addition, transient Wolbachia infections in Anopheles gambiae significantly reduce Plasmodium levels. Here we review the prospects of using a Wolbachia-based approach to reduce human malaria transmission through transinfection of Anopheles mosquitoes.

Evolution of Anopheles gambiae in Relation to Humans and Malaria

The closely related and morphologically indistinguishable mosquito species in the Afrotropical Anopheles gambiae complex differ dramatically in their contribution to malaria transmission, ranging from major vectors through minor or locally important vectors and nonvectors. Radiation of the A. gambiae complex and ongoing diversification within its nominal species appears to be a product of recent and rapid adaptation to environmental heterogeneities, notably those of anthropogenic origin. Polytene chromosome and genomic analyses suggest that paracentric chromosomal inversions and possibly other low-recombination regions have played instrumental roles in this process by facilitating ecotypic differentiation both within and across semipermeable species boundaries. Forthcoming complete genome sequences from several members of the A. gambiae complex will provide powerful tools to accelerate ongoing investigation of how genetic diversification of populations and species has shaped behavioral and physiological traits, such as vector competence, that bear on vectorial importance.

Morphological differentiation may mediate mate-choice between incipient species of Anopheles gambiae s.s

The M and S molecular forms of Anopheles gambiae s.s. have been considered incipient species for more than ten years, yet the mechanism underlying assortative mating of these incipient species has remained elusive. The discovery of the importance of harmonic convergence of wing beat frequency in mosquito mating and its relation to wing size have laid the foundation for exploring phenotypic divergence in wing size of wild populations of the two forms. In this study, wings from field collected mosquitoes were measured for wing length and wing width from two parts of the sympatric distribution, which differ with respect to the strength of assortative mating. In Mali, where assortative mating is strong, as evidenced by low rates of hybridization, mean wing lengths and wing widths were significantly larger than those from Guinea-Bissau. In addition, mean wing widths in Mali were significantly different between molecular forms. In Guinea-Bissau, assortative mating appears comparatively reduced and wing lengths and widths did not differ significantly between molecular forms. The data presented in this study support the hypothesis that wing beat frequency may mediate assortative mating in the incipient species of A. gambiae and represent the first documentation of a morphological difference between the M and S molecular forms.

Confinement of gene drive systems to local populations: a comparative analysis

Mosquito-borne diseases such as malaria and dengue fever pose a major health problem through much of the world. One approach to disease prevention involves the use of selfish genetic elements to drive disease-refractory genes into wild mosquito populations. Recently engineered synthetic drive systems have provided encouragement for this strategy; but at the same time have been greeted with caution over the concern that transgenes may spread into countries and communities without their consent. Consequently, there is also interest in gene drive systems that, while strong enough to bring about local population replacement, are unable to establish themselves beyond a partially isolated release site, at least during the testing phase. Here, we develop simple deterministic and stochastic models to compare the confinement properties of a variety of gene drive systems. Our results highlight several systems with desirable features for confinement-a high migration rate required to become established in neighboring populations, and low-frequency persistence in neighboring populations for moderate migration rates. Single-allele underdominance and single-locus engineered underdominance have the strongest confinement properties, but are difficult to engineer and require a high introduction frequency, respectively. Toxin-antidote systems such as Semele, Merea and two-locus engineered underdominance show promising confinement properties and require lower introduction frequencies. Killer-rescue is self-limiting in time, but is able to disperse to significant levels in neighboring populations. We discuss the significance of these results in the context of a phased release of transgenic mosquitoes, and the need for characterization of local ecology prior to a release.

Improving the population genetics toolbox for the study of the African malaria vector Anopheles nili: microsatellite mapping to chromosomes

BACKGROUND

Anopheles nili is a major vector of malaria in the humid savannas and forested areas of sub-Saharan Africa. Understanding the population genetic structure and evolutionary dynamics of this species is important for the development of an adequate and targeted malaria control strategy in Africa. Chromosomal inversions and microsatellite markers are commonly used for studying the population structure of malaria mosquitoes. Physical mapping of these markers onto the chromosomes further improves the toolbox, and allows inference on the demographic and evolutionary history of the target species.

RESULTS

Availability of polytene chromosomes allowed us to develop a map of microsatellite markers and to study polymorphism of chromosomal inversions. Nine microsatellite markers were mapped to unique locations on all five chromosomal arms of An. nili using fluorescent in situ hybridization (FISH). Probes were obtained from 300-483 bp-long inserts of plasmid clones and from 506-559 bp-long fragments amplified with primers designed using the An. nili genome assembly generated on an Illumina platform. Two additional loci were assigned to specific chromosome arms of An. nili based on in silico sequence similarity and chromosome synteny with Anopheles gambiae. Three microsatellites were mapped inside or in the vicinity of the polymorphic chromosomal inversions 2Rb and 2Rc. A statistically significant departure from Hardy-Weinberg equilibrium, due to a deficit in heterozygotes at the 2Rb inversion, and highly significant linkage disequilibrium between the two inversions, were detected in natural An. nili populations collected from Burkina Faso.

CONCLUSIONS

Our study demonstrated that next-generation sequencing can be used to improve FISH for microsatellite mapping in species with no reference genome sequence. Physical mapping of microsatellite markers in An. nili showed that their cytological locations spanned the entire five-arm complement, allowing genome-wide inferences. The knowledge about polymorphic inversions and chromosomal locations of microsatellite markers has been useful for explaining differences in genetic variability across loci and significant differentiation observed among natural populations of An. nili.

Population genetic structure of Anopheles arabiensis and Anopheles gambiae in a malaria endemic region of southern Tanzania

BACKGROUND

Genetic diversity is a key factor that enables adaptation and persistence of natural populations towards environmental conditions. It is influenced by the interaction of a natural population's dynamics and the environment it inhabits. Anopheles gambiae s.s. and Anopheles arabiensis are the two major and widespread malaria vectors in sub-Saharan Africa. Several studies have examined the ecology and population dynamics of these vectors. Ecological conditions along the Kilombero valley in Tanzania influence the distribution and population density of these two vector species. It remains unclear whether the ecological diversity within the Kilombero valley has affected the population structure of An. gambiae s.l. populations. The goal of this study was to characterise the genetic structure of sympatric An. gambiae s.s and An. arabiensis populations along the Kilombero valley.

METHODOLOGY

Mosquitoes were collected from seven locations in Tanzania: six from the Kilombero valley and one outside the valley (-700 km away) as an out-group. To archive a genome-wide coverage, 13 microsatellite markers from chromosomes X, 2 and 3 were used.

RESULTS

High levels of genetic differentiation among An. arabiensis populations was observed, as opposed to An. gambiae s.s., which was genetically undifferentiated across the 6,650 km2 of the Kilombero valley landscape. It appears that genetic differentiation is not attributed to physical barriers or distance, but possibly by ecological diversification within the Kilombero valley. Genetic divergence among An. arabiensis populations (FST = 0.066) was higher than that of the well-known M and S forms of An. gambiae s. s. in West and Central Africa (FST = 0.035), suggesting that these populations are maintained by some level of reproductive isolation.

CONCLUSION

It was hypothesized that ecological diversification across the valley may be a driving force for observed An. arabiensis genetic divergence. The impact of the observed An. arabiensis substructure to the prospects for new vector control approaches is discussed.

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.

An integrated chromosome map of microsatellite markers and inversion breakpoints for an Asian malaria mosquito, Anopheles stephensi

Anopheles stephensi is one of the major vectors of malaria in the Middle East and Indo-Pakistan subcontinent. Understanding the population genetic structure of malaria mosquitoes is important for developing adequate and successful vector control strategies. Commonly used markers for inferring anopheline taxonomic and population status include microsatellites and chromosomal inversions. Knowledge about chromosomal locations of microsatellite markers with respect to polymorphic inversions could be useful for better understanding a genetic structure of natural populations. However, fragments with microsatellites used in population genetic studies are usually too short for successful labeling and hybridization with chromosomes. We designed new primers for amplification of microsatellite loci identified in the A. stephensi genome sequenced with next-generation technologies. Twelve microsatellites were mapped to polytene chromosomes from ovarian nurse cells of A. stephensi using fluorescent in situ hybridization. All microsatellites hybridized to unique locations on autosomes, and 7 of them localized to the largest arm 2R. Ten microsatellites were mapped inside the previously described polymorphic chromosomal inversions, including 4 loci located inside the widespread inversion 2Rb. We analyzed microsatellite-based population genetic data available for A. stephensi in light of our mapping results. This study demonstrates that the chromosomal position of microsatellites may affect estimates of population genetic parameters and highlights the importance of developing physical maps for nonmodel organisms.

Genetic structure of Anopheles (Nyssorhynchus) marajoara (Diptera: Culicidae) in Colombia

Five Anopheles marajoara Galvão and Damasceno populations, representing diverse ecological conditions, were sampled throughout Colombia and analyzed using nine hypervariable DNA microsatellite loci. The overall genetic diversity (H = 0.58) was lower than that determined for some Brazilian populations using the same markers. The Caquetá population (Colombia) had the lowest gene diversity (H = 0.48), and it was the only population at Hardy-Weinberg equilibrium. Hardy-Weinberg disequilibrium in the remaining four populations was probably caused by the Wahlund effect. The assignment analyses showed two incompletely isolated gene pools separated by the Eastern Andean cordillera. However, other possible geographical barriers (rivers and other mountains) did not play any role in the moderate genetic heterogeneity found among these populations (F(ST) = 0.069). These results are noteworthy, because this species is a putative malaria vector in Colombia.

Population genetic structure of Anopheles gambiae and Anopheles arabiensis in Niger

The increasing usage of long-lasting insecticide-treated nets allows protection of millions of people from malaria infection. Monitoring studies should be planned during any wide-scale malaria control program integrating insecticide-treated materials, to evaluate their effects and effectiveness on epidemiologically relevant parameters. Such operational control interventions may be challenged by insecticide resistance spread within vector populations, as a result of wide insecticide pressure. A nationwide distribution of long-lasting insecticidal nets was implemented throughout Niger in 2005. We studied the population genetic structure of major malaria vectors across Nigerien Sahel, and investigated potential effects of this large malaria control intervention. Wild-caught Anopheles gambiae sensu lato females from seven villages and two wet seasons were genotyped at 12 microsatellite loci. The genetic diversity within both species appeared homogenous between villages and years. The estimated genetic differentiation among samples was very low within both species, indicating high gene flow across the area. An absence of differentiation was also found between 2005 and 2006 wet seasons, for all samples but one, showing that the net distribution did not impact significantly the genetic diversity and structure of vector populations in a single year. We provide valuable results participating to document effects of large malaria control programs, to maximize the efficiency of available tools in future interventions.

Effects of latitude and longitude on the population structure of Culex pipiens s.l., vectors of West Nile virus in North America

We assessed the structure and latitudinal selection that might result in sensitivities to critical day-lengths that trigger diapause between Culex pipiens populations distributed along North-South and East-West axes in eastern North America. Strong population structure between Cx. p. pipiens and Cx. p. quinquefasciatus existed. Among Cx. p. pipiens, a 100-km increase in the latitudinal change resulted in an increased square root of F(ST) by 0.002. A 100-km increase in the longitudinal change caused an increased square root of F(ST) by 0.035. A lack of latitudinal influence on the structure between Cx. p. pipiens populations suggests a uniform signal using the 12 microsatellite markers, which might increase the risk of West Nile virus (WNV) transmission toward northern areas because of longer breeding season, extend host-seeking period, and larger population size. Northern Cx. p. pipiens may have undergone additional generations before diapause is triggered, magnifying population size when WNV amplification is peaking.

Living at the edge: biogeographic patterns of habitat segregation conform to speciation by niche expansion in Anopheles gambiae

BACKGROUND

Ongoing lineage splitting within the African malaria mosquito Anopheles gambiae is compatible with ecological speciation, the evolution of reproductive isolation by divergent natural selection acting on two populations exploiting alternative resources. Divergence between two molecular forms (M and S) identified by fixed differences in rDNA, and characterized by marked, although incomplete, reproductive isolation is occurring in West and Central Africa. To elucidate the role that ecology and geography play in speciation, we carried out a countrywide analysis of An. gambiae M and S habitat requirements, and that of their chromosomal variants, across Burkina Faso.

RESULTS

Maps of relative abundance by geostatistical interpolators produced a distinct pattern of distribution: the M-form dominated in the northernmost arid zones, the S-form in the more humid southern regions. Maps of habitat suitability, quantified by Ecological Niche Factor Analysis based on 15 eco-geographical variables revealed less contrast among forms. M was peculiar as it occurred proportionally more in habitat of marginal quality. Measures of ecological niche breadth and overlap confirmed the mismatch between the fundamental and realized patterns of habitat occupation: forms segregated more than expected from the extent of divergence of their environmental envelope--a signature of niche expansion. Classification of chromosomal arm 2R karyotypes by multilocus genetic clustering identified two clusters loosely corresponding to molecular forms, with 'mismatches' representing admixed individuals due to shared ancestral polymorphism and/or residual hybridization. In multivariate ordination space, these karyotypes plotted in habitat of more marginal quality compared to non-admixed, 'typical', karyotypes. The distribution of 'typical' karyotypes along the main eco-climatic gradient followed a consistent pattern within and between forms, indicating an adaptive role of inversions at this geographical scale.

CONCLUSION

Ecological segregation between M and S is consistent with niche expansion into marginal habitats by chromosomal inversion variants during early lineage divergence; presumably, this process is promoted by inter-karyotype competition in the higher-quality core habitat. We propose that the appearance of favourable allelic combinations in other regions of suppressed recombination (e.g. pericentromeric portions defining speciation islands in An. gambiae) fosters development of reproductive isolation to protect linkage between separate chromosomal regions.

Ecological and genetic relationships of the Forest-M form among chromosomal and molecular forms of the malaria vector Anopheles gambiae sensu stricto

BACKGROUND

Anopheles gambiae sensu stricto, one of the principal vectors of malaria, has been divided into two subspecific groups, known as the M and S molecular forms. Recent studies suggest that the M form found in Cameroon is genetically distinct from the M form found in Mali and elsewhere in West Africa, suggesting further subdivision within that form.

METHODS

Chromosomal, microsatellite and geographic/ecological evidence are synthesized to identify sources of genetic polymorphism among chromosomal and molecular forms of the malaria vector Anopheles gambiae s.s.

RESULTS

Cytogenetically the Forest M form is characterized as carrying the standard chromosome arrangement for six major chromosomal inversions, namely 2La, 2Rj, 2Rb, 2Rc, 2Rd, and 2Ru. Bayesian clustering analysis based on molecular form and chromosome inversion polymorphisms as well as microsatellites describe the Forest M form as a distinct population relative to the West African M form (Mopti-M form) and the S form. The Forest-M form was the most highly diverged of the An. gambiae s.s. groups based on microsatellite markers. The prevalence of the Forest M form was highly correlated with precipitation, suggesting that this form prefers much wetter environments than the Mopti-M form.

CONCLUSION

Chromosome inversions, microsatellite allele frequencies and habitat preference all indicate that the Forest M form of An. gambiae is genetically distinct from the other recognized forms within the taxon Anopheles gambiae sensu stricto. Since this study covers limited regions of Cameroon, the possibility of gene flow between the Forest-M form and Mopti-M form cannot be rejected. However, association studies of important phenotypes, such as insecticide resistance and refractoriness against malaria parasites, should take into consideration this complex population structure.

Ace-1 duplication in Anopheles gambiae: a challenge for malaria control

Background

Insecticide resistance is a rapid and recent evolutionary phenomenon with serious economic and public health implications. In the mosquito Anopheles gambiae s.s., main vector of malaria, resistance to organophosphates and carbamates is mainly due to a single amino-acid substitution in acetylcholinesterase 1 (AChE1). This mutation entails a large fitness cost. However, a resistant duplicated allele of the gene encoding AChE1 (ace-1), potentially associated to a lower fitness cost, recently appeared in An. gambiae.

Methods

Using molecular phenotype data collected from natural populations from West Africa, the frequency of this duplicated allele was investigated by statistical inference. This method is based on the departure from Hardy-Weinberg phenotypic frequency equilibrium caused by the presence of this new allele.

Results

The duplicated allele, Ag-ace-1^D, reaches a frequency up to 0.65 in Ivory Coast and Burkina Faso, and is potentially present in Benin. A previous study showed that Ag-ace-1^D, present in both M and S molecular forms in different West Africa countries, was generated by a single genetic event. This single origin and its present distribution suggest that this new allele is currently spreading.

Conclusion

The spread of this less costly resistance allele could represent a major threat to public health, as it may impede An. gambiae control strategies, and thus increases the risk of malaria outbreaks.

Population structure and geographical subdivision of the Leishmania major vector Phlebotomus papatasi as revealed by microsatellite variation

Multi-locus microsatellite typing (MLMT) has been employed to infer the population structure of Phlebotomus papatasi (Scopoli) (Diptera: Psychodidae) sandflies and assign individuals to populations. Phlebotomus papatasi sandflies were collected from 35 sites in 15 countries. A total of 188 P. papatasi individuals were typed using five microsatellite loci, resulting in 113 different genotypes. Unique microsatellite signatures were observed for some of the populations analysed. Comparable results were obtained when the data were analysed with Bayesian model and distance-based methods. Bayesian statistic-based analyses split the dataset into two distinct genetic clusters, A and B, with further substructuring within each. Population A consisted of five subpopulations representing large numbers of alleles that were correlated with the geographical origins of the sandflies. Cluster B comprised individuals collected in the Middle East and the northern Mediterranean area. The subpopulations B1 and B2 did not, however, show any further correlation to geographical origin. The genetic differentiation between subpopulations was supported by F statistics showing statistically significant (Bonferroni-corrected P < 0.005) values of 0.221 between B2 and B1 and 0.816 between A5 and A4. Identification of the genetic structure of P. papatasi populations is important for understanding the patterns of dispersal of this species and to developing strategies for sandfly control.

The effect of gene drive on containment of transgenic mosquitoes

Mosquito-borne diseases such as malaria and dengue fever continue to be a major health problem through much of the world. Several new potential approaches to disease control utilize gene drive to spread anti-pathogen genes into the mosquito population. Prior to a release, these projects will require trials in outdoor cages from which transgenic mosquitoes may escape, albeit in small numbers. Most genes introduced in small numbers are very likely to be lost from the environment; however, gene drive mechanisms enhance the invasiveness of introduced genes. Consequently, introduced transgenes may be more likely to persist than ordinary genes following an accidental release. Here, we develop stochastic models to analyze the loss probabilities for several gene drive mechanisms, including homing endonuclease genes, transposable elements, Medea elements, the intracellular bacterium Wolbachia, engineered underdominance genes, and meiotic drive. We find that Medea and Wolbachia present the best compromise between invasiveness and containment for the six gene drive systems currently being considered for the control of mosquito-borne disease.

High levels of hybridization between molecular forms of Anopheles gambiae from Guinea Bissau

In the malaria vector Anopheles gambiae Giles sensu stricto, two molecular forms denoted M and S are considered units of incipient speciation within this species. Very low hybrid frequencies and significant genetic differentiation have been found in sympatric M- and S-form populations. We studied the molecular form composition and the degree of genetic differentiation at 15 microsatellites in two samples of An. gambiae collected in two consecutive years from Bissau, Guinea Bissau. High frequencies of M/S hybrids (19-24%) were found in this area. Coincidently, very low levels of genetic differentiation were detected between forms when analysis involved microsatellites mapped at chromosome-3 (mean Fst, 0.000-0.002). The single exception was the X-linked AGXH678, for which high differentiation was measured (Fst, 0.158-0.301). This locus maps near the centromere of chromosome X, a low recombination region in which selection is likely to promote divergence between M and S forms. These results strongly suggest that the degree of isolation between M and S forms, considered the units of incipient speciation within An. gambiae, is not homogenous throughout the species distribution range.

Insertion polymorphisms of SINE200 retrotransposons within speciation islands of Anopheles gambiae molecular forms

Background

SINEs (Short INterspersed Elements) are homoplasy-free and co-dominant genetic markers which are considered to represent useful tools for population genetic studies, and could help clarifying the speciation processes ongoing within the major malaria vector in Africa, Anopheles gambiae s.s. Here, we report the results of the analysis of the insertion polymorphism of a nearly 200 bp-long SINE (SINE200) within genome areas of high differentiation (i.e. "speciation islands") of M and S A. gambiae molecular forms.

Methods

A SINE-PCR approach was carried out on thirteen SINE200 insertions in M and S females collected along the whole range of distribution of A. gambiae s.s. in sub-Saharan Africa. Ten specimens each for Anopheles arabiensis, Anopheles melas, Anopheles quadriannulatus A and 15 M/S hybrids from laboratory crosses were also analysed.

Results

Eight loci were successfully amplified and were found to be specific for A. gambiae s.s.: 5 on 2L chromosome and one on X chromosome resulted monomorphic, while two loci positioned respectively on 2R (i.e. S200 2R12D) and X (i.e. S200 X6.1) chromosomes were found to be polymorphic. S200 2R12D was homozygote for the insertion in most S-form samples, while intermediate levels of polymorphism were shown in M-form, resulting in an overall high degree of genetic differentiation between molecular forms (Fst = 0.46 p < 0.001) and within M-form (Fst = 0.46 p < 0.001). The insertion of S200 X6.1 was found to be fixed in all M- and absent in all S-specimens. This led to develop a novel easy-to-use PCR approach to straightforwardly identify A. gambiae molecular forms. This novel approach allows to overcome the constraints associated with markers on the rDNA region commonly used for M and S identification. In fact, it is based on a single copy and irreversible SINE200 insertion and, thus, is not subjected to peculiar evolutionary patterns affecting rDNA markers, e.g. incomplete homogenization of the arrays through concerted evolution and/or mixtures of M and S IGS-sequences among the arrays of single chromatids.

Conclusion

The approach utilized allowed to develop new easy-to-use co-dominant markers for the analysis of genetic differentiation between M and S-forms and opens new perspectives in the study of the speciation process ongoing within A. gambiae.

Assessment of alternative mating strategies in Anopheles gambiae: Does mating occur indoors?

Mating in Anopheles gambiae has been observed only in outdoor swarms. Here we evaluate whether mating also occurs indoors. Mark-release-recapture of virgin males and females in natural houses showed that mating occurred over a single day even when mosquitoes can leave the house through exit traps and without adaptation to laboratory conditions. In these experiments, insemination rate in the M molecular form of An. gambiae (and An. arabiensis) was higher than that of the S form (15 versus 6%). Under these conditions, smaller females of the M form mated more frequently than larger females of that form. Sampling mosquitoes throughout the day showed that both sexes enter houses around sunrise and leave around sunset, staying indoors together from dawn to dusk. In an area dominated by the M form, the daily rate of insemination in samples from exit traps was approximately 5% higher than in those from entry traps, implying that mating occurred indoors. Importantly, frequency of cross mating between the molecular forms was as high as that between members of the same form, indicating that, indoors, assortative mating breaks down. Altogether, these results suggest that indoor mating is an alternative mating strategy of the M molecular form of An. gambiae. Because naturally occurring mating couples have not yet been observed indoors, this conclusion awaits validation.

The molecular forms of Anopheles gambiae: a phenotypic perspective

The African malaria mosquito Anopheles gambiae is undergoing speciation, being split into the M and S molecular forms. Speciation is the main process promoting biological diversity, thus, new vector species might complicate disease transmission. Genetic differentiation between the molecular forms has been extensively studied, but phenotypic differences between them, the evolutionary forces that generated divergence, and the mechanisms that maintain their genetic isolation have only recently been addressed. Here, we review recent studies suggesting that selection mediated by larval predation and competition promoted divergence between temporary and permanent freshwater habitats. These differences explain the sharp discontinuity in distribution of the molecular forms between rice fields and surrounding savanna, but they can also explain the concurrent cline between humid and arid environments due to the dependence on permanent habitats in the latter. Although less pronounced, differences in adult body size, reproductive output, and longevity also suggest that the molecular forms have adapted to distinct niches. Reproductive isolation between the molecular forms is achieved by spatial swarm segregation, although within-swarm mate recognition appears to play a role in certain locations. The implications of these results to disease transmission and control are discussed and many of the gaps in our understanding are highlighted.

Microsatellite data suggest significant population structure and differentiation within the malaria vector Anopheles darlingi in Central and South America

Background

Anopheles darlingi is the most important malaria vector in the Neotropics. An understanding of A. darlingi's population structure and contemporary gene flow patterns is necessary if vector populations are to be successfully controlled. We assessed population genetic structure and levels of differentiation based on 1,376 samples from 31 localities throughout the Peruvian and Brazilian Amazon and Central America using 5–8 microsatellite loci.

Results

We found high levels of polymorphism for all of the Amazonian populations (mean R_S = 7.62, mean H_O = 0.742), and low levels for the Belize and Guatemalan populations (mean R_S = 4.3, mean H_O = 0.457). The Bayesian clustering analysis revealed five population clusters: northeastern Amazonian Brazil, southeastern and central Amazonian Brazil, western and central Amazonian Brazil, Peruvian Amazon, and the Central American populations. Within Central America there was low non-significant differentiation, except for between the populations separated by the Maya Mountains. Within Amazonia there was a moderate level of significant differentiation attributed to isolation by distance. Within Peru there was no significant population structure and low differentiation, and some evidence of a population expansion. The pairwise estimates of genetic differentiation between Central America and Amazonian populations were all very high and highly significant (F_ST = 0.1859 – 0.3901, P < 0.05). Both the D_A and F_ST distance-based trees illustrated the main division to be between Central America and Amazonia.

Conclusion

We detected a large amount of population structure in Amazonia, with three population clusters within Brazil and one including the Peru populations. The considerable differences in N_e among the populations may have contributed to the observed genetic differentiation. All of the data suggest that the primary division within A. darlingi corresponds to two white gene genotypes between Amazonia (genotype 1) and Central America, parts of Colombia and Venezuela (genotype 2), and are in agreement with previously published mitochondrial COI gene sequences interpreted as incipient species. Overall, it appears that two main factors have contributed to the genetic differentiation between the population clusters: physical distance between the populations and the differences in effective population sizes among the subpopulations.