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

Phenotypic evidence of deltamethrin resistance and identification of selective sweeps in Anopheles mosquitoes on the Bijagós Archipelago, Guinea-Bissau

Vector control in the Bijagós Archipelago of Guinea-Bissau currently relies on pyrethroid insecticide-treated nets. However, data on insecticide resistance in Guinea-Bissau is limited. This study identified deltamethrin resistance in the Anopheles gambiae sensu lato complex on Bubaque island using WHO tube tests in November 2022. Whole genome sequencing of An. gambiae sensu stricto mosquitoes identified six single nucleotide polymorphisms (SNPs) previously associated with, or putatively associated with, insecticide resistance: T791M, L995F, N1570Y, A1746S and P1874L in the vgsc gene, and L119V in the gste2 gene. Twenty additional non-synonymous SNPs were identified in insecticide-resistance associated genes. Four of these SNPs were present at frequencies over 5% in the population: T154S, I126F and G26S in the vgsc gene and A65S in ace1. Genome wide selection scans using Garud's H12 statistic identified two selective sweeps: one in chromosome X and one in chromosome 2R. Both selective sweeps overlap with metabolic genes previously associated with insecticide resistance, including cyp9k1 and the cyp6aa/cyp6p gene cluster. This study presents the first phenotypic testing for deltamethrin resistance and the first whole genome sequence data for Anophelesgambiae mosquitoes from the Bijagós, contributing data of significance for vector control policy in this region.

Population genomic evidence of a putative 'far-west' African cryptic taxon in the Anopheles gambiae complex

The two main Afrotropical malaria vectors - Anopheles coluzzii and An. gambiae - are genetically distinct and reproductively isolated across West Africa. However, populations at the western extreme of their range are assigned as "intermediate" between the two species by whole genome sequence (WGS) data, and as hybrid forms by conventional molecular diagnostics. By exploiting WGS data from 1190 specimens collected across west Africa via the Anopheles gambiae 1000 Genomes network, we identified a putative taxon in the far-west (provisionally named Bissau molecular form), which did not arise by admixture but rather may have originated at the same time as the split between An. coluzzii and An. gambiae. Intriguingly, this taxon lacks insecticide resistance mechanisms commonly observed in the two main species. These findings lead to a change of perspective on malaria vector species in the far-west region with potential for epidemiological implications, and a new challenge for genetic-based mosquito control approaches.

Mitochondrial Variation in Anopheles gambiae and Anopheles coluzzii: Phylogeographic Legacy and Mitonuclear Associations With Metabolic Resistance to Pathogens and Insecticides

Mitochondrial DNA has been a popular marker in phylogeography, phylogeny, and molecular ecology, but its complex evolution is increasingly recognized. Here, we investigated mitochondrial DNA variation in Anopheles gambiae and Anopheles coluzzii, in relation to other species in the Anopheles gambiae complex, by assembling the mitogenomes of 1,219 mosquitoes across Africa. The mitochondrial DNA phylogeny of the Anopheles gambiae complex was consistent with previously reported highly reticulated evolutionary history, revealing important discordances with the species tree. The three most widespread species (An. gambiae, An. coluzzii, and Anopheles arabiensis), known for extensive historical introgression, could not be discriminated based on mitogenomes. Furthermore, a monophyletic clustering of the three saltwater-tolerant species (Anopheles merus, Anopheles melas, and Anopheles bwambae) in the Anopheles gambiae complex also suggested that introgression and possibly selection shaped mitochondrial DNA evolution. Mitochondrial DNA variation in An. gambiae and An. coluzzii across Africa revealed significant partitioning among populations and species. A peculiar mitochondrial DNA lineage found predominantly in An. coluzzii and in the hybrid taxon of the African "far-west" exhibited divergence comparable to the interspecies divergence in the Anopheles gambiae complex, with a geographic distribution matching closely An. coluzzii's geographic range. This phylogeographic relict of the An. coluzzii and An. gambiae split was associated with population and species structure, but not with the rare Wolbachia occurrence. The lineage was significantly associated with single nucleotide polymorphisms in the nuclear genome, particularly in genes associated with pathogen and insecticide resistance. These findings underline potential mitonuclear coevolution history and the role played by mitochondria in shaping metabolic responses to pathogens and insecticides in Anopheles.

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

