Standing genetic variation and chromosome differences drove rapid ecotype formation in a major malaria mosquito.

Species distributed across heterogeneous environments often evolve locally adapted ecotypes, but understanding of the genetic mechanisms involved in their formation and maintenance in the face of gene flow is incomplete. In Burkina Faso, the major African malaria mosquito Anopheles funestus comprises two strictly sympatric and morphologically indistinguishable yet karyotypically differentiated forms reported to differ in ecology and behavior. However, knowledge of the genetic basis and environmental determinants of An. funestus diversification was impeded by lack of modern genomic resources. Here, we applied deep whole-genome sequencing and analysis to test the hypothesis that these two forms are ecotypes differentially adapted to breeding in natural swamps versus irrigated rice fields. We demonstrate genome-wide differentiation despite extensive microsympatry, synchronicity, and ongoing hybridization. Demographic inference supports a split only ~1,300 y ago, closely following the massive expansion of domesticated African rice cultivation ~1,850 y ago. Regions of highest divergence, concentrated in chromosomal inversions, were under selection during lineage splitting, consistent with local adaptation. The origin of nearly all variations implicated in adaptation, including chromosomal inversions, substantially predates the ecotype split, suggesting that rapid adaptation was fueled mainly by standing genetic variation. Sharp inversion frequency differences likely facilitated adaptive divergence between ecotypes by suppressing recombination between opposing chromosomal orientations of the two ecotypes, while permitting free recombination within the structurally monomorphic rice ecotype. Our results align with growing evidence from diverse taxa that rapid ecological diversification can arise from evolutionarily old structural genetic variants that modify genetic recombination.

Diurnal biting of malaria mosquitoes in the Central African Republic indicates residual transmission may be "out of control".

Malaria control interventions target nocturnal feeding of the Anopheles vectors indoors to reduce parasite transmission. Mass deployment of insecticidal bed nets and indoor residual spraying with insecticides, however, may induce mosquitoes to blood-feed at places and at times when humans are not protected. These changes can set a ceiling to the efficacy of these control interventions, resulting in residual malaria transmission. Despite its relevance for disease transmission, the daily rhythmicity of Anopheles biting behavior is poorly documented, most investigations focusing on crepuscular hours and nighttime. By performing mosquito collections 48-h around the clock, both indoors and outdoors, and by modeling biting events using circular statistics, we evaluated the full daily rhythmicity of biting in urban Bangui, Central African Republic. While the bulk of biting by Anopheles gambiae, Anopheles coluzzii, Anopheles funestus, and Anopheles pharoensis occurred from sunset to sunrise outdoors, unexpectedly ∼20 to 30% of indoor biting occurred during daytime. As biting events did not fully conform to any family of circular distributions, we fitted mixtures of von Mises distributions and found that observations were consistent with three compartments, corresponding indoors to populations of early-night, late-night, and daytime-biting events. It is not known whether these populations of biting events correspond to spatiotemporal heterogeneities or also to distinct mosquito genotypes/phenotypes belonging consistently to each compartment. Prevalence of Plasmodium falciparum in nighttime- and daytime-biting mosquitoes was the same. As >50% of biting occurs in Bangui when people are unprotected, malaria control interventions outside the domiciliary environment should be envisaged.

Variation in seed traits among Mediterranean oaks in Tunisia and their ecological significance.

