The ecology of human-caused mortality for a protected large carnivore.

Mitigating human-caused mortality for large carnivores is a pressing global challenge for wildlife conservation. However, mortality is almost exclusively studied at local (within-population) scales creating a mismatch between our understanding of risk and the spatial extent most relevant to conservation and management of wide-ranging species. Here, we quantified mortality for 590 radio-collared mountain lions statewide across their distribution in California to identify drivers of human-caused mortality and investigate whether human-caused mortality is additive or compensatory. Human-caused mortality, primarily from conflict management and vehicles, exceeded natural mortality despite mountain lions being protected from hunting. Our data indicate that human-caused mortality is additive to natural mortality as population-level survival decreased as a function of increasing human-caused mortality and natural mortality did not decrease with increased human-caused mortality. Mortality risk increased for mountain lions closer to rural development and decreased in areas with higher proportions of citizens voting to support environmental initiatives. Thus, the presence of human infrastructure and variation in the mindset of humans sharing landscapes with mountain lions appear to be primary drivers of risk. We show that human-caused mortality can reduce population-level survival of large carnivores across large spatial scales, even when they are protected from hunting.

A circovirus and cycloviruses identified in feces of bobcats (Lynx rufus) in California.

Viruses in the family Circoviridae have small circular single-stranded DNA (ssDNA) genomes. Circoviruses are known to infect a wide variety of animals, with notable disease pathology in psittacine (psittacine beak and feather disease) and porcine (postweaning multisystemic wasting syndrome) species. There is still a dearth of research investigating circoviruses associated with felid species. In six fecal samples collected from bobcats (Lynx rufus) in California from 2010 to 2011, we identified six viruses belonging to the genera Circovirus (n = 1) and Cyclovirus (n = 5), using a high-throughput-sequencing-based approach. Of these, the virus in the genus Circovirus represents a new species, as it shares only 54-60% genome-wide sequence identity with the other members of this genus. The five viruses in the genus Cyclovirus represent three new species, sharing <73% genome-wide sequence identity with all other cycloviruses. Three of the cycloviruses belong to a single putative species and were obtained from the feces of three individual bobcats, sharing 95.7-99.9% sequence identity, whereas the other two unique cycloviruses were identified in a single fecal sample. At present, it is unknown whether the identified viruses infect bobcats, their prey, or their gut parasites.

Novel polyomaviruses identified in fecal samples from four carnivore species.

Polyomaviruses are oncogenic viruses that are generally thought to have co-evolved with their hosts. While primate and rodent polyomaviruses are increasingly well-studied, less is known about polyomaviruses that infect other mammals. In an effort to gain insight into polyomaviruses associated with carnivores, we surveyed fecal samples collected in the USA from bobcats (Lynx rufus), pumas (Puma concolor), Canada lynxes (Lynx canadensis), and grizzly bears (Ursus arctos). Using a viral metagenomic approach, we identified six novel polyomavirus genomes. Surprisingly, four of the six genomes showed a phylogenetic relationship to polyomaviruses found in prey animals. These included a putative rabbit polyomavirus from a bobcat fecal sample and two possible deer-trophic polyomaviruses from Canada lynx feces. One polyomavirus found in a grizzly bear sample was found to be phylogenetically distant from previously identified polyomaviruses. Further analysis of the grizzly bear fecal sample showed that it contained anelloviruses that are known to infect pigs, suggesting that the bear might have preyed on a wild or domestic pig. Interestingly, a polyomavirus genome identified in a puma fecal sample was found to be closely related both to raccoon polyomavirus 1 and to Lyon-IARC polyomavirus, the latter of which was originally identified in human saliva and skin swab specimens but has since been found in samples from domestic cats (Felis catus).

Genomic analyses of gray fox lineages suggest ancient divergence and secondary contact in the southern Great Plains.

