Complete Anopheles funestus mitogenomes reveal an ancient history of mitochondrial lineages and their distribution in southern and central Africa

Distribution and insecticide resistance profile of the major malaria vector Anopheles funestus group across the African continent

There has been significant progress in malaria control in the last 2 decades, with a decline in mortality and morbidity. However, these gains are jeopardised by insecticide resistance, which negatively impacts the core interventions, such as insecticide-treated nets (ITN) and indoor residual spraying (IRS). While most malaria control and research efforts are still focused on Anopheles gambiae complex mosquitoes, Anopheles funestus remains an important vector in many countries and, in some cases, contributes to most of the local transmission. As countries move towards malaria elimination, it is important to ensure that all dominant vector species, including An. funestus, an important vector in some countries, are targeted. The objective of this review is to compile and discuss information related to A. funestus populations' resistance to insecticides and the mechanisms involved across Africa, emphasising the sibling species and their resistance profiles in relation to malaria elimination goals. Data on insecticide resistance in An. funestus malaria vectors in Africa were extracted from published studies. Online bibliographic databases, including Google Scholar and PubMed, were used to search for relevant studies. Articles published between 2000 and May 2023 reporting resistance of An. funestus to insecticides and associated mechanisms were included. Those reporting only bionomics were excluded. Spatial variation in species distribution and resistance to insecticides was recorded from 174 articles that met the selection criteria. It was found that An. funestus was increasingly resistant to the four classes of insecticides recommended by the World Health Organisation for malaria vector control; however, this varied by country. Insecticide resistance appears to reduce the effectiveness of vector control methods, particularly IRS and ITN. Biochemical resistance due to detoxification enzymes (P450s and glutathione-S-transferases [GSTs]) in An. funestus was widely recorded. However, An. funestus in Africa remains susceptible to other insecticide classes, such as organophosphates and neonicotinoids. This review highlights the increasing insecticide resistance of An. funestus mosquitoes, which are important malaria vectors in Africa, posing a significant challenge to malaria control efforts. While An. funestus has shown resistance to the recommended insecticide classes, notably pyrethroids and, in some cases, organochlorides and carbamates, it remains susceptible to other classes of insecticides such as organophosphates and neonicotinoids, providing potential alternative options for vector control strategies. The study underscores the need for targeted interventions that consider the population structure and geographical distribution of An. funestus, including its sibling species and their insecticide resistance profiles, to effectively achieve malaria elimination goals.

Population genetic structure and wing geometric morphometrics of the filarial vector Armigeres subalbatus (Diptera: Culicidae) in Thailand

Armigeres subalbatus (Diptera: Culicidae) is a mosquito species of significant medical and veterinary importance. It is widely distributed across Southeast and East Asia and is commonly found throughout Thailand. This study assessed the genetic diversity and population structure of Ar. subalbatus in Thailand using the cytochrome c oxidase subunit I (COI) gene sequences. Additionally, wing shape variations among these populations were examined using geometric morphometrics (GM). Our results demonstrated that the overall haplotype diversity (Hd) was 0.634, and the nucleotide diversity (π) was 0.0019. Significant negative values in neutrality tests (p < 0.05) indicate that the Ar. subalbatus populations in Thailand are undergoing a phase of expansion following a bottleneck event. The mismatch distribution test suggests that the populations may have started expanding approximately 16,678 years ago. Pairwise genetic differentiation among the 12 populations based on Fst revealed significant differences in 32 pairs (p < 0.05), with the degree of differentiation ranging from 0.000 to 0.419. The GM analysis of wing shape also indicated significant differences in nearly all pairs (p < 0.05), except for between populations from Nakhon Pathom and Samut Songkhram, and between those from Chiang Mai and Mae Hong Son, suggesting no significant difference due to their similar environmental settings. These findings enhance our understanding of the population structure and phenotypic adaptations of mosquito vectors, providing vital insights for the formulation of more efficacious vector control strategies.