Speciation within the Anopheles gambiae complex: high-throughput whole genome sequencing reveals evidence of a putative new cryptic taxon in 'far-west' Africa

The two main Afrotropical malaria vectors - Anopheles coluzzii and An. gambiae - are genetically distinct and reproductively isolated across West Africa. However, populations at the western extreme of their range are assigned as "intermediate" between the two species by whole genome sequence (WGS) data, and as hybrid forms by conventional molecular diagnostics. By exploiting WGS data from 1,190 specimens collected across west Africa via the Anopheles gambiae 1000 Genomes network, we identify a novel putative taxon in the far-west (provisionally named Bissau molecular form), which did not arise by admixture but rather originated at the same time as the split between An. coluzzii and An. gambiae. Intriguingly, these populations lack insecticide resistance mechanisms commonly observed in the two main species. These findings lead to a change of perspective on malaria vector species in the far-west region with potential for epidemiological implications, and a new challenge for genetic-based mosquito control approaches.

Genomic surveillance of Anopheles mosquitoes on the Bijagós Archipelago using custom targeted amplicon sequencing identifies mutations associated with insecticide resistance

BACKGROUND

Insecticide resistance is reducing the efficacy of vector control interventions, consequently threatening efforts to control vector-borne diseases, including malaria. Investigating the prevalence of molecular markers of resistance is a useful tool for monitoring the spread of insecticide resistance in disease vectors. The Bijagós Archipelago (Bijagós) in Guinea-Bissau is a region of stable malaria transmission where insecticide-treated nets are the mainstay for malaria control. However, the prevalence of molecular markers of insecticide resistance in malaria vectors is not well understood.

METHODS

A total of 214 Anopheles mosquitoes were analysed from 13 islands across the Bijagós. These mosquitoes were collected using CDC light traps in November 2019, during the peak malaria transmission season. High-throughput multiplex amplicon sequencing was used to investigate the prevalence of 17 different molecular markers associated with insecticide resistance in four genes: vgsc, rdl, ace1 and gste2.

RESULTS

Of the 17 screened mutations, four were identified in mosquitoes from the Bijagós: vgsc L995F (12.2%), N1570Y (6.2%) and A1746S (0.7%) and rdl A269G (1.1%). This study is the first to report the L995F knock-down resistance (kdr)-west allele in Anopheles melas on the Archipelago. An additional eight non-synonymous single-nucleotide polymorphisms were identified across the four genes which have not been described previously. The prevalences of the vgsc L995F and N1570Y mutations were higher on Bubaque Island than on the other islands in this study; Bubaque is the most populous island in the archipelago, with the greatest population mobility and connection to continental Guinea-Bissau.

CONCLUSIONS

This study provides the first surveillance data for genetic markers present in malaria vectors from islands across the Bijagós Archipelago. Overall prevalence of insecticide resistance mutations was found to be low. However, the identification of the vgsc L995F and N1570Y mutations associated with pyrethroid resistance warrants further monitoring. This is particularly important as the mainstay of malaria control on the islands is the use of pyrethroid insecticide-treated nets.

Conceptual risk assessment of mosquito population modification gene-drive systems to control malaria transmission: preliminary hazards list workshops