The gray fox (Urocyon cinereoargenteus) lineage diverged from all other extant canids at their most basal node and is restricted to the Americas. Previous mitochondrial analysis from coastal populations identified deeply divergent (up to 1 Mya) eastern and western lineages that predate most intraspecific splits in carnivores. We conducted genotyping by sequencing and mitochondrial analysis on gray foxes sampled across North America to determine geographic concordance between nuclear and mitochondrial contact zones and divergence times. We also estimated the admixture within the contact zone between eastern and western gray foxes based on nuclear DNA. Both datasets confirmed that eastern and western lineages met in the southern Great Plains (i.e. Texas and Oklahoma), where they maintained high differentiation. Admixture was generally low, with the majority of admixed individuals carrying <10% ancestry from the other lineage. Divergence times confirmed a mid-Pleistocene split, similar to the mitochondrial estimates. Taken together, findings suggest gray fox lineages represent an ancient divergence event, far older than most intraspecific divergences in North American carnivores. Low admixture may reflect a relatively recent time since secondary contact (e.g. post-Pleistocene) or, alternatively, ecological or reproductive barriers between lineages. Though further research is needed to disentangle these factors, our genomic investigation suggests species-level divergence exists between eastern and western gray fox lineages.

The influence of species life history and distribution characteristics on species responses to habitat fragmentation in an urban landscape.

Fragmentation within urbanized environments often leads to a loss of native species diversity; however, variation exists in responses among-species and among-populations within species. We aimed to identify patterns in species biogeography in an urbanized landscape to understand anthropogenic effects on vertebrate communities and identify species that are more sensitive or resilient to landscape change. We investigated patterns in species richness and species responses to fragmentation in southern Californian small vertebrate communities using multispecies occupancy models and determined factors associated with overall commonness and sensitivity to patch size for 45 small vertebrate species both among and within remaining non-developed patches. In general, smaller patches had fewer species, with amphibian species richness being particularly sensitive to patch size effects. Mammals were generally more common, occurring both in a greater proportion of patches and a higher proportion of the sites within occupied patches. Alternatively, amphibians were generally restricted to larger patches but were more ubiquitous within smaller patches when occupied. Species range size was positively correlated with how common a species was across and within patches, even when controlling for only patches that fell within a species' range. We found sensitivity to patch size was greater for more fecund species and depended on where the patch occurred within a species' range. While all taxa were more likely to occur in patches in the warmer portions of their ranges, amphibians and mammals were more sensitive to fragmentation in these warmer areas as compared to the rest of their ranges. Similarly, amphibians occurred at a smaller proportion of sites within patches in drier portions of their ranges. Mammals occurred at a higher proportion of sites that were also in drier portions of their range while reptiles did not differ in their sensitivity to patch size by range position. We demonstrate that taxonomy, life history, range size and range position can predict commonness and sensitivity of species across this highly fragmented yet biodiverse landscape. The impacts of fragmentation on species communities within an urban landscape depend on scale, with differences emerging among and within species and populations.

Urbanization reduces genetic connectivity in bobcats (Lynx rufus) at both intra- and interpopulation spatial scales.

Urbanization is a major factor driving habitat fragmentation and connectivity loss in wildlife. However, the impacts of urbanization on connectivity can vary among species and even populations due to differences in local landscape characteristics, and our ability to detect these relationships may depend on the spatial scale at which they are measured. Bobcats (Lynx rufus) are relatively sensitive to urbanization and the status of bobcat populations is an important indicator of connectivity in urban coastal southern California. We genotyped 271 bobcats at 13,520 SNP loci to conduct a replicated landscape resistance analysis in five genetically distinct populations. We tested urban and natural factors potentially influencing individual connectivity in each population separately, as well as study-wide. Overall, landscape genomic effects were most frequently detected at the study-wide spatial scale, with urban land cover (measured as impervious surface) having negative effects and topographic roughness having positive effects on gene flow. The negative effect of urban land cover on connectivity was also evident when populations were analyzed separately despite varying substantially in spatial area and the proportion of urban development, confirming a pervasive impact of urbanization largely independent of spatial scale. The effect of urban development was strongest in one population where stream habitat had been lost to development, suggesting that riparian corridors may help mitigate reduced connectivity in urbanizing areas. Our results demonstrate the importance of replicating landscape genetic analyses across populations and considering how landscape genetic effects may vary with spatial scale and local landscape structure.