Anopheles gambiae on remote islands in the Indian Ocean: origins and prospects for malaria elimination by genetic modification of extant populations

The mosquito Anopheles gambiae s.s. is a primary malaria vector throughout sub-Saharan Africa including the islands of the Comoros archipelago (Anjouan, Grande Comore, Mayotte and Mohéli). These islands are located at the northern end of the Mozambique Channel in eastern Africa. Previous studies have shown a relatively high degree of genetic isolation between the Comoros islands and mainland populations of A. gambiae, but the origin of the island populations remains unclear. Here, we analyzed phylogenetic relationships among island and mainland populations using complete mitochondrial genome sequences of individual A. gambiae specimens. This work augments earlier studies based on analysis of the nuclear genome. We investigated the source population of A. gambiae for each island, estimated the number of introductions, when they occurred and explored evidence for contemporary gene flow between island and mainland populations. These studies are relevant to understanding historical patterns in the dispersal of this important malaria vector and provide information critical to assessing their potential for the exploration of genetic-based vector control methods to eliminate this disease. Phylogenetic analysis and haplotype networks were constructed from mitogenome sequences of 258 A. gambiae from the four islands. In addition, 112 individuals from seven countries across sub-Saharan Africa and Madagascar were included to identify potential source populations. Our results suggest that introduction events of A. gambiae into the Comoros archipelago were rare and recent events and support earlier claims that gene flow between the mainland and these islands is limited. This study is concordant with earlier work suggesting the suitability of these oceanic islands as appropriate sites for conducting field trial releases of genetically engineered mosquitoes (GEMs).

Spatial-temporal vector abundance and malaria transmission dynamics in Nchelenge and Lake Mweru islands, a region with a high burden of malaria in northern Zambia

BACKGROUND

Over a decade of vector control by indoor residual spraying (IRS) and long-lasting insecticidal nets (LLINs) distribution on the mainland, and only LLINs on islands had a minimal impact on disease burden in Nchelenge district, northern Zambia. Anopheles funestus and Anopheles gambiae are vectors known only from the mainland. Understanding vector bionomics in the district is necessary for planning and targeting effective vector control. This study aimed to provide information on abundance, seasonality, and Plasmodium falciparum sporozoite infectivity of malaria vectors in Nchelenge, including islands.

METHODS

Mosquitoes were collected in 192 CDC indoor light traps set in 56 households between January 2015 and January 2016. Morphological and molecular species identifications and P. falciparum circumsporoites by ELISA were performed. Mosquito counts and relative abundances from the islands and mainland were compared, and household factors associated with vector counts were determined.

RESULTS

A total of 5888 anophelines were collected during the study. Of these, 5,704 were female Anopheles funestus sensu lato (s.l.) and 248 female An. gambiae s.l. The highest proportion of An. funestus (n = 4090) was from Chisenga Island and An. gambiae (n = 174) was from Kilwa Island. The highest estimated counts per trap for An. funestus s.l. were from Chisenga island, (89.9, p < 0.001) and from the dry season (78.6, p < 001). For An. gambiae the highest counts per trap were from Kilwa island (3.1, p < 0.001) and the rainy season (2.5, p = 0.007). The highest estimated annual entomological inoculation rate was from Chisenga Island with 91.62 ib/p/y followed by Kilwa Island with 29.77 ib/p/yr, and then Mainland with 19.97 ib/p/yr.

CONCLUSIONS

There was varied species abundance and malaria transmission risk across sites and seasons. The risk of malaria transmission was perennial and higher on the islands. The minimal impact of vector control efforts on the mainland was evident, but limited overall. Vector control intervention coverage with effective tools needs to be extended to the islands to effectively control malaria transmission in Nchelenge district.

Advances in the genetic characterization of the malaria vector, Anopheles funestus, and implications for improved surveillance and control

