Skip to main content

Aedes species (Diptera: Culicidae) ecological and host feeding patterns in the north-eastern parts of South Africa, 2014–2018

Abstract

Background

There is a paucity of recent data and knowledge on mosquito diversity and potential vectors of arboviruses in South Africa, with most of the available data dating back to the 1950s–1970s. Aedes and Culex species are the major vectors of some of the principal arboviruses which have emerged and re-emerged in the past few decades.

Methods

In this study we used entomological surveillance in selected areas in the north-eastern parts of South Africa from 2014 to 2018 to assess mosquito diversity, with special emphasis on the Aedes species. The impact of trap types and environmental conditions was also investigated. Identification of the blood meal sources of engorged females collected during the study period was carried out, and DNA barcodes were generated for selected species.

Results

Overall, 18.5% of the total Culicidae mosquitoes collected belonged to the genus Aedes, with 14 species recognised or suspected vectors of arboviruses. Species belonging to the Neomelaniconion subgenus were commonly collected in the Bushveld savanna at conservation areas, especially Aedes mcintoshi and Aedes circumluteolus. Aedes aegypti was present in all sites, albeit in low numbers. Temperature was a limiting factor for the Aedes population, and they were almost exclusively collected at temperatures between 18 °C and 27 °C. The cytochrome oxidase subunit I (COI) barcode fragment was amplified for 21 Aedes species, and for nine of these species it was the first sequence information uploaded on GenBank.

Conclusion

This study provides a better understanding of the diversity and relative abundance of Aedes species in the north-east of South Africa. The information provided here will contribute to future arboviral research and implementation of efficient vector control and prevention strategies.

Graphical abstract

Background

Multiple mosquito species transmit pathogens of medical and veterinary importance to humans and animals, and mosquitoes are therefore considered one of the most important arthropod groups in the public health field [1, 2].

Vector-borne diseases are being reported with greater frequency due to the global movement of humans, animals, and goods, in combination with climate change [3] and the impact of land use and urbanisation [4]; this is especially evident with respect to Aedes-related arboviruses. Information such as the distribution, abundance, and seasonality of vectors, combined with an understanding of their relationship with the environment, is required to develop and implement successful vector control programs [5].

Knowledge of vector host preferences is important for understanding arboviral circulation between vectors, animals, and humans. The feeding behaviour is affected by intrinsic and extrinsic factors [6], and host selection by mosquitoes can be opportunistic and affected by the presence and/or abundance of the various vertebrate species [7].

Another important step for effective monitoring of vectors is the correct identification of the species. Morphological identification uses published keys and relies on external features of adult mosquitoes. For that, the specimens must be handled and stored carefully, a process which is time-consuming and requires specific expertise. Moreover, it does not identify genetic variations and phenotypic plasticity which can be manifested at the species level and in vector competence [8]. DNA barcoding is becoming more popular as an efficient methodology to complement morphological taxonomy. However, available studies using DNA barcoding in southern Africa are mostly focused on Anophelinae mosquitoes [9, 10].

Historically, a large number of studies were published in South Africa during the 1950s and 1970s associated with outbreaks in livestock and in regional surveys on arboviruses affecting humans [11,12,13,14]. Recently, Cornel et al. [15] investigated the diversity and abundance of mosquitoes in southern Africa, and Gorsich et al. [16] compared different trapping methodologies to assess adult mosquito populations in South Africa. However, there is a paucity of recent data about the diversity, density, seasonality, biology, and molecular identification of mosquito species, especially for Culicinae mosquitoes, and Aedes species in particular. The landscape of most of southern Africa has been altered as a result of human activity. Recent data are helpful for understanding how mosquito diversity has changed since the 1950s–1970s and how the risk of arboviral circulation differs between humans and animals.

The aim of this study was to assess the broad patterns of Aedes mosquito species diversity and abundance in different habitats across selected sites in the north-eastern parts of South Africa. Additionally, the influence of climatic features in the faunistic composition of Aedes species within each region was analysed and the relative efficiency of different trap types was compared. The blood meals from engorged females collected during the study period were also determined using molecular methods, and DNA barcodes were sequenced to identify and characterise Aedes species in South Africa. The collected data will assist with risk assessment for disease outbreaks in the region.

Methods

Study areas

The study was conducted in five provinces of South Africa, namely, Gauteng, Limpopo, North West, Mpumalanga, and KwaZulu-Natal (Fig. 1). Mosquitoes were sampled monthly from January 2014 to May 2017 from sentinel sites and from supplementary collections in ad hoc sites following the detection of arboviral occurrence in animal or human hosts. Additional collections were performed from March to April 2017 in and around Kruger National Park (KNP) [16, 17], and mosquitoes caught in this sampling were used only for molecular identification and blood meal analyses. In 2018 only one collection per site was performed from January to May. Conservation areas were sampled in order to collect sylvatic species in those areas where there is less interference from humans and livestock. Sylvatic species can also be important as arboviral vectors among wildlife species. Peri-urban areas were sampled because they are located between the urban and usually more rural farmland areas or between the urban and conservation areas. The ecological features associated with the sites sampled are Highveld (regions vary between 1200 and 1800 m in elevation), Middleveld Bushveld (varying between 600 and 1200 m), and Lowveld Bushveld (varying between 150 and 600 m). A summary of the collections per site are given in Table 1.

Fig. 1
figure 1

Mosquito collection sites and SAWS stations, South Africa, surveyed from January 2014 to May 2018. Squares indicate sentinel sites, triangles indicate ad hoc sites, and stars indicate SAWS stations

Table 1 Number of nights, trapping events, traps set, and mosquitoes collected per genera in the sentinel sites from January 2014 to June 2018

Sampling was performed in different land use types: urban, peri-urban (horse farms), rural, and conservation areas (Additional file 1: Table S1). Trapping was carried out from 15:30–16:00 to 5:00–8:00, and sampling was conducted for 1–3 consecutive nights per site using multiple types of carbon dioxide (CO2)-baited traps: mosquito net, CDC miniature light traps (Centers for Disease Control and Prevention, USA), and BG-Sentinel traps (Biogents AG, Regensburg, Germany). Mosquitoes were collected from the net traps using aspirators and transferred to mesh-topped polystyrene cups. CDC miniature light traps were hung at least 1.5 m from the ground, baited with CO2. BG-Sentinel traps were added from 2017 and were additionally baited with a non-toxic lure. All the traps were placed at least 80 m apart to reduce trap interference.

Collected mosquitoes were immediately euthanised by freezing and were morphologically identified to species using published keys and descriptions [18,19,20,21]. Engorged females were separated by species, collection site, and date, and individually preserved at −80 °C until further analysis. Representative specimens were pinned as reference material, and 1–3 legs were individually preserved for further molecular work. The classification of Wilkerson et al. [22] for Aedini mosquitoes was adopted. Morphologically similar or indistinct species were recorded as belonging to a group as follows: Aedes dentatus, Aedes leesoni, Aedes dentatus/leesoni, or Aedes vexans. Aedes specimens that were too badly damaged to identify morphologically to species were identified to the genus or in a species group following the groups mentioned above. Some adult female species were not distinguishable morphologically and these were recorded as either of the two possible species they could be. The pinned reference material was also compared with the specimens deposited at the National Institute for Communicable Diseases of South Africa (NICD).

Climatic conditions