BACKGROUND AND AIMS: Oaks are the foundation and dominant tree species of most Mediterranean forests. As climate models predict dramatic changes in the Mediterranean basin, a better understanding of the ecophysiology of seed persistence and germination in oaks could help define their regeneration niches. Tunisian oaks occupy distinct geographical areas, which differ in their rainfall and temperature regimes, and are thus a valuable model to investigate relationships between seed traits and species ecological requirements. METHODS: Seed morphological traits, desiccation sensitivity level, lethal freezing temperature, embryonic axis and cotyledon sugar and lipid composition, and seed and acorn germination rates at various constant temperatures were measured in Quercus canariensis, Q. coccifera, Q. ilex and Q. suber, using seeds sampled in 22 Tunisian woodlands. KEY RESULTS: Only faint differences were observed for desiccation sensitivity in the oak species studied. By contrast, the species differed significantly in sensitivity to freezing, germination rates at low temperature and base temperature. Quercus ilex and Q. canariensis, which occur at high elevations where frost events are frequent, showed the lowest freezing sensitivity. A significant correlation was found between hexose contents in the embryonic axis and freezing tolerance. Significant interspecific differences in the time for seeds to germinate and the time for the radicle to pierce the pericarp were observed. The ratio of pericarp mass to acorn mass differed significantly among the species and was negatively correlated with the acorn germination rate. Quercus coccifera, which is frequent in warm and arid environments, showed the highest acorn germination rate and synchrony. CONCLUSIONS: Seed lethal temperature, seed germination time at low temperatures, the ratio of pericarp mass to acorn mass and the embryonic axis hexose content appeared to be key functional traits that may influence the geographical ranges and ecological requirements of Mediterranean oaks in Tunisia.

Association mapping desiccation resistance within chromosomal inversions in the African malaria vector Anopheles gambiae.

Inversion polymorphisms are responsible for many ecologically important phenotypes and are often found under balancing selection. However, the same features that ensure their large role in local adaptation-especially reduced recombination between alternate arrangements-mean that uncovering the precise loci within inversions that control these phenotypes is unachievable using standard mapping approaches. Here, we take advantage of long-term balancing selection on a pair of inversions in the mosquito Anopheles gambiae to map desiccation tolerance via pool-GWAS. Two polymorphic inversions on chromosome 2 of this species (denoted 2La and 2Rb) are associated with arid and hot conditions in Africa and are maintained in spatially and temporally heterogeneous environments. After measuring thousands of wild-caught individuals for survival under desiccation stress, we used phenotypically extreme individuals homozygous for alternative arrangements at the 2La inversion to construct pools for whole-genome sequencing. Genomewide association mapping using these pools revealed dozens of significant SNPs within both 2La and 2Rb, many of which neighboured genes controlling ion channels or related functions. Our results point to the promise of similar approaches in systems with inversions maintained by balancing selection and provide a list of candidate genes underlying the specific phenotypes controlled by the two inversions studied here.

Investigation of mechanisms of bendiocarb resistance in Anopheles gambiae populations from the city of Yaoundé, Cameroon.

BACKGROUND: Resistance to the carbamate insecticide bendiocarb is emerging in Anopheles gambiae populations from the city of Yaoundé in Cameroon. However, the molecular basis of this resistance remains uncharacterized. The present study objective is to investigate mechanisms promoting resistance to bendiocarb in An. gambiae populations from Yaoundé. METHODS: The level of susceptibility of An. gambiae s.l. to bendiocarb 0.1 % was assessed from 2010 to 2013 using bioassays. Mosquitoes resistant to bendiocarb, unexposed and susceptible mosquitoes were screened for the presence of the Ace-1(R) mutation using TaqMan assays. Microarray analyses were performed to assess the pattern of genes differentially expressed between resistant, unexposed and susceptible. RESULTS: Bendiocarb resistance was more prevalent in mosquitoes originating from cultivated sites compared to those from polluted and unpolluted sites. Both An. gambiae and Anopheles coluzzii were found to display resistance to bendiocarb. No G119S mutation was detected suggesting that resistance was mainly metabolic. Microarray analysis revealed the over-expression of several cytochrome P450 s genes including cyp6z3, cyp6z1, cyp12f2, cyp6m3 and cyp6p4. Gene ontology (GO) enrichment analysis supported the detoxification role of cytochrome P450 s with several GO terms associated with P450 activity significantly enriched in resistant samples. Other detoxification genes included UDP-glucosyl transferases, glutathione-S transferases and ABC transporters. CONCLUSION: The study highlights the probable implication of metabolic mechanisms in bendiocarb resistance in An. gambiae populations from Yaoundé and stresses the need for further studies leading to functional validation of detoxification genes involved in this resistance.