Extinction vortex dynamics of top predators isolated by urbanization.

Extinction risk is elevated in small, isolated populations due to demographic and genetic interactions. Therefore, it is critical to model these processes realistically in population viability analyses (PVA) to inform local management and contribute to a greater understanding of mechanisms within the extinction vortex. We conducted PVA's for two small mountain lion populations isolated by urbanization in southern California to predict population growth, extinction probability, and loss of genetic diversity with empirical data. Specifically, we (1) provide the first PVA for isolated mountain lions in the Santa Ana Mountains (SAM) that considers both demographic and genetic risk factors and (2) test the hypothesis that variation in abundance and mortality between the SAM and Santa Monica Mountains (SMM) result in differences in population growth, loss of heterozygosity, and extinction probability. Our models predicted 16-21% probability of local extinction in the SAM due purely to demographic processes over 50 yr with current low levels or no immigration. Our models also predicted that genetic diversity will further erode in the SAM such that concern regarding inbreeding depression is warranted unless gene flow is increased, and that if inbreeding depression occurs, rapid local extinction will be highly likely. Dynamics of the two populations were broadly similar, but they also exhibited differences driven by larger population size and higher mortality in the SAM. Density-independent scenarios predicted a rapidly increasing population in the SMM, whereas growth potential did not differ from a stable trend in the SAM. Demographic extinction probability and loss of heterozygosity were greater in the SMM for density-dependent scenarios without immigration. However, higher levels of immigration had stronger, positive influences on both demographic viability and retention of genetic diversity in the SMM driven by lower abundance and higher adult survival. Our results elucidate demographic and genetic threats to small populations within the extinction vortex, and how these vary relative to demographic structure. Importantly, simulating seemingly attainable increases in connectivity was sufficient to greatly reduce extinction probability. Our work highlights that conservation of large carnivores is achievable within urbanized landscapes, but requires land protection, connectivity, and strategies to promote coexistence with humans.

Novel smacoviruses identified in the faeces of two wild felids: North American bobcat and African lion.

Smacoviruses are small circular single-stranded DNA viruses that appear to be prevalent in faeces of a range of animals and have also been found in a few insect species. In this study, we report the first viral genomes from faeces of free-roaming wild felids on two continents. Two smacoviruses were recovered from the faeces of two North American bobcats (Lynx rufus), and one was recovered from an African lion (Panthera leo). All three genomes are genetically different, sharing 59-69% genome-wide sequence identity to other smacoviruses. These are the first full smacovirus genome sequences associated with a large top-end feline predator, and their presence in these samples suggests that feline faeces are a natural niche for the organisms that these viruses infect.

Urbanization and anticoagulant poisons promote immune dysfunction in bobcats.

Understanding how human activities influence immune response to environmental stressors can support biodiversity conservation across increasingly urbanizing landscapes. We studied a bobcat (Lynx rufus) population in urban southern California that experienced a rapid population decline from 2002-2005 due to notoedric mange. Because anticoagulant rodenticide (AR) exposure was an underlying complication in mange deaths, we aimed to understand sublethal contributions of urbanization and ARs on 65 biochemical markers of immune and organ function. Variance in immunological variables was primarily associated with AR exposure and secondarily with urbanization. Use of urban habitat and AR exposure has pervasive, complex and predictable effects on biochemical markers of immune and organ function in free-ranging bobcats that include impacts on neutrophil, lymphocyte and cytokine populations, total bilirubin and phosphorus. We find evidence of both inflammatory response and immune suppression associated with urban land use and rat poison exposure that could influence susceptibility to opportunistic infections. Consequently, AR exposure may influence mortality and has population-level effects, as previous work in the focal population has revealed substantial mortality caused by mange infection. The secondary effects of anticoagulant exposure may be a worldwide, largely unrecognized problem affecting a variety of vertebrate species in human-dominated environments.

