Phylogeny of Anophelinae using mitochondrial protein coding genes

Genome skimming is a low-cost and robust strategy to assemble complete mitochondrial genomes from ethanol preserved specimens in biodiversity studies

Global loss of biodiversity is an ongoing process that concerns both local and global authorities. Studies of biodiversity mainly involve traditional methods using morphological characters and molecular protocols. However, conventional methods are a time consuming and resource demanding task. The development of high-throughput sequencing (HTS) techniques has reshaped the way we explore biodiversity and opened a path to new questions and novel empirical approaches. With the emergence of HTS, sequencing the complete mitochondrial genome became more accessible, and the number of genome sequences published has increased exponentially during the last decades. Despite the current state of knowledge about the potential of mitogenomics in phylogenetics, this is still a relatively under-explored area for a multitude of taxonomic groups, especially for those without commercial relevance, non-models organisms and with preserved DNA. Here we take the first step to assemble and annotate the genomes from HTS data using a new protocol of genome skimming which will offer an opportunity to extend the field of mitogenomics to under-studied organisms. We extracted genomic DNA from specimens preserved in ethanol. We used Nextera XT DNA to prepare indexed paired-end libraries since it is a powerful tool for working with diverse samples, requiring a low amount of input DNA. We sequenced the samples in two different Illumina platform (MiSeq or NextSeq 550). We trimmed raw reads, filtered and had their quality tested accordingly. We performed the assembly using a baiting and iterative mapping strategy, and the annotated the putative mitochondrion through a semi-automatic procedure. We applied the contiguity index to access the completeness of each new mitogenome. Our results reveal the efficiency of the proposed method to recover the whole mitogenomes of preserved DNA from non-model organisms even if there are gene rearrangement in the specimens. Our findings suggest the potential of combining the adequate platform and library to the genome skimming as an innovative approach, which opens a new range of possibilities of its use to obtain molecular data from organisms with different levels of preservation.

Comparison of malaria incidence rates and socioeconomic-environmental factors between the states of Acre and Rondônia: a spatio-temporal modelling study

BACKGROUND

Plasmodium falciparum malaria is a threat to public health, but Plasmodium vivax malaria is most prevalent in Latin America, where the incidence rate has been increasing since 2016, particularly in Venezuela and Brazil. The Brazilian Amazon reported 193,000 cases in 2017, which were mostly confirmed as P. vivax (~ 90%). Herein, the relationships among malaria incidence rates and the proportion of accumulated deforestation were contrasted using data from the states of Acre and Rondônia in the south-western Brazilian Amazon. The main purpose is to test the hypothesis that the observed difference in incidence rates is associated with the proportion of accumulated deforestation.

METHODS

An ecological study using spatial and temporal models for mapping and modelling malaria risk was performed. The municipalities of Acre and Rondônia were the spatial units of analysis, whereas month and year were the temporal units. The number of reported malaria cases from 2009 until 2015 were used to calculate the incidence rate per 1000 people at risk. Accumulated deforestation was calculated using publicly available satellite images. Geographically weighted regression was applied to provide a local model of the spatial heterogeneity of incidence rates. Time-series dynamic regression was applied to test the correlation of incidence rates and accumulated deforestation, adjusted by climate and socioeconomic factors.

RESULTS

The malaria incidence rate declined in Rondônia but remained stable in Acre. There was a high and positive correlation between the decline in malaria and higher proportions of accumulated deforestation in Rondônia. Geographically weighted regression showed a complex relationship. As deforestation increased, malaria incidence also increased in Acre, while as deforestation increased, malaria incidence decreased in Rondônia. Time-series dynamic regression showed a positive association between malaria incidence and precipitation and accumulated deforestation, whereas the association was negative with the human development index in the westernmost areas of Acre.

CONCLUSION

Landscape modification caused by accumulated deforestation is an important driver of malaria incidence in the Brazilian Amazon. However, this relationship is not linearly correlated because it depends on the overall proportion of the land covered by forest. For regions that are partially degraded, forest cover becomes a less representative component in the landscape, causing the abovementioned non-linear relationship. In such a scenario, accumulated deforestation can lead to a decline in malaria incidence.

Higher risk of malaria transmission outdoors than indoors by Nyssorhynchus darlingi in riverine communities in the Peruvian Amazon

BACKGROUND

