- Open Access
Infection risk varies within urbanized landscapes: the case of coyotes and heartworm
Parasites & Vectors volume 14, Article number: 464 (2021)
Urbanization can have profound effects on ecological interactions. For host–pathogen interactions, differences have been detected between urban and non-urban landscapes. However, host–pathogen interactions may also differ within highly heterogeneous, urbanized landscapes.
We investigated differences in infection risk (i.e., probability of infection) within urbanized landscapes using the coyote (Canis latrans) and mosquito-borne nematode, Dirofilaria immitis (the causative agent for canine heartworm), as a case study. We focused on a coyote population in Chicago for which extensive behavioral and heartworm infection data has been collected between 2001 and 2016. Our objectives were to: (i) determine how onset and duration of the heartworm transmission season varied over the 16-year period and across the urban–suburban gradient; and (ii) investigate how heartworm infection risk in coyotes varied over the years, across the urban–suburban gradient, by coyote characteristics (e.g., age, sex, resident status), and coyote use of the urbanized landscape (e.g., use of urban areas, mosquito habitats).
While onset of the heartworm transmission season differed neither by year nor across the urban–suburban gradient, it was longer closer to the core of Chicago. Of the 315 coyotes sampled, 31.1% were infected with D. immitis. Older coyotes and coyotes sampled in later years (i.e., 2012–2016) were more likely to have heartworm. While coyote location in the urban–suburban gradient was not a significant predictor of infection, the proportion of urban land in coyote home ranges was. Importantly, the size and direction of this association varied by age class. For adults and pups, infection risk declined with urbanization, whereas for subadults it increased. Further, models had a higher predictive power when focusing on resident coyotes (and excluding transient coyotes). The proportion of mosquito habitat in coyote home ranges was not a significant predictor of infection.
Our findings suggest that urbanization may affect host exposure to vectors of D. immitis, that risk of infection can vary within urbanized landscapes, and that urbanization–wildlife infection associations may only be detected for animals with certain characteristics (e.g., age class and resident status).
Urbanization causes a shift in climatic conditions and landscape structure and composition . Temperatures tend to increase with urbanization due to pollution and impervious surfaces [2–4]. Vegetation becomes subdivided into patches surrounded by urban and suburban blocks. This shift in environmental context can have profound effects on processes unfolding in wildlife communities [5, 6]. For example, urban-induced fragmentation of the landscape and/or warming can alter animal behavior [7, 8] and species composition and abundance [9–12], which can in turn influence ecological relationships, such as predation and competition [13, 14]. Host–pathogen interactions can also be affected by urbanization, although it may vary with context, and in many cases, depends on pathogen transmission mode (reviewed in [15–19]).
Vector-borne pathogens are prone to be affected by urbanization because of the insect vector’s dependence on appropriate habitat, warm temperatures, and competent hosts [20, 21]. Whether we should expect an increase or decrease in vector-borne diseases with urbanization is a topic of intense debate, as opposing trends have been detected . For example, the prevalence of both ticks and avian malaria was found to be higher in rural common blackbirds (Turdus merula) than in urban ones . In contrast, West Nile virus seroprevalence was greater in urban than non-urban birds (e.g., [23–25]). Importantly, in addition to contrasting large-scale outcomes across landscapes (e.g., urban vs. rural), differences in infection risk (i.e., probability of infection) can also occur within urbanized landscapes [23, 26]. For instance, differences in the number and types of mosquito habitats (e.g., wetlands, artificial containers; ) across neighborhoods may lead to fine-scale differences in infection outcomes . While differences in infection risk within urbanized landscapes have been detected in vector species [26, 28], whether we should expect to observe similar patterns in urban host populations is less well understood.
The mosquito-borne nematode Dirofilaria immitis is the causative agent for canine heartworm, which is one of the most important parasitic diseases of domestic dogs in North America [29, 30]. Successful development and transmission of D. immitis is dependent on warm, humid conditions [31, 32], along with the presence of competent vectors and hosts [33, 34]. Warm, humid conditions are important for both the onset and duration of the heartworm transmission season and mosquito survival and reproduction. While over 60 mosquito species are susceptible to D. immitis, only nine act as competent vectors  and differ in their ability to adapt to urban settings . Competent hosts include domestic dogs and wild canids, in particular coyotes (Canis latrans). Coyotes in rural or natural areas are considered one of the primary reservoir hosts for D. immitis [36–38]. D. immitis prevalence can be as high as 37% in coyotes sampled in some rural areas [36, 39]. Despite an increasing presence of coyotes in many urban settings [40, 41], the distribution and prevalence of D. immitis in urban coyote populations is relatively unknown .
Here, we investigated how urbanization influences coyote risk of infection with D. immitis. To do this, we leveraged animal behavior and infection data from a well-studied coyote population sampled between 2001 and 2016 in the northwestern portion of the Chicago metropolitan area, which included both urban and suburban regions. Historically, coyotes were rare in the Chicago area, but they increased dramatically during the 1990s and are now common throughout the metropolitan area . While D. immitis prevalence in Chicago is unknown, the number of domestic dog cases reported by veterinary clinics in the Chicago area has increased by over four-fold in the past decade . Our objectives were to (i) determine how onset and duration of the heartworm transmission season varied over the 16-year period and across the urban–suburban gradient; and (ii) investigate how infection risk in coyotes varied over the years, across the urban–suburban gradient, by coyote characteristics (e.g., age, sex), and coyote use of the urbanized landscape (e.g., use of urban areas, mosquito habitats). Because the location and size of resident coyote home ranges vary less than those of transient coyotes [42, 44], we explored coyote use of the urbanized landscape for both resident and transient coyotes and for resident coyotes only.
The Chicago metropolitan area, with a human population of > 9 million people, extends across six counties (i.e., Cook, DuPage, Kane, Lake, McHenry, Will) in northeastern Illinois, USA (41.88° N, 87.63° W). Chicago has a temperate climate, with mean summer and winter temperatures ranging from 26 to 33 °C and from −1 to 3 °C, respectively, and rainfall averaging ~ 845 mm per year . Land cover in the region includes urban, suburban, natural areas, and agriculture. Landscapes in natural and urbanized areas include deciduous and coniferous forest, prairie, floodplain, wetland, open water, and managed green spaces (e.g., parks, greenways, golf courses).
The core of the metropolitan area, downtown Chicago, is situated on the edge of Lake Michigan. While proximity to a major water body can create a cooling effect and cause urban heat to shift westward , Chicago generally experiences only mild cooling that is most pronounced near the lakeshore during the summer  (although see ). This is in part due to Lake Michigan’s downwind location from the southwest winds as well as warm water temperatures in late summer . Further, any cooling effect of Lake Michigan in the core of Chicago is apparently counteracted by densely populated buildings, industrial zones, and train stations . Thus, Chicago’s heat island most likely occurs in the core of Chicago like traditional urban centers [1, 49], although low summer lake temperatures may push the heat island westerly.
Heartworm transmission season
Once infected, the temperature of the mosquito dictates the development of the microfilaria to the infective third stage [50–52]. Microfilarial development occurs above a threshold temperature of 14 °C, the progress of which can be tracked through the accumulation of heartworm development units (HDUs; ) such that 1 HDU is equal to 1 day with an average temperature 1 °C above 14 °C . Infective third-stage larvae will pass from the mosquito to the new host during blood meals only after enough HDUs have been accumulated throughout a given period . The period during which infective larvae are transmitted is called the heartworm transmission season and can be constructed for any region given sufficient climatological data [50–52]. The heartworm transmission season is said to have begun when 30-day HDUs surpass 130 °C and ends when 30-day HDUs drop below 130 °C .
To investigate the influence of urbanization on the onset and duration of the heartworm transmission season, we created four zones, each 15 km wide, to characterize a gradient from urban to suburban landscapes. Zone 1 had an average housing density of ~ 3250/km2, zone 2 of ~ 850/km2, zone 3 of ~ 530/km2, and zone 4 of ~ 400/km2 (Fig. 1) based on the 2010 SILVIS housing density data set (SILVIS Lab Spatial Analysis for Conservation and Sustainability). We obtained daily mean temperature data for the years 2000–2015 from the PRISM Climate Group (PRISM Climate Group, Oregon State University, http://prism.oregonstate.edu). Daily mean temperatures are available as spatial grids of 4-km spatial resolution which are calculated by interpolating climate data obtained from weather monitoring networks. For each zone, we used three 4-km grids that were evenly distributed across the zone (Additional file 1: Figure S1). For each 4-km grid in each year, daily HDUs were calculated by subtracting the threshold temperature of 14 °C from the daily mean temperature . Thirty-day HDUs were constructed by summing each daily HDU with the daily HDUs from the previous 29 days [52, 53,54]. The duration of each heartworm transmission season was determined by the number of months between the initiation and termination of the heartworm transmission season.
Coyotes were captured between February 2001 and December 2016 in the central and northwestern portion of the Chicago metropolitan area (Fig. 1). Most coyotes were captured in forest preserves, golf courses, small woodland parks, or in abandoned lots. Captures were performed opportunistically throughout the year, but primarily during winter and early spring. Coyotes were live-trapped with padded foothold traps and cable restraint devices . With the exception of 19 coyotes, capture locations were recorded using a handheld GPS. For the 19 coyotes without specific location coordinates, we used the coordinates of the center of the park where trapping occurred. All captured individuals were transported in a metal dog carrier to a research laboratory. Each coyote was sedated while in the dog carrier with 2.5 mg/kg Telazol (Fort Dodge Animal Health, Fort Dodge, IA, USA), which was administered intramuscularly in the hind limb. After immobilization, each coyote was sexed, and aged based on reproductive condition and tooth wear [ 56, 57]. Pups were 6–12 months, subadults were between 1 and 2 years, and adults were > 2 years. Pups less than 6 months were excluded since the prepatent period for heartworm is ~ 6 months . Approximately 3 ml of blood was collected from each coyote and poured into serum separator tubes. Tubes were left for ~ 30 min in an upright position and allowed to clot before being centrifuged for 15 min at 1790×g. Serum was extracted and stored in cryovial tubes at −80 °C. Serum was left in the freezer for periods ranging from 6 months to 5 years prior to further analyses. After blood collection, all coyotes were ear tagged and fitted with very-high-frequency radio collars (Advanced Telemetry Systems, Isanti, MN, USA). After recovering from immobilization, animals were released at the capture locations.
Serum samples were submitted to the Veterinary Diagnostic Laboratory at the University of Illinois Urbana-Champaign for heartworm screening. Presence of heartworm antigen was assessed using a membrane-bound ELISA test (SNAP® 4Dx® Plus Test, IDEXX Laboratories Inc.) following the manufacturer’s instructions. SNAP® 4Dx® Plus Test detects proteins produced within the reproductive tract of adult female heartworms and has a sensitivity of 97.5% (95% CI 94.26–99.18) and specificity of 94.0% (95% CI 83.45–98.75) .
Home range analysis and resident status classification
Relocations of all radio-collared coyotes were recorded 2–3 times a week during the day, and once a week during the night . Relocations were estimated using triangulation (with program LOCATE II; Pacer, Truro, Nova Scotia, Canada) with a truck-mounted antenna, or visual sightings. Relocations were used to estimate annual home ranges. We restricted annual home range estimates to individuals with a minimum of 30 relocations during at least 6 consecutive months. All relocations recorded beyond the 12-month period post capture were excluded. We assumed that the land cover types used by each coyote 6–12 months post capture were similar to those used when exposed to vectors of D. immitis.
We used two nonparametric methods to obtain home range estimates for each coyote: (1) we calculated and plotted 95% minimum convex polygons (MCPs); and (2) we used the adaptive local convex hull (a-LoCoH) method [59, 61]. For a-LoCoH, we calculated 95% contours and obtained the value of the adaptive sphere of influence “a” by calculating the maximum distance between two points . We used MCP and a-LoCoH over other home range estimators (e.g., the kernel density estimator) because MCP is most frequently used for very-high-frequency data, and a-LoCoH minimizes the extent to which home ranges cross hard boundaries (e.g., highways, rivers) [42, 60]. Since MCP can overestimate home range size  and a-LoCoH underestimate home range size (particularly if the sample size of locations is relatively small) , we explored heartworm–land cover associations using both methods. Because results were similar across methods (see “Results”), we presented MCP results in the main text and a-LoCoH results in the supplementary materials. All home range analyses were performed using the “adehabitatHR” package  in the statistical program R version 4.0.2 .
Home ranges were imported into ArcGIS version 10.3  and linked to land cover data to estimate the proportion of each land cover type present within each coyote’s home range. We used the 2011 National Land Cover Database (https://www.mrlc.gov) (spatial resolution: 30 m) to subdivide the landscape into different land cover types. Fourteen land cover types were present in the area. We combined eleven of these into two categories: (1) “mosquito habitat” (open water, woody wetlands, and emergent herbaceous wetlands); and (2) “green spaces” (open developed space, mixed forest, evergreen forest, deciduous forest, cultivated crops, pasture/hay, grassland/herbaceous, shrub/scrub). The other three were urban land cover types (high, medium, and low developed urban land) and were examined as independent variables since we were interested in the degree of urbanization. High, medium, and low developed urban land indicates that 80–100%, 50–79%, and 20–49% of the land is impervious surface, respectively. Because all three urban variables were highly correlated with the “green spaces” variable and the “high developed urban land” variable with the “medium developed urban land” variable (r2 > 0.5), we excluded the “green spaces” and "high developed urban land” variables from statistical analyses and focused on “mosquito habitat,” “developed medium urban land,” and “developed low urban land” variables.
We used multiple characteristics to discriminate resident and transient coyotes. Residents repeatedly used an explicit territory across two or more seasons, and transients shifted use areas across seasons and had larger home ranges that overlapped multiple territories [42, 65, 66]. Further, residents were often seen traveling with other coyotes, whereas transients did not, or residents shared the same territories with reproductive pairs occupying a territory.
We ran a linear regression to determine how the duration of the heartworm transmission season varied by year and urban zone (Table 1). The outcome variable in this model was duration of the heartworm transmission season (in months), and predictor variables were year (2001–2015) and urban zone (1–4). We also included latitude to account for any variation associated with collecting temperature data from three 4-km grids located at different latitudes within each zone (Additional file 1: Figure S1). Variation in the onset of the heartworm transmission season with year, urban zone, and latitude was not explored because there was little variation (i.e., onset of the heartworm transmission season occurred in June 97.9% of the time).
We ran a generalized linear model (GLM) to investigate how heartworm infection risk in coyotes varied across years (2001–2016), the urban–suburban gradient (urban zone 1–4), with coyote characteristics (e.g., age, sex), and coyote use of the urbanized landscape (i.e., mosquito habitat, medium developed urban land, and low developed urban land). However, we split the analysis into two because of differences in sample size for some of the fixed effects. One analysis included all coyotes tested for heartworm (n = 315 coyotes with 16 tested more than once) and the other included coyotes for which enough relocation data were obtained to estimate annual home ranges (n = 146 coyotes, a subset of the 315 coyotes). The outcome variable was heartworm infection (yes/no), thus we ran binomial GLMs with logit link functions using the “lme4” R package . Additionally, because we expected infection to increase with age [38, 39, 68], we also included the proportion of adults tested each year as an offset in both analyses (Table 1).
In the analysis that included 315 coyotes, urban zone, age class, and sex were included as categorical fixed effects. Year was included as a continuous fixed effect, and nonlinear relationships with heartworm infection were examined using basis splines (Table 1) using the “splines” R package. Since 16 coyotes were tested more than once (Additional file 1: Table S1), we include “animal ID” as a random intercept and ran a generalized linear mixed model (GLMM) instead of a GLM using the “lme4” package. Further, we grouped observations into a “site” random intercept to account for any significant spatial autocorrelation in model residuals (Moran’s I statistic after including site as a random effect: z = 0.52, P = 0.3). Coyotes were included in the site that was closest to their capture location except if there were man-made barriers (e.g., highways).
In the analysis that included 146 coyotes, we ran four models: (1) for residents and transients using MCP; (2) for residents only using MCP; (3) for residents and transients using a-LoCoH; and (4) for residents only using a-LoCoH. Fixed effects included the proportion of mosquito habitat, medium developed urban land, and low developed urban land in coyote home ranges. Age class and year were included as fixed effects because they were significant predictors in the analysis with 315 coyotes (see “Results”). Resident status was also included in models that examined both residents and transients. Since we expected age class to be an important predictor of infection, we also evaluated whether the association between infection and the proportion of mosquito habitat, medium developed urban land, and low developed urban land in coyote home ranges varied by age class by including interactions between these variables. We included site as a random intercept to account for any significant spatial autocorrelation (Moran’s I statistic after including site as a random effect for residents and transients using MCP: z = 0.08, P = 0.47; for residents only using MCP: z = −0.24, P = 0.59; for residents and transients using a-LoCoH: z = −0.02, P = 0.51; for residents only using a-LoCoH: z = −0.36, P = 0.64). Animal ID was not included as a random intercept because none of the coyotes in this second analysis were resampled more than once.
For all models, the most parsimonious model was identified using an information theory approach, comparing models with different variable combinations, and used the Akaike information criterion corrected for small sample size (AICc) to rank models [69, 70] using the “MuMIn” R package . If at least one model was within 2 ΔAICc values of the top-ranking model, model averaging was used to obtain mean effect sizes and 95% confidence intervals . All continuous predictors were centered and standardized to facilitate interpretation of main effects and to perform model averaging . Multicollinearity among continuous predictors was assessed using the variance inflation factor . Scaled residuals of each model were examined for uniformity using the “DHARMa” package . Model fit was assessed by calculating the marginal and conditional coefficients of determination (rm2 and rc2, respectively) . rm2 is the variance explained by the fixed effects, and rc2 the variance explained by the fixed and random effects .
Onset and duration of the heartworm transmission season
Across all urban zones (1–4) and years (2000–2015), the heartworm transmission season most often began in June (97.9% of the time) and lasted for a period of 2–5 months (mean of 3.56 months). For the duration of the heartworm transmission season, urban zone appeared in all three of the top-ranking models, and year and latitude in one (Additional file 1: Table S2). Model averaging of the top three models showed that urban zone was a significant predictor of infection (Table 2), while year and latitude were not (Table 2). The heartworm transmission season was significantly longer in zone 1 compared to zone 3 and 4 (Fig. 2). Pairwise comparisons revealed that the same was true for zone 2 compared to zone 3 and 4 (zone 2 vs. zone 3: z = 2.88, P = 0.02; zone 2 vs. zone 4: z = 2.89, P = 0.02).
Urbanization and coyote infection risk
Three hundred and fifteen coyotes were captured and tested for heartworm between 2001 and 2016 (16 were captured more than once; Additional file 1: Table S1). The number of animals captured and tested each year ranged from 5 in 2001 to 51 in 2014 (mean = 22.5 per year). Heartworm tests were performed on 94 pups (52 females and 42 males), 108 subadults (53 females and 55 males), and 113 adults (38 females and 75 males). Ninety-eight coyotes were positive for heartworm (31.1%). Prevalence ranged from 7.7% in 2011 (n = 13) to 66.7% in 2016 (n = 21) (Fig. 3a).
When examining infection risk for all captured coyotes (n = 315), the best fit model contained only age class and year (Additional file 1: Table S3). Urban zone and sex were not important predictors of infection because they did not appear in the top-ranking model (Additional file 1: Table S3). A quadratic relation better explained the relationship between heartworm infection and year than a linear relationship (Table 3). Infection risk was lowest in 2008–2011 and increased in 2012–2016 (Fig. 3a). For age class, adults had a higher risk of infection than pups and subadults (Table 3; Fig. 3b).
Of the 315 coyotes tested for heartworm, 146 had enough relocations to estimate annual home ranges (mean number of relocations per animal = 163, range = 45–585). The 146 individuals comprised 37 pups (21 females and 16 males), 46 subadults (23 females and 23 males), and 63 adults (25 females and 38 males). In terms of resident status, this amounted to 107 residents and 39 transients. The years with the lowest number of coyotes tracked were 2005 and 2010 (n = 3), and the years with the highest number of coyotes tracked were 2012 and 2013 (n = 22) and 2014 (n = 23). Association between heartworm infection and age class, year, resident status, and coyote use of the urbanized landscape tended to be similar across all four models (i.e., resident and transient coyotes using MCP, resident coyotes only using MCP, resident and transient coyotes using a-LoCoH, and resident coyotes only using a-LoCoH). However, models with the greatest predictive power used MCP instead of a-LoCoH and focused on resident coyotes only. For top-ranking models using MCP (ΔAICc < 2), the largest r2 value was 0.68 for the resident coyotes only analysis and 0.49 for the resident and transient coyote analysis. For top-ranking models using a-LoCoH, the largest r2 value was 0.47 for the resident coyotes-only analysis and 0.5 for the resident and transient coyote analysis (Additional file 1: Table S3). Results using MCP are summarized in Table 4, Fig. 4, and Additional file 1: Table S4, and results using a-LoCoH are summarized in Additional file 1: Tables S4, S5, and Figure S2.
Age class, proportion of medium developed urban land in coyote home ranges, and the interaction between these two predictors appeared in all top-ranking models (Additional file 1: Table S4, except for resident coyotes only using MCP in which two out of three models had the two predictors and interaction) and thus were the most important predictors of heartworm infection. Further, for the resident and transient models, proportion of mosquito habitat and year were second most important, followed by proportion of low developed urban land in coyote home ranges and resident status (Additional file 1: Table S4). For the resident only models, year appeared in none of the top-ranking models (Additional file 1: Table S4). The interaction between age class and proportion of low developed urban land, and mosquito habitat in coyote home ranges were the least important as they did not appear in any of the top-ranking models (Additional file 1: Table S4).
Model averaging revealed that when resident and transient coyotes were examined, pups had a lower infection risk than adults (Table 4 and Additional file 1: Table S5). Further, infection risk tended to decline as the proportion of medium developed urban land in coyote home ranges increased. However, this association varied by age class, where infection risk declined as the proportion of medium developed urban land in home ranges increased for adults but increased for subadults (Table 4, Fig. 4, and Additional file 1: Table S5 and Fig. S2). Further, changing the model reference level to pups revealed that, for the models that included both resident and transient coyotes, the relationship between infection and proportion of medium developed urban land in home ranges was significantly different between pups and subadults (P < 0.05 for all four model types), where pups, like adults, had a lower risk of infection with more medium developed urban land in their home ranges (Fig. 4 and Additional file 1: Table S2). No significant difference was detected between pups and subadults when focusing on resident coyotes only.
Urbanization can have contrasting effects on host–pathogen interactions [15, 18 19]. Here, we found that urbanization influenced the duration of the heartworm transmission season and infection risk in coyotes. The heartworm transmission season was longer closer to the core of Chicago. Heartworm prevalence in coyotes increased during the study period and with coyote age. Further, the proportion of medium developed urban land in coyote home ranges was an important predictor of infection, but direction and size of the effect varied by age class and models had a higher predictive power when examining resident coyotes only. For adults and pups, infection risk declined with urbanization, whereas for subadults, it increased.
The Chicago coyote population had a heartworm prevalence of 31.1%. Previous studies performed in Madison, Wisconsin and Tucson, Arizona have found heartworm prevalence in urban coyotes to be 35.7% (n = 14)  and 0% (n = 22) , respectively. Importantly, we found that heartworm prevalence fluctuated yearly, with prevalence being as low as 8% in some years and as high as 67% in others, suggesting that there can be notable differences in prevalence across years. Further, heartworm prevalence tended to increase over the 16-year period, a finding that is in line with national trends observed in domestic dogs suggesting that heartworm prevalence is increasing over time across the USA . For northern US states like Illinois, this increase may be associated with an increase in the number and density of mosquito vectors [78, 79], possibly due to the combined effect of shifting climate conditions and few mosquito abatement programs . For Chicago, an increase in coyote numbers over the years  could also be an important factor.
The proportion of coyotes sampled closer to the core of Chicago also increased over the years, which may have contributed to an increase in heartworm prevalence over time. Twenty-six coyotes were sampled in the urban zone 1, of which 24 (92%) were sampled in 2013–2016. While urban zone was not a significant predictor of infection, coyote proximity to the core of Chicago may play a role because the heartworm transmission season tended to be longer closer to the core of Chicago. Additionally, heartworm prevalence in rural coyotes in Illinois is 16% , suggesting that urban coyotes might be at a higher risk of infection than non-urban coyotes. Further, mosquitoes sampled in urban areas can have a higher heartworm prevalence than rural mosquitoes . This could be because one of the main vectors of D. immitis, Aedes albopictus, tends to thrive in urbanized landscapes owing to warmer conditions and the presence of natural and artificial water bodies and containers [35, 80, 81]. A nonsignificant effect of urban zone may be due to a smaller sample size closer to the core Chicago (i.e., only 29 coyotes were tested in zone 1 and 2 combined: 26 in zone 1 and three in zone 2).
Another potential reason for not detecting a significant effect of urban zone on coyote infection could be that measuring proximity of coyotes to the urban core simplifies or underestimates complex patterns occurring within urban patches. For example, mosquito abundance and richness, as well as infection, can vary across distances as small as neighborhoods , and the degree of landscape heterogeneity can influence mosquito diversity [82, 10]. Exploring land cover composition of coyote home ranges provided greater insight for associations with infection risk than proximity to the core of Chicago. The proportion of urban land in coyote home ranges was an important predictor of infection, but only when quantified as medium developed urban land and not low developed urban land. Since impervious surfaces account for 50–79% of total land cover for medium developed urban land and 20–49% for low developed urban land, impervious surfaces and built-up land may explain the observed association.
It was surprising that the proportion of mosquito habitats in coyote home ranges was a nonsignificant predictor of infection risk. In the case of domestic dogs, proximity to mosquito-bearing waters can be an important predictor of infection with D. immitis . One reason for not detecting an association in our system could be that there were other unaccounted for water bodies (e.g., man-made or temporary water bodies such as artificial containers, puddles, tires, trash cans) [84, 85]. Exploring whether the presence of man-made and temporary water bodies versus permanent/vegetated water bodies in coyote home ranges influences heartworm infection would be an important next step to take.
It is noteworthy that the relationship between infection risk and the proportion of medium developed urban land in coyote home ranges varied both in directionality and size by age class. The negative association detected for adults and pups is likely associated with built-up areas having lower mosquito abundance and mosquito species richness than green spaces (e.g., parks, forest preserves; [86, 87]). There are nine competent vectors of D. immitis , most of which breed in wetland and woodland areas . For example, Ae. vexans, a floodwater mosquito, is most frequently found in riparian zones, roadside ditches, and wetlands [35, 89], and is therefore perhaps more commonly found in urban green spaces. Adults and pups may also be at a lower risk of infection with more medium developed urban land in their home range because there is perhaps more mosquito control than in green spaces. The small effect size detected for pups is probably because pups had a lower prevalence than adults. Indeed, heartworm infection risk tended to increase with age, a pattern that is consistent with previous work [38, 39, 68]. The positive association detected for subadults could be because most subadults are transitioning between their natal and new territories. During this dispersal period, subadults may be exposed to a broader range of microhabitats than adults and pups, and thus could have more opportunities to encounter environments with more mosquitoes. An important next step that would help disentangle the importance of these various potential explanations would be to explore whether there is a relationship between infection risk and habitat use within home ranges.
The fact that the directionality of the infection–urbanization association by age class remained the same regardless of the home range estimator used (MCP vs. a-LoCoH) highlights the strength of these associations. However, it is interesting that model predictive power increased when focusing on resident coyotes and excluding transients. Because transients tended to have larger and more complex home ranges than residents, we suspect that the infection–urbanization association might differ, or might be less apparent, if a larger number of transients were included in the analysis (and/or that transients were examined separately). More research is needed to determine whether transient coyotes likely would have the same infection–urbanization association as resident adult coyotes or as subadult coyotes.
We cannot infer that adult coyotes, particularly transients, were likely infected in the areas where they were sampled; however, we can for pups, as coyotes only tend to leave their natal territories as subadults . The fact that the direction of the infection–urbanization association was the same in pups and adults suggests that coyotes from these two age classes may have been infected in the area (or similar environment) where they were sampled. This lends support for the notion that urban wildlife reflect their local environment [25, 91] and highlights the need to carefully consider which types of individuals should be examined to effectively capture any associations with local environments (e.g., pups and resident adults in our case study).
While this study provides insight on how urbanization might influence wildlife infection risk, there were a number of limitations. Firstly, coyote sampling varied across years and urban zones, which limited result interpretations in some cases. That said, to the best of our knowledge, this is one of the first + 15-year urban wildlife disease investigations and provides unique evidence that wildlife disease risk can vary over time in urbanized settings. Another broader limitation worth noting was the inability to account for the time lag between human-derived changes to the landscape, and vector and wildlife response to this change [16, 33]. This could be especially important for urban and suburban land covers, which tend to increase over time. Future work should focus on developing approaches that can approximate time lags between land-use change and host and pathogen response to such changes . Another important limitation of this study was not accounting for socioeconomic factors. Recent studies suggest that a number of socioeconomic factors can influence the distribution of wildlife diseases in urbanized areas (e.g., household income) [92, 93]. Future vector-borne and wildlife disease studies should quantify socioeconomic factors alongside of structural and abiotic components when exploring effects of urbanization on infection risk .
Measuring the effects of urbanization on host–pathogen interactions is becoming an important area of research [15, 19], particularly as urbanized areas continue to expand . Recent work has found differences in disease risk between urban and non-urban wildlife populations (reviewed in ). Here, we found that coyote infection with the vector-borne pathogen D. immitis can vary within urban and suburban areas, and that effects may only be detected for certain age classes, and when using certain metrics of urbanization. While we were not able to make comparisons with rural or wildland coyotes, the fact that we detected differences in infection risk among coyotes residing in different urban and suburban areas highlights the complex way by which vector-borne diseases are transmitted in urbanized landscapes.
Availability of data and materials
Data are available on figshare (https://doi.org/10.6084/m9.figshare.15169146.v1).
Heartworm development unit
Minimum convex polygon
Adaptive local convex hull
Generalized linear model
Generalized linear mixed model
Akaike information criterion corrected for small sample size
Berry BJL. Urbanization. In: Urban ecology. Boston: Springer; 2008.
Arnfield AJ. Two decades of urban climate research: a review of turbulence, exchanges of energy and water, and the urban heat island. Int J Climatol. 2003;23:1–26.
Mishra V, Ganguly AR, Nijssen B, Lettenmaier DP. Changes in observed climate extremes in global urban areas. Environ Res Lett. 2015;10:024005.
Oke TR. The energetic basis of the urban heat island. Q J R Meteorol Soc. 1982;108:1–24.
Johnson MTJ, Munshi-South J. Evolution of life in urban environments. Science. 2017. https://doi.org/10.1126/science.aam8327.
McKinney ML. Urbanization, biodiversity, and conservation. Bioscience. 2002;52:883–90.
Ditchkoff SS, Saalfeld ST, Gibson CJ. Animal behavior in urban ecosystems: modifications due to human-induced stress. Urban Ecosyst. 2006;9:5–12.
Lowry H, Lill A, Wong BBM. Behavioural responses of wildlife to urban environments. Biol Rev. 2013;88:537–49.
Grimm NB, Faeth SH, Golubiewski NE, Redman CL, Wu J, Bai X, et al. Global change and the ecology of cities. Science. 2008;319:756–60.
Chaves LF, Hamer GL, Walker ED, Brown WM, Ruiz MO, Kitron UD. Climatic variability and landscape heterogeneity impact urban mosquito diversity and vector abundance and infection. Ecosphere. 2011;2(6):1–21.
Aronson MFJ, La SFA, Nilon CH, Katti M, Goddard MA, Lepczyk CA, et al. A global analysis of the impacts of urbanization on bird and plant diversity reveals key anthropogenic drivers. Proc R Soc B Biol Sci. 2014;281:20133330.
Battles AC, Kolbe JJ. Miami heat: urban heat islands influence the thermal suitability of habitats for ectotherms. Glob Chang Biol. 2019;25:562–76.
Faeth SH, Warren PS, Shochat E, Marussich WA. Trophic dynamics in urban communities. Bioscience. 2005;55:399–407.
Fischer JD, Cleeton SH, Lyons TP, Miller JR. Urbanization and the predation paradox: the role of trophic dynamics in structuring vertebrate communities. Bioscience. 2012;62:809–18.
Bradley CA, Altizer S. Urbanization and the ecology of wildlife diseases. Trends Ecol Evol. 2007;22:95–102.
Brearley G, Rhodes J, Bradley A, Baxter G, Seabrook L, Lunney D, et al. Wildlife disease prevalence in human-modified landscapes. Biol Rev. 2013;88:427–42.
Mackenstedt U, Jenkins D, Romig T. The role of wildlife in the transmission of parasitic zoonoses in peri-urban and urban areas. Int J Parasitol Parasites Wildl. 2015;4:71–9.
Murray MH, Sánchez CA, Becker DJ, Byers KA, Worsley-Tonks KEL, Craft ME. City sicker? A meta-analysis of wildlife health and urbanization. Front Ecol Environ. 2019;17:575–83.
Hassell JM, Begon M, Ward MJ, Fèvre EM. Urbanization and disease emergence: dynamics at the wildlife–livestock–human interface. Trends Ecol Evol. 2017;32:55–67.
Friggens MM, Beier P. Anthropogenic disturbance and the risk of flea-borne disease transmission. Oecologia. 2010;164:809–20.
Vora N. Impact of anthropogenic environmental alterations on vector-borne diseases. Medscape J Med. 2008;10:238.
Evans KL, Gaston KJ, Sharp SP, McGowan A, Simeoni M, Hatchwell BJ. Effects of urbanisation on disease prevalence and age structure in blackbird Turdus merula populations. Oikos. 2009;118:774–82.
Bradley CA, Gibbs SEJ, Altizer S. Urban land use predicts west Nile virus exposure in songbirds. Ecol Appl. 2008;18:1083–92.
Gibbs SEJ, Wimberly MC, Madden M, Masour J, Yabsley MJ, Stallknecht DE. Factors affecting the geographic distribution of West Nile Virus in Georgia, USA: 2002–2004. Vector-Borne Zoonotic Dis. 2006;6:73–82.
Hamer SA, Lehrer E, Magle SB. Wild birds as sentinels for multiple zoonotic pathogens along an urban to rural gradient in Greater Chicago, Illinois. Zoonoses Public Health. 2012;59:355–64.
LaDeau SL, Allan BF, Leisnham PT, Levy MZ. The ecological foundations of transmission potential and vector-borne disease in urban landscapes. Funct Ecol. 2015;29:889–901.
Wilke ABB, Chase C, Vasquez C, Carvajal A, Medina J, Petrie WD, et al. Urbanization creates diverse aquatic habitats for immature mosquitoes in urban areas. Sci Rep. 2019;9:1–11.
Reisen WK. Landscape epidemiology of vector-borne diseases. Annu Rev Entomol. 2010;55:461–83.
Bowman DD, Atkins CE. Heartworm biology, treatment, and control. Vet Clin North Am Small Anim Pract. 2009;39:1127–58.
Self SW, Pulaski CN, Mcmahan CS, Brown DA, Yabsley MJ, Gettings JR. Regional and local temporal trends in the prevalence of canine heartworm infection in the contiguous United States: 2012–2018. Parasites Vectors. 2019;12:380.
Genchi C, Rinaldi L, Mortarino M, Genchi M, Cringoli G. Climate and Dirofilaria infection in Europe. Vet Parasitol. 2009;163:286–92.
Ledesma N, Harrington L. Fine-scale temperature fluctuation and modulation of (Dirofilaria immitis) larval development in (Aedes aegypti). Vet Parasitol. 2015;209:93–100.
Brown HE, Harrington LC, Kaufman PE, Mckay T, Bowman DD, Nelson CT, et al. Key factors influencing canine heartworm, (Dirofilaria immitis), in the United States. Parasites Vectors. 2012;31:1–9.
Wang D, Bowman DD, Brown HE, Harrington LC, Kaufman PE, McKay T, et al. Factors influencing US canine heartworm (Dirofilaria immitis) prevalence. Parasites Vectors. 2014;7:1–18.
Ledesma N, Harrington L. Mosquito vectors of dog heartworm in the United States: vector status and factors influencing transmission efficiency. Top Companion Anim Med. 2011;26:178–85.
Weinmann CJ, Garcia R. Coyotes and canine heartworm in California. J Wildl Dis. 1980;16:217–22.
Sacks BN. Increasing prevalence of canine heartworm in coyotes from California. J Wildl Dis. 1998;34:386–9.
Nelson TA, Gregory DG, Laursen JR. Canine heartworms in coyotes in Illinois. J Wildl Dis. 2003;39:593–9.
Aher A, Caudill D, Caudill G, Butryn RS, Wolf D, Fox M, et al. Prevalence, genetic analyses, and risk factors associated with heartworm (Dirofilaria immitis) in wild coyotes (Canis latrans) from Florida, USA. J Wildl Dis. 2016;52:785–92.
Bateman PW, Fleming PA. Big city life: carnivores in urban environments. J Zool. 2012;287:1–23.
Gehrt SD, Riley SPD. Coyotes (Canis latrans). In: Urban carnivores: ecology, conflict, and conservation. Baltimore: The Johns Hopkins University Press; 2010. p. 79–95.
Gehrt SD, Anchor C, White LA. Home range and landscape use of coyotes in a metropolitan landscape: conflict or coexistence? J Mammal. 2009;90:1045–57.
American Heartworm Society. Incidence maps. 2021. https://www.heartwormsociety.org/pet-owner-resources/incidence-maps.
Hinton JW, Chamberlain MJ, Van Manen FT. Long-distance movements of transient coyotes in Eastern North Carolina. Am Midl Nat. 2012;168:281–8.
United States Climate Data. Climate Chicago—Illinois, your weather service-world climate. 2020. https://www.usclimatedata.com/climate/chicago/illinois/united-states/usil0225.
Oswald EM, Rood RB, Zhang K, Gronlund CJ, O’neill MS, White-Newsome JL, et al. An investigation into the spatial variability of near-surface air temperatures in the Detroit, Michigan, metropolitan region. J Appl Meteorol Climatol. 2012;51:1290–304.
Coseo P, Larsen L. How factors of land use/land cover, building configuration, and adjacent heat sources and sinks explain Urban Heat Islands in Chicago. Landsc Urban Plan. 2014;125:117–29.
Ackerman B. Temporal march of the Chicago heat island. J Clim Appl Meteorol. 1985;24:547–54.
Gala TS, Alfraihat R, Mulugeta G, Gala TS. Ecological evaluation of urban heat island in Chicago city, USA hydrological application of remote sensing view project ecological evaluation of urban heat island in Chicago city, USA. J Atmos Pollut. 2016;4:23–9.
Kutz FW, Dobson RC. Effects of temperature on the development of Dirofilaria immitis (Leidy) in Anopheles quadrimaculatus and on vector mortality resulting from this development 2,3. Ann Entomol Soc Am. 1974;67:325–31.
Christensen BM, Hollander AL. Effect of temperature on vector–parasite relationships of Aedes trivittatus and Dirofilaria immitis. Proc Helminthol Soc Wash. 1978;45:115–9.
Fortin JF, Slocombe JOD. Temperature requirements for the development of Dirofilaria immitis in Aedes triseriatus and Ae. vexans. Mosq News. 1981;41:625–33.
Slocombe JOD, Surgeoner GA, Srivastava B. Determination of heartworm transmission period and its use in diagnosis and control. In: Proc Hear Symp. 1990. p. 19–26.
Sacks BN, Woodward DL, Colwell AE. A long-term study of non-native-heartworm transmission among coyotes in a Mediterranean ecosystem. Oikos. 2003;102:478–90.
Gehrt SD, Brown JL, Anchor C. Is the urban coyote a misantrhopic synanthrope? The case from Chicago. Cities Environ. 2011;4:3.
Nellis CH, Wetmore SP, Keith LB. Age-related characteristics of coyote canines. J Wildl Manag. 1978;42:680–3.
Scrivner JH, Johnson CA, Sego CA. Use of cementum annuli and eye-lens weight for aging coyotes. Wildl Soc Bull. 2014;38:874–7.
McCall JW, Genchi C, Kramer LH, Guerrero J, Venco L. Chapter 4 Heartworm disease in animals and humans. In: Advances in parasitology. London: Academic Press; 2008. p. 193–285.
Henry LG, Brunson KJ, Walden HS, Wenzlow N, Beachboard SE, Barr KL, et al. Comparison of six commercial antigen kits for detection of (Dirofilaria immitis) infections in canines with necropsy-confirmed heartworm status. Vet Parasitol. 2018;254:178–82.
Getz WM, Wilmers CC. A local nearest-neighbor convex-hull construction of home ranges and utilization distributions. Ecography. 2004;27:489–505.
Getz WM, Fortmann-Roe S, Cross PC, Lyons AJ, Ryan SJ, Wilmers CC. LoCoH: nonparameteric kernel methods for constructing home ranges and utilization distributions. PLoS ONE. 2007;2:e207.
Calenge C. Home range estimation in R: the adehabitatHR package. 2015.
R Development Core Team. R: a language and environment for statistical computing. Vienna, Austria; 2020.
ESRI ESRI. ArcGIS professional GIS for the desktop, version 10.3.1, United States. 2015.
Gese EM, Morey PS, Gehrt SD. Influence of the urban matrix on space use of coyotes in the Chicago metropolitan area. J Ethol. 2012;30(3):413–25.
Poessel SA, Breck SW, Gese EM. Spatial ecology of coyotes in the Denver metropolitan area: influence of the urban matrix. J Mammal. 2016;97:1414–27.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. arXiv Prepr arXiv. 2014;1406.5823.
Chitwood MC, Swingen MB, Lashley MA, Flowers JR, Palamar MB, Apperson CS, et al. Parasitology and serology of free-ranging coyotes (Canis latrans) in North Carolina, USA. J Wildl Dis. 2015;51:664–9.
Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer; 2002.
Johnson JB, Omland KS. Model selection in ecology and evolution. Trends Ecol Evol. 2004;19:101–8.
Barton K. Package ‘MuMIn’. R package version 1. 2013.
Grueber CE, Nakagawa S, Laws RJ, Jamieson IG. Multimodel inference in ecology and evolution: challenges and solutions. J Evol Biol. 2011;24:699–711.
Fox J, Weisberg S. Car: an R companion to applied regression, 2nd edition. R package. 2011.
Hartig F. Package ‘DHARMa’ residual diagnostics for hierarchical (multi-level/mixed) regression models. 2017.
Nakagawa S, Schielzeth H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol Evol. 2013;4:133–42.
Pluemer M, Dubay S, Drake D, Crimmins S, Veverka T, Hovanec H, et al. Red foxes (Vulpes vulpes) and coyotes (Canis latrans) in an urban landscape: prevalence and risk factors for disease. J Urban Ecol. 2019;5:1–9.
Grinder M, Krausman PR. Morbidity-mortality factors and survival of an urban coyote population in Arizona. J Wildl Dis. 2001;37:312–7.
Hahn MB, Eisen L, McAllister J, Savage HM, Mutebi JP, Eisen RJ. Updated reported distribution of Aedes (Stegomyia) aegypti and Aedes (Stegomyia) albopictus (Diptera: Culicidae) in the United States, 1995–2016. J Med Entomol. 2017;54:1420–4.
Proestos Y, Christophides GK, Ergüler K, Tanarhte M, Waldock J, Lelieveld J. Present and future projections of habitat suitability of the Asian tiger mosquito, a vector of viral pathogens, from global climate simulation. Philos Trans R Soc B Biol Sci. 2015;370:1–16.
Paras KL, O’Brien VA, Reiskind MH. Comparison of the vector potential of different mosquito species for the transmission of heartworm, (Dirofilaria immitis), in rural and urban areas in and surrounding Stillwater, Oklahoma, USA. Med Vet Entomol. 2014;28:60–7.
Li Y, Kamara F, Zhou G, Puthiyakunnon S, Li C, Liu Y, et al. Urbanization increases (Aedes albopictus) larval habitats and accelerates mosquito development and survivorship. PLoS Negl Trop Dis. 2014;8:e3301.
Ruiz MO, Chaves LF, Hamer GL, Sun T, Brown WM, Walker ED, et al. Local impact of temperature and precipitation on West Nile virus infection in Culex species mosquitoes in northeast Illinois, USA. Parasites Vectors. 2010;3:19.
Trájer A, Rengei A, Farkas-Iványi K, Bede-Fazekas Á. Impacts of urbanisation level and distance from potential natural mosquito breeding habitats on the abundance of canine Dirofilariosis. Acta Vet Hung. 2016;64:340–59.
Crocker W, Maute K, Webb C, French K. Mosquito assemblages associated with urban water bodies; implications for pest and public health threats. Landsc Urban Plan. 2017;162:115–25.
Bartlett-Healy K, Unlu I, Obenauer P, Hughes T, Healy S, Crepeau T, et al. Larval mosquito habitat utilization and community dynamics of Aedes albopictus and Aedes japonicus (Diptera: Culicidae). J Med Entomol. 2012;49:813–24.
Ferraguti M, Martínez-De La Puente J, Roiz D, Ruiz S, Soriguer R, Figuerola J. Effects of landscape anthropization on mosquito community composition and abundance. Sci Rep. 2016;6:1–9.
Medeiros-Sousa AR, Fernandes A, Ceretti-Junior W, Barreto A, Wilke B, Toledo MM. Mosquitoes in urban green spaces: using an island biogeographic approach to identify drivers of species richness and composition. Sci Rep. 2017;7:17826.
Rey J, Walton W, Wolfe R, Connelly C, O’Connell S, Berg J, et al. North American wetlands and mosquito control. Int J Environ Res Public Health. 2012;9:4537–605.
Goddard J, Varnado WC, Harrison BA. An annotated list of the mosquitoes (Diptera: Culicidae) of Mississippi. J Vector Ecol. 2010;35:79–88.
Gese EM, Rongstad OJ, Mytton WR. Home range and habitat use of coyotes in southeastern Colorado. J Wildl Manag. 1988;52:640–6.
Frantz A, Pottier MA, Karimi B, Corbel H, Aubry E, Haussy C, et al. Contrasting levels of heavy metals in the feathers of urban pigeons from close habitats suggest limited movements at a restricted scale. Environ Pollut. 2012;168:23–8.
Ayral F, Artois J, Zilber AL, Widén F, Pounder KC, Aubert D, et al. The relationship between socioeconomic indices and potentially zoonotic pathogens carried by wild Norway rats: a survey in Rhône, France (2010–2012). Epidemiol Infect. 2015;143:586–99.
Murray MH, Fidino M, Fyffe R, Byers KA, Pettengill JB, Sondgeroth KS, et al. City sanitation and socioeconomics predict rat zoonotic infection across diverse neighbourhoods. Zoonoses Public Health. 2020;67:673–83.
Moll RJ, Cepek JD, Lorch PD, Dennis PM, Tans E, et al. What does urbanization actually mean? A framework for urban metrics in wildlife research. J Appl Ecol. 2019;56:1289.
Seto KC, Fragkias M, GüNeralp B, Reilly MK. A meta-analysis of global urban land expansion. PLoS ONE. 2011;6:23777.
We extend many thanks to technicians and graduate students in the Gehrt lab for field and technical assistance. We also thank the anonymous reviewers for their helpful comments.
Funding was provided by Donna Alexander from the Cook County Animal and Rabies Control, the Max McGraw Wildlife Foundation, the Forest Preserve District of Cook County, CVM Research Office UMN Ag Experiment Station General Ag Research Funds (MIN-62-098), and UMN AHC Seed Grant.
Ethics approval and consent to participate
Coyote captures were approved by the Ohio State University (IACUC ID: 2013A00000012-89 R1) and by Illinois Department of Natural Resources (permit: IDNR W17.0122).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Map of the Chicago metropolitan area. Black circular lines delineate the four urban zones. Red squares are grids where temperature data were gathered. Temperature data were obtained from the PRISM Climate Group (PRISM Climate Group, Oregon State University, http://prism.oregonstate.edu). Figure S2. Relationship between heartworm infection and the proportion of medium developed urban land in coyote home ranges by age class. Panel a includes both resident and transient coyotes (n = 146), and panel b resident coyotes only (n = 107). The lines are mean estimates of the probability of heartworm infection by proportion urban land in coyote home ranges. The shaded bands are 95% confidence intervals. Coyote home ranges were estimated using 95% adaptive local convex hulls (95% a-LoCoH). Table S1. Number of recaptured coyotes and years of captures (n = 16). Coyotes were grouped based on whether they tested positive or negative on the first and second occasion. Table S2. Linear regression models predicting the duration of the heartworm transmission season (n = 192). Models are ranked based on ΔAICc. Table S3. Generalized linear mixed models (GLMMs) predicting heartworm infection (n = 315). GLMMs are ranked based on ΔAICc. Table S4. Top twenty generalized linear mixed models predicting heartworm infection. Models within ΔAICc < 2 from the best fit model were included in model averaging. Table S5. Model averaging results from binomial generalized linear mixed models of the probability of heartworm infection in coyotes (n = 146) using the adaptive local convex hull (a-LoCoH). Predictors were obtained from the top-ranking models (ΔAICc < 2; Table S4).
About this article
Cite this article
Worsley-Tonks, K.E.L., Gehrt, S.D., Anchor, C. et al. Infection risk varies within urbanized landscapes: the case of coyotes and heartworm. Parasites Vectors 14, 464 (2021). https://doi.org/10.1186/s13071-021-04958-1
- Home range