Gene expression divergence between malaria vector sibling species Anopheles gambiae and An. coluzzii from rural and urban Yaoundé Cameroon.

Divergent selection based on aquatic larval ecology is a likely factor in the recent isolation of two broadly sympatric and morphologically identical African mosquito species, the malaria vectors Anopheles gambiae and An. coluzzii. Population-based genome scans have revealed numerous candidate regions of recent positive selection, but have provided few clues as to the genetic mechanisms underlying behavioural and physiological divergence between the two species, phenotypes which themselves remain obscure. To uncover possible genetic mechanisms, we compared global transcriptional profiles of natural and experimental populations using gene-based microarrays. Larvae were sampled as second and fourth instars from natural populations in and around the city of Yaoundé, capital of Cameroon, where the two species segregate along a gradient of urbanization. Functional enrichment analysis of differentially expressed genes revealed that An. coluzzii--the species that breeds in more stable, biotically complex and potentially polluted urban water bodies--overexpresses genes implicated in detoxification and immunity relative to An. gambiae, which breeds in more ephemeral and relatively depauperate pools and puddles in suburbs and rural areas. Moreover, our data suggest that such overexpression by An. coluzzii is not a transient result of induction by xenobiotics in the larval habitat, but an inherent and presumably adaptive response to repeatedly encountered environmental stressors. Finally, we find no significant overlap between the differentially expressed loci and previously identified genomic regions of recent positive selection, suggesting that transcriptome divergence is regulated by trans-acting factors rather than cis-acting elements.

Behavioural divergence of sympatric Anopheles funestus populations in Burkina Faso.

BACKGROUND: In Burkina Faso, two chromosomal forms of the malaria vector Anopheles funestus, Folonzo and Kiribina, are distinguished by contrasting frequencies of shared polymorphic chromosomal inversions. Sympatric and synchronous populations of Folonzo and Kiribina mate assortatively, as indicated by a significant deficit of heterokaryotypes, and genetic associations among inversions on independently segregating chromosome arms. The present study aimed to assess, by intensive longitudinal sampling, whether sympatric Folonzo and Kiribina populations are characterized by behavioural differences in key malaria vectorial parameters. METHODS: The study was conducted in two adjacent villages near Ouagadougou, in the dry savanna of central Burkina Faso. Mosquito adult resting behaviour of both forms was compared based on parallel indoor/outdoor collections across six breeding seasons; 8,235 fully karyotyped samples of half-gravid females were analysed in total. Additionally, indoor/outdoor human biting behaviour, host selection, and Plasmodium falciparum sporozoite rate was assessed and compared between chromosomal forms. RESULTS: The Kiribina form was numerically predominant in the area. However, the Folonzo form was significantly over-represented in indoor resting collections and showed stronger post-prandial endophily, while Kiribina predominated outdoors. Neither form was statistically distinguishable in human biting behaviour, and both were more likely to seek human blood meals indoors than outside. The human blood index and sporozoite rate were comparably high in both chromosomal forms in indoor collections (>89% and >8%, respectively). CONCLUSIONS: Both Kiribina and Folonzo chromosomal forms are formidable malaria vectors in Burkina Faso. However, the significantly greater tendency for the Kiribina form to rest outdoors despite its pronounced anthropophily suggests that uniform exposure of the overall An. funestus population to indoor-based vector control tools cannot be expected; Kiribina is more likely to evade indoor interventions and escape unharmed outdoors, reducing the efficacy of malaria control. Accordingly, more efficient methods to detect Kiribina and Folonzo, and a more complete understanding of their distribution and behaviour in Africa are advocated.

Exposure to disinfectants (soap or hydrogen peroxide) increases tolerance to permethrin in Anopheles gambiae populations from the city of Yaoundé, Cameroon.