Feline Immunodeficiency Virus Cross-Species Transmission: Implications for Emergence of New Lentiviral Infections.

Owing to a complex history of host-parasite coevolution, lentiviruses exhibit a high degree of species specificity. Given the well-documented viral archeology of human immunodeficiency virus (HIV) emergence following human exposures to simian immunodeficiency virus (SIV), an understanding of processes that promote successful cross-species lentiviral transmissions is highly relevant. We previously reported natural cross-species transmission of a subtype of feline immunodeficiency virus, puma lentivirus A (PLVA), between bobcats (Lynx rufus) and mountain lions (Puma concolor) for a small number of animals in California and Florida. In this study, we investigate host-specific selection pressures, within-host viral fitness, and inter- versus intraspecies transmission patterns among a larger collection of PLV isolates from free-ranging bobcats and mountain lions. Analyses of proviral and viral RNA levels demonstrate that PLVA fitness is severely restricted in mountain lions compared to that in bobcats. We document evidence of diversifying selection in three of six PLVA genomes from mountain lions, but we did not detect selection among 20 PLVA isolates from bobcats. These findings support the hypothesis that PLVA is a bobcat-adapted virus which is less fit in mountain lions and under intense selection pressure in the novel host. Ancestral reconstruction of transmission events reveals that intraspecific PLVA transmission has occurred among panthers (Puma concolor coryi) in Florida following the initial cross-species infection from bobcats. In contrast, interspecific transmission from bobcats to mountain lions predominates in California. These findings document outcomes of cross-species lentiviral transmission events among felids that compare to the emergence of HIV from nonhuman primates.IMPORTANCE Cross-species transmission episodes can be singular, dead-end events or can result in viral replication and spread in the new species. The factors that determine which outcome will occur are complex, and the risk of new virus emergence is therefore difficult to predict. We used molecular techniques to evaluate the transmission, fitness, and adaptation of puma lentivirus A (PLVA) between bobcats and mountain lions in two geographic regions. Our findings illustrate that mountain lion exposure to PLVA is relatively common but does not routinely result in communicable infections in the new host. This is attributed to efficient species barriers that largely prevent lentiviral adaptation. However, the evolutionary capacity for lentiviruses to adapt to novel environments may ultimately overcome host restriction mechanisms over time and under certain ecological circumstances. This phenomenon provides a unique opportunity to examine cross-species transmission events leading to new lentiviral emergence.

Genome-wide expression reveals multiple systemic effects associated with detection of anticoagulant poisons in bobcats (Lynx rufus).

Anticoagulant rodenticides (ARs) are indiscriminate toxicants that threaten nontarget predatory and scavenger species through secondary poisoning. Accumulating evidence suggests that AR exposure may have disruptive sublethal consequences on individuals that can affect fitness. We evaluated AR-related effects on genome-wide expression patterns in a population of bobcats in southern California. We identify differential expression of genes involved in xenobiotic metabolism, endoplasmic reticulum stress response, epithelial integrity and both adaptive and innate immune function. Further, we find that differential expression of immune-related genes may be attributable to AR-related effects on leucocyte differentiation. Collectively, our results provide an unprecedented understanding of the sublethal effects of AR exposure on a wild carnivore. These findings highlight potential detrimental effects of ARs on a wide variety of species worldwide that may consume poisoned rodents and indicate the need to investigate gene expression effects of other toxicants added to natural environments by humans.

Determining the drivers of population structure in a highly urbanized landscape to inform conservation planning.

