Next Article in Journal
The Effect of Feeding on Briareum violacea Growth, Survival and Larval Development under Temperature and Salinity Stress
Next Article in Special Issue
Modeling of Valeriana wallichii Habitat Suitability and Niche Dynamics in the Himalayan Region under Anticipated Climate Change
Previous Article in Journal
Immuno-PET Imaging of Atherosclerotic Plaques with [89Zr]Zr-Anti-CD40 mAb—Proof of Concept
Previous Article in Special Issue
Spatiotemporal Variations of Plague Risk in the Tibetan Plateau from 1954–2016
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 11
Issue 3
10.3390/biology11030409
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Phylodynamic and Spread of the Invasive Asian Malaria Vectors, Anopheles stephensi, in Sudan

by
Mustafa Abubakr
1,
Hamza Sami
1,
Isam Mahdi
1,
Omnia Altahir
2,
Hanadi Abdelbagi
2,
Nouh Saad Mohamed
2 and
Ayman Ahmed
1,2,3,4,5,*
1
Directorate of the Integrated Vector Management (IVM), Federal Ministry of Health, Khartoum 11111, Sudan
2
Molecular Biology Unit, Sirius Training and Research Centre, Khartoum 11111, Sudan
3
Institute of Endemic Diseases, University of Khartoum, Khartoum 11111, Sudan
4
Swiss Tropical and Public Health Institute (Swiss TPH), 4123 Allschwil, Switzerland
5
Faculty of Science, University of Basel, Petersplatz 1, 4001 Basel, Switzerland
*
Author to whom correspondence should be addressed.
Biology 2022, 11(3), 409; https://doi.org/10.3390/biology11030409
Submission received: 22 January 2022 / Revised: 14 February 2022 / Accepted: 17 February 2022 / Published: 7 March 2022
(This article belongs to the Special Issue The Impact of Climate Change on Ecosystems and Human Health)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Climate change increases the suitability of some environments for the establishment of newly introduced insects and is a major driver for the spread of mosquitoes that transmit diseases. Anopheles stephensi is a disease vector mosquito that transmits malaria and is naturally endemic in Asia. This vector newly emerged in Africa (first report from Djibouti in 2012), where annual malaria infections and deaths are the highest worldwide. This mosquito has different ecology and behavior from previously known malaria vectors in Africa, which makes control difficult for local under-resourced health systems. Considering the capacity of this vector to transmit at least two malaria-causing parasites (Plasmodium falciparum and Plasmodium vivax), we investigated its distribution and population structure in Sudan and assessed the potential risk of its further spread into neighboring countries. Using morphological and genomic sequencing techniques, we confirmed the presence of Anopheles stephensi along the borders of six countries previously assumed as free, including Chad, Egypt, Eritrea, Libya, Republic of Central Africa, and South Sudan. African countries need to enhance vector surveillance and control services and utilize genomics tools for tracking the dynamics of invasive disease vectors.

Abstract

Anopheles stephensi is an invasive Asian malaria vector that initially emerged in Africa in 2012 and was reported in Sudan in 2019. We investigated the distribution and population structure of An. stephensi throughout Sudan by using sequencing and molecular tools. We confirmed the presence of An. stephensi in eight border-states, identifying both natural and human-made breeding sites. Our analysis revealed the presence of 20 haplotypes with different distributions per state. This study revealed a countrywide spread of An. stephensi in Sudan, with confirmed presence in borders states with Chad, Egypt, Eritrea, Ethiopia, Libya, Republic of Central Africa, and South Sudan. Detection of An. stephensi at points of entry with these countries, particularly Chad, Libya, and South Sudan, indicates the rapid previously undetected spread of this invasive vector. Our phylogenetic and haplotype analysis suggested local establishment and evolutionary adaptation of the vector to different ecological and environmental conditions in Sudan. Urgent engagement of the global community is essential to control and prevent further spread into Africa.

Graphical Abstract

1. Introduction

Malaria is a serious life-threatening illness that is a public health risk for humans mainly in rural tropical areas worldwide. In 2019, up to 252 million malaria cases and 460,000 related deaths were estimated globally [1]. Africa suffers the highest burden of malaria, with over 94% of malaria cases and related deaths reported from African countries [1]. This proportionally high risk of malaria in Africa is attributed to the existence of several competent vectors, including different species of An. gambiae and An. funestus complexes [2,3]. More importantly, pregnant women and children under 5 years of age are the most vulnerable groups affected by the disease, with around 67% malaria reported deaths among children under the age of 5 years [1]. Malaria control programs worldwide rely on two approaches for the prevention and control of the disease: chemotherapy for case management and insecticide-based vector control tools (insecticide-treated bed nets (ITNs), and indoor residual spraying (IRS)) [4]. However, both approaches are challenged by rapidly growing resistance to antimalaria drugs [5] and insecticide resistance among disease vectors [6]. The recently approved malaria vaccine might offer a better alternative for protecting people at high risk of malaria infection [7]. The recent World Malaria Report 2021 highlighted an increase in malaria cases (14 million) and deaths (47,000), mainly attributed to the ongoing COVID-19 pandemic [8].
Malaria is an important threat to public health throughout Sudan, with 1% of the world malaria cases and deaths reported there [2]. According to the 2019 World Malaria Report, around 4.4 million malaria cases were estimated in Sudan, contributing more than 45% of the total number of cases reported in the World Health Organization (WHO) Eastern Mediterranean Region (EMRO). About 13,000 malaria-related deaths were estimated in the country [2]. Unfortunately, morbidity and mortality of malaria in Sudan increased by nearly 50% between 2015 and 2019 [2]. Although the country’s entire population lives with the risk of malaria, endemicity and burden are heterogeneously distributed [9]. This variation in malaria transmission can be attributed to several risk factors, including varied distribution and composition of competent vectors of malaria per different State [3,10], ecology, environment, climate change [3,11,12], conflicts, and human population displacement [9,11]. The composition of malaria vectors in Sudan includes different members of the An. gambiae and An. funestus complexes, An. pharoensis, and An. Arabiensis [2,3,10]. The distribution and spread of disease and disease vectors could be significantly influenced by climate change [12].
Recent emergence of the malaria vector An. stephensi in Sudan carries a serious threat of malaria epidemics and may be an additional risk factor contributing to the current increase in malaria morbidity and mortality in Sudan [9,13,14,15]. Emergence of this invasive malaria vector in Djibouti was associated with an epidemic [16,17]. After initial reports on emergence of An. stephensi in Sudan in 2019 [14,15], unusual increases in malaria cases were reported in different regions, suggesting the involvement of An. stephensi, similarly to Djibouti [1,2,9]. This study was prompted by the alarming WHO alert encouraging countries in the Horn of African to enhance their vector surveillance to delineate the spread of this invasive malaria vector in the region [18]. We report our findings related to the spread of An. stephensi in Sudan and its presence and absence in regional representative States.

2. Materials and Methods

2.1. Mosquito Collection

We followed two strategies to collect mosquito samples, (1) active surveys targeting An. stephensi mosquitoes and (2) routine vector surveillance, which aim to collect and detect vectors for all diseases. To collect Anopheles mosquitoes, we used aspirators and BG-Sentinel traps for adults and dippers for aquatic stages (larvae and pupae).