BACKGROUND: The rapid expansion of insecticide resistance is limiting the efficiency of malaria vector control interventions. However, current knowledge of factors inducing pyrethroid resistance remains incomplete. In the present study, the role of selection at the larval stage by disinfectants, such as soap and hydrogen peroxide (H2O2), on adult mosquito resistance to permethrin was investigated. METHODS: Field Anopheles gambiae sensu lato larvae, were exposed to variable concentrations of soap and H2O2. Larvae surviving to acute toxicity assays after 24 hours were reared to the adult stage and exposed to permethrin. The susceptibility level of adults was compared to the untreated control group. The effect of soap or hydrogen peroxide selection on the length of larval development and emergence rate was assessed. RESULT: Larval bioassays analysis showed a more acute effect of hydrogen peroxide on mosquito larvae compared to soap. The regression lines describing the dose mortality profile showed higher mean and variance to hydrogen peroxide than to soap. The duration of larval development (<5 days) and adults emergence rates (1 to 77%) were shorter and lower compare to control. Anopheles gambiae s.l. larvae surviving to selection with either soap or hydrogen peroxide or both, produced adults who were up to eight-times more resistant to permethrin than mosquitoes from the untreated control group. CONCLUSION: The present study shows that selective pressure exerted by non-insecticidal compounds such as soap and hydrogen peroxide affect An. gambiae s.l. tolerance to pyrethroids. This requires further studies with regard to the adaptation of An. gambiae s.l. to polluted habitats across sub-Saharan Africa cities.

Adaptive divergence between incipient species of Anopheles gambiae increases resistance to Plasmodium.

The African malaria mosquito Anopheles gambiae is diversifying into ecotypes known as M and S forms. This process is thought to be promoted by adaptation to different larval habitats, but its genetic underpinnings remain elusive. To identify candidate targets of divergent natural selection in M and S, we performed genomewide scanning in paired population samples from Mali, followed by resequencing and genotyping from five locations in West, Central, and East Africa. Genome scans revealed a significant peak of M-S divergence on chromosome 3L, overlapping five known or suspected immune response genes. Resequencing implicated a selective target at or near the TEP1 gene, whose complement C3-like product has antiparasitic and antibacterial activity. Sequencing and allele-specific genotyping showed that an allelic variant of TEP1 has been swept to fixation in M samples from Mali and Burkina Faso and is spreading into neighboring Ghana, but is absent from M sampled in Cameroon, and from all sampled S populations. Sequence comparison demonstrates that this allele is related to, but distinct from, TEP1 alleles of known resistance phenotype. Experimental parasite infections of advanced mosquito intercrosses demonstrated a strong association between this TEP1 variant and resistance to both rodent malaria and the native human malaria parasite Plasmodium falciparum. Although malaria parasites may not be direct agents of pathogen-mediated selection at TEP1 in nature--where larvae may be the more vulnerable life stage--the process of adaptive divergence between M and S has potential consequences for malaria transmission.

Physiological correlates of ecological divergence along an urbanization gradient: differential tolerance to ammonia among molecular forms of the malaria mosquito Anopheles gambiae.