Malaria remains an important public health problem in Peru where incidence has been increasing since 2011. Of over 55,000 cases reported in 2017, Plasmodium vivax was the predominant species (76%), with P. falciparum responsible for the remaining 24%. Nyssorhynchus darlingi (previously Anopheles darlingi) is the main vector in Amazonian Peru, where hyperendemic Plasmodium transmission pockets have been found. Mazán district has pronounced spatial heterogeneity of P. vivax malaria. However, little is known about behavior, ecology or seasonal dynamics of Ny. darlingi in Mazán. This study aimed to gather baseline information about bionomics of malaria vectors and transmission risk factors in a hyperendemic malaria area of Amazonian Peru.

METHODS

To assess vector biology metrics, five surveys (two in the dry and three in the rainy season), including collection of sociodemographic information, were conducted in four communities in 2016-2017 on the Napo (Urco Miraño, URC; Salvador, SAL) and Mazán Rivers (Visto Bueno, VIB; Libertad, LIB). Human-biting rate (HBR), entomological inoculation rate (EIR) and human blood index (HBI) were measured to test the hypothesis of differences in entomological indices of Ny. darlingi between watersheds. A generalized linear mixed effect model (GLMM) was constructed to model the relationship between household risk factors and the EIR.

RESULTS

Nyssorhynchus darlingi comprised 95% of 7117 Anophelinae collected and its abundance was significantly higher along the Mazán River. The highest EIRs (3.03-4.54) were detected in March and June in URC, LIB and VIB, and significantly more Ny. darlingi were infected outdoors than indoors. Multivariate analysis indicated that the EIR was >12 times higher in URC compared with SAL. The HBI ranged from 0.42-0.75; humans were the most common blood source, followed by Galliformes and cows. There were dramatic differences in peak biting time and malaria incidence with similar bednet coverage in the villages.

CONCLUSIONS

Nyssorhynchus darlingi is the predominant contributor to malaria transmission in the Mazán District, Peru. Malaria risk in these villages is higher in the peridomestic area, with pronounced heterogeneities between and within villages on the Mazán and the Napo Rivers. Spatiotemporal identification and quantification of the prevailing malaria transmission would provide new evidence to orient specific control measures for vulnerable or at high risk populations.

Genetic diversity of Nyssorhynchus (Anopheles) darlingi related to biting behavior in western Amazon

Background

In the Amazon Basin, Nyssorhynchus (Anopheles) darlingi is the most aggressive and effective malaria vector. In endemic areas, behavioral aspects of anopheline vectors such as host preference, biting time and resting location post blood meal have a key impact on malaria transmission dynamics and vector control interventions. Nyssorhynchus darlingi presents a range of feeding and resting behaviors throughout its broad distribution.

Methods

To investigate the genetic diversity related to biting behavior, we collected host-seeking Ny. darlingi in two settlement types in Acre, Brazil: Granada (~ 20-year-old, more established, better access by road, few malaria cases) and Remansinho (~ 8-year-old, active logging, poor road access, high numbers malaria cases). Mosquitoes were classified by the location of collection (indoors or outdoors) and time (dusk or dawn).

Results

Genome-wide SNPs, used to assess the degree of genetic divergence and population structure, identified non-random distributions of individuals in the PCA for both location and time analyses. Although genetic diversity related to behavior was confirmed by non-model-based analyses and F_ST values, model-based STRUCTURE detected considerable admixture of these populations.

Conclusions

To our knowledge, this is the first study to detect genetic markers associated with biting behavior in Ny. darlingi. Additional ecological and genomic studies may help to understand the genetic basis of mosquito behavior and address appropriate surveillance and vector control.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3498-4) contains supplementary material, which is available to authorized users.

Novel molecular evidence of population structure in Anopheles (Kerteszia) bellator from Brazilian Atlantic Forest

Anopheles bellator is a primary malaria vector in the Atlantic Forest. Partial sequences of timeless and Clock genes were used to assess the genetic differentiation of five Brazilian populations, which showed strong population structure (e.g. high F_ST values and fixed differences) in all pairwise comparisons between Bahia sample and the others from Paraná, São Paulo and Rio de Janeiro states. Also, the resulting phylogenetic trees clearly grouped the sequences from Bahia in a different cluster with high bootstrap values. Among southern and southeastern populations low levels of genetic differentiation were found suggesting a general stability of the genetic structure.

Malaria vector species in Amazonian Peru co-occur in larval habitats but have distinct larval microbial communities