The field-testing and eventual adoption of genetically-engineered mosquitoes (GEMs) to control vector-borne pathogen transmission will require them meeting safety criteria specified by regulatory authorities in regions where the technology is being considered for use and other locales that might be impacted. Preliminary risk considerations by researchers and developers may be useful for planning the baseline data collection and field research used to address the anticipated safety concerns. Part of this process is to identify potential hazards (defined as the inherent ability of an entity to cause harm) and their harms, and then chart the pathways to harm and evaluate their probability as part of a risk assessment. The University of California Malaria Initiative (UCMI) participated in a series of workshops held to identify potential hazards specific to mosquito population modification strains carrying gene-drive systems coupled to anti-parasite effector genes and their use in a hypothetical island field trial. The hazards identified were placed within the broader context of previous efforts discussed in the scientific literature. Five risk areas were considered i) pathogens, infections and diseases, and the impacts of GEMs on human and animal health, ii) invasiveness and persistence of GEMs, and interactions of GEMs with target organisms, iii) interactions of GEMs with non-target organisms including horizontal gene transfer, iv) impacts of techniques used for the management of GEMs and v) evolutionary and stability considerations. A preliminary hazards list (PHL) was developed and is made available here. This PHL is useful for internal project risk evaluation and is available to regulators at prospective field sites. UCMI project scientists affirm that the subsequent processes associated with the comprehensive risk assessment for the application of this technology should be driven by the stakeholders at the proposed field site and areas that could be affected by this intervention strategy.

RNA-Seq-Pop: Exploiting the sequence in RNA sequencing-A Snakemake workflow reveals patterns of insecticide resistance in the malaria vector Anopheles gambiae

We provide a reproducible and scalable Snakemake workflow, called RNA-Seq-Pop, which provides end-to-end analysis of RNA sequencing data sets. The workflow allows the user to perform quality control, perform differential expression analyses and call genomic variants. Additional options include the calculation of allele frequencies of variants of interest, summaries of genetic variation and population structure, and genome-wide selection scans, together with clear visualizations. RNA-Seq-Pop is applicable to any organism, and we demonstrate the utility of the workflow by investigating pyrethroid resistance in selected strains of the major malaria mosquito, Anopheles gambiae. The workflow provides additional modules specifically for An. gambiae, including estimating recent ancestry and determining the karyotype of common chromosomal inversions. The Busia laboratory colony used for selections was collected in Busia, Uganda, in November 2018. We performed a comparative analysis of three groups: a parental G24 Busia strain; its deltamethrin-selected G28 offspring; and the susceptible reference strain Kisumu. Measures of genetic diversity reveal patterns consistent with that of laboratory colonization and selection, with the parental Busia strain exhibiting the highest nucleotide diversity, followed by the selected Busia offspring, and finally, Kisumu. Differential expression and variant analyses reveal that the selected Busia colony exhibits a number of distinct mechanisms of pyrethroid resistance, including the Vgsc-995S target-site mutation, upregulation of SAP genes, P450s and a cluster of carboxylesterases. During deltamethrin selections, the 2La chromosomal inversion rose in frequency (from 33% to 86%), supporting a previous link with pyrethroid resistance. RNA-Seq-Pop is hosted at: github.com/sanjaynagi/rna-seq-pop. We anticipate that the workflow will provide a useful tool to facilitate reproducible, transcriptomic studies in An. gambiae and other taxa.

Climate and Environmental Changes and Their Potential Effects on the Dynamics of Chagas Disease: Hybridization in Rhodniini (Hemiptera, Triatominae)

Chagas disease affects about eight million people. In view of the issues related to the influence of anthropogenic changes in the dynamics of the distribution and reproductive interaction of triatomines, we performed experimental crosses between species of the Rhodniini tribe in order to evaluate interspecific reproductive interactions and hybrid production capacity. Reciprocal crossing experiments were conducted among Rhodnius brethesi × R. pictipes, R. colombiensis × R. ecuadoriensis, R. neivai × R. prolixus, R. robustus × R. prolixus, R. montenegrensis × R. marabaensis; R. montenegrensis × R. robustus, R. prolixus × R. nasutus and R. neglectus × R. milesi. With the exception of crosses between R. pictipes ♀ × R. brethesi ♂, R. ecuadoriensis ♀ × R. colombiensis ♂ and R. prolixus ♀ × R. neivai ♂, all experimental crosses resulted in hybrids. Our results demonstrate that both allopatric and sympatric species produce hybrids, which can generate concern for public health agencies in the face of current anthropogenic events. Thus, we demonstrate that species of the Rhodniini tribe are capable of producing hybrids under laboratory conditions. These results are of great epidemiological importance and raise an important discussion about the influence of climatic and environmental interactions on Chagas disease dynamics.