BACKGROUND: Limitations in the ability of organisms to tolerate environmental stressors affect their fundamental ecological niche and constrain their distribution to specific habitats. Evolution of tolerance, therefore, can engender ecological niche dynamics. Forest populations of the afro-tropical malaria mosquito Anopheles gambiae have been shown to adapt to historically unsuitable larval habitats polluted with decaying organic matter that are found in densely populated urban agglomerates of Cameroon. This process has resulted in niche expansion from rural to urban environments that is associated with cryptic speciation and ecological divergence of two evolutionarily significant units within this taxon, the molecular forms M and S, among which reproductive isolation is significant but still incomplete. Habitat segregation between the two forms results in a mosaic distribution of clinally parapatric patches, with the M form predominating in the centre of urban agglomerates and the S form in the surrounding rural localities. We hypothesized that development of tolerance to nitrogenous pollutants derived from the decomposition of organic matter, among which ammonia is the most toxic to aquatic organisms, may affect this pattern of distribution and process of niche expansion by the M form. RESULTS: Acute toxicity bioassays indicated that populations of the two molecular forms occurring at the extremes of an urbanization gradient in Yaounde, the capital of Cameroon, differed in their response to ammonia. The regression lines best describing the dose-mortality profile differed in the scale of the explanatory variable (ammonia concentration log-transformed for the S form and linear for the M form), and in slope (steeper for the S form and shallower for the M form). These features reflected differences in the frequency distribution of individual tolerance thresholds in the two populations as assessed by probit analysis, with the M form exhibiting a greater mean and variance compared to the S form. CONCLUSIONS: In agreement with expectations based on the pattern of habitat partitioning and exposure to ammonia in larval habitats in Yaounde, the M form showed greater tolerance to ammonia compared to the S form. This trait may be part of the physiological machinery allowing forest populations of the M form to colonize polluted larval habitats, which is at the heart of its niche expansion in densely populated human settlements in Cameroon.

Towards the invasion of wild and rural forested areas in Gabon (Central Africa) by the Asian tiger mosquito Aedes albopictus: Potential risks from the one health perspective.

BACKGROUND: Since its first record in urban areas of Central-Africa in the 2000s, the invasive mosquito, Aedes albopictus, has spread throughout the region, including in remote villages in forested areas, causing outbreaks of Aedes-borne diseases, such as dengue and chikungunya. Such invasion might enhance Ae. albopictus interactions with wild animals in forest ecosystems and favor the spillover of zoonotic arboviruses to humans. The aim of this study was to monitor Ae. albopictus spread in the wildlife reserve of La Lopé National Park (Gabon), and evaluate the magnitude of the rainforest ecosystem colonization. METHODOLOGY: From 2014 to 2018, we used ovitraps, larval surveys, BG-Sentinel traps, and human landing catches along an anthropization gradient from La Lopé village to the natural forest in the Park. CONCLUSIONS: We detected Ae. albopictus in gallery forest up to 15 km away from La Lopé village. However, Ae. albopictus was significantly more abundant at anthropogenic sites than in less anthropized areas. The number of eggs laid by Ae. albopictus decreased progressively with the distance from the forest fringe up to 200m inside the forest. Our results suggested that in forest ecosystems, high Ae. albopictus density is mainly observed at interfaces between anthropized and natural forested environments. Additionally, our data suggested that Ae. albopictus may act as a bridge vector of zoonotic pathogens between wild and anthropogenic compartments.

The role of transthoracic echocardiography in the diagnosis and management of acute type A aortic syndrome.

BACKGROUND: Transthoracic echocardiography (TTE) has been traditionally considered inadequate for the diagnosis of acute type A aortic syndrome (AAAS). In the last decade, high-resolution probes and harmonic imaging have been implemented in new echocardiographic systems. However, studies assessing the diagnostic accuracy of TTE for the identification of AAAS in large populations using modern ultrasound technology are lacking. METHODS: The diagnostic value of harmonic imaging TTE was assessed in 270 consecutive patients with suspected AAAS in whom TTE was the initial diagnostic test. RESULTS: Acute type A aortic syndrome was diagnosed in 67 patients and excluded in 203 patients (disease prevalence 25%). Sixty-two patients had a classic acute type A aortic dissection, and 5, an acute type A intramural hematoma. Image quality achieved was considered optimal in 244 patients (90%). In the whole study population, TTE had sensitivity, specificity, positive predictive value, and negative predictive value for the diagnosis of AAAS of 87%, 91%, 75%, and 95%, respectively. When evaluating only patients with optimal image quality, these values increased to 97%, 100%, 100%, and 99%, respectively. Forty-seven patients with clear-cut evidence of AAAS were transferred immediately to the operative room, where transesophageal echocardiography confirmed the diagnosis obtained by TTE in all patients. CONCLUSIONS: Transthoracic echocardiography is a useful imaging modality for the diagnosis of classic acute type A aortic dissection. It cannot be used as the sole screening technique for detecting AAAS, but in the light of the predictive values observed, patients with optimal image quality and clear-cut diagnosis of AAAS should proceed to the operative room, whereas in patients with negative or indeterminate studies, other imaging techniques are needed to refine the diagnosis.