Anopheles mosquitoes present a major public health challenge in sub-Saharan Africa; notably, as vectors of malaria that kill over half a million people annually. In parts of the east and southern Africa region, one species in the Funestus group, Anopheles funestus, has established itself as an exceptionally dominant vector in some areas, it is responsible for more than 90% of all malaria transmission events. However, compared to other malaria vectors, the species is far less studied, partly due to difficulties in laboratory colonization and the unresolved aspects of its taxonomy and systematics. Control of An. funestus is also increasingly difficult because it has developed widespread resistance to public health insecticides. Fortunately, recent advances in molecular techniques are enabling greater insights into species identity, gene flow patterns, population structure, and the spread of resistance in mosquitoes. These advances and their potential applications are reviewed with a focus on four research themes relevant to the biology and control of An. funestus in Africa, namely: (i) the taxonomic characterization of different vector species within the Funestus group and their role in malaria transmission; (ii) insecticide resistance profile; (iii) population genetic diversity and gene flow, and (iv) applications of genetic technologies for surveillance and control. The research gaps and opportunities identified in this review will provide a basis for improving the surveillance and control of An. funestus and malaria transmission in Africa.

Phylogenetic Analysis of Some Species of the Anopheles hyrcanus Group (Diptera: Culicidae) in China Based on Complete Mitochondrial Genomes

Some species of the Hyrcanus group are vectors of malaria in China. However, the member species are difficult to identify accurately by morphology. The development of sequencing technologies offers the possibility of further studies based on the complete mitochondrial genome. In this study, samples of mosquitoes of the Hyrcanus group were collected in China between 1997 and 2015. The mitochondrial genomes of ten species of the Hyrcanus group were analyzed, including the structure and base composition, codon usage, secondary structure of tRNA, and base difference sites in protein coding regions. Phylogenetic analyses using maximum-likelihood and Bayesian inference were performed based on mitochondrial genes and complete mitochondrial genomes The mitochondrial genome of 10 Hyrcanus group members ranged from 15,403 bp to 15,475 bp, with an average 78.23% (A + T) content, comprising of 13 PCGs (protein coding genes), 22 tRNAs, and 2 rRNAs. Site differences between some closely related species in the PCGs were small. There were only 36 variable sites between Anopheles sinensis and Anopheles belenrae for a variation ratio of 0.32% in all PCGs. The pairwise interspecies distance based on 13 PCGs was low, with an average of 0.04. A phylogenetic tree constructed with the 13 PCGs was consistent with the known evolutionary relationships. Some phylogenetic trees constructed by single coding regions (such as COI or ND4) or combined coding regions (COI + ND2 + ND4 + ND5 or ND2 + ND4) were consistent with the phylogenetic tree constructed using the 13 PCGs. The phylogenetic trees constructed using some coding genes (COII, ND5, tRNAs, 12S rRNA, and 16S rRNA) differed from the phylogenetic tree constructed using PCGs. The difference in mitochondrial genome sequences between An. sinensis and An. belenrae was very small, corresponding to intraspecies difference, suggesting that the species was in the process of differentiation. The combination of all 13 PCG sequences was demonstrated to be optimal for phylogenetic analysis in closely related species.

Expanded geographic distribution and host preference of Anopheles gibbinsi (Anopheles species 6) in northern Zambia

Background

Nchelenge District in northern Zambia suffers from holoendemic malaria transmission despite a decade of yearly indoor residual spraying (IRS) and insecticide-treated net (ITN) distributions. One hypothesis for this lack of impact is that some vectors in the area may forage in the early evening or outdoors. Anopheles gibbinsi specimens were identified in early evening mosquito collections performed in this study area, and further insight was gleaned into this taxon, including characterizing its genetic identity, feeding preferences, and potential role as a malaria vector.

Methods

Mosquitoes were collected in July and August 2019 by CDC light traps in Nchelenge District in indoor sitting rooms, outdoor gathering spaces, and animal pens from 16:00–22:00. Host detection by PCR, COI and ITS2 PCR, and circumsporozoite (CSP) ELISA were performed on all samples morphologically identified as An. gibbinsi, and a subset of specimens were selected for COI and ITS2 sequencing. To determine risk factors for increased abundance of An. gibbinsi, a negative binomial generalized linear mixed-effects model was performed with household-level variables of interest.

