Urbanization favors the proliferation of Aedesaegypti and Culexquinquefasciatus in urban areas of Miami-Dade County, Florida

Urbanization processes are increasing globally. Anthropogenic alterations in the environment have profound effects on biodiversity. Decreased biodiversity due to biotic homogenization processes as a consequence of urbanization often result in increased levels of mosquito vector species and vector-borne pathogen transmission. Understanding how anthropogenic alterations in the environment will affect the abundance, richness, and composition of vector mosquito species is crucial for the implementation of effective and targeted mosquito control strategies. We hypothesized that anthropogenic alterations in the environment are responsible for increasing the abundance of mosquito species that are adapted to urban environments such as Aedes aegypti and Culex quinquefasciatus. Therefore, our objective was to survey mosquito relative abundance, richness, and community composition in Miami-Dade County, Florida, in areas with different levels of urbanization. We selected 24 areas, 16 remote areas comprised of natural and rural areas, and 8 urban areas comprised of residential and touristic areas in Miami-Dade County, Florida. Mosquitoes were collected weekly in each area for 24 h for 5 consecutive weeks from August to October 2020 using BG-Sentinel traps baited with dry ice. A total of 36,645 mosquitoes were collected, from which 34,048 were collected in the remote areas and 2,597 in the urban areas. Our results show a clear and well-defined pattern of abundance, richness, and community composition according to anthropogenic modifications in land use and land cover. The more urbanized a given area the fewer species were found and those were primary vectors of arboviruses, Ae. aegypti and Cx. quinquefasciatus.

There are over one billion cases of mosquito-borne diseases reported worldwide every year 1 . The epidemiology of mosquito-borne diseases is affected by how modifications of the natural environment alter the interactions between vector mosquito species, hosts, and pathogens 2 . Anthropogenic changes in the environment in the form of climate change, urbanization, and biodiversity loss favor the proliferation of some vector mosquito species, such as Aedes aegypti and Aedes albopictus, and are important drivers for arbovirus transmission in urban areas [3][4][5][6][7][8][9][10][11] .
The epidemiology of mosquito-borne diseases continues to be uncertain since the extent to which mosquito species are dispersing across spatiotemporal scales and how urbanization processes are affecting their presence and abundance remains unclear 2,12,13 . Recent estimates have shown that 129 of the world's countries and territories have conducive environments for the proliferation of mosquito vector species and are at risk of dengue transmission 1,14 .
Urbanization processes are ubiquitous in the contemporary world. Land use and land cover transformation to better suit the needs of the human population have led to the degradation of natural environments 15,16 . The human footprint can be found in most areas of the planet, including uninhabited areas [17][18][19] . The human