Understanding the environmental contributors to population structure is of paramount importance for conservation in urbanized environments. We used spatially explicit models to determine genetic population structure under current and future environmental conditions across a highly fragmented, human-dominated environment in Southern California to assess the effects of natural ecological variation and urbanization. We focused on 7 common species with diverse habitat requirements, home-range sizes, and dispersal abilities. We quantified the relative roles of potential barriers, including natural environmental characteristics and an anthropogenic barrier created by a major highway, in shaping genetic variation. The ability to predict genetic variation in our models differed among species: 11-81% of intraspecific genetic variation was explained by environmental variables. Although an anthropogenically induced barrier (a major highway) severely restricted gene flow and movement at broad scales for some species, genetic variation seemed to be primarily driven by natural environmental heterogeneity at a local level. Our results show how assessing environmentally associated variation for multiple species under current and future climate conditions can help identify priority regions for maximizing population persistence under environmental change in urbanized regions.

Interactions between demography, genetics, and landscape connectivity increase extinction probability for a small population of large carnivores in a major metropolitan area.

The extinction vortex is a theoretical model describing the process by which extinction risk is elevated in small, isolated populations owing to interactions between environmental, demographic, and genetic factors. However, empirical demonstrations of these interactions have been elusive. We modelled the dynamics of a small mountain lion population isolated by anthropogenic barriers in greater Los Angeles, California, to evaluate the influence of demographic, genetic, and landscape factors on extinction probability. The population exhibited strong survival and reproduction, and the model predicted stable median population growth and a 15% probability of extinction over 50 years in the absence of inbreeding depression. However, our model also predicted the population will lose 40-57% of its heterozygosity in 50 years. When we reduced demographic parameters proportional to reductions documented in another wild population of mountain lions that experienced inbreeding depression, extinction probability rose to 99.7%. Simulating greater landscape connectivity by increasing immigration to greater than or equal to one migrant per generation appears sufficient to largely maintain genetic diversity and reduce extinction probability. We provide empirical support for the central tenet of the extinction vortex as interactions between genetics and demography greatly increased extinction probability relative to the risk from demographic and environmental stochasticity alone. Our modelling approach realistically integrates demographic and genetic data to provide a comprehensive assessment of factors threatening small populations.

Pathogen exposure varies widely among sympatric populations of wild and domestic felids across the United States.

Understanding how landscape, host, and pathogen traits contribute to disease exposure requires systematic evaluations of pathogens within and among host species and geographic regions. The relative importance of these attributes is critical for management of wildlife and mitigating domestic animal and human disease, particularly given rapid ecological changes, such as urbanization. We screened > 1000 samples from sympatric populations of puma (Puma concolor), bobcat (Lynx rufus), and domestic cat (Felis catus) across urban gradients in six sites, representing three regions, in North America for exposure to a representative suite of bacterial, protozoal, and viral pathogens (Bartonella sp., Toxoplasma gondii, feline herpesvirus-1, feline panleukopenea virus, feline calicivirus, and feline immunodeficiency virus). We evaluated prevalence within each species, and examined host trait and land cover determinants of exposure; providing an unprecedented analysis of factors relating to potential for infections in domesticated and wild felids. Prevalence differed among host species (highest for puma and lowest for domestic cat) and was greater for indirectly transmitted pathogens. Sex was inconsistently predictive of exposure to directly transmitted pathogens only, and age infrequently predictive of both direct and indirectly transmitted pathogens. Determinants of pathogen exposure were widely divergent between the wild felid species. For puma, suburban land use predicted increased exposure to Bartonella sp. in southern California, and FHV-1 exposure increased near urban edges in Florida. This may suggest interspecific transmission with domestic cats via flea vectors (California) and direct contact (Florida) around urban boundaries. Bobcats captured near urban areas had increased exposure to T. gondii in Florida, suggesting an urban source of prey Bobcats captured near urban areas in Colorado and Florida had higher FIV exposure, possibly suggesting increased intraspecific interactions through pile-up of home ranges. Beyond these regional and pathogen specific relationships, proximity to the wildland-urban interface did not generally increase the probability of disease exposure in wild or domestic felids, empha- sizing the importance of local ecological determinants. Indeed, pathogen exposure was often negatively associated with the wildland-urban interface for all felids. Our analyses suggest cross-species pathogen transmission events around this interface may be infrequent, but followed by self-sustaining propagation within the new host species. virus; puma (Puma concolor); Toxoplasma gondii; urbanization.

