The last bastion? X chromosome genotyping of Anopheles gambiae species pair males from a hybrid zone reveals complex recombination within the major candidate 'genomic island of speciation'.

Speciation with gene flow may be aided by reduced recombination helping to build linkage between genes involved in the early stages of reproductive isolation. Reduced recombination on chromosome X has been implicated in speciation within the Anopheles gambiae complex, species of which represent the major Afrotropical malaria vectors. The most recently diverged, morphologically indistinguishable, species pair, A. gambiae and Anopheles coluzzii, ubiquitously displays a 'genomic island of divergence' spanning over 4 Mb from chromosome X centromere, which represents a particularly promising candidate region for reproductive isolation genes, in addition to containing the diagnostic markers used to distinguish the species. Very low recombination makes the island intractable for experimental recombination studies, but an extreme hybrid zone in Guinea Bissau offers the opportunity for natural investigation of X-island recombination. SNP analysis of chromosome X hemizygous males revealed: (i) strong divergence in the X-island despite a lack of autosomal divergence; (ii) individuals with multiple-recombinant genotypes, including likely double crossovers and localized gene conversion; (iii) recombination-driven discontinuity both within and between the molecular species markers, suggesting that the utility of the diagnostics is undermined under high hybridization. The largely, but incompletely protected nature of the X centromeric genomic island is consistent with a primary candidate area for accumulation of adaptive variants driving speciation with gene flow, while permitting some selective shuffling and removal of genetic variation.

Evaluation of a protocol for remote identification of mosquito vector species reveals BG-Sentinel trap as an efficient tool for Anopheles gambiae outdoor collection in Burkina Faso.

BACKGROUND: Feasibility and costs of monitoring efforts aimed to monitor mosquito species are strictly dependent on the presence of skilled entomologists directly in the field. However, in several contexts this is not possible or easy to organize, thus limiting the possibility to obtain crucial information on presence/abundance of potential disease vectors and of new invasive species. Digital imaging approaches could be extremely useful in the frame of medical entomology to overcome this limit. This work describes a surveillance approach to collect and morphologically identify host-seeking malaria vectors based on remote transmission of digital images of specimens collected with ad hoc modified traps. METHODS: CDC light trap (CDC) and the BG-Sentinel trap (BG), both baited with BG-lure and CO2, were modified in order to have collected mosquitoes immobilized on a bi-dimensional surface. The performance of the two traps in the field was comparatively tested by Latin-square experiments in two villages of Burkina Faso under low and high mosquito densities. The efficiency of identifications based the inspection of digital images versus microscopic identifications of collected specimens was compared. RESULTS: A total of 1,519 mosquitoes belonging to 16 species were collected, of which 88.5% were microscopically identified as Anopheles gambiae s.l. (mainly Anopheles coluzzii, 85.7%). During dry season BG collected 15 times more females than CDC outdoors, whereas indoors the BG collected 0.4 times less than CDC. During rainy season the ratio BG/CDC was 6.4 and 0.7 outdoors and indoors, respectively. The efficiency of digital images versus microscopic identifications of collected specimens was 97.9%, 95.6% and 81.5% for Culicidae, Anophelinae and An. gambiae s.l., respectively. CONCLUSIONS: Results strongly encourage the use of BG-trap for collecting host-seeking An. gambiae particularly in the outdoor environment, providing new perspectives to the challenge of collecting this fraction of the biting population, whose epidemiological relevance is increasing due to the success of large-scale implementation of indoor malaria vector control strategies. Moreover, results show that the transmission of digital images of specimens collected by the ad hoc modified host-seeking traps efficiently allows identification of malaria vectors, thus opening the perspective to easily carry out mosquito monitoring also in the absence of entomologists directly in the field.

Remarkable diversity of intron-1 of the para voltage-gated sodium channel gene in an Anopheles gambiae/Anopheles coluzzii hybrid zone.