Results
Mosquitoes were collected in 24 collection sites across Miami-Dade County, from which 16 were located in remote areas and 8 in urban areas. A total of 36,645 mosquitoes were collected, from which 34,048 were collected in the remote areas and 2,597 in the urban areas. A total of 26 species were collected in the remote sites and 20 in the urban sites. The species richness per collection site ranged from 8 to 19 species among the remote sites with an average of 11 species per site, whereas among the urban sites the species richness ranged from 3 to 18 with an average of 6 species per site. However, the urban site U8 is a new urban development and borders natural areas. If the urban site U8 had been removed from the analyses the species richness among the urban sites would have ranged from 3 to 7 species with an average of 5 species per site.
All species collected in the urban sites were also collected in the remote sites. However, the opposite was not true, Aedes bahamensis, Anopheles atropos, Coquillettidia perturbans, Psorophora ferox, Uranotaenia lowii, and Uranotaenia sapphirina were collected in the remote sites but not on the urban sites. The most abundant species in the remote areas were Cx. nigripalpus (13,338 females and 3 males), followed by Anopheles crucians (6,509 females and 7 males) and Culex erraticus (6,086 females and 18 males). On the other hand, the most abundant species in the urban areas were Ae. aegypti (762 females and 399 males), followed by Cx. nigripalpus (547 females and 7 males) and Cx. quinquefasciatus (247 females and 185 males) ( Table 1).
The mean number of mosquitoes collected in the remote and urban collection sites revealed that most of the species were more commonly found in the remote areas. From the 26 species collected, only 5 species had a mean number higher than 20% in the urban areas, Ae. aegypti, Ae. triseriatus, Cx. quinquefasciatus, Mansonia titillans, and Wyeomyia vanduzeei. Furthermore, Ae. aegypti was the most commonly found species in urban areas with a mean value higher than 90%. On the other hand, Ae. albopictus, An. crucians, Anopheles quadrimaculatus, and Cx. erraticus were almost exclusively collected in the remote areas (Fig. 1).  www.nature.com/scientificreports/ Our results indicate a robust pattern of mosquito relative abundance and species richness and composition between remote and urban collection sites. The remote collection sites had a higher species richness with many abundant species that were not commonly found in the urban collection sites such as Aedes atlanticus, An. quadrimaculatus, and Cx. erraticus. On the other hand, Ae. aegypti and Cx. quinquefasciatus were more commonly found in the urban collection sites. Furthermore, the species richness and abundance were higher in areas less impacted by urbanization, being inversely proportional to the presence of buildings and proximity to roads (Fig. 2).
The analysis of the percent composition of vector mosquito species in the remote and urban collection sites revealed that the urban collection sites had fewer and more dominant mosquito vector species than the remote collection sites. Aedes aegypti and Cx. quinquefasciatus comprised most of the mosquitoes collected in the urban sites with Ae. aegypti being the most abundant species collected in these sites except collection site U8, which was located in a natural-urban transition zone. On the other hand, Ae. aegypti and Cx. quinquefasciatus were not collected in 6 of the remote sites and apart from the remote sites R4 and R5 were collected in negligible numbers in the remote collection sites. On the remote collection sites, Cx. nigripalpus was the dominant species but not as dominant as Ae. aegypti and Cx. quinquefasciatus in the urban collection sites. These results showed a clear pattern of species richness and relative abundance in both the remote and urban collection sites according to the environment and respective resource availability. The mosquito community composition was substantially different in remote and urban collection sites disregarding geographical proximity as observed in the comparison between the collection sites U1 and R2 and R3, U2 and R15 and R16, and U7 and R14 (Fig. 3).
The Permutational Multivariate Analysis of Variance (PERMANOVA) yielded significant results for the comparison between mosquito species composition in remote and urban areas (F = 8.707; P = 0.0001). The PERMANOVA analysis also revealed significant differences in the mosquito species composition in areas with different land uses (F = 3.205; P < 0.0003). The subsequent SIMPER (Similarity Percentage) analysis comparing remote and urban areas showed that Cx. nigripalpus, An. crucians, Cx. erraticus, and Ae. aegypti contributed the most to the observed differences, whereas in the analysis considering different land uses Cx. nigripalpus, An. crucians, Cx. erraticus, and Ae. atlanticus contributed the most to the observed differences (Supplementary Table 2).
Results from the Moran's Index of Spatial Autocorrelation were not statistically significant failing to reject the null hypothesis, therefore, indicating that autocorrelation is not present in the samples analyzed here (Observed: − 0.002; Expected: − 0.043; Standard Deviation: 0.050; and P-value: 0.420). The Generalized Linear Mixed Methods (GLMM) regression results showed a statistically significant association between the Normalized Difference www.nature.com/scientificreports/ Vegetation Index (NDVI) and the mean distance from major roads and their interaction with species richness ( Table 2). Species richness and relative abundance served as the dependent variables and were analyzed using a normal distribution for each urbanization measurement (Normalized Difference Vegetation Index (NDVI), 2-D Building Footprint, and the mean distance from major roads) and their interaction as independent variables, and week as a random effect with collection sites nested in time (weeks).