Genetic isolation within the malaria mosquito Anopheles melas.

Anopheles melas is a brackish water-breeding member of the Anopheles gambiae complex that is distributed along the coast of West Africa and is a major malaria vector within its range. Because little is known about the population structure of this species, we analysed 15 microsatellite markers and 1161 bp of mtDNA in 11 A. melas populations collected throughout its range. Compared with its sibling species A. gambiae, A. melas populations have a high level of genetic differentiation between them, representing its patchy distribution due to its fragmented larval habitat that is associated with mangroves and salt marsh grass. Populations clustered into three distinct groups representing Western Africa, Southern Africa and Bioko Island populations that appear to be mostly isolated. Fixed differences in the mtDNA are present between all three clusters, and a Bayesian clustering analysis of the microsatellite data found no evidence for migration from mainland to Bioko Island populations, and little migration was evident between the Southern to the Western cluster. Surprisingly, mtDNA divergence between the three A. melas clusters is on par with levels of divergence between other species of the A. gambiae complex, and no support for monophyly was observed in a maximum-likelihood phylogenetic analysis. Finally, an approximate Bayesian analysis of microsatellite data indicates that Bioko Island A. melas populations were connected to the mainland populations in the past, but became isolated, presumably when sea levels rose after the last glaciation period (≥10 000-11 000 bp). This study has exposed species-level genetic divergence within A. melas and also has implications for control of this malaria vector.

Anopheles gambiae distribution and insecticide resistance in the cities of Douala and Yaoundé (Cameroon): influence of urban agriculture and pollution.

BACKGROUND: Urban malaria is becoming a major health priority across Africa. A study was undertaken to assess the importance of urban pollution and agriculture practice on the distribution and susceptibility to insecticide of malaria vectors in the two main cities in Cameroon. METHODS: Anopheline larval breeding sites were surveyed and water samples analysed monthly from October 2009 to December 2010. Parameters analysed included turbidity, pH, temperature, conductivity, sulfates, phosphates, nitrates, nitrites, ammonia, aluminium, alkalinity, iron, potassium, manganese, magnesium, magnesium hardness and total hardness. Characteristics of water bodies in urban areas were compared to rural areas and between urban sites. The level of susceptibility of Anopheles gambiae to 4% DDT, 0.75% permethrin, 0.05% deltamethrin, 0.1% bendiocarb and 5% malathion were compared between mosquitoes collected from polluted, non polluted and cultivated areas. RESULTS: A total of 1,546 breeding sites, 690 in Yaoundé and 856 in Douala, were sampled in the course of the study. Almost all measured parameters had a concentration of 2- to 100-fold higher in urban compare to rural breeding sites. No resistance to malathion was detected, but bendiocarb resistance was present in Yaounde. Very low mortality rates were observed following DDT or permethrin exposure, associated with high kdr frequencies. Mosquitoes collected in cultivated areas, exhibited the highest resistant levels. There was little difference in insecticide resistance or kdr allele frequency in mosquitoes collected from polluted versus non-polluted sites. CONCLUSION: The data confirm high selection pressure on mosquitoes originating from urban areas and suggest urban agriculture rather than pollution as the major factor driving resistance to insecticide.

Effects of insemination and blood-feeding on locomotor activity of wild-derived females of the malaria mosquito Anopheles coluzzii.