Introgression of a synthetic sex ratio distortion transgene into different genetic backgrounds of Anopheles coluzzii

The development of genetically modified mosquitoes (GMM) and their subsequent field release offers innovative approaches for vector control of malaria. A non-gene drive self-limiting male-bias Ag(PMB)1 strain has been developed in a 47-year-old laboratory G3 strain of Anopheles gambiae s.l. When Ag(PMB)1 males are crossed to wild-type females, expression of the endonuclease I-PpoI during spermatogenesis causes the meiotic cleavage of the X chromosome in sperm cells, leading to fertile offspring with a 95% male bias. However, World Health Organization states that the functionality of the transgene could differ when inserted in different genetic backgrounds of Anopheles coluzzii which is currently a predominant species in several West-African countries and thus a likely recipient for a potential release of self-limiting GMMs. In this study, we introgressed the transgene from the donor Ag(PMB)1 by six serial backcrosses into two recipient colonies of An. coluzzii that had been isolated in Mali and Burkina Faso. Scans of informative Single Nucleotide Polymorphism (SNP) markers and whole-genome sequencing analysis revealed a nearly complete introgression of chromosomes 3 and X, but a remarkable genomic divergence in a large region of chromosome 2 between the later backcrossed (BC6) transgenic offspring and the recipient paternal strains. These findings suggested to extend the backcrossing breeding strategy beyond BC6 generation and increasing the introgression efficiency of critical regions that have ecological and epidemiological implications through the targeted selection of specific markers. Disregarding differential introgression efficiency, we concluded that the phenotype of the sex ratio distorter is stable in the BC6 introgressed An. coluzzii strains.

High Levels of Admixture in Anopheles gambiae Populations from Côte d'Ivoire Revealed by Multilocus Genotyping

Anopheles coluzzii and An. gambiae-the two most recently radiated species of the An. gambiae complex and the major Afrotropical malaria vector species-are identified by markers in the X-centromeric IGS rDNA region. Putative IGS-hybrids are rarely found in the field, except in restricted areas where genomic studies have led to the hypothesis that the observed IGS-patterns are due to cryptic taxa rather than to hybridization between the two species. We investigated the genome-wide levels of admixture in two villages in Côte d'Ivoire where high levels of IGS-hybrids have been detected, confirming unparalleled high frequencies in the coastal village. Genotyping of 24 Ancestry Informative Markers (AIMs) along the three chromosomes produced discordant results between the IGS-marker and the multilocus genotype obtained for AIMs across the whole genome (29%) as well as AIMs on chromosome-X (considered to be fundamental for species reproductive isolation) only (21%). Results highlight a complicated pattern of admixture that deserves deeper genomic analyses to understand better possible underlying causes (from extensive processes of hybridization to the existence of different cryptic taxa), and stress the need of developing advanced diagnostics for An. coluzzii, An. gambiae and putative new taxa, instrumental for assessing taxon-specific epidemiological characters.

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.

Modelling spatiotemporal trends in the frequency of genetic mutations conferring insecticide target-site resistance in African mosquito malaria vector species

BACKGROUND

Resistance in malaria vectors to pyrethroids, the most widely used class of insecticides for malaria vector control, threatens the continued efficacy of vector control tools. Target-site resistance is an important genetic resistance mechanism caused by mutations in the voltage-gated sodium channel (Vgsc) gene that encodes the pyrethroid target-site. Understanding the geographic distribution of target-site resistance, and temporal trends across different vector species, can inform strategic deployment of vector control tools.

RESULTS