Climate data were obtained from 11 permanent stations of the South African Weather Service (SAWS) closest to the sampling sites [23] (Fig. 1). Data obtained from SAWS included total hourly rainfall, hourly average wind speed, hourly average air temperature, and hourly average humidity. For this analysis, calculations for the climatic variables were as follows: air temperature (°C) (average daily, average nightly from 16:00 to 7:00, average prior 48 h, average prior 15 days, average prior 30 days), rainfall (mm) (total daily, total nightly from 16:00 to 7:00, total prior 48 h, total prior 15 days, total prior 30 days), wind speed (m/s) (average daily, average nightly from 16:00 to 7:00, average prior 48 h), and humidity (%) (average daily, average nightly from 16:00 to 7:00, average prior 48 h, average prior 15 days, average prior 30 days).

DNA extraction, DNA barcoding, and blood meal analyses

DNA was extracted using the DNeasy® Blood & Tissue Kit (Qiagen, Valencia, CA, USA), according to the manufacturer’s guidelines. DNA from the leg (1 up to 3 legs) of each mosquito and DNA from freshly engorged females was extracted.

The barcode region of mtDNA of subunit I of the cytochrome oxidase (COI) gene was amplified using universal primers [24]. The 50 µl polymerase chain reaction (PCR) consisted of 5 μl of the extracted DNA, 1 μl of 10 mM deoxynucleotide triphosphates (dNTPs), 10 μl of buffer, 0.5 μl Phusion® High-Fidelity DNA Polymerase (Thermo Fisher Scientific™), and 1 μl of 20 μM of each primer. PCR conditions were as follow: 98 °C for 30 s followed by 35 cycles of 98 °C for 30 s, 52 °C for 45 s, and 72 °C for 30 s, with a final extension of 72 °C for 5 min using a thermocycler (Applied Biosystems™).

Molecular identification of the blood meal was performed targeting cytochrome b (CYTB). Published primers targeting mammal [16, 17] and avian [16, 18] genetic markers were selected for amplification, and combined with Phusion® High-Fidelity DNA Polymerase (Thermo Fisher Scientific™). The PCR was performed in a final volume of 50 µl which consisted of 20 µM of each primer, 10 mM dNTPs, and 10 µl DNA. The reaction mix was then subjected to an initial incubation of 95 °C for 5 min followed by 36 cycles of 94 °C for 30 s, 55 °C (mammal) or 60 °C (avian) for 45 s, and 72 °C for 90 s, with a final extension of 72 °C for 7 min using a thermocycler (Applied Biosystems™).

Gel electrophoresis and sequencing

PCR products were viewed by 2.0% agarose gel electrophoresis containing ethidium bromide. Amplicons of the correct size were excised from the gel and purified using a Zymoclean™ Gel DNA Recovery Kit (Zymo Research, CA, USA) according to the manufacturer’s instructions. Purified amplicons were bidirectionally sequenced using the BigDye® Direct Cycle Sequencing Kit (Thermo Fisher Scientific™, MA, USA) and sent to the University of Pretoria DNA sequencing facility or Inqaba Biotec™ (Pretoria, South Africa) for Sanger sequencing.

Data analysis and phylogenetic analysis

Analyses and interpretation were performed separately for sentinel collections and ad hoc sites due to the variation in the number and method of collection events. Individual site adult mosquito densities were calculated as mean mosquitoes per trap-night by dividing the total collection at each site by the number of collection nights and number of traps utilised. There were 21 trap failure events as a result of battery failure or equipment damage from environmental factors or animals, and these collections were excluded prior to analysis.

Species richness and species diversity were calculated for each site and for each trap type. Species richness is reported as the number of mosquito species collected at each site. Species diversity was estimated by calculating the Simpson diversity index. The presumption for normality was tested using the Shapiro–Wilk test. To investigate whether there were significant differences between Aedes species composition within trap types, we applied the Wilcoxon rank-sum test (n = 2, two independent variables) or Kruskal–Wallis H tests (n > 2, three or more independent variables). Significant values were adjusted by the Bonferroni correction for multiple tests. The trap analyses were performed using data from collections at sentinel and ad hoc sites in 2017 and 2018. CDC light traps were considered the same trap type independently of the type of light used. The effect of environmental factors was tested using Spearman correlation analysis to examine the relation between the different variables and log10 abundance of Aedes mosquitoes per trap-night. Regression analysis was also performed to clarify the relationship and R-squared values were calculated. Differences were considered statistically significant at a value of p ≤ 0.05. Species richness and species diversity were calculated using the ‘entropart’ [25] package in R version 3.6.1 [26]. All statistical analyses were completed using IBM SPSS version 26 software (IBM Corp., Armonk, NY, USA). Graphics were generated using Microsoft Excel (version 2010) [27].

The resulting sequences were edited and analysed using the QIAGEN CLC Main Workbench version 8.0.1 [28]. Sequences were compared with the databases available in GenBank and the Barcode of Life Data System (BOLD) accordingly with the gene fragment, using the Basic Local Alignment Search Tool (BLAST). Multiple sequence alignments were compiled using the online version of Multiple Alignment using Fast Fourier Transform (MAFFT) with default parameters. The model of best fit was identified, and a maximum likelihood tree was produced using Molecular Evolutionary Genetics Analysis (MEGA) version 7.0 [29]. For barcode analyses, sequence statistics and pairwise sequence divergence were calculated based on the Kimura 2-parameter (K2P) model.

Results

Mosquito collection

A total of 64,603 adult mosquitoes belonging to 11 genera were collected across sentinel and ad hoc sites. The most common genus collected was Culex (38.90%, N = 25,131), followed by Anopheles (33.27%, N = 21,494), Aedes (18.63%, N = 12,037), Mansonia (6.17%, N = 3987), and other genera combined (3.03%, N = 1954, Uranotaenia, Aedeomyia, Ficalbia, Coquillettidia, Mimomyia, Culiseta, and Eretmapodites). Culex, Aedes, Anopheles, and Mansonia were the most important genera collected with respect to medical and veterinary significance. A total of 12,037 adult mosquitoes belonging to the Aedes genus were collected in all the sites from January 2014 to May 2018 (Table 1). At the sentinel sites, a total of 6702 aedine species were collected (10.37% of the total Culicidae collected), while a total of 5335 (8.26% of the total) were collected at the ad hoc sites.

While the overall number of Aedes mosquitoes collected in Kyalami and Boschkop was relatively low, these two sites had the highest biological diversity index. On average, species diversity was highest in Kyalami (Simpson’s Diversity Index 1− D = 0.88), followed by Boschkop (0.78), Lapalala (0.74), and Marakele (0.30) (Additional file 1: Table S2). Lapalala had the highest Aedes species richness, but it was not consistently the most biologically diverse site. This is most likely because of Aedes (Neomelaniconion) mcintoshi (42.62%) comprising a large proportion of the collection from this site, thus skewing the diversity index due to the abundance of one species. A similar result was observed in Marakele: even though the species richness was high, a large proportion of the collection were from Ae. mcintoshi (82.81%), making the biological diversity index exceptionally low (0.31).