2.2. Active Surveys Targeting An. stephensi

We implemented active An. stephensi-targeted surveys in States considered at high risks for emergence. The risk of An. stephensi emergence in each State was assessed based on three factors: (a) previous reports about vector presence in the State (Gedaref, Kassala, and the Red Sea States), (b) sharing a border with a State or country with confirmed vector presence (the Blue Nile and Sennar states), and (c) unusual increases in malaria cases (North Darfur State) (Figure 1).

2.3. Routine Vector Surveillance

In addition to active surveys, the national department of Integrated Vector Management (IVM) conducted a series of extensive trainings for national public health officers on collection, identification, and control of disease vectors throughout the 18 States. Routine vector surveillance operating throughout the country captured the presence of An. stephensi in Sennar, Gezira, Khartoum, and River Nile States.

2.4. The Morphological Identification of Mosquito Samples

We transported the mosquito larvae that we have collected to the insectary at the National Public Health laboratory where we reared them to adults. We have used the standard morphological keys for Afrotropical Anopheles mosquitoes for the morphological identification of the mosquito samples to species level [19].

2.5. DNA Extraction from Mosquito

Following manufacturer instructions, we extracted the total DNA from the mosquito samples using QiaAmp tissue extraction kits (Qiagen, Hilden, Germany). We checked the quality of the DNA using a nanodrop spectrophotometer (ND1000, Houston, TX, USA) following the manufacturer’ guidelines and then preserved them at −20 °C until molecular examination.

2.6. Polymerase Chain Reaction (PCR)

We have used Folmer primers (LCO1490 and HCO2198) to amplify the cytochrome oxidase 1 (CO1) region of the mitochondrial DNA of the mosquito genomes in a thermocycler PCR machine (Applied Biosystems, ThermoFisher Scientific, Budapest, Hungary) [20]. We have run PCR by taking 2 µL of the extracted DNA and adding to a 4 µL PCR master mix (Solis Biodyne, Tartu, Estonia), which consisted of 1 U DNA polymerase, 12.5 mM MgCl2, and 4 mM dNTPs. PCR cycling conditions were as follows: initial denaturation at 95 °C for 5 min, 35 cycles of denaturation at 95 °C for 30 s, annealing at 58 °C for 30 s, and extension at 72 °C for 30 s, with a final extension step of 72 °C for 10 min. Following amplification, we visualized PCR products on gel electrophoresis (Major Sciences, Saratoga, CA, USA.) by loading PCR products on 2% agarose gel and placed them onto 150 V and 75 A for 1 h. The amplified PCR products were sequenced after being checked for band quality, including band sharpness and thickness intensity under UV-light using the Sanger Deoxyribonucleic acid sequencing method by 3730XL DNA analyzer (Applied Biosystems, Waltham, MA, USA) through Macrogen company (Macrogen Inc, Amsterdam, The Netherlands).

2.7. Sequences’ Identity Confirmation

We have checked the identity of our sequences of the cytochrome c oxidase 1 region of the mitochondrial DNA by comparing their similarity with the worldwide-published sequences of An. stephensi, using the online BLAST nucleotide algorithm available in the NCBI GenBank database (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 1 October 2021). We have deposited our sequences of An. stephensi that we have obtained during this study into the GenBank database at NCBI under accession numbers OK216338 to OK216737.

2.8. Bioinformatics Analysis