In Amazonian Peru, the primary malaria vector, Nyssorhynchus darlingi (formerly Anopheles darlingi), is difficult to target using standard vector control methods because it mainly feeds and rests outdoors. Larval source management could be a useful supplementary intervention, but to determine its feasibility, more detailed studies on the larval ecology of Ny. darlingi are essential. We conducted a multi-level study of the larval ecology of Anophelinae mosquitoes in the peri-Iquitos region of Amazonian Peru, examining the environmental characteristics of the larval habitats of four species, comparing the larval microbiota among species and habitats, and placing Ny. darlingi larval habitats in the context of spatial heterogeneity in human malaria transmission. We collected Ny. darlingi, Nyssorhynchus rangeli (formerly Anopheles rangeli), Nyssorhynchus triannulatus s.l. (formerly Anopheles triannulatus s.l.), and Nyssorhynchus sp. nr. konderi (formerly Anopheles sp. nr. konderi) from natural and artificial water bodies throughout the rainy and dry seasons. We found that, consistent with previous studies in this region and in Brazil, the presence of Ny. darlingi was significantly associated with water bodies in landscapes with more recent deforestation and lower light intensity. Nyssorhynchus darlingi presence was also significantly associated with a lower vegetation index, other Anophelinae species, and emergent vegetation. Though they were collected in the same water bodies, the microbial communities of Ny. darlingi larvae were distinct from those of Ny. rangeli and Ny. triannulatus s.l., providing evidence either for a species-specific larval microbiome or for segregation of these species in distinct microhabitats within each water body. We demonstrated that houses with more reported malaria cases were located closer to Ny. darlingi larval habitats; thus, targeted control of these sites could help ameliorate malaria risk. The co-occurrence of Ny. darlingi larvae in water bodies with other putative malaria vectors increases the potential impact of larval source management in this region.

Venezuela's humanitarian crisis, resurgence of vector-borne diseases, and implications for spillover in the region

In the past 5-10 years, Venezuela has faced a severe economic crisis, precipitated by political instability and declining oil revenue. Public health provision has been affected particularly. In this Review, we assess the impact of Venezuela's health-care crisis on vector-borne diseases, and the spillover into neighbouring countries. Between 2000 and 2015, Venezuela witnessed a 359% increase in malaria cases, followed by a 71% increase in 2017 (411 586 cases) compared with 2016 (240 613). Neighbouring countries, such as Brazil, have reported an escalating trend of imported malaria cases from Venezuela, from 1538 in 2014 to 3129 in 2017. In Venezuela, active Chagas disease transmission has been reported, with seroprevalence in children (<10 years), estimated to be as high as 12·5% in one community tested (n=64). Dengue incidence increased by more than four times between 1990 and 2016. The estimated incidence of chikungunya during its epidemic peak is 6975 cases per 100 000 people and that of Zika virus is 2057 cases per 100 000 people. The re-emergence of many vector-borne diseases represents a public health crisis in Venezuela and has the possibility of severely undermining regional disease elimination efforts. National, regional, and global authorities must take action to address these worsening epidemics and prevent their expansion beyond Venezuelan borders.

Regional variation in life history traits and plastic responses to temperature of the major malaria vector Nyssorhynchus darlingi in Brazil

The primary Brazilian malaria vector, Nyssorhynchus darlingi (formerly Anopheles darlingi), ranges from 0°S-23°S across three biomes (Amazonia, Cerrado, Mata Atlântica). Rising temperatures will increase mosquito developmental rates, and models predict future malaria transmission by Ny. darlingi in Brazil will shift southward. We reared F1 Ny. darlingi (progeny of field-collected females from 4 state populations across Brazil) at three temperatures (20, 24, 28 °C) and measured key life-history traits. Our results reveal geographic variation due to both genetic differences among localities and plastic responses to temperature differences. Temperature significantly altered all traits: faster larval development, shorter adult life and overall lifespan, and smaller body sizes were seen at 28 °C versus 20 °C. Low-latitude Amazonia mosquitoes had the fastest larval development at all temperatures, but at 28 °C, average development rate of high-latitude Mata Atlântica mosquitoes was accelerated and equivalent to low-latitude Amazonia. Body size of adult mosquitoes from the Mata Atlântica remained larger at all temperatures. We detected genetic variation in the plastic responses among mosquitoes from different localities, with implications for malaria transmission under climate change. Faster development combined with larger body size, without a tradeoff in adult longevity, suggests vectorial capacities of some Mata Atlântica populations may significantly increase under warming climates.