For the ad hoc sites, Mnisi (25 species from 937 Aedes mosquitoes identified) and Jozini (17 species from 2348 Aedes mosquitoes identified) had the highest Aedes species richness. Aedes mcintoshi (45.99%, n = 431) and Aedes (Stegomyia) aegypti (22.30%, n = 209) were the most common species in Mnisi. Jozini exhibited a dominant species of Aedes (Aedimorphus) durbanensis (66.26%, n = 1556). The species diversity indexes were highest in KNP Shingwedzi (Simpson’s Diversity Index = 0.82), followed by Mnisi (0.70), KNP Skukuza (0.65), Roodeplaat (0.64), Vulpro (0.62), Benoni (0.52), Jozini (0.51), Matikwane (0.32), and Pretoria North (0.07) (Additional file 1: Table S3).

In this survey, 14 Aedes species which are recognised or suspected vectors of mosquito-borne viruses in southern Africa were collected. The most abundant potential vector collected throughout the study was the floodwater species Ae. mcintoshi, which occurred in high numbers in Middleveld locations such as Marakele and Lapalala (Limpopo Province at conservation sites) (Additional file 1: Fig. S1). Aedes fowleri and Ae. vittatus also appeared in high numbers specifically in Lapalala. Aedes aegypti occurred in all sentinel sites but in low numbers. Species from the dentatus group including Ae. dentatus, Ae. cumminsii, and Ae. pachyurus, the leesoni group, Ae. juppi, and Ae. unidentatus were common in the Highveld (Gauteng province at peri-urban sites) (Additional file 1: Fig. S1).

Effect of trap type on Aedes spp. abundance

There was a significant difference between the abundance of Aedes mosquitoes captured by different traps in sentinel sites (Kruskal–Wallis test: p-value < 0.05, Additional file 1: Table S4). In Marakele and Lapalala there was a significant difference in aedine captured between BG-Sentinel traps and net traps (Wilcoxon rank-sum test, V = 18.5—Marakele and V = 0—Lapalala, p < 0.05), and in Boschkop we found a significant difference between BG-Sentinel trap and CDC light traps (CDC-LT) (Wilcoxon rank-sum test, V = 10, p < 0.05). Net traps exhibited the highest Aedes species richness, with 24 species collected, followed by CDC-LT with 22 and BG-Sentinel with 20 aedine identified. There were two species that were collected exclusively in BG-Sentinel traps (Ae. vexans group and Aedes breedensis), one species only found in net traps (Aedes mixtus/microstictus), and one species found only in CDC-LT (Aedes aerarius). The Aedes species diversity index was highest at the CDC-LT (Simpson’s Diversity Index = 0.83), followed by BG-Sentinel (0.75) and net traps (0.56).

The ad hoc site collections yielded no significant “between-trap” difference in Aedes species captured (Kruskal–Wallis test: p = 0.376). The two sites with the largest Aedes collection, Mnisi and Jozini, showed a significant difference between BG-Sentinel and CDC-LT (Wilcoxon rank-sum test, V = 12, p < 0.05), and BG-Sentinel and net traps (Wilcoxon rank-sum test, V = 36, p < 0.05) at Mnisi, and a significant difference between BG-Sentinel and CDC traps (the only trap type used at this site) (Wilcoxon rank-sum test, V = 17, p < 0.05) in Jozini.

Overall, looking at the most abundant recognised or potential vectors collected, Ae. aegypti was predominantly collected in BG-Sentinel traps (56.87%, n = 484) and Ae. mcintoshi in net traps (93.01%, n = 1278).

Climatic conditions

The population structure and size of Aedes species fluctuated with season during the study period (Fig. 2). An overall trend was evident in the mean number of Aedes mosquitoes per trap-night, which was lower in the drier months (June to October) than during the wetter months (November to mid-April). The population peaks appeared to correlate with rainfall events and the highest mean temperatures. Aedes mcintoshi, the most abundant aedine within the study, exhibited higher population numbers during the rainy season. The population peaked in January to March 2014, and during the same period in 2017, coinciding with the highest rainfall (Additional file 1: Fig. S2). Aedes (Fredwardsius) vittatus appeared shortly after the first episodes of rain each year, and their population increased from November to January (Additional file 1: Fig. S2). Interestingly, Aedes fowleri mosquitoes were only collected in the early wet season (November and December) of 2016 and not in the other years.

Fig. 2
figure 2

Aedes mosquito abundance per trap-night from January 2014 to June 2018 at the sentinel sites

A positive correlation was observed between the abundance of Aedes per trap-night and the environmental conditions including average temperature 30 days prior to collection, average rainfall 30 days prior to the collection event, and average humidity 15 days prior to collection, with Spearman rank correlations of 0.533 (p < 0.001), 0.498 (p < 0.001), and 0.388 (p < 0.001), respectively (Additional file 1: Table S5). There was no significant correlation with average daily rainfall. The effect of wind on breeding, feeding, and flight was also considered as an important factor, particularly during the collection event. However, there was a negative correlation between the abundance of Aedes per trap-night and wind variables tested (Additional file 1: Table S3). The effect of air temperature on the abundance of aedine mosquitoes showed a relationship with aedine abundance during temperatures between 18 °C and 27 °C. The regression relationship (R2) between the abundance of Aedes per trap-night and average temperature 30 and 15 days prior to collection, average rainfall 30 and 15 days prior to the collection event, and average humidity 30 and 15 prior to collection was quite low, ranging between 0.031 and 0.157, respectively, which is shown by a scatter plot (Additional file 1: Fig. S3).

Molecular identification

A total of a 52 COI sequences were generated in this study from 21 Aedes species belonging to nine subgenera: Stegomyia (4 Ae. aegypti, 2 Ae. simpsoni, 1 Ae. ledgeri), Aedimorphus (1 Ae. hirsutus, 2 Ae. vexans, 3 Ae. fowleri, 4 Ae. cumminsii, 1 Ae. pachyurus, 5 Ae. eritreae, 3 Ae. durbanensis, 1 Ae. quasiunivittatus), Catageiomyia (5 Ae. microstictus), Neomelaniconion (2 Ae. aurovenatus, 4 Ae. mcintoshi, 1 Ae. circumluteolus, 3 Ae. unidentatus), Fredwardsius (2 Ae. vittatus), Ochlerotatus (3 Ae. juppi), Mucidus (2 Ae. sudanensis), Albuginosus (2 Ae. haworthi), and Diceromyia (1 Ae. furcifer).

DNA barcode sequences consisting of 507 bp produced in this study and sequences retrieved from NCBI GenBank and BOLD reported from Africa and worldwide were aligned and used to build a maximum likelihood tree (Fig. 3) as confirmation of morphological mosquito identification. The results show the evolutionary distances using the Tamura-Nei model with a bootstrap tested with 1000 replicates (30). The phylogenetic tree of COI showed Lutzomyia longipalpis species separated as outgroup from the aedine tested. The rest of the aedine taxa were divided into two clusters: one was a large cluster, and the other was a well-supported cluster containing only isolates belonging to the subgenus Mucidus (Ae. sudanensis and Aedes alternans). The larger cluster was further divided into a small cluster comprising Ochlerotatus species and a larger cluster of the remaining subgenera. Species belonging to Aedimorphus subgenera were recovered in different clusters, suggesting that this subgenus is not monophyletic in this study. Sequences produced here from mosquitoes which were identified morphologically as Ae. cumminsii (Additional file 1: Fig. S4) clustered with Ae. pachyurus (of the same group) and Ae. quasiunivittatus (of another group within the same subgenus) and did not cluster together with sequences from Ae. cumminsii from Kenya (KU187000.1, KU187001.1, and MK300225), Guinea (MN552300.1), or Senegal (MG242484). Species belonging to the Neomelaniconion subgenus were divided into two clusters: a small one containing only Ae. aurovenatus, and a larger one with no separation between Ae. unidentatus, Ae. mcintoshi, and Ae. circumluteolus.