We have used GENtle software (v1.9.4) to check the correctness of sequencing peaks and nucleotides base-calling errors that could occur during sequencing [21]. We trimmed the primer’s sequences to reduce sequencing mismatching that occurs at the start of the sequencing process [22]. We analysed the sequences to identify novel cox1 gene sequence polymorphism using MEGA7 software by aligning the sequences with the reference sequence KT899888.1; it belongs to the An. stephensi strain Hor that was originally collected from India and reared in the Third Military Medical University, China [22,23]. We considered the nucleotide substitution model with the lowest Bayesian Information Criterion (BIC) scores as the best-fit model; Tamura 3-parameter model [24]. We have modelled the non-uniformity of evolutionary rates among sites using a discrete Gamma distribution for creating the phylogenetic tree [24]. We have obtained the sequence diversity parameters including the number of haplotypes (Hap), haplotypes diversity (Hapd), segregating sites (S), and average number of nucleotide differences between two sequences (Pi) using the software DnaSP v5.10 [25]. We have constructed haplotypes’ network by developing the median-joining network using popART software (v4.8) (http://popart.otago.ac.nz, accessed on 1 October 2021). To test the natural selection theory and the population expansion of An. stephensi, we performed the FuFs statistics and Tajima’s D tests. We have estimated the degree of the An. stephensi genetic diversity by comparing our sequences to previously published sequences and calculated the pairwise fixation index (Fst) using DnaSP [25].

3. Results

Our findings confirm the widespread invasive malaria vector, An. Stephensi, in Sudan. Additionally, we identified 20 haplotypes among the An. stephensi populations that we have collected from different states (Figure 1).

3.1. Surveillance Data

We collected mosquito samples by using active surveys from six States distributed throughout Sudan, namely Gedaref, Kassala, Red Sea, River Nile, North Darfur, and Sennar States. We preserved and shipped the collected samples to Khartoum, where we have identified the presence of An. stephensi morphologically and molecularly. We confirmed the identity of our sequence by performing a BLAST search on the NCBI database. We have collected four different species of Anopheles mosquitoes from the different states, namely An. arabiensis, An. macupalpis, An. rufipes, and An. stephensi (Table 1). However, the proportion of An. stephensi from the total collected Anopheles mosquito samples differs per State, with the highest percentage in the Red Sea States (100%), 94% in Gedarif, and only 2% in North Darfur. We did not detect An. stephensi in the Blue Nile State (Table 1) and (Figure 1). Additionally, the presence of An. stephensi was later confirmed through routine surveillance in four new States, including Gezira (1 sample), Khartoum (five samples), Northern (three samples), North Kordofan (one sample), and South Darfur (seven samples) (Figure 1).

3.2. Breeding Sites of An. stephensi

Several different types of breeding sites were positive for the presence of An. stephensi during the larvae collection, as shown in Figure 2. Cement, clay, and plastic containers were the most common and productive containers in the eastern region, while iron containers (e.g., barrels) were the most productive containers in the western region of the country (Figure 2).

3.3. Bioinformatic Analysis

3.3.1. Haplotype Analysis

Results of the sequence alignment compared to the reference sequence KT899888.1 showed the presence of different nucleotides polymorphisms at different sites in the CO1 gene. However, all detected polymorphisms and nucleotides substitutions were nonsynonymous substitutions (Figure 3). We have identified 20 haplotypes (named as Hap 01-20) among the collected Sudanese populations of An. Stephensi, with a Hapd of 0.6091 ± 0.00077. The number of segregating sites (S) detected was 21, and the average number of nucleotide differences between two sequences (Pi) was 0.00486. The most frequent Hap detected was Hap01, which was detected among 245 (61.3%) samples, followed by Hap02 at 38 (9.5%) and Hap03 at 23 (5.8%). The remaining haplotypes constituted a low frequency (Figure 3).

3.3.2. Phylogenetic Analysis

We have constructed the phylogenetic tree based on the Tamura 3-parameter model by adding reference sequences to investigate sequence relations based on origins. Sequences added for the analysis are presented in Supplementary Files (Table S1). The phylogenetic tree showed that Hap01, 02, 04, 05, 06, and 13 were all clustered with sequences reported from Pakistan, Sri Lanka, India, and the previously deposited sequences from Sudan. Hap03 clustered with sequences from Pakistan and Sudan. Hap14, 15, 16, 17, 18, and 19 were clustered separately, while Hap07, 08, 09, 10, 11, and 12 were clustered with sequences reported only from Saudi Arabia (Figure 4)

3.3.3. Worldwide Haplotypes Analysis

We trimmed the ends of the obtained sequences to standardize the length of sequences before performing sequence alignment, which has resulted in 19 different haplotypes. The Hapd of the worldwide sequences was 0.5850 ± 0.00064. The average number of pairwise nucleotide differences was 1.338, and nucleotide diversity was 0.00385. Further, we grouped sequences based on origin and analyzed them to investigate genetic diversity and the natural selection theory (Table 2). Countries with the highest number of haplotypes detected were Pakistan and Saudi Arabia, with four haplotypes each. However, study sequences showed a high number of haplotypes among the Red Sea, Kassal, and Al Gedarif states, with 11 Haps. Their Hapd was 0.539 ± 0.00169, 0.579 ± 0.00339, and 0.64 ± 0.00279, respectively. Although we have identified 11 different haplotypes, only a few diversities were present, indicating that the population in these three states is partially conserved or their divergence has started recently. However, we observed the complete divergence between the haplotypes for haplotypes present in North Darfur state where Hapd was 1.0 ± 0.00926. The natural selection theory tests, Tajima D, and FuFs tests revealed that the only significance of the FuFs test was for Pakistan (Table 2). For most of the analyzed sequences, Tajima D and FuFs tests were negatively insignificant (p value > 0.05) in North Darfur, Red Sea, Kassala, and Gedarif (Table 2).
We computed the pairwise Fst test to estimate the degree of gene flow among the different populations. The values of Fst test among the populations ranged from 0.004 to 0.927 and were all statistically significant (p value < 0.05) (Table S2. in Supplementary Files). Values of Fst test were not very high between Sudanese An. stephensi sequences. A high value of Fst test indicates that most of the populations were genetically differentiated (Table 3).
We clustered all sequences belonging to An. stephensi from Sudan in a single population to reduce the bias in population genetic differentiation. The values of Fst test for the genetic differentiation were also statistically significant, p value < 0.05 (Table S3 in Supplementary Files).
The constructed haplotype network shows the existence of a shared haplotype, Hap01 (Figure 5). Hap01 was a major haplotype in all the populations and was consistently placed in the center of the haplotype networks. Only sequences from Saudi Arabia, Iran, and Ethiopia did not include Hap01 (Table S4 in Supplementary Files). The haplotype network exhibited a typical star-like expansion from the main founder Hap01, with the presence of several unique haplotypes presented in certain populations (Figure 5).
Analyzing the distribution of the different haplotypes that we have identified in this study revealed the predominant distribution of Hap01 throughout most of the country. Hap17 and Hap18 were confined to the River Nile state; meanwhile, Hap15 and Hap16 were confined to the North Darfur state (Figure 6).
In Supplementary Files Table S4, we show the distribution of the different haplotypes among different regions.

4. Discussion

We report the large spread of the invasive malaria vector An. stephensi in Sudan and the high haplotype diversity, indicating vector local establishment and adaptation to the different environments in the country (Figure 6). The emergence and spread of An. stephensi in Africa, including Sudan, is of global importance because it threatens regional public health [14,15,16,26]. This was further indicated by the prompt release of vector alert by WHO in 2019 [18]. The spread of this invasive vector carries the threat of malaria epidemics in urban settings of Africa, particularly in densely populated cities such as Khartoum, Kadugli, and Al Fashir [9,13,14]. Documentation of emergence of An. stephensi in Africa is recent, first in Djibouti in 2012 [16] and then in Ethiopia in 2016 [26].
Since the first report of the emergence in country in 2019 [14,15], widespread An. stephensi across Sudan was suggested by epidemiological reports that highlighted recent increases in malaria cases and epidemics in some regions, particularly in the southeast and southwest [9,11]. This is similar to the unusual increase in malaria cases in Djibouti associated with the first emergence of An. stephensi [16,17].
Climate change is a major driver for the spread of invasive disease vectors as it increases suitability of new environments for novel vectors to establish and eventually accelerate the invasion [27,28]. The influence of climate change on the emergence and distribution of vector-borne diseases and their outbreaks in Sudan is rapidly growing [29]. This is further underscored by the emergence of several arboviruses and their vectors including dengue, Crimean–Congo Hemorrhagic Fever (CCHF), Chikungunya, and Rift Valley fever in different regions of the country [11,12,30,31,32,33,34]. Therefore, the current rapid spread of An. stephensi throughout Sudan could be attributed to climate change. In particular, we report the presence of this invasive vector in areas located at more than 1000 km from the nearest point that is predicted to be environmentally suitable for the spread of An. stephensi [35]. More importantly, local and international studies have reported that climate change is a major driver for the increase in the transmission of malaria and other vector-borne diseases [36,37]. As a sub-Saharan country where most of the environment throughout the country constitutes deserts and semi-deserts, the survival of diseases vectors is mainly limited to microenvironments surrounding water bodies with grassy covers such as banks of rivers [3]. In such environments, the transmission of vector-borne diseases including malaria is climate dependent, with factors such as rainfall, maximum temperature, relative humidity, and the level of river and flooding having direct and significant influences on malaria transmission and epidemiology [36]. Unfortunately, during recent years, Sudan has suffered from extreme weather events, including heavy rains, severe flooding, and rainstorms, throughout the country [12]. The direct impacts of these climate change phenomena on malaria transmission include but are not limited to the substantial increase in suitable breeding sites, increasing the environmental suitability for vector survival beyond the parasite’s incubation period inside the vectors (the extrinsic incubation period) long enough to breed [38,39,40,41]. Furthermore, heavy rains and flooding are hindering public health services and interventions by limiting accessibility to implementation sites and direct interfering with interventions, such as washing out insecticides [41]. This is particularly alarming when reviewed with our finding of An. stephensi aquatic stages (larvae and pupae) in rainwater ponds on the ground, confirming the adaptability of this vector to locally available environments (Figure 2).
Studies in Djibouti and Ethiopia indicated the role of An. stephensi in changing malaria transmission and epidemiology in the region, which is further underscored by several malaria epidemics in the area [9,13,17,18]. The risk of malaria epidemics in the urban settings in Africa is predicted if An. Stephensi spread into these settings [13,42]. However, several epidemics of malaria have occurred in Sudan in association with extreme weather events, but limited or no entomological investigations were carried out. An epidemic of malaria occurred in the Gezira state in 2013 following heavy rains and flash flood [38]. In correspondence with unprecedented heavy rains and flooding throughout Sudan, malaria cases in the country reached the epidemic threshold; however, the majority of cases were reported among the fragile populations living in humanitarian areas in the war-torn regions of Kordofan and Darfur [9,13,43]. Over 110, 103, and 45 thousand cases of malaria were reported from South, North, and East Darfur states, respectively [43].
Field reports and experimental studies confirmed the competency of An. stephensi populations in Africa to transmit both Plasmodium falciparum and P. vivax [18,44,45]. However, little is known about the dynamics and routes of vector introduction in Sudan and other countries in the Horn of Africa [14,15], although a recent study has suggested that An. stephensi in Sudan was introduced from Ethiopia [15]. However, Hap01 has never been reported from Ethiopia, yet it is the major haplotype in Sudan from which other haplotypes might have originated. This suggests that An. stephensi in Sudan has not been introduced from Saudi Arabia, Iran, nor Ethiopia due to the lack of Hap01 in their populations of An. stephensi (Table S4). Our median-joining haplotype network indicates that An. stephensi in Sudan was most likely introduced from countries that share the common haplotype Hap01, namely Pakistan, India, and Sri Lanka (Figure 5). However, we exclude Sri Lanka as a potential source for the origin of this vector because it has been reported there in 2017 as an invasive vector rather than endemic, while evidence confirmed its presence in Sudan since 2016 [15,44]. The extreme gaps in knowledge about bionomics, actual distribution, and susceptibility of An. stephensi populations with respect to locally available vector control tools in Africa urge the need for the implementation of research to generate evidence from the field in order to guide intervention and inform policymakers [14,15,16,42]. High coordination and timely data sharing are essential for improving the implementation of vector control strategies and response plans [45]. It is critical to enhance the implementation of International Health Regulations (IHRs 2005), particularly at points of entry/exit between current An. stephensi free and already invaded areas [18,45]. Considering the large size of the country and the borders with seven countries, including Chad, Central African Republic, Egypt, Eritrea, Ethiopia, Libya, and South Sudan, the current spread of An. stephensi in the country is very alarming nationally and regionally. We confirm the presence of An. stephensi in eight States with open international borders (Figure 1). The local establishment of An. stephensi in Sudan is further highlighted by the high haplotypic variation (20 haplotypes) detected among the collected populations of An. stephensi (Figure 3, Figure 4 and Figure 5).
Previous studies suggested that An. stephensi was only breeding in artificial water containers in African countries [46,47]. Our study also found An. stephensi larvae in ground breeding sites resulting from leaking water supply pipelines, in agreement with reports from Iran [48]. Variations in the nature of the breeding sites might be attributed to vector adaptation to local ecological, climate, and environmental conditions. Similar preference variation in breeding sites was observed for Aedes aegypti. In East Sudan, Ae. aegypti larvae were commonly found in high density in clay and cement containers with lower densities in iron containers, while the reverse was found for Ae. aegypti in West Sudan (Ahmed et al., unpublished data). This could be attributed to the fact that during the rainy/transmission season in West Sudan, the weather is relatively colder. Heat-absorbent iron containers with warmer water might be more attractive for egg-laying female mosquitoes (influencing breeding site selection) or might increase the hatchability of eggs and eventually the productivity of the breeding site [49,50]. Such an adaptation to on-ground breeding could be how the spread of An. stephensi in Sudan is influenced by climate change, specifically heavy rain and flooding [12,27]. In our study, on-ground concrete tanks for water storage were the most productive breeding-site (Figure 3).
Recent malaria epidemics have occurred in areas of Sudan with relatively high coverage of the vector control tools, indoor residual house spraying, and long-lasting insecticidal nets [4,9,13]. The biological and/or behavioral resistance of An. stephensi populations in Africa to the currently used insecticides was reported by WHO and recent studies [8,51]. Therefore, the Ministries of Health in Africa need to restructure their vector surveillance and control systems to become more capable in preventing, detecting, and controlling diseases vectors and proactively response to any change in their species composition [52]. The national malaria control programmes in the area should modify their vector surveillance systems to be more vigilant for the emergence and spread of invasive disease vectors and include house inspection for vectors breeding in human-made water containers into their package of routine surveillance [42,53]. Vector surveillance systems should incorporate genomics sequencing and molecular tools into their vector identification protocol to avoid former delay in the detection of novel vectors due to morphological similarity [15,42]. While IVM programmes could implement larvae source management to control An. stephensi and the co-existing Ae. Aegypti, however, the co-effectiveness of this approach should be systematically evaluated in the field. Governments and local and international partners of the health cluster including WHO and Global Fund should support the implementation of more cost-effective and environmentally friendly (non-insecticide based) vector control interventions [52]. Particular attention should be provided for protecting women, children under five years old, and poor and fragile communities, especially those living in humanitarian settings [9,12,13,52].
As part of the action plan and response strategy to the emergence of An. stephensi in Sudan (Vector Control Strategic and Response Plan for Sudan 2021–2025), the IVM department and the Federal Ministry of Health have mobilized their local resources to implement active vector surveillance to confirm the presence or absence of An. stephensi in specific States of high risk (The Red Sea, Kassala, Gedaref, River Nile, Sinnar, and North Darfur States) (Figure 1). This risk assessment was carried out based on vector presence in/nearby the specific area and is based on the epidemiological recording of malaria cases (Sudan Strategic Plan for the Surveillance and Control of An. stephensi in Sudan, 2021–2025).
Interestingly, all areas where active surveys were implemented were positive for the presence of An. Stephensi, except for the Blue Nile State (Figure 1). In parallel, the IVM department invested in expanding vector surveillance and control by training national surveillance officers on the surveillance and control of different disease vectors instead of only focusing on An. arabiensis. Unfortunately, it seems that the previous lack of implementing integrated vector surveillance allowed this vector to spread throughout the country undetected and become established [14,15].
Alarmingly, an experimental study has shown that the infection of An. stephensi with malaria parasites facilitate Rift Valley fever viruses [54]. This could justify the recently reported changes and emergence events of Rift Valley fever in Sudan [30,33,34]. Similarly, another independent experimental study revealed the capacity of An. stephensi to transmit Chikungunya virus [55], and this is further supported by the correspondence of the massive outbreak of Chikungunya in Sudan in the areas heavily manifested by An. stephensi [30]. More importantly, if An. stephensi was confirmed to transmit the Chikungunya virus in field conditions, the recent large outbreak of Chikungunya (27,540 cases reported) might be an indicator for the spread of this vector into Chad [56]. Therefore, entomological investigations are essential during outbreaks of arboviral diseases in the area to determine the role of An. stephensi in the transmission of these deadly and economically devastating diseases.
Although vector-borne diseases constitute a major threat to public health in Sudan, surveillance and control of disease vector services in the country are limited by several persistent issues, including poverty, poor sanitation, unplanned urbanization, the underfunded and weak health system, and the high turnover of trained personnel [1,53]. Furthermore, limited resources, relatively high rate of conflicts, humanitarian crisis, massive human population displacement, and climate change are the major drivers for the emergence and outbreaks of vector-borne diseases in Sudan [11,30,31,32,33,34].

5. Conclusions

The invasive malaria vector An. stephensi is well spread throughout Sudan and has potentially been introduced into neighboring assumedly An. Stephensi-free countries already. The population structure, haplotypic variation, and change in breeding site preferences of An. stephensi in Sudan indicate adaptation and establishment and suggest that the invasion occurred earlier than 2016, possibly through multiple introductions most likely from Pakistan and/or India. The rapid expansion of the geographical distribution of this vector emphasizes the need for global engagement with all relevant stakeholders, including donors and technical supporters such as different UN agencies, national programs, and the Ministry of Health to invest early in the control and elimination of this invasive vector before it fully adapts and become endemic. Furthermore, countries in Central and West Africa and their health partners should be alarmed by the potential risk of introducing this invasive vector into their areas, which further necessitates strict implementation of IHRs to control disease vectors, particularly at points of entry/exit. Additionally, vector surveillance and control systems in Africa should shift from doing business as usual surveillance to be more vigilant for the spread of invasive diseases vectors through the incorporation of genomics and molecular tools into their vector surveillance. Moreover, health stakeholders, mainly WHO and Global Fund, should invest more and support the implementation of the environmentally friendly non-insecticides-based vector control interventions such as larvae source management, which might be a more cost-effective tool for integrated vector management.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology11030409/s1, Table S1: The accessions number of sequences of An. stephensi that are currently available in the GenBank and were included in the analysis and comparison of our sequences. Table S2: p-values of the pairwise Fst test was computed to estimate the degree of gene flow among the different populations. Table S3: Pairwise fixation index (Fst values) between An. stephensi populations calculated from the nucleotide data set of the mitochondrial cytochrome c oxidase subunit 1 (cox1) gene. Table S4: The distribution of the detected haplotypes among the different study sites compared to previously published An. stephensi cytochrome c oxidase 1 (CO1) sequences.

Author Contributions

Conceptualization, M.A. and A.A.; methodology, M.A., N.S.M. and A.A.; software, N.S.M.; formal analysis, N.S.M. and A.A.; investigation, M.A., H.S., I.M., O.A., H.A. and N.S.M.; resources, M.A.; data curation, N.S.M.; visualization, M.A., N.S.M. and A.A.; supervision, A.A.; project administration, M.A.; funding acquisition, M.A.; writing—original draft preparation, M.A. and A.A.; writing—review and editing, M.A., H.S., I.M., O.A., H.A., N.S.M. and A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated during this study are included this article and the Supplementary Materials.

Acknowledgments

We are grateful for our colleague Lisa Crump at the Swiss TPH for her revision. Moreover, we are thankful for our colleagues at the National Public Health Laboratory, Tropical Medicine Research Institute, and Federal and State Ministries of Health for their support in the collection of samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. The World Health Organization. World Malaria Report 2020. Available online: https://cms.who.int/publications-detail-redirect/9789240015791 (accessed on 16 July 2021).
  2. The World Health Organization. World Malaria Report 2019. Available online: https://www.who.int/publications/i/item/world-malaria-report-2019 (accessed on 9 October 2019).
  3. Ageep, T.B.; Damiens, D.; Alsharif, B.; Ahmed, A.; Salih, E.H.; Ahmed, F.T.; Diabaté, A.; Lees, R.S.; Gilles, J.R.; El Sayed, B.B. Participation of Irradiated Anopheles Arabiensis Males in Swarms Following Field Release in Sudan. Malar. J. 2014, 13, 484. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Bhatt, S.; Weiss, D.J.; Cameron, E.; Bisanzio, D.; Mappin, B.; Dalrymple, U.; Battle, K.E.; Moyes, C.L.; Henry, A.; Eckhoff, P.A.; et al. The Effect of Malaria Control on Plasmodium Falciparum in Africa between 2000 and 2015. Nature 2015, 526, 207–211. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Mohamed, N.S.; Abdelbagi, H.; Osman, H.A.; Ahmed, A.E.; Yousif, A.M.; Edris, Y.B.; Osman, E.Y.; Elsadig, A.R.; Siddig, E.E.; Mustafa, M.; et al. A Snapshot of Plasmodium Falciparum Malaria Drug Resistance Markers in Sudan: A Pilot Study. BMC Res. Notes 2020, 13, 512. [Google Scholar] [CrossRef] [PubMed]
  6. Abuelmaali, S.; Ahmed, A.; AbdHamid, M.; Elaagip, A.; Basheer, M.; Weetman, D. Poster of the Impacts of the Agricultural Schemes on Insecticide Resistance. In Proceedings of the Roll Back Malaria (RBM) Vector Control Working Group (VCWG-11) Annual Meeting 2016, Geneva, Switzerland, 10 December 2013. [Google Scholar] [CrossRef]
  7. Mohamed, N.S.; AbdElbagi, H.; Elsadig, A.R.; Ahmed, A.E.; Mohammed, Y.O.; Elssir, L.T.; Elnour, M.-A.B.; Ali, Y.; Ali, M.S.; Altahir, O.; et al. Assessment of Genetic Diversity of Plasmodium Falciparum Circumsporozoite Protein in Sudan: The RTS, S Leading Malaria Vaccine Candidate. Malar. J. 2021, 20, 436. [Google Scholar] [CrossRef] [PubMed]
  8. The World Health Organization. World Malaria Report 2021. Available online: https://www.who.int/publications/i/item/9789240040496 (accessed on 8 December 2021).
  9. Mohamed, N.S.; Ali, Y.; Muneer, M.S.; Siddig, E.E.; Sibley, C.H.; Ahmed, A. Malaria Epidemic in Humanitarian Crisis Settings the Case of South Kordofan State, Sudan. J. Infect. Dev. Ctries 2021, 15, 168–171. [Google Scholar] [CrossRef]
  10. Abdelwhab, O.F.; Elaagip, A.; Albsheer, M.M.; Ahmed, A.; Paganotti, G.M.; Abdel Hamid, M.M. Molecular and Morphological Identification of Suspected Plasmodium Vivax Vectors in Central and Eastern Sudan. Malar. J. 2021, 20, 132. [Google Scholar] [CrossRef]
  11. Ahmed, A.; Eldigail, M.; Elduma, A.; Breima, T.; Dietrich, I.; Ali, Y.; Weaver, S.C. First Report of Epidemic Dengue Fever and Malaria Co-Infections among Internally Displaced Persons in Humanitarian Camps of North Darfur, Sudan. Int. J. Infect. Dis. 2021, 108, 513–516. [Google Scholar] [CrossRef] [PubMed]
  12. Ahmed, A.; Mohamed, N.S.; Siddig, E.E.; Algaily, T.; Sulaiman, S.; Ali, Y. The Impacts of Climate Change on Displaced Populations: A Call for Actions. J. Clim. Change Health 2021, 3, 100057. [Google Scholar] [CrossRef]
  13. Takken, W.; Lindsay, S. Increased Threat of Urban Malaria from Anopheles Stephensi Mosquitoes, Africa. Emerg. Infect. Dis. 2019, 25, 1431–1433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ahmed, A.; Khogali, R.; Elnour, M.-A.B.; Nakao, R.; Salim, B. Emergence of the Invasive Malaria Vector Anopheles Stephensi in Khartoum State, Central Sudan. Parasites Vectors 2021, 14, 511. [Google Scholar] [CrossRef]
  15. Ahmed, A.; Pignatelli, P.; Elaagip, A.; Hamid, M.M.A.; Alrahman, O.F.; Weetman, D. Invasive Malaria Vector Anopheles Stephensi Mosquitoes in Sudan, 2016–2018. EID 2021, 27, 2952–2954. [Google Scholar] [CrossRef] [PubMed]
  16. Faulde, M.K.; Rueda, L.M.; Khaireh, B.A. First Record of the Asian Malaria Vector Anopheles Stephensi and Its Possible Role in the Resurgence of Malaria in Djibouti, Horn of Africa. Acta Trop. 2014, 139, 39–43. [Google Scholar] [CrossRef] [PubMed]
  17. de Santi, V.P.; Khaireh, B.A.; Chiniard, T.; Pradines, B.; Taudon, N.; Larréché, S.; Mohamed, A.B.; de Laval, F.; Berger, F.; Gala, F.; et al. Role of Anopheles Stephensi Mosquitoes in Malaria Outbreak, Djibouti, 2019. Emerg. Infect. Dis. 2021, 27, 1697–1700. [Google Scholar] [CrossRef]
  18. WHO. Vector Alert: Anopheles Stephensi Invasion and Spread. Available online: https://www.who.int/news-room/detail/26-08-2019-vector-alert-anopheles-stephensi-invasion-and-spread (accessed on 31 May 2020).
  19. Coetzee, M. Key to the Females of Afrotropical Anopheles Mosquitoes (Diptera: Culicidae). Malar. J. 2020, 19, 70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Folmer, O.; Black, M.; Hoeh, W.; Lutz, R.; Vrijenhoek, R. DNA Primers for Amplification of Mitochondrial Cytochrome c Oxidase Subunit I from Diverse Metazoan Invertebrates. Mol. Marine Biol. Biotechnol. 1994, 3, 294–299. [Google Scholar]
  21. Manske, M. GENtle, a Free Multi-Purpose Molecular Biology Tool. Ph.D. Thesis, Universität zu Köln, Köln, Germany, 2006. [Google Scholar]
  22. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  23. Hao, Y.-J.; Zou, Y.-L.; Ding, Y.-R.; Xu, W.-Y.; Yan, Z.-T.; Li, X.-D.; Fu, W.-B.; Li, T.-J.; Chen, B. Complete Mitochondrial Genomes of Anopheles Stephensi and An. Dirus and Comparative Evolutionary Mitochondriomics of 50 Mosquitoes. Sci. Rep. 2017, 7, 7666. [Google Scholar] [CrossRef]
  24. Tamura, K. Estimation of the Number of Nucleotide Substitutions When There Are Strong Transition-Transversion and G+C-Content Biases. Mol. Biol. Evol. 1992, 9, 678–687. [Google Scholar] [CrossRef] [Green Version]
  25. Rozas, J.; Sánchez-DelBarrio, J.C.; Messeguer, X.; Rozas, R. DnaSP, DNA Polymorphism Analyses by the Coalescent and Other Methods. Bioinformatics 2003, 19, 2496–2497. [Google Scholar] [CrossRef]
  26. Carter, T.E.; Yared, S.; Gebresilassie, A.; Bonnell, V.; Damodaran, L.; Lopez, K.; Ibrahim, M.; Mohammed, S.; Janies, D. First Detection of Anopheles Stephensi Liston, 1901 (Diptera: Culicidae) in Ethiopia Using Molecular and Morphological Approaches. Acta Trop. 2018, 188, 180–186. [Google Scholar] [CrossRef]
  27. Iwamura, T.; Guzman-Holst, A.; Murray, K.A. Accelerating Invasion Potential of Disease Vector Aedes Aegypti under Climate Change. Nat. Commun. 2020, 11, 2130. [Google Scholar] [CrossRef]
  28. Crowl, T.A.; Crist, T.O.; Parmenter, R.R.; Belovsky, G.; Lugo, A.E. The Spread of Invasive Species and Infectious Disease as Drivers of Ecosystem Change. Front. Ecol. Environ. 2008, 6, 238–246. [Google Scholar] [CrossRef]
  29. Ahmed, A.; Dietrich, I.; LaBeaud, A.D.; Lindsay, S.W.; Musa, A.; Weaver, S.C. Risks and Challenges of Arboviral Diseases in Sudan: The Urgent Need for Actions. Viruses 2020, 12, 81. [Google Scholar] [CrossRef] [Green Version]
  30. Ahmed, A.; Ali, Y.; Mohamed, N.S. Arboviral Diseases: The Emergence of a Major yet Ignored Public Health Threat in Africa. Lancet Planet. Health 2020, 4, e555. [Google Scholar] [CrossRef]
  31. Ahmed, A.; Ali, Y.; Elmagboul, B.; Mohamed, O.; Elduma, A.; Bashab, H.; Mahamoud, A.; Khogali, H.; Elaagip, A.; Higazi, T. Dengue Fever in the Darfur Area, Western Sudan. Emerg. Infect. Dis. 2019, 25, 2126. [Google Scholar] [CrossRef]
  32. Ahmed, A.; Elduma, A.; Magboul, B.; Higazi, T.; Ali, Y. The First Outbreak of Dengue Fever in Greater Darfur, Western Sudan. Trop. Med. Infect. Dis. 2019, 4, 43. [Google Scholar] [CrossRef] [Green Version]
  33. Ahmed, A.; Ali, Y.; Elduma, A.; Eldigail, M.H.; Mhmoud, R.A.; Mohamed, N.S.; Ksiazek, T.G.; Dietrich, I.; Weaver, S.C. Unique Outbreak of Rift Valley Fever in Sudan, 2019. Emerg. Infect. Dis. 2020, 26, 3030–3033. [Google Scholar] [CrossRef] [PubMed]
  34. Ahmed, A.; Mahmoud, I.; Eldigail, M.; Elhassan, R.M.; Weaver, S.C. The Emergence of Rift Valley Fever in Gedaref State Urges the Need for a Cross-Border One Health Strategy and Enforcement of the International Health Regulations. Pathogens 2021, 10, 885. [Google Scholar] [CrossRef] [PubMed]
  35. Sinka, M.E.; Pironon, S.; Massey, N.C.; Longbottom, J.; Hemingway, J.; Moyes, C.L.; Willis, K.J. A New Malaria Vector in Africa: Predicting the Expansion Range of Anopheles Stephensi and Identifying the Urban Populations at Risk. Proc. Natl. Acad. Sci. USA 2020, 117, 24900–24908. [Google Scholar] [CrossRef]
  36. Hussien, H.H. Malaria’s Association with Climatic Variables and an Epidemic Early Warning System Using Historical Data from Gezira State, Sudan. Heliyon 2019, 5, e01375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Tian, H.; Li, N.; Li, Y.; Kraemer, M.U.G.; Tan, H.; Liu, Y.; Li, Y.; Wang, B.; Wu, P.; Cazelles, B.; et al. Malaria Elimination on Hainan Island despite Climate Change. Commun. Med. 2022, 2, 12. [Google Scholar] [CrossRef]
  38. Elsanousi, Y.E.A.; Elmahi, A.S.; Pereira, I.; Debacker, M. Impact of the 2013 Floods on the Incidence of Malaria in Almanagil Locality, Gezira State, Sudan. PLoS Curr. 2018, 10. [Google Scholar] [CrossRef]
  39. Franklinos, L.H.V.; Jones, K.E.; Redding, D.W.; Abubakar, I. The Effect of Global Change on Mosquito-Borne Disease. Lancet Infect. Dis. 2019, 19, e302–e312. [Google Scholar] [CrossRef]
  40. Colón-González, F.J.; Sewe, M.O.; Tompkins, A.M.; Sjödin, H.; Casallas, A.; Rocklöv, J.; Caminade, C.; Lowe, R. Projecting the Risk of Mosquito-Borne Diseases in a Warmer and More Populated World: A Multi-Model, Multi-Scenario Intercomparison Modelling Study. Lancet Planet. Health 2021, 5, e404–e414. [Google Scholar] [CrossRef]
  41. Watts, N.; Amann, M.; Arnell, N.; Ayeb-Karlsson, S.; Belesova, K.; Boykoff, M.; Byass, P.; Cai, W.; Campbell-Lendrum, D.; Capstick, S.; et al. The 2019 Report of The Lancet Countdown on Health and Climate Change: Ensuring That the Health of a Child Born Today Is Not Defined by a Changing Climate. Lancet 2019, 394, 1836–1878. [Google Scholar] [CrossRef] [Green Version]
  42. Ahmed, A. The Threat of the Invasive An. Stephensi in the African Region and Lessons for Vector Control Policy. In Proceedings of the ISNTD Bites 2021 Vector Control & Vector-borne Diseases, London, UK, 27 October 2021. [Google Scholar] [CrossRef]
  43. United Nations Office for the Coordination of Humanitarian Affairs in 2019 Malaria Breached the Epidemic Threshold in Sudan. Available online: https://reports.unocha.org/en/country/sudan/card/1Av57iny1Q/ (accessed on 17 January 2020).
  44. Gayan Dharmasiri, A.G.; Perera, A.Y.; Harishchandra, J.; Herath, H.; Aravindan, K.; Jayasooriya, H.T.R.; Ranawaka, G.R.; Hewavitharane, M. First Record of Anopheles Stephensi in Sri Lanka: A Potential Challenge for Prevention of Malaria Reintroduction. Malar. J. 2017, 16, 326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Ahmed, A. Urgent Call for a Global Enforcement of the Public Sharing of Health Emergencies Data: Lesson Learned from Serious Arboviral Disease Epidemics in Sudan. Int. Health 2020, 12, 238–240. [Google Scholar] [CrossRef] [PubMed]
  46. Seyfarth, M.; Khaireh, B.A.; Abdi, A.A.; Bouh, S.M.; Faulde, M.K. Five Years Following First Detection of Anopheles Stephensi (Diptera: Culicidae) in Djibouti, Horn of Africa: Populations Established-Malaria Emerging. Parasitol. Res. 2019, 118, 725–732. [Google Scholar] [CrossRef] [PubMed]
  47. Tadesse, F.G.; Ashine, T.; Teka, H.; Esayas, E.; Messenger, L.A.; Chali, W.; Meerstein-Kessel, L.; Walker, T.; Wolde Behaksra, S.; Lanke, K.; et al. Anopheles Stephensi Mosquitoes as Vectors of Plasmodium Vivax and Falciparum, Horn of Africa, 2019. Emerg. Infect. Dis. 2021, 27, 603–607. [Google Scholar] [CrossRef] [PubMed]
  48. Mehravaran, A.; Vatandoost, H.; Oshaghi, M.A.; Abai, M.R.; Edalat, H.; Javadian, E.; Mashayekhi, M.; Piazak, N.; Hanafi-Bojd, A.A. Ecology of Anopheles Stephensi in a Malarious Area, Southeast of Iran. Acta Med. Iran. 2012, 50, 61–65. [Google Scholar]
  49. Munga, S.; Minakawa, N.; Zhou, G.; Barrack, O.-O.J.; Githeko, A.K.; Yan, G. Oviposition Site Preference and Egg Hatchability of Anopheles Gambiae: Effects of Land Cover Types. J. Med. Entomol. 2005, 42, 993–997. [Google Scholar] [CrossRef] [Green Version]
  50. Yang, H.M.; Macoris, M.L.G.; Galvani, K.C.; Andrighetti, M.T.M.; Wanderley, D.M.V. Assessing the Effects of Temperature on the Population of Aedes Aegypti, the Vector of Dengue. Epidemiol. Infect. 2009, 137, 1188–1202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Yared, S.; Gebressielasie, A.; Damodaran, L.; Bonnell, V.; Lopez, K.; Janies, D.; Carter, T.E. Insecticide Resistance in Anopheles Stephensi in Somali Region, Eastern Ethiopia. Malar. J. 2020, 19, 180. [Google Scholar] [CrossRef] [PubMed]
  52. Ahmed, A.; Abubakr, M.; Ali, Y.; Siddig, E.E.; Mohamed, N.S. Vector Control Strategy for Anopheles Stephensi in Africa. Lancet Microbe 2022. [Google Scholar] [CrossRef]
  53. Ahmed, A. Current Status of Mosquito-Borne Arboviruses in Sudan, and Challenges of Surveillance and Responses. In Proceedings of the Mosquito-Borne Arboviruses: The Rising Global Threat, Malaria Consortium Webinar, Internet, 10 February 2021. [Google Scholar]
  54. Vaughan, J.A.; Turell, M.J. Facilitation of Rift Valley Fever Virus Transmission by Plasmodium Berghei Sporozoites in Anopheles Stephensi Mosquitoes. Am. J. Trop. Med. Hyg. 1996, 55, 407–409. [Google Scholar] [CrossRef]
  55. Yadav, P.; Gokhale, M.D.; Barde, P.V.; Singh, D.K.; Mishra, A.C.; Mourya, D.T. Experimental Transmission of Chikungunya Virus by Anopheles Stephensi Mosquitoes. Acta Virol. 2003, 47, 45–47. [Google Scholar] [PubMed]
  56. The World Health Organization Chikungunya Outbreak in Chad. Available online: https://www.who.int/emergencies/disease-outbreak-news/item/chikungunya-chad (accessed on 8 February 2022).
Biology 11 00409 g001 550
Figure 1. Sudanese States (highlighted in light blue) where active surveys were implemented to detect An. stephensi mosquitoes. The red circles indicate the number of An. stephensi mosquitoes collected.
Figure 1. Sudanese States (highlighted in light blue) where active surveys were implemented to detect An. stephensi mosquitoes. The red circles indicate the number of An. stephensi mosquitoes collected.
Biology 11 00409 g001
Biology 11 00409 g002 550
Figure 2. Representative breeding sites that were positive for presence of aquatic stages of An. stephensi (larvae and pupae). (A) Plastic cistern, (B) plastic barrel, (C) mud pot, (D) clay pot (Zeir/Jar), (E) ground water-basin, (F) plastic jerrycan, (G) iron barrel, (H,I) leakage of broken water supply, and (J) rainwater pond trapped in a rocky valley.
Figure 2. Representative breeding sites that were positive for presence of aquatic stages of An. stephensi (larvae and pupae). (A) Plastic cistern, (B) plastic barrel, (C) mud pot, (D) clay pot (Zeir/Jar), (E) ground water-basin, (F) plastic jerrycan, (G) iron barrel, (H,I) leakage of broken water supply, and (J) rainwater pond trapped in a rocky valley.
Biology 11 00409 g002
Biology 11 00409 g003 550
Figure 3. Sequence alignment of the 20 Sudanese An. stephensi haplotypes. Substitutions were indicted with their nucleotide codes; no deletion nor insertions were present. The dots (.) indicate identical nucleotides at the specified position in comparison with the reference sequence KT899888.1 (ref). * Nucleotide substitution positions, based on the start of the complete mitochondrial cytochrome c oxidase 1 (CO1) gene, are read vertically.
Figure 3. Sequence alignment of the 20 Sudanese An. stephensi haplotypes. Substitutions were indicted with their nucleotide codes; no deletion nor insertions were present. The dots (.) indicate identical nucleotides at the specified position in comparison with the reference sequence KT899888.1 (ref). * Nucleotide substitution positions, based on the start of the complete mitochondrial cytochrome c oxidase 1 (CO1) gene, are read vertically.
Biology 11 00409 g003
Biology 11 00409 g004 550
Figure 4. Phylogenetic tree showing the relationship between the Sudanese An. stephensi haplotypes with 15 reference sequences. Sudanese haplotypes (Sudan Hap01–Sudan Hap20) are in bold. The reference sequences along with their accession numbers and origin of isolate were included for each. Drosophila melanogaster was used as an outgroup taxon.
Figure 4. Phylogenetic tree showing the relationship between the Sudanese An. stephensi haplotypes with 15 reference sequences. Sudanese haplotypes (Sudan Hap01–Sudan Hap20) are in bold. The reference sequences along with their accession numbers and origin of isolate were included for each. Drosophila melanogaster was used as an outgroup taxon.
Biology 11 00409 g004
Biology 11 00409 g005 550
Figure 5. Median-joining haplotype network of the 19 Sudanese An. stephensi haplotypes and worldwide An. stephensi sequences. Haplotypes of each region are presented in color code. Black dashes between the haplotype lines represent the number of substitutions.
Figure 5. Median-joining haplotype network of the 19 Sudanese An. stephensi haplotypes and worldwide An. stephensi sequences. Haplotypes of each region are presented in color code. Black dashes between the haplotype lines represent the number of substitutions.
Biology 11 00409 g005
Biology 11 00409 g006 550
Figure 6. Pie charts inserted to show the haplotypic composition of An. stephensi population across the states of Sudan. Haplotypes of each region are presented in color code.
Figure 6. Pie charts inserted to show the haplotypic composition of An. stephensi population across the states of Sudan. Haplotypes of each region are presented in color code.
Biology 11 00409 g006
Table
Table 1. Identity and percentage of Anopheles mosquito samples collected by conducting active entomological surveys from different States of Sudan in 2021.
Table 1. Identity and percentage of Anopheles mosquito samples collected by conducting active entomological surveys from different States of Sudan in 2021.
No.StateAn. stephensiAn. arabiensisAn. macupalpisAn. rufipes% of An. stephensiTotal
1Blue Nile012400-124
2Gedarif163100094.2%173
3Kassala106800057.1%186
4North Darfur62196002.1%285
5Red Sea309000100%309
6River Nile11254004.2%265
7Sinnar2923202510.1%286
Total6309216025 1628
Table
Table 2. Diversity and neutrality indices for An. stephensi populations calculated from the nucleotide data set of the mitochondrial cytochrome c oxidase 1 (CO1) gene.
Table 2. Diversity and neutrality indices for An. stephensi populations calculated from the nucleotide data set of the mitochondrial cytochrome c oxidase 1 (CO1) gene.
PopulationNSHapHapd ± VarHapdPiTajimaDFuFs
North Darfur6861.0 ± 0.009260.009−0.6231−3.178
River Nile11660.873 ± 0.004990.006170.1874−1.375
Khartoum5440.9 ± 0.025920.005750.2734−1.012
Red Sea18711110.539 ± 0.001690.00268−1.2208−4.434
Kassala9411110.579 ± 0.003390.00357−1.1133−3.962
El Gedarif9611110.64 ± 0.002790.00523−0.396−1.963
Al Gezira1n.d.1n.d.n.d.n.d.n.d.
India27230.274 ± 0.011430.00082−0.9543−1.052
Iran201n.d.0n.d.n.d.
Pakistan29340.2 ± 0.009550.00059−1.7326−3.324 *
Saudi Arabia8340.643 ± 0.03390.00216−1.4475−1.832
Sri Lanka18120.471 ± 0.006780.001351.16621.215
United Arab Emirates1n.d.1n.d.n.d.n.d.n.d.
Djibouti1n.d.1n.d.n.d.n.d.n.d.
Ethiopia1n.d.1n.d.n.d.n.d.n.d.
Sudan ¥3120.667 ± 0.098770.00192n.d.0.201
N: Number of sequences; S: number of segregating sites; Hap: number of haplotypes; Hapd ± VarHapd: haplotype diversity ± variance of haplotype diversity; Pi: nucleotides diversity; n.d.: not determined; ¥: indicates previously published sequences from Sudan; *: significance level < 0.05.
Table
Table 3. Pairwise fixation index (Fst test values) between An. stephensi populations calculated from the nucleotide data set of CO1 gene.
Table 3. Pairwise fixation index (Fst test values) between An. stephensi populations calculated from the nucleotide data set of CO1 gene.
Populations *North DarfurRiver NileKhartoumRed SeaKassalaGedarifIndiaIranPakistanSaudi ArabiaSri Lanka
River Nile0.058----------
Khartoum0.0550.036---------
Red Sea0.4230.2030.415--------
Kassala0.3940.1710.3720.004-------
El Gedarif0.3600.1550.3240.0540.014------
India0.4840.2820.5130.0300.0580.136-----
Iran0.6240.5790.6880.6780.6210.5450.876----
Pakistan0.4840.2730.5090.0480.0640.1350.0230.906---
Saudi Arabia0.6520.6720.6820.7960.7450.6410.8790.9270.885--
Sri Lanka0.4850.3040.5120.0650.1060.1680.1260.8240.2240.863-
Sudan0.4000.1250.2590.2260.1700.1560.4160.8000.3970.8050.431
* Populations consisting of one sequence were not included in the comparison.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abubakr, M.; Sami, H.; Mahdi, I.; Altahir, O.; Abdelbagi, H.; Mohamed, N.S.; Ahmed, A. The Phylodynamic and Spread of the Invasive Asian Malaria Vectors, Anopheles stephensi, in Sudan. Biology 2022, 11, 409. https://doi.org/10.3390/biology11030409

AMA Style

Abubakr M, Sami H, Mahdi I, Altahir O, Abdelbagi H, Mohamed NS, Ahmed A. The Phylodynamic and Spread of the Invasive Asian Malaria Vectors, Anopheles stephensi, in Sudan. Biology. 2022; 11(3):409. https://doi.org/10.3390/biology11030409

Chicago/Turabian Style

Abubakr, Mustafa, Hamza Sami, Isam Mahdi, Omnia Altahir, Hanadi Abdelbagi, Nouh Saad Mohamed, and Ayman Ahmed. 2022. "The Phylodynamic and Spread of the Invasive Asian Malaria Vectors, Anopheles stephensi, in Sudan" Biology 11, no. 3: 409. https://doi.org/10.3390/biology11030409

APA Style

Abubakr, M., Sami, H., Mahdi, I., Altahir, O., Abdelbagi, H., Mohamed, N. S., & Ahmed, A. (2022). The Phylodynamic and Spread of the Invasive Asian Malaria Vectors, Anopheles stephensi, in Sudan. Biology, 11(3), 409. https://doi.org/10.3390/biology11030409

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop