Assessing connectivity despite high diversity in island populations of a malaria mosquito

Fine-scale spatial distribution of deltamethrin resistance and population structure of Anopheles funestus and Anopheles arabiensis populations in Southern Mozambique

BACKGROUND

Insecticide resistance in malaria vectors can be spatially highly heterogeneous, yet population structure analyses frequently find relatively high levels of gene flow among mosquito populations. Few studies have contemporaneously assessed phenotypic, genotypic and population structure analysis on mosquito populations and none at fine geographical scales. In this study, genetic diversity, population structure, and insecticide resistance profiles of Anopheles funestus and Anopheles arabiensis were examined across mosquito populations from and within neighbouring villages.

METHODS

Mosquitoes were collected from 11 towns in southern Mozambique, as well as from different neighbourhoods within the town of Palmeira, during the peak malaria transmission season in 2016. CDC bottle bioassay and PCR assays were performed with Anopheles mosquitoes at each site to determine phenotypic and molecular insecticide resistance profiles, respectively. Microsatellite analysis was conducted on a subsample of mosquitoes to estimate genetic diversity and population structure.

RESULTS

Phenotypic insecticide resistance to deltamethrin was observed in An. funestus sensu stricto (s.s.) throughout the area, though a high level of mortality variation was seen. However, 98% of An. funestus s.s. were CYP6P9a homozygous resistant. An. arabiensis was phenotypically susceptible to deltamethrin and 99% were kdr homozygous susceptible. Both Anopheles species exhibited high allelic richness and heterozygosity. Significant deviations from Hardy-Weinberg equilibrium were observed, and high linkage disequilibrium was seen for An. funestus s.s., supporting population subdivision. However, the FST values were low for both anophelines (- 0.00457 to 0.04213), Nm values were high (9.4-71.8 migrants per generation), AMOVA results showed almost 100% genetic variation among and within individuals, and Structure analysis showed no clustering of An. funestus s.s. and An. arabiensis populations. These results suggest high gene flow among mosquito populations.

CONCLUSION

Despite a relatively high level of phenotypic variation in the An. funestus population, molecular analysis shows the population is admixed. These data indicate that CYP6P9a resistance markers do not capture all phenotypic variation in the area, but also that resistance genes of high impact are likely to easily spread in the area. Conversely, other strategies, such as transgenic mosquito release programmes will likely not face challenges in this locality.

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.

Selection of sites for field trials of genetically engineered mosquitoes with gene drive

Novel malaria control strategies using genetically engineered mosquitoes (GEMs) are on the horizon. Population modification is one approach wherein mosquitoes are engineered with genes rendering them refractory to the malaria parasite, Plasmodium falciparum, coupled with a low-threshold, Cas9-based gene drive. When released into a wild vector population, GEMs preferentially transmit these parasite-blocking genes to their offspring, ultimately modifying a vector population into a nonvector one. Deploying this technology awaits ecologically contained field trial evaluations. Here, we consider a process for site selection, the first critical step in designing a trial. Our goal is to identify a site that maximizes prospects for success, minimizes risk, and serves as a fair, valid, and convincing test of efficacy and impacts of a GEM product intended for large-scale deployment in Africa. We base site selection on geographic, geological, and biological, rather than social or legal, criteria. We recognize the latter as critically important but not as a first step in selecting a site. We propose physical islands as being the best candidates for a GEM field trial and present an evaluation of 22 African islands. We consider geographic and genetic isolation, biological complexity, island size, and topography and identify two island groups that satisfy key criteria for ideal GEM field trial sites.

The origin of island populations of the African malaria mosquito, Anopheles coluzzii

Anopheles coluzzii is a major malaria vector throughout its distribution in west-central Africa. Here we present a whole-genome study of 142 specimens from nine countries in continental Africa and three islands in the Gulf of Guinea. This sample set covers a large part of this species' geographic range. Our population genomic analyses included a description of the structure of mainland populations, island populations, and connectivity between them. Three genetic clusters are identified among mainland populations and genetic distances (F_ST) fits an isolation-by-distance model. Genomic analyses are applied to estimate the demographic history and ancestry for each island. Taken together with the unique biogeography and history of human occupation for each island, they present a coherent explanation underlying levels of genetic isolation between mainland and island populations. We discuss the relationship of our findings to the suitability of São Tomé and Príncipe islands as candidate sites for potential field trials of genetic-based malaria control strategies.

Genome variation and population structure among 1142 mosquitoes of the African malaria vector species Anopheles gambiae and Anopheles coluzzii

Mosquito control remains a central pillar of efforts to reduce malaria burden in sub-Saharan Africa. However, insecticide resistance is entrenched in malaria vector populations, and countries with a high malaria burden face a daunting challenge to sustain malaria control with a limited set of surveillance and intervention tools. Here we report on the second phase of a project to build an open resource of high-quality data on genome variation among natural populations of the major African malaria vector species Anopheles gambiae and Anopheles coluzzii We analyzed whole genomes of 1142 individual mosquitoes sampled from the wild in 13 African countries, as well as a further 234 individuals comprising parents and progeny of 11 laboratory crosses. The data resource includes high-confidence single-nucleotide polymorphism (SNP) calls at 57 million variable sites, genome-wide copy number variation (CNV) calls, and haplotypes phased at biallelic SNPs. We use these data to analyze genetic population structure and characterize genetic diversity within and between populations. We illustrate the utility of these data by investigating species differences in isolation by distance, genetic variation within proposed gene drive target sequences, and patterns of resistance to pyrethroid insecticides. This data resource provides a foundation for developing new operational systems for molecular surveillance and for accelerating research and development of new vector control tools. It also provides a unique resource for the study of population genomics and evolutionary biology in eukaryotic species with high levels of genetic diversity under strong anthropogenic evolutionary pressures.

Recent advances and future perspectives in vector-omics

We have reviewed recent progress and the remaining challenges in vector-omics. We have highlighted several technologies and applications that facilitate novel biological insights beyond achieving a reference-quality genome assembly. Among other topics, we have discussed the applications of chromatin conformation capture, chromatin accessibility assays, optical mapping, full-length RNA sequencing, single cell RNA analysis, proteomics, and population genomics. We anticipate that we will witness a great expansion in vector-omics research not only in its application in a broad range of species, but also its ability to uncover novel genetic elements and tackle previously inaccessible regions of the genome. It is our hope that the continued innovation in device portability, cost reduction, and informatics support will in the foreseeable future bring vector-omics to every vector laboratory and field station in the world, which will have an unparalleled impact on basic research and the control of vector-borne infectious diseases.

Transcontinental dispersal of Anopheles gambiae occurred from West African origin via serial founder events

The mosquito Anopheles gambiae s.s. is distributed across most of sub-Saharan Africa and is of major scientific and public health interest for being an African malaria vector. Here we present population genomic analyses of 111 specimens sampled from west to east Africa, including the first whole genome sequences from oceanic islands, the Comoros. Genetic distances between populations of A. gambiae are discordant with geographic distances but are consistent with a stepwise migration scenario in which the species increases its range from west to east Africa through consecutive founder events over the last ~200,000 years. Geological barriers like the Congo River basin and the East African rift seem to play an important role in shaping this process. Moreover, we find a high degree of genetic isolation of populations on the Comoros, confirming the potential of these islands as candidate sites for potential field trials of genetically engineered mosquitoes for malaria control.