Results

Comparison of COI and ITS2 An. gibbinsi reference sequences to the NCBI database revealed > 99% identity to “Anopheles sp. 6” from Kenya. More than 97% of specimens were morphologically and molecularly consistent with An. gibbinsi. Specimens were primarily collected in animal pen traps (59.2%), followed by traps outdoors near where humans gather (24.3%), and traps set indoors (16.5%). Host DNA detection revealed a high propensity for goats, but 5% of specimens with detected host DNA had fed on humans. No specimens were positive for Plasmodium falciparum sporozoites. Animal pens and inland households > 3 km from Lake Mweru were both associated with increased An. gibbinsi abundance.

Conclusions

This is the first report of An. gibbinsi in Nchelenge District, Zambia. This study provided a species identity for unknown “An. sp. 6” in the NCBI database, which has been implicated in malaria transmission in Kenya. Composite data suggest that this species is largely zoophilic and exophilic, but comes into contact with humans and the malaria parasites they carry. This species should continue to be monitored in Zambia and neighbouring countries as a potential malaria vector.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04231-5.

Mitogenome Analyses Reveal Limited Introduction of Anopheles coluzzii Into the Central African Islands of São Tomé and Príncipe

Islands possess physical characteristics that make them uniquely well-suited for initial field trials of new genetic-based technologies applied to African malaria vectors. This has led to efforts to characterize the degree of isolation of island mosquito populations. São Tomé and Príncipe (STP) is a country composed of two small islands in the Gulf of Guinea (Central Africa) where Anopheles coluzzii is the primary malaria vector. Several studies have shown a relatively high degree of genetic isolation between A. coluzzii populations in STP and the mainland compared with pairs of mainland populations separated by equivalent distances. Here, we analyzed complete mitochondrial genomes of individual A. coluzzii specimens from STP and neighboring mainland countries. The objectives are to describe the history of A. coluzzii establishment in STP, specifically to address several questions germane to their suitability as sites for a field trial release of genetically engineered mosquitoes (GEMs). These questions include: (i) What are the origins of A. coluzzii populations in STP?; (ii) How many introductions occurred?; (iii) When was A. coluzzii introduced into STP? and (iv) Is there ongoing, contemporary gene flow into STP from mainland populations? Phylogenetic analysis and haplotype networks were constructed from sequences of 345 A. coluzzii from STP, and 107 individuals from 10 countries on or near the west coast of Africa. Analysis of these data suggest that there have been two introductions of A. coluzzii onto the island of São Tomé that occurred roughly 500 years ago and that these originated from mainland West Africa. It appears that A. coluzzii has never been introduced into Príncipe Island directly from mainland Africa, but there have been at least four introductions originating from São Tomé. Our findings provide further support for the notion that contemporary populations of A. coluzzii on São Tomé and Príncipe are genetically isolated from mainland populations of this mosquito species.

A Minimally Morphologically Destructive Approach for DNA Retrieval and Whole-Genome Shotgun Sequencing of Pinned Historic Dipteran Vector Species

Museum collections contain enormous quantities of insect specimens collected over the past century, covering a period of increased and varied insecticide usage. These historic collections are therefore incredibly valuable as genomic snapshots of organisms before, during, and after exposure to novel selective pressures. However, these samples come with their own challenges compared with present-day collections, as they are fragile and retrievable DNA is low yield and fragmented. In this article, we tested several DNA extraction procedures across pinned historic Diptera specimens from four disease vector genera: Anopheles, Aedes, Culex, and Glossina. We identify an approach that minimizes morphological damage while maximizing DNA retrieval for Illumina library preparation and sequencing that can accommodate the fragmented and low yield nature of historic DNA. We identify several key points in retrieving sufficient DNA while keeping morphological damage to a minimum: an initial rehydration step, a short incubation without agitation in a modified low salt Proteinase K buffer (referred to as "lysis buffer C" throughout), and critical point drying of samples post-extraction to prevent tissue collapse caused by air drying. The suggested method presented here provides a solid foundation for exploring the genomes and morphology of historic Diptera collections.