We develop a Bayesian statistical spatiotemporal model to interpret species-specific trends in the frequency of the most common resistance mutations, Vgsc-995S and Vgsc-995F, in three major malaria vector species Anopheles gambiae, An. coluzzii, and An. arabiensis over the period 2005-2017. The models are informed by 2418 observations of the frequency of each mutation in field sampled mosquitoes collected from 27 countries spanning western and eastern regions of Africa. For nine selected countries, we develop annual predictive maps which reveal geographically structured patterns of spread of each mutation at regional and continental scales. The results show associations, as well as stark differences, in spread dynamics of the two mutations across the three vector species. The coverage of ITNs was an influential predictor of Vgsc allele frequencies, with modelled relationships between ITN coverage and allele frequencies varying across species and geographic regions. We found that our mapped Vgsc allele frequencies are a significant partial predictor of phenotypic resistance to the pyrethroid deltamethrin in An. gambiae complex populations.

CONCLUSIONS

Our predictive maps show how spatiotemporal trends in insecticide target-site resistance mechanisms in African An. gambiae vary across individual vector species and geographic regions. Molecular surveillance of resistance mechanisms will help to predict resistance phenotypes and track their spread.

Entomological Survey Confirms Changes in Mosquito Composition and Abundance in Senegal and Reveals Discrepancies among Results by Different Host-Seeking Female Traps

Mosquitoes-borne diseases are major public health issues particularly in Africa. Vector control interventions and human-made environmental/climatic changes significantly affect the distribution and abundance of vector species. We carried out an entomological survey targeting host-seeking mosquitos in two different ecological contexts-coastal and inland-in Senegal, by CDC-light and BG-sentinel traps. Results show high predominance of Culex quinquefasciatus (90%) and of Anopheles arabiensis within malaria vectors (46%), with mean numbers of females/trap/nights =8 and <1, respectively, reinforcing previous evidence of changes in species composition and abundance, highlighting thus increasing risk of transmission of filariasis and emerging arboviruses in the Senegambia region. From the methodological perspective, results show a higher specificity of BG traps for Cx. quinquefasciatus and of CDC traps for An. gambiae s.l. and highlight that, despite both traps target the host-seeking fraction of the population, they provide different patterns of species abundance, temporal dynamics and host-seeking activity, leading to possible misinterpretation of the species bionomics. This draws attention to the need of taking into account trapping performance, in order to provide realistic quantification of the number of mosquitoes per units of space and time, the crucial parameter for evaluating vector-human contact, and estimating risk of pathogen transmission.

Genotyping of Ascaris spp. infecting humans and pigs in Italy, Slovakia and Colombia

BACKGROUND

The systematics and taxonomy of Ascaris lumbricoides and Ascaris suum, two of the world's most widespread nematodes, still represent a highly debated scientific issue. Two different transmission scenarios have been described according to endemicity: separated host-specific transmission cycles in endemic regions, and a single pool of infection shared by humans and pigs in non-endemic regions. The swine roundworm A. suum is now recognized as an important cause of human ascariasis also in endemic areas such as China, where cross-infections and hybridization have also been reported, as well as in non-endemic regions like Italy. This study aimed to investigate the molecular epidemiology of human and pig ascariasis in three countries representing different epidemiological scenarios: Italy as a non-endemic country, Colombia as an endemic country, and Slovakia as a non-endemic country, but with a poor socio-economic context linked to some focal populations of Roma settlements.

MATERIALS AND METHODS

A total of 237 nematodes were analysed: 46 from Colombia (13 from humans, 33 from pigs), 114 from Slovakia (20 from humans, 94 from pigs) and 77 from Italy (17 from humans and 60 from pigs). Genotyping by PCR-RFLP of nuclear (ITS) and sequencing of mitochondrial (cox1) target regions were performed. ITS genotypes were used to estimate the Hardy-Weinberg (HW) equilibrium according to hosts and country of origin. The partial cox1 sequences were used to analyse genetic polymorphisms according to hosts and country of origin, as well as to infer the network of haplotypes, their evolutionary relationships and geographical distribution.

RESULTS

110 quality cox1 sequences were obtained. Haplotype network revealed three main groups corresponding to clade A, B and C. Clade C included most of the human cases from Italy, while those from Slovakia and Colombia were grouped in clade B. Ascaris from Italian and Colombian pigs showed HW equilibrium at the ITS marker, while disequilibrium was found in A. lumbricoides from Slovak pigs, which suggest a high unexpected amount of roundworms of human origin circulating also in pigs.