High-accuracy detection of malaria vector larval habitats using drone-based multispectral imagery

Interest in larval source management (LSM) as an adjunct intervention to control and eliminate malaria transmission has recently increased mainly because long-lasting insecticidal nets (LLINs) and indoor residual spray (IRS) are ineffective against exophagic and exophilic mosquitoes. In Amazonian Peru, the identification of the most productive, positive water bodies would increase the impact of targeted mosquito control on aquatic life stages. The present study explores the use of unmanned aerial vehicles (drones) for identifying Nyssorhynchus darlingi (formerly Anopheles darlingi) breeding sites with high-resolution imagery (~0.02m/pixel) and their multispectral profile in Amazonian Peru. Our results show that high-resolution multispectral imagery can discriminate a profile of water bodies where Ny. darlingi is most likely to breed (overall accuracy 86.73%- 96.98%) with a moderate differentiation of spectral bands. This work provides proof-of-concept of the use of high-resolution images to detect malaria vector breeding sites in Amazonian Peru and such innovative methodology could be crucial for LSM malaria integrated interventions.

The most efficient malaria vector in the Latin American region is Nyssorhynchus darlingi (formerly Anopheles darlingi). In Amazonian Peru, where malaria is endemic, Ny. darlingi feeds both indoors and outdoors (endophagy, exophagy), depending on the local environment, and rests outdoors (exophily). LLINs and IRS, the most common tools employed for vector control, target endophagic and endophilic mosquitoes. Thus, they are only partially effective against Ny. darlingi. Control of the aquatic stages of vector mosquitoes, larval source management (LSM), targets the most productive breeding sites nearest to human habitation. In four riverine communities, we used drones with high-resolution imagery as a key initial step to analyze water bodies within the estimated flight range of Ny. darlingi, ~ 1 km. We found distinctive spectral profiles for water bodies that were positive versus negative for Ny. darlingi. The methodology and analysis reported here provide the basis for testing whether LSM can be combined successfully with LLINs and IRS to contribute to the elimination of transmission in malaria hotspots in the Amazon.

Nyssorhynchus dunhami: bionomics and natural infection by Plasmodium falciparum and P. vivax in the Peruvian Amazon

BACKGROUND Nyssorhynchus dunhami, a member of the Nuneztovari Complex, has been collected in Brazil, Colombia, and Peru and described as zoophilic. Although to date Ny. dunhami has not been documented to be naturally infected by Plasmodium, it is frequently misidentified as other Oswaldoi subgroup species that are local or regional malaria vectors. OBJECTIVES The current study seeks to verify the morphological identification of Nuneztovari Complex species collected in the peri-Iquitos region of Amazonian Peru, to determine their Plasmodium infection status, and to describe ecological characteristics of their larval habitats. METHODS We collected Ny. nuneztovari s.l. adults in 2011-2012, and Ny. nuneztovari s.l. larvae and adults in 2016-2017. When possible, samples were identified molecularly using cytochrome c oxidase subunit I (COI) barcode sequencing. Adult Ny. nuneztovari s.l. from 2011-2012 were tested for Plasmodium using real-time PCR. Environmental characteristics associated with Ny. nuneztovari s.l. larvae-positive water bodies were evaluated. FINDINGS We collected 590 Ny. nuneztovari s.l. adults and 116 larvae from eight villages in peri-Iquitos. Of these, 191 adults and 111 larvae were identified by COI sequencing; all were Ny. dunhami. Three Ny. dunhami were infected with P. falciparum, and one with P. vivax, all collected from one village on one night. Ny. dunhami larvae were collected from natural and artificial water bodies, and their presence was positively associated with other Anophelinae larvae and amphibians, and negatively associated with people living within 250m. MAIN CONCLUSIONS Of Nuneztovari Complex species, we identified only Ny. dunhami across multiple years in eight peri-Iquitos localities. This study is, to our knowledge, the first report of natural infection of molecularly identified Ny. dunhami with Plasmodium. We advocate the use of molecular identification methods in this region to monitor Ny. dunhami and other putative secondary malaria vectors to more precisely evaluate their importance in malaria transmission.

A Multi-Gene Analysis and Potential Spatial Distribution of Species of the Strodei Subgroup of the Genus Nyssorhynchus (Diptera: Culicidae)