Fig. 3
figure 3

Phylogenetic tree of sequences from Aedes mosquitoes based on 111 sequences and 507 base pairs of COI gene. The tree was constructed with MEGA 7, using the maximum likelihood method and the Tamura-Nei model with 1000 bootstrap replicates. The tree with the highest log likelihood (−6035.72) is shown. GenBank accession numbers are indicated. Numbers on internal branches indicate bootstrap values. Samples which are part of this study are marked with a black triangle

The sequences analysed showed a nucleotide diversity of 0.097. Total nucleotide composition of the COI fragment varied slightly across the tested samples, and the sequences were highly adenine/thymine (AT)-rich, ranging from 66.27 to 71.40%, while guanine/cytosine (GC) content ranged from 28.60 to 33.73%. The distance matrix calculated for the COI sequence for all the samples showed an overall K2P average distance of 0.105. Data analysis for sequence divergence revealed intraspecies divergence ranging from 0.0 to 2.83%; the higher value was used as a cut-off to define the limit for species identification, and interspecies divergence ranged from 3.24 to 16.27%. The highest genetic distance was for Ae. sudanensis collected in KNP Shingwedzi (from 9.91 to 15.29%) to the other samples. The K2P interspecific distance between the members of Neomelaniconion subgenera ranged from 0.20 to 9.92%, with the highest genetic distance between Aedes aurovenatus (from 8.12 to 9.92%) and the other members of Neomelaniconion species. The K2P interspecific distance ranged from 0.20 to 1.20% between Ae. mcintoshi, Ae. circumluteolus, and Ae. unidentatus.

Blood meal analyses

Engorged Aedes females comprised < 1% of the adults collected. That was expected due to characteristics of the trap methodology used for collection, which is more attractive to host-seeking adult females and not to recently blood-fed individuals. Despite this trap limitation, 112 freshly blood fed-females were caught between January 2014 and May 2018. The specimens tested represented 16 Aedes species. Sampling conducted in Kruger National Park (Shingwedzi area) and Marakele National Park had the greatest number of engorged females (n = 19/112 per site, 16.96%), followed by Lapalala (n = 18/112, 16.07%). Of the individuals tested, 64/112 (54.46%) were successfully identified for blood meal origin, which showed 62/64 (96.88%) mammal and 2/64 (3.12%) avian host species (Table 2).

Table 2 Aedes species, collection site, and host identified through sequencing of partial cytochrome b gene

The results revealed 16 different mammal species, of which buffalo (Syncerus caffer) was the predominant species detected (17%), followed by domestic cow (Bos taurus, 14%) and humans (Homo sapiens, 13%). Two different avian species, Gyps coprotheres (vulture) and Dryocopus gambensis (bushshrike), were detected. In species with multiple engorged females collected and successfully tested, host blood meal differed by up to five different hosts across the Aedes species (Additional file 1: Table S6).

Discussion

To assess aedine mosquito community diversity in the north-eastern parts of South Africa, the present study collected adult mosquitoes based on adult trapping methods commonly used in southern Africa [15, 16, 31].

Aedes species were overall found in lower numbers compared to Culex, Anopheles, and other genera. In total, across all trap types and collection sites, 30 species belonging to the genus Aedes were collected. However, aedine species are considered generally day-biting mosquitoes, and as this study used traps to collect mosquitoes from sunset to sunrise, this would likely have contributed to bias against the capture of Aedes species. Larval sampling is an additional method that could further complement adult collection in future studies to understand the diversity of an Aedes community [32] but was not used in this study. The trap types used here were also known to be biased towards host-seeking females due to the use of CO2, which could justify the low numbers of relative freshly engorged females collected during the study period. A possible method that could be incorporated in future studies to increase the number of freshly blood-fed females collected would be to sample mosquitoes from resting sites [33, 34]. However, despite the low number of engorged females tested, a better understanding of host selection by Aedes species was possible. All examined specimens primarily or exclusively fed on mammals, which confirms prior literature on Aedes species’ host selection [20, 35, 36].

The comparison of the host preference suggests that Aedes are not specialist but rather more opportunistic mammophilic blood feeders. The specimens identified feeding on cattle were collected in rural (Mnisi and Jozini) or peri-urban (Benoni) areas. At rural sites, cattle were the primary livestock present, while in Benoni, a peri-urban area of small holdings outside of Johannesburg and Pretoria, there are many horses but also other animals such as cattle and domestic dogs. Consequently, it is not surprising that all the species found to feed on cattle were collected from these sites. Species that were identified feeding on goats and sheep were collected in sites localised in peri-urban areas (Boschkop and Kyalami) on horse farms where no cattle were present, although goats and sheep were commonly found. The Aedes species that were detected feeding on wildlife were mainly collected in conservation areas (Marakele, Lapalala, and KNP).

Aedes mcintoshi/circumluteolus was the most abundant aedine sampled in the conservation regions. They are recognised as floodwater species that feed on livestock such as sheep, goats, and cattle, but also feed on humans [20]. Here, Ae. mcintoshi exhibited the highest feeding diversity and was identified feeding on cattle, hippopotamus, white rhinoceros, humans, kudu, bushbuck, blue wildebeest, and vultures, while Ae. circumluteolus was identified feeding on cattle. The variety of blood-feeding hosts can be attributed to the composition and abundance/availability of the hosts in each collection site. Additionally, these species are known as potential arboviral vectors, particularly in the Bushveld savanna (Marakele and Lapalala), and they were abundantly collected in the net traps. The species are widely distributed in South Africa and have been collected in the Karoo, Highveld, Lowveld, and the coast of Kwazulu-Natal Province [11, 37, 38]. They play a major role in Rift Valley fever (RVF) occurrence in South Africa, principally in the Karoo and Highveld (1950s and 1970s), during years that experienced much higher than average rainfall [51]. In such conditions, Ae. mcintoshi appeared to be the most important vector in the inland regions, while Ae. circumluteolus were more infected with RVF virus in KwaZulu-Natal Province [39]. More recently, in 2010–2011, an epidemic occurred with more than 14,000 animal cases recorded in eight out of nine provinces in South Africa [40], and in 2018, RVF was detected in a single farm in the Free State Province [41]. Furthermore, beyond their role in the maintenance of RVF, other arboviruses have been isolated from wild-caught mosquitoes in South Africa, including Wesselsbron virus (WSLB) and Middelburg virus (MIDV) in Ae. mcintoshi/luridus, Ae. circumluteolus, and Ae. juppi/caballus [42, 43].