Discussion
Our results show a clear and well-defined pattern of abundance, richness, and community composition according to anthropogenic modifications in land use and land cover in Miami-Dade. The more urbanized a given area the fewer species were found and Ae. aegypti and Cx. quinquefasciatus were the most dominant species. Our study shows that land use and land cover transformation of natural areas into urban areas of South Florida can affect mosquito relative abundance, richness, and community composition by favoring mosquito vector species that are adapted to thrive in urban environments. As a consequence, human populations will have increased contact with mosquito vector species, especially Ae. aegypti and Cx. quinquefasciatus, increasing the risk of vector-borne disease transmission. Figure 2. Impact of urbanization on species richness and mosquito relative abundance. Left: Species richness relationship with NDVI, building footprint, and mean distance from roads (in meters); Right: Relative abundance relationship with NDVI, building footprint, and mean distance from roads (in meters). www.nature.com/scientificreports/ Our results also show that the community composition, species richness, and relative abundance were not related to the geographic proximity of the remote and urban collection sites. Even adjacent remote and urban collection sites had completely distinct mosquito community composition, species richness, and relative abundance. These results support the hypothesis that the availability of resources at the micro-geographic scale (i. e., neighborhood-level) is a major driver for the proliferation or decrease of populations of vector mosquito species according to their specific ecology and behavior 22,58 . For example, both Ae. aegypti and Cx. quinquefasciatus were collected in relatively high numbers in the remote areas R4 and R5. The remote area R4 is a conservation  www.nature.com/scientificreports/ zone closed to the general public, however, an abandoned checkpoint flooded with rainwater was a potentially conducive habitat for the proliferation of Ae. aegypti and Cx. quinquefasciatus. A similar scenario was observed in the remote area R5, in which the presence of a human-made ditch to divert rainwater also made available suitable habitats for these vector species. The transformation of natural environments into urban spaces and the increased density of urban land cover surfaces (e.g., asphalt, concrete, etc.) within existing urbanized locations is a major driver for biodiversity loss 6,59,60 . As a result, the steep increase in urbanization in the last decades has led to biodiversity loss on a global scale 15 .
In this context, the mosquito community composition is greatly affected by urbanization, in which mosquito abundance and species richness are significantly affected and tend to decrease proportionally with urbanization levels [61][62][63][64] . However, some mosquito vector species can thrive in urban environments and greatly benefit from the resources made available by the rise of transient underserved populations, absence of natural predators, human-made aquatic habitats, increased human population density, and warmer temperatures due to global warming. Aedes aegypti and Cx. quinquefasciatus are among those species that can thrive in urban areas, and coincidentally are the primary vectors of chikungunya, dengue, yellow fever, West Nile, and Zika viruses 48,[65][66][67][68][69] . This is in agreement with the results from the GLMM analysis, in which NDVI and mean distance from roads were significantly associated with reduced species richness.
Our results also revealed the presence of many mosquito vector species in the remote areas of Miami-Dade. Anopheles quadrimaculatus, the primary vector of human malaria in the southeast United States and Cx. erraticus, the bridge vector of Eastern Equine Encephalitis in the southern United States, were almost exclusively found in the remote areas and were only collected in high numbers in the urban area U8, a new urban development bordering natural areas. Aedes albopictus is not commonly found in urban areas in Miami-Dade County 9,22 , and is relegated to restricted and well-defined areas such as cemeteries 70 . Therefore, it was not unexpected that approximately 90% of Ae. albopictus specimens were collected in the remote areas. However, despite the epidemiological importance of these species, their distant relationship with humans and inability of invading and thriving in urban areas of Miami-Dade greatly decreases their relevance in disease transmission when compared with mosquito vectors of anthroponotic vector-borne diseases. In this context, as our results revealed, Ae. aegypti is considerably more abundant in urban areas and is the primary vector of arboviruses that use humans as amplification hosts without the need for bridge vectors and primary hosts, such as dengue and Zika, representing a much greater public health threat 71 .
Urbanization has a major impact on the epidemiology of vector-borne disease transmission. It not only provides all the resources necessary for the survival of vector mosquito species such as Ae. aegypti and Cx. quinquefasciatus, but also provides shelter from the elements. Urban features such as tire shops or underground subway stations may allow populations of vector mosquito species to survive scorching and freezing temperatures that normally would kill them 72,73 . Furthermore, these highly productive urban environments that are responsible for the proliferation of vector mosquitoes are often located in populous areas, increasing, even more, the contact between humans and mosquito vectors 28,70 .
The impact of urbanization in the proliferation of vector mosquito species and arbovirus transmission must be considered under the Integrated Vector Management (IVM) framework 74 . Environmental ordinances and good practices, as well as simple modifications in the urban built environment at the early stages of development, can substantially decrease human exposure to mosquito vectors by attenuating social inequities and consequently social determinants of health 75,76 .
Even though several new technologies for controlling vector mosquito species in urban areas are being developed (e.g., genetically modified mosquitoes and Wolbachia-infect mosquitoes), they have not been validated to be used on large scales under real-world conditions 77,78 , and further entomological and epidemiological validation is still needed before they can be included and implemented under the IVM framework 79 . In this context, the World Health Organization (WHO) recommendation to control vector mosquito populations relies on the removal of aquatic habitats for immature mosquitoes and targeted insecticide application when needed. However, the current levels of proliferation of vector mosquito species in urban areas mediated by an overabundance of resources 3,9,22 make it virtually impossible to achieve the desired results of a safe mosquito abundance threshold to avoid arbovirus transmission 80 .
This study is not without limitations. We did not collect data across all weather and season variations that would have brought further insight into the natural variation in the mosquito community composition and abundance. We have collected mosquitoes using BG-Sentinel traps, which are the gold standard for collecting Aedes Stegomyia species. Even though BG-Sentinel traps have been proven effective to collect other less anthropophilic mosquito species and the fact that we have collected 26 mosquito species during this study, including many species that are notably not anthropophilic, we may have underestimated the presence and relative abundance of mosquito species that are less attracted by the BG-Sentinel traps.

Conclusion
The relationships between mosquito vectors, human hosts, and pathogens are driven by environmental conditions. Vector-borne diseases are conditioned by the environment, and anthropogenic changes have a direct influence on their epidemiology 7 . Climate change and urbanization increase the risk of arbovirus transmission by increasing the presence and abundance of mosquito vector species, therefore, increasing their contact with human populations. The findings of this study shed light on the effect of urbanization on the community composition of mosquitoes by reducing species richness and increasing the abundance of Ae. aegypti and Cx. quinquefasciatus in a non-random process of biotic homogenization. Large urban areas hold diverse socio-ecological conditions that can be highly conducive to both mosquito vector proliferation and local arbovirus transmission. Miami-Dade is www.nature.com/scientificreports/ one of the most critical entry points into the United States, with an elevated influx of people arriving and departing from endemic areas. Miami-Dade is, therefore, a sentinel or harbinger of what other cities in the contiguous United States will experience this century with climate change, population growth, regional trade, and human movement. The findings from this study highlight the importance of understanding how anthropogenic changes in the environment create an overabundance of resources that are responsible for sustaining the invasion, spread, and colonization of urban areas by vector mosquito species.