BACKGROUND: Behavioural shifts in the canonical location and timing of biting have been reported in natural populations of anopheline malaria vectors following the implementation of insecticide-based indoor vector control interventions. These modifications increase the likelihood of human-vector contact and allow mosquitoes to avoid insecticides, both conditions being favourable to residual transmission of the malarial parasites. The biting behaviour of mosquitoes follows rhythms that are under the control of biological clocks and environmental conditions, modulated by physiological states. In this work we explore modifications of spontaneous locomotor activity expressed by mosquitoes in different physiological states to highlight phenotypic variability associated to circadian control that may contribute to explain residual transmission in the field. METHODS: The F10 generation progeny of field-collected Anopheles coluzzii from southwestern Burkina Faso was tested using an automated recording apparatus (Locomotor Activity Monitor, TriKinetics Inc.) under LD 12:12 or DD light regimens in laboratory-controlled conditions. Activity recordings of each test were carried out for a week with 6-day-old females belonging to four experimental treatments, representing factorial combinations of two physiological variables: insemination status (virgin vs inseminated) and gonotrophic status (glucose fed vs blood fed). Chronobiological features of rhythmicity in locomotor activity were explored using periodograms, diversity indices, and generalized linear mixed modelling. RESULTS: The average strength of activity, onset of activity, and acrophase were modulated by both nutritional and insemination status as well as by the light regimen. Inseminated females showed a significant excess of arrhythmic activity under DD. When rhythmicity was observed in DD, females displayed sustained activity also during the subjective day. CONCLUSIONS: Insemination and gonotrophic status influence the underlying light and circadian control of chronobiological features of locomotor activity. Overrepresentation of arrhythmic chronotypes as well as the sustained activity of inseminated females during the subjective day under DD conditions suggests potential activity of natural populations of A. coluzzii during daytime under dim conditions, with implications for residual transmission of malarial parasites.

A PCR-RFLP method for genotyping of inversion 2Rc in Anopheles coluzzii.

BACKGROUND: Genotyping of polymorphic chromosomal inversions in malaria vectors such as An. coluzzii Coetzee & Wilkerson is important, both because they cause cryptic population structure that can mislead vector analysis and control and because they influence epidemiologically relevant eco-phenotypes. The conventional cytogenetic method of genotyping is an impediment because it is labor intensive, requires specialized training, and can be applied only to one gender and developmental stage. Here, we circumvent these limitations by developing a simple and rapid molecular method of genotyping inversion 2Rc in An. coluzzii that is both economical and field-friendly. This inversion is strongly implicated in temporal and spatial adaptations to climatic and ecological variation, particularly aridity. METHODS: Using a set of tag single-nucleotide polymorphisms (SNPs) strongly correlated with inversion orientation, we identified those that overlapped restriction enzyme recognition sites and developed four polymerase chain reaction (PCR) restriction fragment length polymorphism (RFLP) assays that distinguish alternative allelic states at the tag SNPs. We assessed the performance of these assays using mosquito population samples from Burkina Faso that had been cytogenetically karyotyped as well as genotyped, using two complementary high-throughput molecular methods based on tag SNPs. Further validation was performed using mosquito population samples from additional West African (Benin, Mali, Senegal) and Central African (Cameroon) countries. RESULTS: Of four assays tested, two were concordant with the 2Rc cytogenetic karyotype > 90% of the time in all samples. We recommend that these two assays be employed in tandem for reliable genotyping. By accepting only those genotypic assignments where both assays agree, > 99% of assignments are expected to be accurate. CONCLUSIONS: We have developed tandem PCR-RFLP assays for the accurate genotyping of inversion 2Rc in An. coluzzii. Because this approach is simple, inexpensive, and requires only basic molecular biology equipment, it is widely accessible. These provide a crucial tool for probing the molecular basis of eco-phenotypes relevant to malaria epidemiology and vector control.

Ecological plasticity to ions concentration determines genetic response and dominance of Anopheles coluzzii larvae in urban coastal habitats of Central Africa.