Aedes aegypti was detected in this study across all habitats sampled, although in lower numbers. In urban sites, Matikwane and Pretoria North, Ae. aegypti was the dominant species. In general, this species was predominantly collected in BG-Sentinel traps. This trap is baited with lure as an attractant and is considered the gold standard for collection of Aedes mosquitoes, especially Stegomyia species [44, 45]. Aedes aegypti is the main vector for the four most important arboviruses, dengue (DENV), Zika (ZIKV), chikungunya (CHIKV), and yellow fever (YFV) viruses [46,47,48,49,50]. This species in Africa breeds predominantly in the domestic environment and in the ancestral sylvatic habitat [51]. In South Africa, there is an occurrence of a rural non-anthropophilic population of Ae. aegypti [52], which could explain why conservation areas exhibited lower overall numbers in comparison to urban and peri-urban sites. Previous studies have shown that both forms occur in South Africa: a non-anthropophilic community and anthropophilic populations [52, 53]. Such a finding could explain why some arboviruses such as CHIKV do not spread over large areas in southern Africa and have remained in rural settings [54]. Aedes aegypti was found to feed on domestic cats and on avian species. In other parts of Africa, this species has been detected feeding on humans, monkeys, oxen, goats, cats, dogs, rats, and fowl [55].

Aedes dentatus/leesoni group specimens were mainly collected in the peri-urban Highveld sites in Gauteng Province. Females of this group are often morphologically indistinct, particularly when there is interbreeding or when damage occurs during collection. Aedes dentatus is abundant primarily in the Highveld and feeds readily on humans and larger domestic animals [56]. In this study, the Ae. dentatus group was found feeding on horses in Benoni, a peri-urban site where mosquito sampling was performed in horse farms. Rift Valley, Sindbis (SINV), MIDV, Spondweni (SPOV), and Shokwe (SHOV) viruses were previously isolated from mosquitoes from this group [57,58,59,60,61].

Significant variation was found in the abundance of Aedes caught over the course of this survey. The highest numbers were collected during the wet season, especially between November and mid-April. Many studies have shown the influence of temperature and rainfall on Aedes populations [62,63,64]. The rainy season in most of the northern provinces in southern Africa is from November to April [65]. An increase in rainfall promotes the occurrence of temporary pools, a favourable habitat for many floodwater Aedes species, which include some major arboviral species [15]. The results showed rainfall as being the most influential climatic variable on Aedes populations, as Aedes mcintoshi, Ae. vittatus, and Ae. fowleri peaked soon after the beginning of the rains. This is an expected result, as they are considered floodwater species, with flooding events as a trigger for egg hatching [66]. The apparent drier months during the period from mid-2014 to 2016 when compared to the wetter seasons in 2016/2017 show variations in rainfall within South Africa during the study period [67, 68]. Annual rainfall has been shown to be decreasing in the north-eastern parts of the country [67, 68], while air temperature has been increasing over the past 50 years [75]. In Limpopo Province and parts of northern Mpumalanga Province, the changes are evident, with annual mean rain days translating to nearly 16 days and the temperature increasing by approximately 1 °C over a period of 50 years (1960–2010) [68]. Climate change can alter the geographical distribution of vectors and diseases such as DENV or CHIKV [69,70,71]. Temperature can be a survival-limiting factor for distribution of Aedes species [72]. In our results, Aedes individuals were largely collected at temperatures between 18 °C and 27 °C. Under laboratory conditions, Ae. aegypti has optimal adult longevity at around 21 °C [72], and below 14–15 °C they demonstrate reduced mobility and struggle to take a blood meal [73, 74], which supports the temperature thresholds and collection numbers in this study. Humidity can play a further integral role in mosquito survival, as they desiccate rapidly in dry conditions, thereby decreasing the survival rates [75]. Wind speed can also influence collection events, as high winds affect the flight behaviour of mosquitoes, and a 3 km/h wind speed can significantly reduce the host-seeking flight behaviour of present populations [76]. Humidity and wind speed during the trapping events correlated strongly with the abundance of mosquitoes, supporting such an association with mosquito behaviour and collection.

The genomic sampling in mosquito diversity has been improving recently; however, only 0.8% of the 3556 species recognised by the Mosquito Taxonomic Inventory are represented in NCBI GenBank, and the majority of them belong to Anophelinae [77]. There is also a lack of sequences belonging to African Aedes species in NCBI GenBank and BOLD databases, and this study enabled the addition of COI sequences for nine species: Ae. eritreae, Ae. pachyurus, Ae. ledgeri, Ae. juppi, Ae. microstictus, Ae. aurovenatus, Ae. unidentatus, Ae. durbanensis, and Ae. haworthi. Most of these species were last recorded in the 1950s to 1970s. Aedes eritreae had not been registered since 1971 when McIntosh et al. collected it in a single locality in the former province of Transvaal [56]. That was also true for Ae. pachyurus, where this species had only been collected in the Cape Provinces [78] and KwaZulu-Natal Province [56].

The Neomelaniconion cluster included sequences generated from species of Ae. mcintoshi, Ae. circumluteolus, Ae. unidentatus, and Ae. aurovenatus. With the exception of Ae. aurovenatus, all species clustered together, indicating limitations of the COI gene as a marker to separate these species. The specimens are morphologically similar but can be distinguished based on unguis shapes and wing scale patterns. Remarkably, though, they were very similar genetically. Evidence from Kenya suggests that Ae. mcintoshi is a complex of species, although the results were discordant between the COI and internal transcribed spacer (ITS) [79]. In Madagascar, ITS2 was also insufficient to resolve the species and species complexes in the Neomelaniconion subgenus [80]. Further studies using different markers in combination with morphology are needed to resolve uncertainties in the Neomelaniconion subgenus.

Aedes cumminsii collected in South Africa did not cluster together with sequences from specimens identified as this species in Kenya, Guinea, and Senegal. This species was originally described in Ghana and is widely distributed in Africa. Subspecies based on subtle differences in abdomen scaling have been described, such as ssp. mediopunctatus (Theobald) [18], and it is not surprising that specimens from South Africa are different. Aedes cumminsii likely represents a complex of species which will require further studies to elucidate their taxonomy and importance as vectors of SPOV, MIDV, SINV, and RVF [57, 81,82,83]. McIntosh [56] suggested that only ssp. mediopunctatus occurs in South Africa. However, Jupp [20] used morphological characters from mediopunctatus in his key, naming the species cumminsii. Morphologically, species analysed in this study had a proboscis entirely dark-scaled, scutum indistinctly marked, scutellum with narrow yellow scales, hind tarsal ungues armed, large mosquito with wing length about 5.0 mm, in agreement with descriptions for the typical Ae. cumminsii species. However, they also had a long medium basal band on the abdominal tergites, which according to Edwards [18] is more characteristic of ssp. mediopunctatus. Therefore, the results presented here suggest that ssp. mediopunctatus occurs in north-eastern areas of South Africa. Further studies would be necessary to clarify the status of Ae. cumminsii in South Africa.

Conclusions

There is a paucity of information on Aedes species ecology across southern Africa. Because of the importance of Aedes as disease vectors and the expansion of arboviruses globally over the past two decades due to the increased travel, trade, and climate change, investigations are urgently needed to update the records of sylvatic species. Although the survey area was limited geographically to the north-eastern parts of South Africa, and despite its limited overall trap coverage, the study provides a better understanding of the diversity and relative abundance of Aedes species in these areas where arboviruses are known to be active.