BACKGROUND: Genomic differentiation between Anopheles gambiae and Anopheles coluzzii--the major malaria vectors in sub-Saharan Africa--is localized into large "islands" toward the centromeres of chromosome-X and the two autosomes. Linkage disequilibrium between these genomic islands was first detected between species-specific polymorphisms within ribosomal DNA genes (IGS-rDNA) on the X-chromosome and a single variant at position 702 of intron 1 (Int-1702) of the para Voltage-Gated Sodium Channel (VGSC) gene on chromosome arm 2 L. Intron-1 sequence data from West and Central Africa revealed two clearly distinct and species-specific haplogroups, each characterized by very low polymorphism, which has been attributed to a selective sweep. The aim of this study was to analyse Int-1 sequence diversity in A. gambiae and A. coluzzii populations from the Far-West of their range, in order to assess whether this selective-sweep signature could persist in a zone of high interspecific hybridization. METHODS: A 531 bp region of VGSC Int-1 was sequenced in 21 A. coluzzii, 31 A. gambiae, and 12 hybrids from The Gambia and Guinea Bissau, located within the Far-West geographical region, and in 53 A. gambiae s.l. samples from the rest of the range. RESULTS: Far-West samples exhibit dramatic Int-1 polymorphism, far higher within each country than observed throughout the rest of the species range. Moreover, patterning of haplotypes within A. coluzzii confirms previous evidence of a macro-geographic subdivision into a West and a Central African genetic cluster, and reveals a possible genetic distinction of A. coluzzii populations from the Far-West. CONCLUSIONS: The results suggest a relaxation of selective pressures acting across the VGSC gene region in the hybrid zone. Genetic differentiation in the Far-West could be attributable to a founder effect within A. coluzzii, with subsequent extensive gene flow with secondarily-colonizing A. gambiae, potentially yielding a novel insight on the dynamic processes impacting genetic divergence of these key malaria vectors.

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.

Wide cross-reactivity between Anopheles gambiae and Anopheles funestus SG6 salivary proteins supports exploitation of gSG6 as a marker of human exposure to major malaria vectors in tropical Africa.

BACKGROUND: The Anopheles gambiae gSG6 is an anopheline-specific salivary protein which helps female mosquitoes to efficiently feed on blood. Besides its role in haematophagy, gSG6 is immunogenic and elicits in exposed individuals an IgG response, which may be used as indicator of exposure to the main African malaria vector A. gambiae. However, malaria transmission in tropical Africa is sustained by three main vectors (A. gambiae, Anopheles arabiensis and Anopheles funestus) and a general marker, reflecting exposure to at least these three species, would be especially valuable. The SG6 protein is highly conserved within the A. gambiae species complex whereas the A. funestus homologue, fSG6, is more divergent (80% identity with gSG6). The aim of this study was to evaluate cross-reactivity of human sera to gSG6 and fSG6. METHODS: The A. funestus SG6 protein was expressed/purified and the humoral response to gSG6, fSG6 and a combination of the two antigens was compared in a population from a malaria hyperendemic area of Burkina Faso where both vectors were present, although with a large A. gambiae prevalence (>75%). Sera collected at the beginning and at the end of the high transmission/rainy season, as well as during the following low transmission/dry season, were analysed. RESULTS: According to previous observations, both anti-SG6 IgG level and prevalence decreased during the low transmission/dry season and showed a typical age-dependent pattern. No significant difference in the response to the two antigens was found, although their combined use yielded in most cases higher IgG level. CONCLUSIONS: Comparative analysis of gSG6 and fSG6 immunogenicity to humans suggests the occurrence of a wide cross-reactivity, even though the two proteins carry species-specific epitopes. This study supports the use of gSG6 as reliable indicator of exposure to the three main African malaria vectors, a marker which may be useful to monitor malaria transmission and evaluate vector control measures, especially in conditions of low malaria transmission and/or reduced vector density. The Anopheles stephensi SG6 protein also shares 80% identity with gSG6, suggesting the attractive possibility that the A. gambiae protein may also be useful to assess human exposure to several Asian malaria vectors.

Comparative analyses reveal discrepancies among results of commonly used methods for Anopheles gambiaemolecular form identification.