Evolution of puma lentivirus in bobcats (Lynx rufus) and mountain lions (Puma concolor) in North America.

Mountain lions (Puma concolor) throughout North and South America are infected with puma lentivirus clade B (PLVB). A second, highly divergent lentiviral clade, PLVA, infects mountain lions in southern California and Florida. Bobcats (Lynx rufus) in these two geographic regions are also infected with PLVA, and to date, this is the only strain of lentivirus identified in bobcats. We sequenced full-length PLV genomes in order to characterize the molecular evolution of PLV in bobcats and mountain lions. Low sequence homology (88% average pairwise identity) and frequent recombination (1 recombination breakpoint per 3 isolates analyzed) were observed in both clades. Viral proteins have markedly different patterns of evolution; sequence homology and negative selection were highest in Gag and Pol and lowest in Vif and Env. A total of 1.7% of sites across the PLV genome evolve under positive selection, indicating that host-imposed selection pressure is an important force shaping PLV evolution. PLVA strains are highly spatially structured, reflecting the population dynamics of their primary host, the bobcat. In contrast, the phylogeography of PLVB reflects the highly mobile mountain lion, with diverse PLVB isolates cocirculating in some areas and genetically related viruses being present in populations separated by thousands of kilometers. We conclude that PLVA and PLVB are two different viral species with distinct feline hosts and evolutionary histories. Importance: An understanding of viral evolution in natural host populations is a fundamental goal of virology, molecular biology, and disease ecology. Here we provide a detailed analysis of puma lentivirus (PLV) evolution in two natural carnivore hosts, the bobcat and mountain lion. Our results illustrate that PLV evolution is a dynamic process that results from high rates of viral mutation/recombination and host-imposed selection pressure.

Novel gammaherpesviruses in North American domestic cats, bobcats, and pumas: identification, prevalence, and risk factors.

UNLABELLED: Gammaherpesviruses (GHVs) are a diverse and rapidly expanding group of viruses associated with a variety of disease conditions in humans and animals. To identify felid GHVs, we screened domestic cat (Felis catus), bobcat (Lynx rufus), and puma (Puma concolor) blood cell DNA samples from California, Colorado, and Florida using a degenerate pan-GHV PCR. Additional pan-GHV and long-distance PCRs were used to sequence a contiguous 3.4-kb region of each putative virus species, including partial glycoprotein B and DNA polymerase genes. We identified three novel GHVs, each present predominantly in one felid species: Felis catus GHV 1 (FcaGHV1) in domestic cats, Lynx rufus GHV 1 (LruGHV1) in bobcats, and Puma concolor GHV 1 (PcoGHV1) in pumas. To estimate infection prevalence, we developed real-time quantitative PCR assays for each virus and screened additional DNA samples from all three species (n = 282). FcaGHV1 was detected in 16% of domestic cats across all study sites. LruGHV1 was detected in 47% of bobcats and 13% of pumas across all study sites, suggesting relatively common interspecific transmission. PcoGHV1 was detected in 6% of pumas, all from a specific region of Southern California. The risk of infection for each host varied with geographic location. Age was a positive risk factor for bobcat LruGHV1 infection, and age and being male were risk factors for domestic cat FcaGHV1 infection. Further characterization of these viruses may have significant health implications for domestic cats and may aid studies of free-ranging felid ecology. IMPORTANCE: Gammaherpesviruses (GHVs) establish lifelong infection in many animal species and can cause cancer and other diseases in humans and animals. In this study, we identified the DNA sequences of three GHVs present in the blood of domestic cats (Felis catus), bobcats (Lynx rufus), and pumas (Puma concolor; also known as mountain lions, cougars, and panthers). We found that these viruses were closely related to, but distinct from, other known GHVs of animals and represent the first GHVs identified to be native to these feline species. We developed techniques to rapidly and specifically detect the DNA of these viruses in feline blood and found that the domestic cat and bobcat viruses were widespread across the United States. In contrast, puma virus was found only in a specific region of Southern California. Surprisingly, the bobcat virus was also detected in some pumas, suggesting relatively common virus transmission between these species. Adult domestic cats and bobcats were at greater risk for infection than juveniles. Male domestic cats were at greater risk for infection than females. This study identifies three new viruses that are widespread in three feline species, indicates risk factors for infection that may relate to the route of infection, and demonstrates cross-species transmission between bobcats and pumas. These newly identified viruses may have important effects on feline health and ecology.

