Human practices promote presence and abundance of disease-transmitting mosquito species

Humans alter the environment at unprecedented rates through habitat destruction, nutrient pollution and the application of agrochemicals. This has recently been proposed to act as a potentially significant driver of pathogen-carrying mosquito species (disease vectors) that pose a health risk to humans and livestock. Here, we use a unique set of locations along a large geographical gradient to show that landscapes disturbed by a variety of anthropogenic stressors are consistently associated with vector-dominated mosquito communities for a wide range of human and livestock infections. This strongly suggests that human alterations to the environment promote the presence and abundance of disease vectors across large spatial extents. As such, it warrants further studies aimed at unravelling mechanisms underlying vector prevalence in mosquito communities, and opens up new opportunities for preventative action and predictive modelling of vector borne disease risks in relation to degradation of natural ecosystems.

www.nature.com/scientificreports/ but differ in how they are affected by humans. To obtain an overall assessment of anthropogenic disturbance, we quantified 5 ubiquitous anthropogenic pressures 1 : (1) organophosphate pesticide abundance, (2) eutrophication,

Results and discussion
The five pressures were consistently higher outside vs.  (Table S2). These pressures generate a composite measure of disturbance ( Fig. 1) and confirm that water bodies outside KNP are consistently more impacted by anthropogenic pressures. To quantify whether these higher levels of disturbance are linked to changes in mosquito communities, we trapped 3,918 females of 39 mosquito species and species complexes (Diptera: Culicidae). The paired trapping design ( Fig. 1) allowed us to simultaneously collect mosquitoes from waterbodies inside and outside KNP 31 . Mosquito abundance outside the national park was on average 2.9 times higher (range: 1.5-10 times higher) than inside KNP (F (1,24) = 17.3, p < 0.001; Fig. 2A-D). This pattern resembles mosquito abundances along gradients of human disturbance in the Mediterranean 24 and tropical rainforests in Thailand 29 . More importantly, we observed pronounced shifts in the species composition of mosquitoes inside vs. outside KNP (ANOSIM: r = 0.15; p < 0.001) for each of the four geographical regions ( Fig. 2A-D, Table S3). Despite differences in overall richness (35 species outside vs 30 inside KNP) and abovementioned differences in species composition, we observed no differences between different alpha diversity metrics inside vs. outside KNP (species richness, Shannon's J, Shannon's H'max and Simpson's evenness; Figures S2A-C). Known vectors of human and livestock diseases were more abundant outside vs. inside KNP; these species explained 79% (± SD 3) of the variation in communities between paired regions ( Fig. 2E; SIMPER: r = 0.48; p < 0.001). This pattern was consistent for each of the regions (Malelane 82%, Satara 78%, Punda Maria 76%, Skukuza 77%; Figures S3A-D, S4). It was also largely consistent across vectors and infections: Aedes aegypti, a vector for dengue, chikungunya, yellow fever and Zika 32 , was more common outside KNP. Multiple Culex species were also more common outside, including vectors for West Nile, sindbis, Wesselsbron and filariasis ( Fig. 2E-F). Only two vector species, Ae. vexans and Cx. theileri, vectors from Rift Valley fever 32 , were more abundant inside KNP. These patterns are in line with studies carried out along gradients of deforestation, which have been shown to lead to increased presence of malaria vectors such as An. gambiae and An. darlingi 18,[33][34][35] .
Together, these results show that human disturbances are strongly associated with increased mosquito abundances and shifts in community composition towards known disease vectors, as illustrated in conceptual Fig. 2G. They highlight the experimentally-demonstrated importance of changes in environmental conditions 5,6 by identifying its consistent pattern across large geographic extents and climatic conditions. These results emphasize that anthropogenic disturbances to the landscape and the larval environment represent an important mechanism driving vector distributions. Given the global extent and intensity of the investigated anthropogenic pressures 1 , these results are likely relevant for a wide array of vector-borne pathogens and provides a mechanism for the association between ecosystem degradation and disease 36,37 . This raises important questions on how different human activities drive vector prevalence in mosquito communities, and presents new opportunities for targeted preventative action as well as predictive modelling of vector borne disease risks in relation to ecosystem services.