Methods
Study design. In this study, we identified and selected 24 areas: (i) 16 remote areas with high normalized difference vegetation index (NDVI) values. These areas were comprised of natural and rural areas with low population density or complete absence of humans, and no or minimum urban development allowed; and (ii) 8 urban areas with low NDVI values comprised of a university campus, and residential, and touristic areas with high human population density in Miami-Dade County, Florida (Fig. 4). We quantified the change in the NDVI obtained from Landsat satellite imagery mapped at 30 m spatial resolution. NDVI, which ranges from − 1 to + 1, provides a direct measurement of photosynthetic activity and is positively correlated with moisture availability, evapotranspiration, and vegetation biomass 81 . NDVI values less than or equal to ~ 0.1 are typically associated with urbanized surfaces (e.g., pavement, bare soils, and rooftops) or water bodies. We selected a cloud-free, atmospherically corrected Landsat from 2017 from https:// earth explo rer. usgs. gov/ during the month of November, which coincides with the approximate start of the dry season in South Florida 82,83 . To assess vegetation conditions around the traps, we utilized GIS software to create 1-km buffers around each trap and calculated the mean and standard deviation (SD) of NDVI within 1-km radii, which coincides approximately with the maximum flight distance of many species collected at the traps. The NDVI for each trap buffer provides a measure of spatial variation in vegetation state around each trap. For traps that were proximate to the shore (i.e., within 1 km of a shoreline), we masked out any water portions of the radius so that only NDVI values from land areas were included in NDVI. The masking process reduced the sample area from 3.15 km 2 for inland locations (i.e., greater than 1 km from shorelines) to a minimum of 0.76 km 2 for one island site with an extensive shoreline (Table 3). Then, we used the same 1-km buffers around each trap to calculate the number of pixels with buildings in each buffer around the trap sites to create a 2-D Building Footprint as one measure of urban development. We also calculated the mean distance from major roads in meters for each of the buffers (i.e., mean for all the pixels that comprise a buffer) as a proxy for urbanization. The layers are freely available at the Miami-Dade County's Open Data Hub-https:// gis-mdc. opend ata. arcgis. com/.

Collection of mosquitoes.
Mosquitoes were collected from August to October 2020 using BG-Sentinel traps (Biogents AG, Regensburg, Germany) baited with dry ice 84 . Mosquitoes were collected weekly at each collection site for 24 h for 5 consecutive weeks. Traps were placed in shaded areas that were protected from direct solar radiation, wind, and precipitation to enhance mosquito collections. The collected mosquitoes were transported to the Miami-Dade County Mosquito Control Laboratory and subsequently morphologically identified to species using taxonomic keys 85 .

Statistical analyses.
To compare the mosquito species composition in remote and urban areas as well as in areas with different land uses we performed a PERMANOVA with 9,999 permutations based on Bray-Curtis distances 86,87 . First, we subsetted the data into two groups to compare the mosquito species composition between rural and urban areas, and then we subsequently subsetted the data into 5 groups according to land use: natural, rural, university, residential, and touristic areas. Then we used the SIMPER method for assessing which species has contributed the most to the observed differences between groups of samples 88 . Analyses were done using PAST v3.2 89 .
We assessed spatial autocorrelation between samples using Moran's Index of Spatial Autocorrelation using ArcMap v10.5, in which failing to accept the alternative hypothesis indicates autocorrelation is not present in the data 90 . We checked for multicollinearity between the covariates NDVI, 2-D Building Footprint, and the mean distance from major roads using the variance inflation factor (VIF) values. VIF values below 10.00 indicate the assumptions were met and the covariates were not colinear 91,92 . There was no collinearity between the covariates, as the highest VIF value yielded 1.433 (Supplementary Table 3). We used the Eta squared (η 2 ) as the proportion of variance in the continuous target field explained by an effect to determine effect size 93 (Supplementary Table 4). Then, we performed a Generalized Linear Mixed Methods (GLMM) regression using a normal distribution considering species richness and relative abundance as the dependent variables and each urbanization measurement (Normalized Difference Vegetation Index (NDVI), 2-D Building Footprint, and the mean distance from major roads) and their interaction as independent variables, and week as a random effect with collection sites nested in time (weeks). Statistical analyses were done in SPSS V28.0.