Nyssorhynchus strodei (Root) is an understudied species of potential epidemiological importance, having been found naturally infected in Brazil with Plasmodium falciparum Welch, Plasmodium vivax Grassi & Feletti, and Plasmodium malariae Feletti & Grassi. It belongs to the strodei subgroup that includes Nyssorhynchus albertoi (Unti), Nyssorhynchus arthuri (Unti), Nyssorhynchus rondoni (Neiva & Pinto), Nyssorhynchus striatus (Sant'Ana & Sallum), and three unnamed species, Nyssorhynchus arthuri B, Nyssorhynchus arthuri C, and Nyssorhynchus arthuri D. As the accurate identification of vector species is of fundamental importance for public health entomology, the aim of the study was to provide additional information for the presence of seven species that had been previously misidentified as Ny. strodei. Bayesian phylogenetic analyses using DNA sequences of the genes COI mtDNA, white, CAD and CAT nuclear genes confirmed Ny. albertoi, Ny. arthuri, Ny. strodei, and Ny. striatus as distinct groups within the strodei subgroup and corroborated the presence of three undescribed species under the name Ny. arthuri. Results of the GMYC model analysis corroborated Ny. arthuri B, Ny. arthuri C, and Ny. arthuri D; however, they did not distinguish between Ny. strodei and Ny. albertoi. Predicted distribution of seven species based on maximum entropy in MaxEnt showed that each species has its specific ecological niche suitability.

Complete mitogenome of Anopheles sinensis and mitochondrial insertion segments in the nuclear genomes of 19 mosquito species

Anopheles sinensis is a major malarial vector in China and Southeast Asia. The mitochondria is involved in many important biological functions. Nuclear mitochondrial DNA segments (NUMTs) are common in eukaryotic organisms, but their characteristics are poorly understood. We sequenced and analyzed the complete mitochondrial (mt) genome of An. sinensis. The mt genome is 15,418 bp long and contains 13 protein-coding genes (PCGs), two rRNAs, 22 tRNAs and a large non-coding region. Its gene arrangement is similar to previously published mosquito mt genomes. We identified and analyzed the NUMTs of 19 mosquito species with both nuclear genomes and mt genome sequences. The number, total length and density of NUMTs are significantly correlated with genome size. About half of NUMTs are short (< 200 bp), but larger genomes can house longer NUMTs. NUMTs may help explain genome size expansion in mosquitoes. The expansion due to mitochondrial insertion segments is variable in different insect groups. PCGs are transferred to nuclear genomes at a higher frequency in mosquitoes, but NUMT origination is more different than in mammals. Larger-sized nuclear genomes have longer mt genome sequences in both mosquitoes and mammals. The study provides a foundation for the functional research of mitochondrial genes in An. sinensis and helps us understand the characteristics and origin of NUMTs and the potential contribution to genome expansion.

Exploring malaria vector diversity on the Amazon Frontier

BACKGROUND

Deforestation in the Amazon and the social vulnerability of its settler communities has been associated with increased malaria incidence. The feeding biology of the most important malaria vectors in the region, notably Nyssorhynchus darlingi, compounds efforts to control vectors and reduce transmission of what has become known as "Frontier Malaria". Exploring Anophelinae mosquito diversity is fundamental to understanding the species responsible for transmission and developing appropriate management and intervention strategies for malaria control in the Amazon River basin.

METHODS

This study describes Anophelinae mosquito diversity from settler communities affected by Frontier Malaria in the states of Acre, Amazonas and Rondônia by analysing COI gene data using cluster and tree-based species delimitation approaches.

RESULTS

In total, 270 specimens from collection sites were sequenced and these were combined with 151 reference (GenBank) sequences in the analysis to assist in species identification. Conservative estimates found that the number of species collected at these sites was between 23 (mPTP partition) and 27 (strict ABGD partition) species, up to 13 of which appeared to be new. Nyssorhynchus triannulatus and Nyssorhynchus braziliensis displayed exceptional levels of intraspecific genetic diversity but there was little to no support for putative species complex status.

CONCLUSIONS

This study demonstrates that Anophelinae mosquito diversity continues to be underestimated in poorly sampled areas where frontier malaria is a major public health concern. The findings will help shape future studies of vector incrimination and transmission dynamics in these areas and support efforts to develop more effective vector control and transmission reduction strategies in settler communities in the Amazon River basin.

An Anopheles by Any Other Name …?

The physical characteristics of a nematoceran fly that evokes the image of a mosquito of the genus Anopheles Meigen (Diptera: Culicidae) has existed since the genus was defined in the early part of the twentieth century. Is that image likely to change with the recent proposal to elevate four Neotropical subgenera to generic status based on relationships generated by phylogenetic analyses of mitochondrial protein-coding genes? Reasons for retaining the traditional concept of Anopheles, with its traditionally recognized subgenera, are presented. However, as the ranking of taxa as genera or subgenera is subjective, the choice of rank of a genus-group name may be considered a subjective decision and a matter of user preference.

Mosquitoes (Diptera: Culicidae) From the Northwestern Brazilian Amazon: Araçá River

The mosquito fauna (Diptera: Culicidae) from two remote localities along the Araçá River, within the Municipality of Barcelos, towards the northern border of the Brazilian State of Amazonas, were sampled using CDC, Shannon, Malaise, and Suspended traps, along with net sweeping and immature collections. During June 2010, 111 collections yielded more than 23,500 mosquitoes distributed in 15 genera, representing 119 different species, together with eight morphospecies, which may represent undescribed new taxa. Among the species collected, there is one new distributional record for Brazil and nine new distributional records for the State of Amazonas. With the highest number of species, the genus Culex Linnaeus also had the largest number of individuals followed by Aedes Meigen with the second highest number of species. The most abundant species was Culex (Melanoconion) gnomatos Sallum, Hutchings & Ferreira followed by Culex (Melanoconion) portesi Senevet & Abonnenc, Culex (Culex) mollis Dyar & Knab, Aedes (Ochlerotatus) fulvus (Wiedemann), Culex (Melanoconion) pedroi Sirivanakarn & Belkin, Culex (Melanoconion) crybda Dyar, Aedes (Ochlerotatus) nubilus (Theobald), and Anopheles (Anopheles) peryassui Dyar & Knab. The epidemiological implications of mosquito species found are discussed and are compared with other mosquito inventories from the Amazon region. As the first standardized mosquito inventory of the Araçá River, with the identification of 127 species level taxa, the number of mosquito species which have been collected along the northern tributaries of the middle Rio Negro Basin (i.e., Padauari and Araçá Rivers) increased significantly.

Abundance of impacted forest patches less than 5 km2 is a key driver of the incidence of malaria in Amazonian Brazil

The precise role that deforestation for agricultural settlements and commercial forest products plays in promoting or inhibiting malaria incidence in Amazonian Brazil is controversial. Using publically available databases, we analyzed temporal malaria incidence (2009–2015) in municipalities of nine Amazonian states in relation to ecologically defined variables: (i) deforestation (rate of forest clearing over time); (ii) degraded forest (degree of human disturbance and openness of forest canopy for logging) and (iii) impacted forest (sum of deforested and degraded forest patches). We found that areas affected by one kilometer square of deforestation produced 27 new malaria cases (r² = 0.78; F1,10 = 35.81; P < 0.001). Unexpectedly, we found both a highly significant positive correlation between number of impacted forest patches less than 5 km² and malaria cases, and that these patch sizes accounted for greater than ~95% of all patches in the study area. There was a significantly negative correlation between extraction forestry economic indices and malaria cases. Our results emphasize not only that deforestation promotes malaria incidence, but also that it directly or indirectly results in a low Human Development Index, and favors environmental conditions that promote malaria vector proliferation.

Mitochondrial genomes of Anopheles arabiensis, An. gambiae and An. coluzzii show no clear species division

Here we report the complete mitochondrial sequences of 70 individual field collected mosquito specimens from throughout Sub-Saharan Africa. We generated this dataset to identify species specific markers for the following Anopheles species and chromosomal forms: An. arabiensis, An. coluzzii (The Forest and Mopti chromosomal forms) and An. gambiae (The Bamako and Savannah chromosomal forms).  The raw Illumina sequencing reads were mapped to the NC_002084 reference mitogenome sequence. A total of 783 single nucleotide polymorphisms (SNPs) were detected on the mitochondrial genome, of which 460 are singletons (58.7%). None of these SNPs are suitable as molecular markers to distinguish among An. arabiensis, An. coluzzii and An. gambiae or any of the chromosomal forms. The lack of species or chromosomal form specific markers is also reflected in the constructed phylogenetic tree, which shows no clear division among the operational taxonomic units considered here.