In Central Africa, the malaria vector Anopheles coluzzii is predominant in urban and coastal habitats. However, little is known about the environmental factors that may be involved in this process. Here, we performed an analysis of 28 physicochemical characteristics of 59 breeding sites across 5 urban and rural sites in coastal areas of Central Africa. We then modelled the relative frequency of An. coluzzii larvae to these physicochemical parameters in order to investigate environmental patterns. Then, we assessed the expression variation of 10 candidate genes in An. coluzzii, previously incriminated with insecticide resistance and osmoregulation in urban settings. Our results confirmed the ecological plasticity of An. coluzzii larvae to breed in a large range of aquatic conditions and its predominance in breeding sites rich in ions. Gene expression patterns were comparable between urban and rural habitats, suggesting a broad response to ions concentrations of whatever origin. Altogether, An. coluzzii exhibits a plastic response to occupy both coastal and urban habitats. This entails important consequences for malaria control in the context of the rapid urban expansion in Africa in the coming years.

Genetic association of physically unlinked islands of genomic divergence in incipient species of Anopheles gambiae.

Previous efforts to uncover the genetic underpinnings of ongoing ecological speciation of the M and S forms of the African malaria vector Anopheles gambiae revealed two centromere-proximal islands of genetic divergence on X and chromosome 2. Under the assumption of considerable ongoing gene flow between M and S, these persistently divergent genomic islands were widely considered to be 'speciation islands'. In the course of microarray-based divergence mapping, we discovered a third centromere-associated island of divergence on chromosome 3, which was validated by targeted re-sequencing. To test for genetic association between the divergence islands on all three chromosomes, SNP-based assays were applied in four natural populations of M and S spanning West, Central and East Africa. Genotyping of 517 female M and S mosquitoes revealed nearly complete linkage disequilibrium between the centromeres of the three independently assorting chromosomes. These results suggest that despite the potential for inter-form gene flow through hybridization, actual (realized) gene flow between M and S may be substantially less than commonly assumed and may not explain most shared variation. Moreover, the possibility of very low gene flow calls into question whether diverged pericentromeric regions-characterized by reduced levels of variation and recombination-are in fact instrumental rather than merely incidental to the speciation process.

Breakpoint structure of the Anopheles gambiae 2Rb chromosomal inversion.

BACKGROUND: Alternative arrangements of chromosome 2 inversions in Anopheles gambiae are important sources of population structure, and are associated with adaptation to environmental heterogeneity. The forces responsible for their origin and maintenance are incompletely understood. Molecular characterization of inversion breakpoints provides insight into how they arose, and provides the basis for development of molecular karyotyping methods useful in future studies. METHODS: Sequence comparison of regions near the cytological breakpoints of 2Rb allowed the molecular delineation of breakpoint boundaries. Comparisons were made between the standard 2R+b arrangement in the An. gambiae PEST reference genome and the inverted 2Rb arrangements in the An. gambiae M and S genome assemblies. Sequence differences between alternative 2Rb arrangements were exploited in the design of a PCR diagnostic assay, which was evaluated against the known chromosomal banding pattern of laboratory colonies and field-collected samples from Mali and Cameroon. RESULTS: The breakpoints of the 7.55 Mb 2Rb inversion are flanked by extensive runs of the same short (72 bp) tandemly organized sequence, which was likely responsible for chromosomal breakage and rearrangement. Application of the molecular diagnostic assay suggested that 2Rb has a single common origin in An. gambiae and its sibling species, Anopheles arabiensis, and also that the standard arrangement (2R+b) may have arisen twice through breakpoint reuse. The molecular diagnostic was reliable when applied to laboratory colonies, but its accuracy was lower in natural populations. CONCLUSIONS: The complex repetitive sequence flanking the 2Rb breakpoint region may be prone to structural and sequence-level instability. The 2Rb molecular diagnostic has immediate application in studies based on laboratory colonies, but its usefulness in natural populations awaits development of complementary molecular tools.