Radiation with reticulation marks the origin of a major malaria vector

Advances in genomics have led to an appreciation that introgression is common, but its evolutionary consequences are poorly understood. In recent species radiations the sharing of genetic variation across porous species boundaries can facilitate adaptation to new environments and generate novel phenotypes, which may contribute to further diversification. Most Anopheles mosquito species that are of major importance as human malaria vectors have evolved within recent and rapid radiations of largely nonvector species. Here, we focus on one of the most medically important yet understudied anopheline radiations, the Afrotropical Anopheles funestus complex (AFC), to investigate the role of introgression in its diversification and the possible link between introgression and vector potential. The AFC comprises at least seven morphologically similar species, yet only An. funestus sensu stricto is a highly efficient malaria vector with a pan-African distribution. Based on de novo genome assemblies and additional whole-genome resequencing, we use phylogenomic and population genomic analyses to establish species relationships. We show that extensive interspecific gene flow involving multiple species pairs has shaped the evolutionary history of the AFC since its diversification. The most recent introgression event involved a massive and asymmetrical movement of genes from a distantly related AFC lineage into An. funestus, an event that predated and plausibly facilitated its subsequent dramatic geographic range expansion across most of tropical Africa. We propose that introgression may be a common mechanism facilitating adaptation to new environments and enhancing vectorial capacity in Anopheles mosquitoes.

Genetic Diversity of Anopheles coustani (Diptera: Culicidae) in Malaria Transmission Foci in Southern and Central Africa

Despite ongoing malaria control efforts implemented throughout sub-Saharan Africa, malaria remains an enormous public health concern. Current interventions such as indoor residual spraying with insecticides and use of insecticide-treated bed nets are aimed at targeting the key malaria vectors that are primarily endophagic and endophilic. Anopheles coustani s.l., an understudied vector of malaria, is a species previously thought to exhibit mostly zoophilic behavior. Like many of these understudied species, An. coustani has greater anthropophilic tendencies than previously appreciated, is often both endophagic and exophagic, and carries Plasmodium falciparum sporozoites. The aim of this study was to explore genetic variation of An. coustani mosquitoes and the potential of this species to contribute to malaria parasite transmission in high transmission settings in Zambia and the Democratic Republic of the Congo (DRC). Morphologically identified An. coustani specimens that were trapped outdoors in these study sites were analyzed by PCR and sequencing for species identification and bloodmeal sources, and malaria parasite infection was determined by ELISA and qPCR. Fifty An. coustani s.s. specimens were confirmed by analysis of mitochondrial DNA cytochrome c oxidase subunit I (COI) and ribosomal internal transcribed spacer region 2 (ITS2). Maximum likelihood phylogenetic analysis of COI and ITS2 sequences revealed two distinct phylogenetic groups within this relatively small regional collection. Our findings indicate that both An. coustani groups have anthropophilic and exophagic habits and come into frequent contact with P. falciparum, suggesting that this potential alternative malaria vector might elude current vector control measures in northern Zambia and southern DRC.

Genetic differentiation and population structure of Anopheles funestus from Uganda and the southern African countries of Malawi, Mozambique, Zambia and Zimbabwe

Background

Anopheles funestus (s.s.) is a primary vector of the malaria parasite Plasmodium falciparum in Africa, a human pathogen that causes almost half a million deaths each year. The population structure of An. funestus was examined in samples from Uganda and the southern African countries of Malawi, Mozambique, Zambia and Zimbabwe.

Methods

Twelve microsatellites were used to estimate the genetic diversity and differentiation of An. funestus from 13 representative locations across five countries. These were comprised of four sites from Uganda, three from Malawi and two each from Mozambique, Zambia and Zimbabwe.

Results