The most abundant aedine species caught were those belonging to the Neomelaniconion subgenus, which are recognised as potential arboviral vectors, especially Ae. mcintoshi, enabling the propagation of arboviral infection in animals and humans in the study area. Aedes aegypti was present in almost all sites sampled, predominantly in urban areas, albeit in low numbers. This study also provided initial data on the feeding behaviour of Aedes mosquitoes in the collection sites, reflecting potential reservoirs for the viruses they transmit in Africa and risk of transmission to human and domestic species. We provided information about the feeding behaviour of sylvatic species including Ae. eritreae/karooensis, Ae. leesoni/alboventralis, Ae. pachyurus, and Ae. aurovenatus, which was not previously available in the literature. The results presented here show that molecular characterisation was able to be successfully employed for species identification of Aedes mosquitoes in South Africa, and contribute to expanding the DNA barcode references available for African Aedes. Apart from confirming the presence and genetic data of several common vector species across the study region, this study added COI sequences for nine additional species to the open-access online databases, namely Ae. eritreae, Ae. pachyurus, Ae. ledgeri, Ae. juppi, Ae. microstictus, Ae. aurovenatus, Ae. unidentatus, Ae. durbanensis, and Ae. haworthi.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its Additional files.

Abbreviations

BG:

Biogents

CHIKV:

Chikungunya virus

CO2 :

Carbon dioxide

COI:

Cytochrome oxidase subunit I

CDC:

Centers for Disease Control and Prevention

DENV:

Dengue virus

ITS:

Internal transcribed spacer

KNP:

Kruger National Park

K2P:

Kimura 2-parameter model

MAFFT:

Multiple alignment using fast Fourier transform

MEGA:

Molecular evolutionary genetics analysis

MIDV:

Middelburg virus

NICD:

National Institute for Communicable Diseases of South Africa

RVF:

Rift Valley fever

SAWS:

South African Weather Service

SINV:

Sindbis virus

SHOV:

Shokwe virus

SPOV:

Spondweni virus

spp:

Species

ssp:

Subspecies

WSLB:

Wesselsbron virus

YFV:

Yellow fever virus

ZIKV:

Zika virus