CONCLUSIONS

This study updates and extends the current understanding of Ascaris species and genotypes circulating in different epidemiological scenarios, with particular attention to the inclusion of human-derived Ascaris in the phylogenetic cluster C. Despite the evidence of HW equilibrium in the ITS in pig-derived Italian samples, the amount of genetic variation seems to support the existence of two closely related species.

Unravelling population structure heterogeneity within the genome of the malaria vector Anopheles gambiae

Background

Whole genome re-sequencing provides powerful data for population genomic studies, allowing robust inferences of population structure, gene flow and evolutionary history. For the major malaria vector in Africa, Anopheles gambiae, other genetic aspects such as selection and adaptation are also important. In the present study, we explore population genetic variation from genome-wide sequencing of 765 An. gambiae and An. coluzzii specimens collected from across Africa. We used t-SNE, a recently popularized dimensionality reduction method, to create a 2D-map of An. gambiae and An. coluzzii genes that reflect their population structure similarities.

Results

The map allows intuitive navigation among genes distributed throughout the so-called “mainland” and numerous surrounding “island-like” gene clusters. These gene clusters of various sizes correspond predominantly to low recombination genomic regions such as inversions and centromeres, and also to recent selective sweeps. Because this mosquito species complex has been studied extensively, we were able to support our interpretations with previously published findings. Several novel observations and hypotheses are also made, including selective sweeps and a multi-locus selection event in Guinea-Bissau, a known intense hybridization zone between An. gambiae and An. coluzzii.

Conclusions

Our results present a rich dataset that could be utilized in functional investigations aiming to shed light onto An. gambiae s.l genome evolution and eventual speciation. In addition, the methodology presented here can be used to further characterize other species not so well studied as An. gambiae, shortening the time required to progress from field sampling to the identification of genes and genomic regions under unique evolutionary processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07722-y.

Evidence for natural hybridization and novel Wolbachia strain superinfections in the Anopheles gambiae complex from Guinea

Wolbachia, a widespread bacterium which can influence mosquito-borne pathogen transmission, has recently been detected within Anopheles (An.) species that are malaria vectors in Sub-Saharan Africa. Although studies have reported Wolbachia strains in the An. gambiae complex, apparent low density and prevalence rates require confirmation. In this study, wild Anopheles mosquitoes collected from two regions of Guinea were investigated. In contrast with previous studies, RNA was extracted from adult females (n = 516) to increase the chances for the detection of actively expressed Wolbachia genes, determine Wolbachia prevalence rates and estimate relative strain densities. Molecular confirmation of mosquito species and Wolbachia multilocus sequence typing (MLST) were carried out to analyse phylogenetic relationships of mosquito hosts and newly discovered Wolbachia strains. Strains were detected in An. melas (prevalence rate of 11.6%-16/138) and hybrids between An. melas and An. gambiae sensu stricto (prevalence rate of 40.0%-6/15) from Senguelen in the Maferinyah region. Furthermore, a novel high-density strain, termed wAnsX, was found in an unclassified Anopheles species. The discovery of novel Wolbachia strains (particularly in members, and hybrids, of the An. gambiae complex) provides further candidate strains that could be used for future Wolbachia-based malaria biocontrol strategies.

Novel genotyping approaches to easily detect genomic admixture between the major Afrotropical malaria vector species, Anopheles coluzzii and An. gambiae