BACKGROUND: Anopheles gambiae M and S molecular forms, the major malaria vectors in the Afro-tropical region, are ongoing a process of ecological diversification and adaptive lineage splitting, which is affecting malaria transmission and vector control strategies in West Africa. These two incipient species are defined on the basis of single nucleotide differences in the IGS and ITS regions of multicopy rDNA located on the X-chromosome. A number of PCR and PCR-RFLP approaches based on form-specific SNPs in the IGS region are used for M and S identification. Moreover, a PCR-method to detect the M-specific insertion of a short interspersed transposable element (SINE200) has recently been introduced as an alternative identification approach. However, a large-scale comparative analysis of four widely used PCR or PCR-RFLP genotyping methods for M and S identification was never carried out to evaluate whether they could be used interchangeably, as commonly assumed. RESULTS: The genotyping of more than 400 A. gambiae specimens from nine African countries, and the sequencing of the IGS-amplicon of 115 of them, highlighted discrepancies among results obtained by the different approaches due to different kinds of biases, which may result in an overestimation of MS putative hybrids, as follows: i) incorrect match of M and S specific primers used in the allele specific-PCR approach; ii) presence of polymorphisms in the recognition sequence of restriction enzymes used in the PCR-RFLP approaches; iii) incomplete cleavage during the restriction reactions; iv) presence of different copy numbers of M and S-specific IGS-arrays in single individuals in areas of secondary contact between the two forms. CONCLUSIONS: The results reveal that the PCR and PCR-RFLP approaches most commonly utilized to identify A. gambiae M and S forms are not fully interchangeable as usually assumed, and highlight limits of the actual definition of the two molecular forms, which might not fully correspond to the two A. gambiae incipient species in their entire geographical range. These limits are discussed and operational suggestions on the choice of the most convenient method for large-scale M- and S-form identification are provided, also taking into consideration technical aspects related to the epidemiological characteristics of different study areas.

Highly specific PCR-RFLP assays for karyotyping the widespread 2Rb inversion in malaria vectors of the Anopheles gambiae complex.

BACKGROUND: Chromosomal inversion polymorphisms play a role in adaptation to heterogeneous environments. Inversion polymorphisms are implicated in the very high ecological flexibility of the three main malaria vector species of the Afrotropical Anopheles gambiae complex, facilitating the exploitation of anthropogenic environmental modifications and promoting a strong association with humans. In addition to extending the species' spatial and temporal distribution, inversions are associated with epidemiologically relevant mosquito behavior and physiology, underscoring their medical importance. We here present novel PCR-RFLP based assays strongly predictive of genotype for the cosmopolitan 2Rb inversion in An. coluzzii and An. gambiae, a development which overcomes the numerous constraints inherent to traditional cytological karyotyping. METHODS: We designed PCR-RFLP genotyping assays based on tag SNPs previously computationally identified as strongly predictive (> 95%) of 2Rb genotype. We targeted those tags whose alternative allelic states destroyed or created the recognition site of a commercially available restriction enzyme, and designed assays with distinctive cleavage profiles for each inversion genotype. The assays were validated on 251 An. coluzzii and 451 An. gambiae cytologically karyotyped specimens from nine countries across Africa and one An. coluzzii laboratory colony. RESULTS: For three tag SNPs, PCR-RFLP assays (denoted DraIII, MspAI, and TatI) reliably produced robust amplicons and clearly distinguishable electrophoretic profiles for all three inversion genotypes. Results obtained with the DraIII assay are ≥ 95% concordant with cytogenetic assignments in both species, while MspAI and TatI assays produce patterns highly concordant with cytogenetic assignments only in An. coluzzii or An. gambiae, respectively. Joint application of species-appropriate pairs of assays increased the concordance levels to > 99% in An. coluzzii and 98% in An. gambiae. Potential sources of discordance (e.g. imperfect association between tag and inversion, allelic dropout, additional polymorphisms in the restriction target site, incomplete or failed restriction digestion) are discussed. CONCLUSIONS: The availability of highly specific, cost effective and accessible molecular assays for genotyping 2Rb in An. gambiae and An. coluzzii allows karyotyping of both sexes and all developmental stages. These novel tools will accelerate deeper investigations into the role of this ecologically and epidemiologically important chromosomal inversion in vector biology.

In Silico Karyotyping of Chromosomally Polymorphic Malaria Mosquitoes in the Anopheles gambiae Complex.