All loci were highly polymorphic across the populations with high allelic richness and heterozygosity. A high genetic diversity was observed with 2–19 alleles per locus and an average number of seven alleles. Overall, expected heterozygosity (H_e) ranged from 0.65 to 0.79. When samples were pooled three of the 12 microsatellite loci showed Hardy–Weinberg equilibrium. Unsupervised Bayesian clustering analysis of microsatellite data revealed two clusters with An. funestus samples from Mozambique, Uganda and Zambia falling into one group and Malawi and Zimbabwe into another. The overall genetic differentiation between the populations was moderate (F_ST = 0.116). Pairwise differentiation between the pairs was low but significant. A weak but significant correlation was established between genetic and geographical distance for most populations.

Conclusions

High genetic diversity revealed by the loci with low to moderate differentiation, identified two clusters among the An. funestus populations. Further research on the population dynamics of An. funestus in east and southern Africa is essential to understand the implications of this structuring and what effect it may have on the efficient implementation of mosquito vector control strategies.

High Plasmodium falciparum genetic diversity and temporal stability despite control efforts in high transmission settings along the international border between Zambia and the Democratic Republic of the Congo

Background

While the utility of parasite genotyping for malaria elimination has been extensively documented in low to moderate transmission settings, it has been less well-characterized in holoendemic regions. High malaria burden settings have received renewed attention acknowledging their critical role in malaria elimination. Defining the role for parasite genomics in driving these high burden settings towards elimination will enhance future control programme planning.

Methods

Amplicon deep sequencing was used to characterize parasite population genetic diversity at polymorphic Plasmodium falciparum loci, Pfama1 and Pfcsp, at two timepoints in June–July 2016 and January–March 2017 in a high transmission region along the international border between Luapula Province, Zambia and Haut-Katanga Province, the Democratic Republic of the Congo (DRC).

Results

High genetic diversity was observed across both seasons and in both countries. No evidence of population structure was observed between parasite populations on either side of the border, suggesting that this region may be one contiguous transmission zone. Despite a decline in parasite prevalence at the sampling locations in Haut-Katanga Province, no genetic signatures of a population bottleneck were detected, suggesting that larger declines in transmission may be required to reduce parasite genetic diversity. Analysing rare variants may be a suitable alternative approach for detecting epidemiologically important genetic signatures in highly diverse populations; however, the challenge is distinguishing true signals from potential artifacts introduced by small sample sizes.

Conclusions

Continuing to explore and document the utility of various parasite genotyping approaches for understanding malaria transmission in holoendemic settings will be valuable to future control and elimination programmes, empowering evidence-based selection of tools and methods to address pertinent questions, thus enabling more efficient resource allocation.

Characterization of the complete mitogenome of Anopheles aquasalis, and phylogenetic divergences among Anopheles from diverse geographic zones

Whole mitogenome sequences (mtDNA) have been exploited for insect ecology studies, using them as molecular markers to reconstruct phylogenies, or to infer phylogeographic relationships and gene flow. Recent Anopheles phylogenomic studies have provided information regarding the time of deep lineage divergences within the genus. Here we report the complete 15,393 bp mtDNA sequences of Anopheles aquasalis, a Neotropical human malaria vector. When comparing its structure and base composition with other relevant and available anopheline mitogenomes, high similarity and conserved genomic features were observed. Furthermore, 22 mtDNA sequences comprising anopheline and Dipteran sibling species were analyzed to reconstruct phylogenies and estimate dates of divergence between taxa. Phylogenetic analysis using complete mtDNA sequences suggests that A. aquasalis diverged from the Anopheles albitarsis complex ~28 million years ago (MYA), and ~38 MYA from Anopheles darlingi. Bayesian analysis suggests that the most recent ancestor of Nyssorhynchus and Anopheles + Cellia was extant ~83 MYA, corroborating current estimates of ~79–100 MYA. Additional sampling and publication of African, Asian, and North American anopheline mitogenomes would improve the resolution of the Anopheles phylogeny and clarify early continental dispersal routes.