The two most efficient and most recently radiated Afrotropical vectors of human malaria - Anopheles coluzzii and An. gambiae - are identified by single-locus diagnostic PCR assays based on species-specific markers in a 4 Mb region on chromosome-X centromere. Inherently, these diagnostic assays cannot detect interspecific autosomal admixture shown to be extensive at the westernmost and easternmost extremes of the species range. The main aim of this study was to develop novel, easy-to-implement tools for genotyping An. coluzzii and An. gambiae-specific ancestral informative markers (AIMs) identified from the Anopheles gambiae 1000 genomes (Ag1000G) project. First, we took advantage of this large set of data in order to develop a multilocus approach to genotype 26 AIMs on all chromosome arms valid across the species range. Second, we tested the multilocus assay on samples from Guinea Bissau, The Gambia and Senegal, three countries spanning the westernmost hybridization zone, where conventional species diagnostic is problematic due to the putative presence of a novel "hybrid form". The multilocus assay was able to capture patterns of admixture reflecting those revealed by the whole set of AIMs and provided new original data on interspecific admixture in the region. Third, we developed an easy-to-use, cost-effective PCR approach for genotyping two AIMs on chromosome-3 among those included in the multilocus approach, opening the possibility for advanced identification of species and of admixed specimens during routine large scale entomological surveys, particularly, but not exclusively, at the extremes of the range, where WGS data highlighted unexpected autosomal admixture.

Evolution of the Insecticide Target Rdl in African Anopheles Is Driven by Interspecific and Interkaryotypic Introgression

The evolution of insecticide resistance mechanisms in natural populations of Anopheles malaria vectors is a major public health concern across Africa. Using genome sequence data, we study the evolution of resistance mutations in the resistance to dieldrin locus (Rdl), a GABA receptor targeted by several insecticides, but most notably by the long-discontinued cyclodiene, dieldrin. The two Rdl resistance mutations (296G and 296S) spread across West and Central African Anopheles via two independent hard selective sweeps that included likely compensatory nearby mutations, and were followed by a rare combination of introgression across species (from A. gambiae and A. arabiensis to A. coluzzii) and across nonconcordant karyotypes of the 2La chromosomal inversion. Rdl resistance evolved in the 1950s as the first known adaptation to a large-scale insecticide-based intervention, but the evolutionary lessons from this system highlight contemporary and future dangers for management strategies designed to combat development of resistance in malaria vectors.

Mapping trends in insecticide resistance phenotypes in African malaria vectors

Mitigating the threat of insecticide resistance in African malaria vector populations requires comprehensive information about where resistance occurs, to what degree, and how this has changed over time. Estimating these trends is complicated by the sparse, heterogeneous distribution of observations of resistance phenotypes in field populations. We use 6,423 observations of the prevalence of resistance to the most important vector control insecticides to inform a Bayesian geostatistical ensemble modelling approach, generating fine-scale predictive maps of resistance phenotypes in mosquitoes from the Anopheles gambiae complex across Africa. Our models are informed by a suite of 111 predictor variables describing potential drivers of selection for resistance. Our maps show alarming increases in the prevalence of resistance to pyrethroids and DDT across sub-Saharan Africa from 2005 to 2017, with mean mortality following insecticide exposure declining from almost 100% to less than 30% in some areas, as well as substantial spatial variation in resistance trends.

Monitoring and molecular profiling of contemporary insecticide resistance status of malaria vectors in Guinea-Bissau

Despite reduction in the prevalence of malaria, Guinea-Bissau (GB) is still widely affected by the disease that is primarily vectored by Anopheles gambiae s.l. mosquitoes. Monitoring mosquito susceptibility and investigating the insecticide resistance status is an integral part of malaria control actions. Here, mosquito populations from five regions of GB: Bafatá, Bissau, Buba, Cacheu and Gabu were monitored for species ID and insecticide resistance, using diagnostic and intensity WHO bioassays, as well as molecular assays. Phenotypic and molecular identification of species showed the presence of An. gambiae s.s. (S form), An. coluzzii (M form) and An. arabiensis, as well as rare An. arabiensis/ An. gambiae hybrids. Resistance to permethrin and deltamethrin was found in all Anopheles populations assayed, with the intensity of resistance for permethrin being moderate to high, as confirmed by bioassays performed at concentration intensities of 5X and 10X. Consistent to these findings, molecular analysis showed a higher frequency of knock-down resistance (kdr) mutations (L1014F, L1014S, reaching > 90% in some areas) compared to previous studies in the same region, as well as detected for the first time the presence of the super kdr mutation (N1575Y) in GB. The "iAche" (G119S) resistance mutation was also found in GB in low frequencies (up to 12.41%). Additionally, the synergistic PBO-permethrin bioassays suggested partial involvement of non target (metabolic and/or reduced penetration) resistance mechanism. Expression analysis of known pyrethroid metabolisers indicated the slight overexpression and possible association of the cytochrome P450s CYP6Z1, CYP4G16 with the pyrethroid resistance phenotype. The findings should guide future evidence-based resistance management strategies in GB.