Chromosomal inversion polymorphisms play an important role in adaptation to environmental heterogeneities. For mosquito species in the Anopheles gambiae complex that are significant vectors of human malaria, paracentric inversion polymorphisms are abundant and are associated with ecologically and epidemiologically important phenotypes. Improved understanding of these traits relies on determining mosquito karyotype, which currently depends upon laborious cytogenetic methods whose application is limited both by the requirement for specialized expertise and for properly preserved adult females at specific gonotrophic stages. To overcome this limitation, we developed sets of tag single nucleotide polymorphisms (SNPs) inside inversions whose biallelic genotype is strongly correlated with inversion genotype. We leveraged 1,347 fully sequenced An. gambiae and Anopheles coluzzii genomes in the Ag1000G database of natural variation. Beginning with principal components analysis (PCA) of population samples, applied to windows of the genome containing individual chromosomal rearrangements, we classified samples into three inversion genotypes, distinguishing homozygous inverted and homozygous uninverted groups by inclusion of the small subset of specimens in Ag1000G that are associated with cytogenetic metadata. We then assessed the correlation between candidate tag SNP genotypes and PCA-based inversion genotypes in our training sets, selecting those candidates with >80% agreement. Our initial tests both in held-back validation samples from Ag1000G and in data independent of Ag1000G suggest that when used for in silico inversion genotyping of sequenced mosquitoes, these tags perform better than traditional cytogenetics, even for specimens where only a small subset of the tag SNPs can be successfully ascertained.

Fine-Mapping Complex Inversion Breakpoints and Investigating Somatic Pairing in the Anopheles gambiae Species Complex Using Proximity-Ligation Sequencing.

Chromosomal inversions are fundamental drivers of genome evolution. In the main Afrotropical malaria vector species, belonging to the Anopheles gambiae species complex, inversions play an important role in local adaptation and have a rich history of cytological study. Despite the importance and ubiquity of some chromosomal inversions across the species complex, inversion breakpoints are often challenging to map molecularly due to the presence of large repetitive regions. Here, we develop an approach that uses Hi-C sequencing data to molecularly fine-map the breakpoints of inversions. We demonstrate that this approach is robust and likely to be widely applicable for both identification and fine-mapping inversion breakpoints in species whose inversions have heretofore been challenging to characterize. We apply our method to interrogate the previously unknown inversion breakpoints of 2Rbc and 2Rd in An. coluzzii We found that inversion breakpoints occur in large repetitive regions, and, strikingly, among three inversions analyzed, two breakpoints appear to be reused in two separate inversions. These breakpoint-adjacent regions are strongly enriched for the presence of a 30 bp satellite repeat sequence. Because low frequency inversion breakpoints are not correlated with genomic regions containing this satellite, we suggest that interrupting this particular repeat may result in arrangements with higher relative fitness. Additionally, we use heterozygous individuals to quantitatively investigate the impacts of somatic pairing in the regions immediately surrounding inversion breakpoints. Finally, we discuss important considerations for possible applications of this approach for inversion breakpoint identification in a range of organisms.

Phenotypic and genotypic pyrethroid resistance of Aedes albopictus, with focus on the 2017 chikungunya outbreak in Italy.

BACKGROUND: The highly invasive mosquito species Aedes albopictus has become a major health concern in temperate areas due to its role as vector of exotic arboviruses. Pyrethroid insecticides represent the main tools for limiting the circulation of such mosquito-borne viruses. The present work aim to extend previous reports on phenotypic pyrethroid-resistance in European Ae. albopictus, to identify its genetic basis and to monitor the geographical distribution of resistant genotypes, with a particular focus on sites experiencing the 2017 chikungunya outbreak in Italy. RESULTS: Bioassays, performed according to World Health Organization protocols, showed full susceptibility to deltamethrin (concentration = 0.05%) and varying levels of resistance to permethrin (0.75%) and/or α-cypermethrin (0.05%) across Italy, with highest levels in the core of the 2017 chikungunya outbreak. Partial genotyping of the VSSC gene revealed widespread distribution of V1016G mutation and confirmed its association with pyrethroid resistance. CONCLUSION: The results obtained show that the condition for the spread of pyrethroid resistance in Ae. albopictus in Europe exists under strong selective pressure due to intensive insecticide spraying to control exotic arbovirus outbreak or high levels of nuisance. The results draw attention to the need for an evidence-based implementation of mosquito nuisance control, taking insecticide resistance management into consideration. © 2019 Society of Chemical Industry.

Chromosomal plasticity and evolutionary potential in the malaria vector Anopheles gambiae sensu stricto: insights from three decades of rare paracentric inversions.

BACKGROUND: In the Anopheles gambiae complex, paracentric chromosomal inversions are non-randomly distributed along the complement: 18/31 (58%) of common polymorphic inversions are on chromosome arm 2R, which represents only approximately 30% of the complement. Moreover, in An. gambiae sensu stricto, 6/7 common polymorphic inversions occur on 2R. Most of these inversions are considered markers of ecological adaptation that increase the fitness of the carriers of alternative karyotypes in contrasting habitats. However, little is known about the evolutionary forces responsible for their origin and subsequent establishment in field populations. RESULTS: Here, we present data on 82 previously undescribed rare chromosomal inversions (RCIs) recorded during extensive field sampling in 16 African countries over a 30 year period, which may shed light on the dynamics of chromosomal plasticity in An. gambiae. We analyzed breakpoint distribution, length, and geographic distribution of RCIs, and compared these measures to those of the common inversions. We found that RCIs, like common inversions, are disproportionately clustered on 2R, which may indicate that this arm is especially prone to breakages. However, contrasting patterns were observed between the geographic distribution of common inversions and RCIs. RCIs were equally frequent across biomes and on both sides of the Great Rift Valley (GRV), whereas common inversions predominated in arid ecological settings and west of the GRV. Moreover, the distribution of RCI lengths followed a random pattern while common inversions were significantly less frequent at shorter lengths. CONCLUSION: Because 17/82 (21%) RCIs were found repeatedly at very low frequencies - at the same sampling location in different years and/or in different sampling locations - we suggest that RCIs are subject mainly to drift under unperturbed ecological conditions. Nevertheless, RCIs may represent an important reservoir of genetic variation for An. gambiae in response to environmental changes, further testifying to the considerable evolutionary potential hidden within this pan-African malaria vector.

Distribution and chromosomal characterization of the Anopheles gambiae complex in Angola.

Mosquitoes of the Anopheles gambiae complex (N = 1,336) were sampled (2001-2005) across Angola to identify taxa, study inversion polymorphisms, and detect the circumsporozoite protein of Plasmodium falciparum. Anopheles gambiae s.s. was found in all sites; it was characterized as M-form in localities of the tropical dry and semi-desertic belts, whereas the S-form was predominant in comparatively more humid and less anthropized sites. Both forms were characterized by low degrees of chromosomal polymorphism based solely on the 2La inversion, a pattern usually associated with An. gambiae populations from forested, humid, and derived savanna areas. Unexpectedly, this pattern was also observed in M-form populations collected in dry/pre-desertic areas, where this form largely predominates over An. arabiensis, which was also detected in central/inland sites. Anopheles melas was found in northern coastal sites. Three of 534 An. gambiae s.s. were positive for P. falciparum CS-protein, whereas none of the 105 An. melas were positive.

Anopheles gambiae complex along The Gambia river, with particular reference to the molecular forms of An. gambiae s.s.

BACKGROUND: The geographic and temporal distribution of M and S molecular forms of the major Afrotropical malaria vector species Anopheles gambiae s.s. at the western extreme of their range of distribution has never been investigated in detail. MATERIALS AND METHODS: Collections of indoor-resting An. gambiae s.l. females were carried out along a ca. 400 km west to east transect following the River Gambia from the western coastal region of The Gambia to south-eastern Senegal during 2005 end of rainy season/early dry season and the 2006 rainy season. Specimens were identified to species and molecular forms by PCR-RFLP and the origin of blood-meal of fed females was determined by ELISA test. RESULTS: Over 4,000 An. gambiae s.l. adult females were collected and identified, 1,041 and 3,038 in 2005 and 2006, respectively. M-form was mainly found in sympatry with Anopheles melas and S-form in the western part of the transect, and with Anopheles arabiensis in the central part. S-form was found to prevail in rural Sudan-Guinean savannah areas of Eastern Senegal, in sympatry with An. arabiensis. Anopheles melas and An. arabiensis relative frequencies were generally lower in the rainy season samples, when An. gambiae s.s. was prevailing. No large seasonal fluctuations were observed for M and S-forms. In areas where both M and S were recorded, the frequency of hybrids between them ranged from to 0.6% to 7%. DISCUSSION: The observed pattern of taxa distribution supports the hypothesis of a better adaptation of M-form to areas characterized by water-retaining alluvial deposits along the Gambia River, characterized by marshy vegetation, mangrove woods and rice cultivations. In contrast, the S-form seems to be better adapted to free-draining soil, covered with open woodland savannah or farmland, rich in temporary larval breeding sites characterizing mainly the eastern part of the transect, where the environmental impact of the Gambia River is much less profound and agricultural activities are mainly rain-dependent. Very interestingly, the observed frequency of hybridization between the molecular forms along the whole transect was much higher than has been reported so far for other areas. CONCLUSION: The results support a bionomic divergence between the M and S-forms, and suggest that the western extreme of An. gambiae s.s. geographical distribution may represent an area of higher-than-expected hybridization between the two molecular forms.