References

  1. Tolle MA. Mosquito-borne diseases. Curr Probl Pediatr Adolesc Health Care. 2009;39(4):97–140.

    Article  PubMed  Google Scholar 

  2. Beerntsen BT, James AA, Christensen BM. Genetics of mosquito vector competence. Microbiol Mol Biol Rev. 2000;64(1):115–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Weaver SC, Reisen WK. Present and future arboviral threats. Antiviral Res. 2010;85(2):328–45.

    Article  CAS  PubMed  Google Scholar 

  4. Leisnham PT, Juliano SA. Impacts of climate, land use, and biological invasion on the ecology of immature Aedes mosquitoes: implications for La Crosse emergence. EcoHealth. 2012;9(2):217–28.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Ndoen E, Wild C, Dale P, Sipe P, Dale M. Relationships between anopheline mosquitoes and topography in West Timor and Java, Indonesia. Malar J. 2010;9.

  6. Takken W, Verhulst NO. Host preferences of blood feeding mosquitoes. Annu Rev Entomol. 2013;58:433–53.

    Article  CAS  PubMed  Google Scholar 

  7. Kilpatrick AM, Kramer LD, Jones MJ, Marra PP, Daszak P. West Nile virus epidemics in North America are driven by shifts in mosquito feeding behavior. PLoS Biol. 2006;4(4):e82.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Hebert PDN, Cywinska A, Ball SL, deWaard JR. Biological identifications through DNA barcodes. Proc R Soc Lond B Biol Sci. 2003;270:313–21.

    Article  CAS  Google Scholar 

  9. Erlank E, Koekemoer LL, Coetzee M. The importance of morphological identification of African anopheline mosquitoes (Diptera: Culicidae) for malaria control programmes. Malar J. 2018;17(1):43.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Norris LC, Norris DE. Phylogeny of anopheline (Diptera: Culicidae) species in southern Africa, based on nuclear and mitochondrial gene. J Vector Ecol. 2015;40:16–27.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Huang YM. A new african species of Aedes (Diptera: Culicidae). Mosq Syst. 1985;17(2).

  12. Muspratt J. Research on South African Culicini (Diptera, Culicidae). Ill. A check-list of the species and their distribution, with noles on taxonomy, bionomics and identification. J Entomol Soc South Afr. 1955;18(2).

  13. Jupp PG. Culex pipens pipens and Culex pipens quinquefasciatus in South Africa: morphological and reproduction evidence in facour of their status as two species. Mosq Syst. 1978;10(4).

  14. Muspratt J. The Stegomyia mosquitoes of South Africa and some neighbouring territories. Memoirs Entomol Soc South Afr. 1956;4:1–38.

    Google Scholar 

  15. Cornel AJ, Lee Y, Almeida APG, Johnson T, Mouatcho J, Venter M, et al. Mosquito community composition in South Africa and some neighboring countries. Parasit Vectors. 2018;11(1):331.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Gorsich EE, Beechler BR, van Bodegom PM, Govender D, Guarido MM, Venter M, et al. A comparative assessment of adult mosquito trapping methods to estimate spatial patterns of abundance and community composition in southern Africa. Parasit Vectors. 2019;12(1):462.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Schrama M, Hunting ER, Beechler B, Guarido MM, Govender D, Nijland W, et al. Human practices promote presence and abundance of disease-transmitting mosquito species. Sci Rep. 2020.

  18. Edwards FW. Mosquitoes of the Ethiopian Region. III-Culicine Adults and Pupae: The Oxford University Press; 1941.

  19. Gillies MT, Coetzee M. A supplement to the Anophelinae of South of the Sahara (Afrotropical Region) Johannesburg: The South African Institute for Medical Research; 1987.

  20. Jupp PG. Mosquitoes os Southern Africa. Culicinae and Toxorhynchitinae: Ekogilde Publishers; 1996.

  21. Service MW. Handbook to the Afrotropical toxorhynchitine and Culicine mosquitoes, excepting Aedes and Culex. London: British Museum (Natural History); 1990.

  22. Wilkerson RC, Linton YM, Fonseca DM, Schultz TR, Price DC, Strickman DA. Making mosquito taxonomy useful: a stable classification of tribe Aedini that balances utility with current knowledge of evolutionary relationships. PLoS ONE. 2015;10(7):e0133602.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Service SAW. About Climate at SAWS: Our Services 2019 [Available from: http://www.weathersa.co.za/.

  24. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrate. Mol Mar Biol Biotechnol. 1994;3(5):294–9.

    CAS  PubMed  Google Scholar 

  25. Marcon E, Hérault B. entropart: AnRPackage to Measure and Partition Diversity. Journal of Statistical Software. 2015;67(8).

  26. R Core Team R: A language and environment for statisticalcomputing. Vienna, Austria: R Foundation for Statistical Computing; 2003 [Available from: https://CRAN.R-project.org/doc/FAQ/R-FAQ.html.

  27. Corporation M. Microsoft Office Excel Redmond, WA: Microsoft Corporation; 2010.

  28. CLC Genomics Workbench 8.0.1: Qiagen Bioinformatics.; [Available from: https://www.qiagenbioinformatics.com/.

  29. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Feselsntein J. Phylogenies and the comparative method. Am Nat. 1985;125:1–15.

    Article  Google Scholar 

  31. Pachka H, Annelise T, Alan K, Power T, Patrick K, Veronique C, et al. Rift Valley fever vector diversity and impact of meteorological and environmental factors on Culex pipiens dynamics in the Okavango Delta, Botswana. Parasit Vectors. 2016;9(1):434.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Rueda LM. Global diversity of mosquitoes (Insecta: Diptera: Culicidae) in freshwater. In: Freshwater animal diversity assessment. Developments in hydrobiology 2008. p. 477–87.

  33. Khomar N, Pollack RJ, Spielman A. A nestable fiber pot for sampling resting mosquitoes. J Am Mosq Control Assoc. 1995;11(4):463–7.

    Google Scholar 

  34. Panella NA, Crockett RJ, Biggerstaff BJ, Komar N. The Centers for Disease Control and Prevention resting trap: a novel device for collecting resting mosquitoes. J Am Mosq Control Assoc. 2011;27(3):323–5.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Schonenberger AC, Wagner S, Tuten HC, Schaffner F, Torgerson P, Furrer S, et al. Host preferences in host-seeking and blood-fed mosquitoes in Switzerland. Med Vet Entomol. 2016;30(1):39–52.

    Article  CAS  PubMed  Google Scholar 

  36. Ponlawat A, Harrington LC. Blood feeding patterns of Aedes aegypti and Aedes albopictus in Thailand. J Med Entomol. 2005;42:844–9.

    Article  PubMed  Google Scholar 

  37. Jupp PG, McIntoshi BM, Nevill M. A survey of the mosquito and Culicoides faunas at two localities in the Karoo region of South Africa with some observations on bionomics. Onderstepoort J Vet Res. 1980;47:1–6.

    CAS  PubMed  Google Scholar 

  38. McIntoshi BM. The aedine subgenus Neomelaniconion Newstead (Culicidae, Diptera) in southern Africa with descriptions of two new specie. J Entomol Soc South Afr. 1971;34(2):319–33.

    Google Scholar 

  39. Jupp PG. Mosquitoes as vectors of human disease in South Africa. South African Family Pract. 2014;47(9):68–72.

    Article  Google Scholar 

  40. National Institute for Communicable Diseases. Viral haemorrhagic fever outbreaks, South Africa, 2011. Communicable Diseases Surveillance Bulletin. 2012;10(1).

  41. National Institute for Communicable Diseases. An update on human cases of Rift Valley fever on a farm in the Jacobsdal area, Free State Province. NICD Communicable Diseases Communiqué. 2018.

  42. Jupp PG, Kemp A. Studies on an outbreak of Wesselsbron virus in the Free State Province, South Africa. J Am Mosq Control Assoc. 1998;14(1):40–5.

    CAS  PubMed  Google Scholar 

  43. Smithburn KC, Kokernot RH, Weinbren MP, de Meillon B. Studies on arthropod-borne viruses of Tongaland. VII. Simbu virus, a hitherto unknown agent isolated from Aedes (Banksinella) circumluteolus Theo. S Afr J Med Sci. 1957;22:113–20.

    CAS  PubMed  Google Scholar 

  44. KrÖCkel U, Rose A, Eiras ÁE, Geier M. New tools for surveillance of adult yellow fever mosquitoes: comparison of trap catches with human landing rates in an urban environment. J Am Mosq Control Assoc. 2006;22(2):229–38.

    Article  PubMed  Google Scholar 

  45. Li Y, Su X, Zhou G, Zhang H, Puthiyakunnon S, Shuai S, et al. Comparative evaluation of the efficiency of the BG-Sentinel trap, CDC light trap and Mosquito-oviposition trap for the surveillance of vector mosquitoes. Parasit Vectors. 2016;9(1):446.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Barrett AD, Higgs S. Yellow fever: a disease that has yet to be conquered. Annu Rev Entomol. 2007;52:209–29.

    Article  CAS  PubMed  Google Scholar 

  47. Burt FJ, Chen W, Miner JJ, Lenschow DJ, Merits A, Schnettler E, et al. Chikungunya virus: an update on the biology and pathogenesis of this emerging pathogen. Lancet Infect Dis. 2017;17(4):e107–17.

    Article  CAS  PubMed  Google Scholar 

  48. Kindhauser MK, Allen T, Frank V, Santhana RS, Dye C. Zika: the origin and spread of a mosquito-borne virus. Bull World Health Organ. 2016;94(9):675–86.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Jansen CC, Beebe NW. The dengue vector Aedes aegypti: what comes next. Microbes Infect. 2010;12(4):272–9.

    Article  PubMed  Google Scholar 

  50. Thangamani S, Huang J, Hart CE, Guzman H, Tesh RB. Vertical transmission of Zika Virus in Aedes aegypti mosquitoes. Am J Trop Med Hyg. 2016;95(5):1169–73.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Weetman D, Kamgang B, Badolo A, Moyes CL, Shearer FM, Coulibaly M, et al. Aedes mosquitoes and Aedes-borne arboviruses in Africa: current and future threats. Int J Environ Res Public Health. 2018;15(2):220.

    Article  PubMed Central  Google Scholar 

  52. Jupp PG, Kemp A, Frangos C. The potencial for dengue in South Africa: morphology and the taxonomic status of Aedes aegypti populations. Mosq Syst. 1991;23(3).

  53. Kemp A, Jupp PG. Potential for dengue in South Africa Mosquito Ecology with particular reference to Aedes aegypti. J Am Mosq Control Assoc. 1991;7(4).

  54. Jupp PG, Kemp A. What is the potential for future outbreaks of chikungunya, dengue and yellow fever in southern Africa? SAMJ (S Afr Med J). 1996;86(1).

  55. Weitz B. Feeding habits of bloodsucking arthropods. Exp Parasitol. 1960;9(1):63–82.

    Article  Google Scholar 

  56. McIntoshi BM. A taxonomic revision of certain Aedes species (Diptera: Culicidae) of the subgenus Aedimorphus in southern Africa. J Entomol Soc S Afr. 1975;38(2):251–87.

    Google Scholar 

  57. McIntoshi BM, Kokernot RH, Paterson HE, De Meillon B. Isolation of Spondweni virus from four species of culicine mosquitoes and a report of two laboratory infections with the virus. SAMJ (S Afr Med J). 1961;35:647–50.

    Google Scholar 

  58. Centers for Disease Control and Prevention. Arbovirus Catalog.

  59. McIntoshi BM. Rift valley fever. 1. Vectors studies in the field. J S Afr Vet Assoc. 1972;43(4):391–5.

    Google Scholar 

  60. McIntoshi BM. Isolation of Semliki Forest virus from Aedes (Aedimorphus) argenteopunctatus (Theobald) collected in portuguese East Africa. Trans R Cosoc Trop Med Hyg. 1961;55(2):192–8.

    Article  Google Scholar 

  61. Worth CB, Paterson HE, De Meillon B. The incidence of arthropod-borne viruses in a population of culicine mosquitoes in Tongaland, Union of South Africa (January, 1956, through April, 1960). Am J Trop Med Hyg. 1961;10:583–92.

    Article  CAS  PubMed  Google Scholar 

  62. Mello PS, Leake CJ. Climatic and geographic influences on arboviral infections and vectors. Rev Sci Tech Off Int Epizoot. 2000;19(1):41–54.

    Article  Google Scholar 

  63. Valdez LD, Sibona GJ, Condat CA. Impact of rainfall on Aedes aegypti populations. Ecol Modell. 2018;385:96–105.

    Article  Google Scholar 

  64. Wong MC, Mok HY, Ma HM, Lee MW, Fok MY. A climate model for predicting the abundance of Aedes mosquitoes in Hong Kong. Meteorol Appl. 2011;18(1):105–10.

    Article  Google Scholar 

  65. White F. The Vegetation of Africa: Unesco; 1983.

  66. Gillett JD. Contribution to the oviposition-cycle by the individual mosquitoes in a population. J Insect Physiol. 1962;8:665–81.

    Article  Google Scholar 

  67. Kruger AC, Nxumalo MP. Historical rainfall trends in South Africa: 1921–2015. Water SA. 2017;43(2):285.

    Article  Google Scholar 

  68. MacKellar N, New M, Jack C. Observed and modelled trends in rainfall and temperature for South Africa: 1960–2010. S Afr J Sci. 2014;110(7/8):1–13.

    Article  Google Scholar 

  69. Fischer D, Thomas SM, Suk JE, Sudre B, Hess A, Tjaden NB, et al. Climate change effects on Chikungunya transmission in Europe: geospatial analysis of vector’s climatic suitability and virus’ temperature requirements. Int J Health Geogr. 2013;12:51.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Leta S, Beyene TJ, De Clercq EM, Amenu K, Kraemer MUG, Revie CW. Global risk mapping for major diseases transmitted by Aedes aegypti and Aedes albopictus. Int J Infect Dis. 2018;67:25–35.

    Article  PubMed  Google Scholar 

  71. Li R, Xu L, Bjørnstad ON, Liu K, Song T, Chen A, et al. Climate-driven variation in mosquito density predicts the spatiotemporal dynamics of dengue. Proc Natl Acad Sci USA. 2019;119(9):3624–9.

    Article  CAS  Google Scholar 

  72. Brady OJ, Johansson MA, Guerra CA, Bhatt S, Golding N, Pigott DM, et al. Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings. Parasites Vectors. 2013;6.

  73. Nawrock SJ, Hawley WA. Estimation of the northern limits of distribution of Aedes albopictus in North America. J Am Mosq Control Assoc. 1987;3(2):314–7.

    Google Scholar 

  74. Yang HM, Macoris ML, Galvani KC, Andrighetti MT, Wanderley DM. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiol Infect. 2009;137(8):1188–202.

    Article  CAS  PubMed  Google Scholar 

  75. Patz JA, Githeko AK, McCarty JP, Hussein S, Confalonieri U, de Wet N. Climate change and infectious diseases. Climate Change and Human Health: Risks and Responses. Geneva: World Health Organization; 2003.

  76. Service MW. Effects of wind on the behaviour and distribution of mosquitoes and blackflies. Int J Biometeorol. 1980;24(4):347–53.

    Article  Google Scholar 

  77. Ruzzante L, Reijnders M, Waterhouse RM. Of genes and genomes: mosquito evolution and diversity. Trends Parasitol. 2019;35(1):32–51.

    Article  CAS  PubMed  Google Scholar 

  78. Muspratt J. Research on South African Culicini (Diptera, Culicidae) II—taxonomy relating to eight species of Aedes. J Entomol Soc South Afr. 1953;16(2):83–93.

    Google Scholar 

  79. Tchouassi DP, Bastos AD, Sole CL, Diallo M, Lutomiah J, Mutisya J, et al. Population genetics of two key mosquito vectors of Rift Valley Fever virus reveals new insights into the changing disease outbreak patterns in Kenya. PLoS Negl Trop Dis. 2014;8(12):e3364.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Kengne P, Le Goff G, Fontenille D. Molecular genetic investigation of morphological species of members of the Neomelaniconion subgenus (Diptera: Culicidae: Aedini) from Madagascar using ribosomal internal transcribed spacer 2. J Med Entomol. 2009;46(2):403–7.

    Article  PubMed  Google Scholar 

  81. Robin Y, Cornet M, Brès P, Héry G, Chateau R. Isolement d'une souche de virus middelburg à partir d'un lot d'Aedes (A.) cumminsi récoltés à Bandia (Sénégal). Bull Soc Pathol Exot. 1969;62:112–118.

  82. McIntosh BM, Jupp PG, De Sousa J. Further isolations of the arboviruses from mosquitoes collected in Tongaland, South Africa, 1960–1968. J Med Entomol. 1972;9:155–9.

    Article  CAS  PubMed  Google Scholar 

  83. Paterson HE, Bronsden P, Levitt J, Worth CB. Some Culicine mosquitoes (Diptera, Culicidae) at Ndumu, Republic of South Africa. A study of their host preferences and host range. Med Proc. 1964;10:188–92.

    Google Scholar 

Download references

Acknowledgements

We are grateful to the staff at Lapalala Wilderness, Marataba Conservation, and South African National Parks for logistical assistance and permission to collect mosquitoes. We acknowledge Danny Govender for her valuable contribution in the collection of the mosquitoes.

Funding

This study was supported by a cooperative agreement (no. 5 NU2GGH001874-02-00) funded by the Centers for Disease Control and Prevention. The findings and conclusions in this manuscript are those of the authors and do not necessarily represent the official position of the funding agencies. This study was also funded in part through scholarships from the South African National Research Foundation and The World Academy of Sciences. A.P.G. Almeida has been a recipient of the Visiting Professor Programme by the University of Pretoria, SA, and acknowledges the Global Health and Tropical Medicine unit. This study was also funded by the Gratama Fund (grant number 2016.08), the Uyttenboogaart-Eliasen Foundation (SUB.2016.12.08), and the RCN-IDEAS travel grant to M. Schrama and E.E. Gorsich.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualisation and design the experiments, MMG, APGA, MV and BDB. The analyses were performed by MMG, MAR, EEG and MS. Field work and identification of mosquitoes were done by MMG, MAR, LEOB, TJ, EEG and MS. MMG wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Marietjie Venter.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. Collection sites for mosquito samples surveyed from January 2014 to May 2018. Table S2. Number of Aedes collected at sentinel sites, relative abundance (%), and diversity indexes, January 2014 to June 2018. Table S3. Number of Aedes collected at ad hoc sites, relative abundance (%), and diversity indexes, January 2014 to June 2018. Table S4. Comparison of the abundance of Aedes captured by the three trap types at the sentinel sites from 2017 and 2018. Table showing the results of the Kruskal–Wallis multiple tests. Table S5. Correlation coefficient between Aedes per trap-night and the environmental conditions. Table S6. Aedes blood meals identified in various mammal and avian species. Figure S1. Relative abundance (%) of Aedes vectors per trap-night from January 2014 to June 2018 at the sentinel sites. Figure S2. Most common potential Aedes vectors abundance per trap-night, rainfall, and temperatures from January 2014 to June 2018 at the sentinel sites. Figure S3. Regression relationship between log10 of the Aedes per trap-night and average temperature 30 and 15 days prior to collection, average rainfall 30 and 15 days prior to the collection event, and average humidity 30 and 15 days prior to collection, January 2014 to June 2018. Figure S4. The morphological characteristics of an adult female of Ae. cumminsii collected in South Africa. A) Dorsal view; B) Lateral view.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Guarido, M.M., Riddin, M.A., Johnson, T. et al. Aedes species (Diptera: Culicidae) ecological and host feeding patterns in the north-eastern parts of South Africa, 2014–2018. Parasites Vectors 14, 339 (2021). https://doi.org/10.1186/s13071-021-04845-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-021-04845-9

Keywords