Genome-wide divergence among invasive populations of Aedes aegypti in California

BACKGROUND

In the summer of 2013, Aedes aegypti Linnaeus was first detected in three cities in central California (Clovis, Madera and Menlo Park). It has now been detected in multiple locations in central and southern CA as far south as San Diego and Imperial Counties. A number of published reports suggest that CA populations have been established from multiple independent introductions.

RESULTS

Here we report the first population genomics analyses of Ae. aegypti based on individual, field collected whole genome sequences. We analyzed 46 Ae. aegypti genomes to establish genetic relationships among populations from sites in California, Florida and South Africa. Based on 4.65 million high quality biallelic SNPs, we identified 3 major genetic clusters within California; one that includes all sample sites in the southern part of the state (South of Tehachapi mountain range) plus the town of Exeter in central California and two additional clusters in central California.

CONCLUSIONS

A lack of concordance between mitochondrial and nuclear genealogies suggests that the three founding populations were polymorphic for two main mitochondrial haplotypes prior to being introduced to California. One of these has been lost in the Clovis populations, possibly by a founder effect. Genome-wide comparisons indicate extensive differentiation between genetic clusters. Our observations support recent introductions of Ae. aegypti into California from multiple, genetically diverged source populations. Our data reveal signs of hybridization among diverged populations within CA. Genetic markers identified in this study will be of great value in pursuing classical population genetic studies which require larger sample sizes.

Multi-population genomic analysis of malaria parasites indicates local selection and differentiation at the gdv1 locus regulating sexual development

Parasites infect hosts in widely varying environments, encountering diverse challenges for adaptation. To identify malaria parasite genes under locally divergent selection across a large endemic region with a wide spectrum of transmission intensity, genome sequences were obtained from 284 clinical Plasmodium falciparum infections from four newly sampled locations in Senegal, The Gambia, Mali and Guinea. Combining these with previous data from seven other sites in West Africa enabled a multi-population analysis to identify discrete loci under varying local selection. A genome-wide scan showed the most exceptional geographical divergence to be at the early gametocyte gene locus gdv1 which is essential for parasite sexual development and transmission. We identified a major structural dimorphism with alternative 1.5 kb and 1.0 kb sequence deletions at different positions of the 3'-intergenic region, in tight linkage disequilibrium with the most highly differentiated single nucleotide polymorphism, one of the alleles being very frequent in Senegal and The Gambia but rare in the other locations. Long non-coding RNA transcripts were previously shown to include the entire antisense of the gdv1 coding sequence and the portion of the intergenic region with allelic deletions, suggesting adaptive regulation of parasite sexual development and transmission in response to local conditions.

Genetic diversity of the African malaria vector Anopheles gambiae

The sustainability of malaria control in Africa is threatened by the rise of insecticide resistance in Anopheles mosquitoes, which transmit the disease. To gain a deeper understanding of how mosquito populations are evolving, here we sequenced the genomes of 765 specimens of Anopheles gambiae and Anopheles coluzzii sampled from 15 locations across Africa, and identified over 50 million single nucleotide polymorphisms within the accessible genome. These data revealed complex population structure and patterns of gene flow, with evidence of ancient expansions, recent bottlenecks, and local variation in effective population size. Strong signals of recent selection were observed in insecticide-resistance genes, with several sweeps spreading over large geographical distances and between species. The design of new tools for mosquito control using gene-drive systems will need to take account of high levels of genetic diversity in natural mosquito populations.

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.