Distribution of knock-down resistance mutations in Anopheles gambiae molecular forms in west and west-central Africa.

BACKGROUND: Knock-down resistance (kdr) to DDT and pyrethroids in the major Afrotropical vector species, Anopheles gambiae sensu stricto, is associated with two alternative point mutations at amino acid position 1014 of the voltage-gated sodium channel gene, resulting in either a leucine-phenylalanine (L1014F), or a leucine-serine (L1014S) substitution. In An. gambiae S-form populations, the former mutation appears to be widespread in west Africa and has been recently reported from Uganda, while the latter, originally recorded in Kenya, has been recently found in Gabon, Cameroon and Equatorial Guinea. In M-form populations surveyed to date, only the L1014F mutation has been found, although less widespread and at lower frequencies than in sympatric S-form populations. METHODS: Anopheles gambiae M- and S-form specimens from 19 sites from 11 west and west-central African countries were identified to molecular form and genotyped at the kdr locus either by Hot Oligonucleotide Ligation Assay (HOLA) or allele-specific PCR (AS-PCR). RESULTS: The kdr genotype was determined for about 1,000 An. gambiae specimens. The L1014F allele was found at frequencies ranging from 6% to 100% in all S-form samples (N = 628), with the exception of two samples from Angola, where it was absent, and coexisted with the L1014S allele in samples from Cameroon, Gabon and north-western Angola. The L1014F allele was present in M-form samples (N = 354) from Benin, Nigeria, and Cameroon, where both M- and S-forms were sympatric. CONCLUSION: The results represent the most comprehensive effort to analyse the overall distribution of the L1014F and L1014S mutations in An. gambiae molecular forms, and will serve as baseline data for resistance monitoring. The overall picture shows that the emergence and spread of kdr alleles in An. gambiae is a dynamic process and that there is marked intra- and inter-form heterogeneity in resistance allele frequencies. Further studies are needed to determine: i) the importance of selection pressure exerted by both agricultural and public health use of pyrethroid insecticides, ii) the phenotypic effects, particularly when the two mutations co-occur; and iii) the epidemiological importance of kdr for both pyrethroid- and DDT-based malaria control operations, particularly if/when the two insecticides are to be used in concert.

Molecular karyotyping of the 2La inversion in Anopheles gambiae.

The African malaria vector Anopheles gambiae is polymorphic for alternative arrangements on the left arm of chromosome 2 (2La and 2L+(a)) that are non-randomly distributed with respect to degree of aridity. Detailed studies on the ecological role of inversion 2La have been hindered by the technical demands of traditional karyotype analysis and by sex- and stage-specific limitations on the availability of polytene chromosomes favorable for analysis. Recent molecular characterization of both inversion breakpoints presented the opportunity to develop a polymerase chain reaction (PCR)-based method for karyotype analysis. Here we report the development of this molecular diagnostic assay and the results of extensive field validation. When tested on 765 An. gambiae specimens sampled across Africa, the molecular approach compared favorably with traditional cytologic methods, correctly scoring > 94% of these specimens. By providing ready access to the 2La karyotype, this tool lays groundwork for future studies of the ecological genomics of this medically important species.

PCR-based karyotyping of Anopheles gambiae inversion 2Rj identifies the BAMAKO chromosomal form.

BACKGROUND: The malaria vector Anopheles gambiae is polymorphic for chromosomal inversions on the right arm of chromosome 2 that segregate nonrandomly between assortatively mating populations in West Africa. One such inversion, 2Rj, is associated with the BAMAKO chromosomal form endemic to southern Mali and northern Guinea Conakry near the Niger River. Although it exploits a unique ecology and both molecular and chromosomal data suggest reduced gene flow between BAMAKO and other A. gambiae populations, no molecular markers exist to identify this form. METHODS: To facilitate study of the BAMAKO form, a PCR assay for molecular karyotyping of 2Rj was developed based on sequences at the breakpoint junctions. The assay was extensively validated using more than 700 field specimens whose karyotypes were determined in parallel by cytogenetic and molecular methods. As inversion 2Rj also occurs in SAVANNA populations outside the geographic range of BAMAKO, samples were tested from Senegal, Cameroon and western Guinea Conakry as well as from Mali. RESULTS: In southern Mali, where 2Rj polymorphism in SAVANNA populations was very low and most of the 2Rj homozygotes were found in BAMAKO karyotypes, the molecular and cytogenetic methods were almost perfectly congruent. Elsewhere agreement between the methods was much poorer, as the molecular assay frequently misclassified 2Rj heterozygotes as 2R+j standard homozygotes. CONCLUSION: Molecular karyotyping of 2Rj is robust and accurate on 2R+j standard and 2Rj inverted homozygotes. Therefore, the proposed approach overcomes the lack of a rapid tool for identifying the BAMAKO form across developmental stages and sexes, and opens new perspectives for the study of BAMAKO ecology and behaviour. On the other hand, the method should not be applied for molecular karyotyping of j-carriers within the SAVANNA chromosomal form.