Treat Yourself: Pilot Testing a New Method to Treat Mange in Wild Carnivores.

Mange is a skin disease caused by mites that parasitize an animal's skin, often yielding inflamed immune responses and hair loss. At a population level, mange may reduce survival and cause population declines. Many forms of mange can be treated quite effectively when an animal is in hand; however, this is not often feasible for many free-ranging wildlife populations. Some animals, particularly territorial carnivores, will rub or roll to scent mark and transmit information about their presence to other individuals. We posited that rub stations comprised, in part, of anthelmintic medication and foreign scents that induce rubbing could be used to remotely treat mange in the wild. We deployed 39 rub stations containing lure and dye in Santa Monica Mountains National Recreation Area, Southern California, USA, October-November 2022. Carnivores rubbed or rolled at >97% of rub stations, with coyotes (Canis latrans), gray foxes (Urocyon cinereoargenteus), and bobcats (Lynx rufus) being the most abundant species. Time to first rub or roll was generally <1 wk. Several sympatric species (e.g., mule deer, Odocoileus hemionus) were detected at rub stations but did not rub. Our pilot test provides strong evidence that treating mange in wild carnivores may be possible using the remote medicinal rub stations we describe. Future efforts to add medicine to rub stations and monitor for a change in mange prevalence are a logical next step.

Vector-Borne Disease in Wild Mammals Impacted by Urban Expansion and Climate Change.

Ecologies of zoonotic vector-borne diseases may shift with climate and land use change. As many urban-adapted mammals can host ectoparasites and pathogens of human and animal health concern, our goal was to compare patterns of arthropod-borne disease among medium-sized mammals across gradients of rural to urban landscapes in multiple regions of California. DNA of Anaplasma phagocytophilum was found in 1-5% of raccoons, coyotes, and San Joaquin kit foxes; Borrelia burgdorferi in one coyote, rickettsiae in two desert kit foxes, and Yersinia pestis in two coyotes. There was serological evidence of rickettsiae in 14-37% of coyotes, Virginia opossums, and foxes; and A. phagocytophilum in 6-40% of coyotes, raccoons, Virginia opossums, and foxes. Of six flea species, one Ctenocephalides felis from a raccoon was positive for Y. pestis, and Ct. felis and Pulex simulans fleas tested positive for Rickettsia felis and R. senegalensis. A Dermacentor similis tick off a San Joaquin kit fox was PCR-positive for A. phagocytophilum. There were three statistically significant risk factors: risk of A. phagocytophilum PCR-positivity was threefold greater in fall vs the other three seasons; hosts adjacent to urban areas had sevenfold increased A. phagocytophilum seropositivity compared with urban and rural areas; and there was a significant spatial cluster of rickettsiae within greater Los Angeles. Animals in areas where urban and rural habitats interconnect can serve as sentinels during times of change in disease risk.