Methods
Mosquito trapping and identification have been described previously 31 ; mosquito sampling occurred across 4 regions and 112 trapping nights, with simultaneous collections for 3-4 consecutive nights at all sampling points in a region (4 inside and 4 outside KNP). Sampling points within a region were selected based on multiple criteria. The primary selection criterion was to sample from waterbodies that were representative of the region, including a diversity of wetlands rather than those with the highest catch rates 28 . Additional criteria stipulated that the water bodies were at least 1 km away from one another to avoid sampling mosquitoes from adjacent water bodies, as mean mosquito dispersal distances range from 35 m to 1.4 km 31 . In this setup, we follow the ecosystem boundaries framework 38 and quantified 5 pressures associated with human impacts in the different ecosystems. The concentration of different persistent organophosphates (POPs, e.g., DDT and breakdown products) were measured at 2 sampling points per region using multi-residue analysis (GC-ECD and GC-MS) by the African Research Council (Roodeplaat, South Africa). We determined phosphate levels (PO 4 3− ) using a photospectrometer (Merck Spectroquant Nova 60) in one litre of water, which was composed of 20 subsamples of 50 ml per sampling point. The fraction of bare and sparsely vegetated area at each sampling point was assessed using satellite data derived from the Sentinel-2 sensor acquired in January 2017. The 2-band Enhanced Vegetation Index (EVI) 39 was calculated from a monthly maximum EVI composite avoiding atmospheric disturbances. To derive the percentage of barren and sparsely vegetated areas EVI thresholds of 0.15 and 0.35 were used respectively. The proportion of barren and sparsely vegetated pixels (10*10 m) within a 150 m radius around each sample point was taken as a representative sample of the vegetation cover for each of the trapping locations. The fraction of bare and sparsely vegetated area was assessed using monthly maximum Enhanced Vegetation Index (EVI; using Sentinel-2data Jan-March 2017). In a radius of 150 m around each sample point, we estimated the percentage of pixels (10*10 m) between EVI thresholds of 0.15 and 0.35, respectively. Animal densities at each sampling point were determined using the gridded livestock of the world map (fao.org/livestock-systems/global-distributions; resolution: 0.05′*0.05′ degrees ~ 5*5 km) and the 2018 KNP African buffalo (Syncerus caffer) counts. All numbers were transferred to a raster file with 0.05′*0.05′ grid cells, after which the average density was calculated in a buffer Scientific RepoRtS | (2020) 10:13543 | https://doi.org/10.1038/s41598-020-69858-3 www.nature.com/scientificreports/ www.nature.com/scientificreports/ ring of 2 km diameter around each sampling point. Human population densities for each of the regions inside and outside KNP were determined using the 2011 population census using data at the scale of the municipality (https ://www.stats sa.gov.za).
To test for differences in mosquito abundance inside vs. outside KNP, we used general linear models comparing the number of mosquitos collected at each site aggregated across each night. Region (Punda Maria, Satara, Skukuza and Malelane), disturbance (inside vs. outside KNP) and their interaction were included as main effects. The number of mosquitoes collected was overdispersed, log-transformed for normality, and assessed using Quantile Quantile-plots and a Levene's test (P = 0.06). Independence assumptions of the regression model were evaluated using plots of model residuals by location and distance between sites.
Differences in community composition inside vs. outside KNP were tested based on the species composition data aggregated across all nights at each site, because multiple trapping nights are needed to capture rare species 26,31 . The data were analysed using a non-parametric analysis of similarities analysis (ANOSIM) and visualized with non-metric multidimensional scaling (NMDS). The ANOSIM analysis is a non-parametric test for differences in mosquito communities among traps that compares the ranks of Bray-Curtis dissimilarity measures from samples collected inside vs. outside KNP 40,41 . A SIMPER analysis was used to assess which taxa are responsible for shifts in community composition 42 . Overall patterns of richness and diversity (Shannon's J, Shannon's H'max and Simpson's evenness) are also reported ( Fig. S2A-C), but the SIMPER analysis provides is more appropriate, because it was developed to identify the species responsible for shifts in community measures 42 . All analyses were conducted in R using the lme4 and vegan packages 43,44 .