Massive introgression drives species radiation at the range limit of Anopheles gambiae.

Impacts of introgressive hybridisation may range from genomic erosion and species collapse to rapid adaptation and speciation but opportunities to study these dynamics are rare. We investigated the extent, causes and consequences of a hybrid zone between Anopheles coluzzii and Anopheles gambiae in Guinea-Bissau, where high hybridisation rates appear to be stable at least since the 1990s. Anopheles gambiae was genetically partitioned into inland and coastal subpopulations, separated by a central region dominated by A. coluzzii. Surprisingly, whole genome sequencing revealed that the coastal region harbours a hybrid form characterised by an A. gambiae-like sex chromosome and massive introgression of A. coluzzii autosomal alleles. Local selection on chromosomal inversions may play a role in this process, suggesting potential for spatiotemporal stability of the coastal hybrid form and providing resilience against introgression of medically-important loci and traits, found to be more prevalent in inland A. gambiae.

Dissecting functional components of reproductive isolation among closely related sympatric species of the Anopheles gambiae complex.

Explaining how and why reproductive isolation evolves and determining which forms of reproductive isolation have the largest impact on the process of population divergence are major goals in the study of speciation. By studying recent adaptive radiations in incompletely isolated taxa, it is possible to identify barriers involved at early divergence before other confounding barriers emerge after speciation is complete. Sibling species of the Anopheles gambiae complex offer opportunities to provide insights into speciation mechanisms. Here, we studied patterns of reproductive isolation among three taxa, Anopheles coluzzii, An. gambiae s.s. and Anopheles arabiensis, to compare its strength at different spatial scales, to dissect the relative contribution of pre- versus postmating isolation, and to infer the involvement of ecological divergence on hybridization. Because F1 hybrids are viable, fertile and not uncommon, understanding the dynamics of hybridization in this trio of major malaria vectors has important implications for how adaptations arise and spread across the group, and in planning studies of the safety and efficacy of gene drive as a means of malaria control. We first performed a systematic review and meta-analysis of published surveys reporting on hybrid prevalence, showing strong reproductive isolation at a continental scale despite geographically restricted exceptions. Second, we exploited our own extensive field data sets collected at a regional scale in two contrasting environmental settings, to assess: (i) levels of premating isolation; (ii) spatio/temporal and frequency-dependent dynamics of hybridization, (iii) relationship between reproductive isolation and ecological divergence and (iv) hybrid viability penalty. Results are in accordance with ecological speciation theory predicting a positive association between the strength of reproductive isolation and degree of ecological divergence, and indicate that postmating isolation does contribute to reproductive isolation among these species. Specifically, only postmating isolation was positively associated with ecological divergence, whereas premating isolation was correlated with phylogenetic distance.

On the distribution and genetic differentiation of Anopheles gambiae s.s. molecular forms.

This paper summarises published and unpublished data on the spatial and temporal distribution, and on the genetic characterisation of molecular forms M and S of Anopheles gambiae s.s. The two forms are characterised by a high level of gene-flow restriction, by a largely overlapping geographical and temporal distribution, and by a low degree of genetic differentiation. Floating paracentric inversions on chromosome-2 are shown to be shared by the two forms, although with very different frequencies of alternative arrangements, confirming that these inversions are most probably involved in ecotypic adaptation, rather than in the building of reproductive barriers. Further studies and tools are needed to throw light on the genetic and biological differentiation of M and S to improve the knowledge of the real composition of the vector system, of its demography, population genetics and dynamics, also in view of the possible consequences on the transmission of human pathogens in sub-Saharan Africa. Preliminary results and perspectives of the use of transposable element insertion sites as markers of genetic differentiation and tools for population genetic studies are discussed.