Gene flow and pathogen transmission among bobcats (Lynx rufus) in a fragmented urban landscape

Habitat connectivity and host relatedness influence virus spread across an urbanising landscape in a fragmentation-sensitive carnivore

Spatially heterogeneous landscape factors such as urbanisation can have substantial effects on the severity and spread of wildlife diseases. However, research linking patterns of pathogen transmission to landscape features remains rare. Using a combination of phylogeographic and machine learning approaches, we tested the influence of landscape and host factors on feline immunodeficiency virus (FIV_Lru) genetic variation and spread among bobcats (Lynx rufus) sampled from coastal southern California. We found evidence for increased rates of FIV_Lru lineage spread through areas of higher vegetation density. Furthermore, single-nucleotide polymorphism (SNP) variation among FIV_Lru sequences was associated with host genetic distances and geographic location, with FIV_Lru genetic discontinuities precisely correlating with known urban barriers to host dispersal. An effect of forest land cover on FIV_Lru SNP variation was likely attributable to host population structure and differences in forest land cover between different populations. Taken together, these results suggest that the spread of FIV_Lru is constrained by large-scale urban barriers to host movement. Although urbanisation at fine spatial scales did not appear to directly influence virus transmission or spread, we found evidence that viruses transmit and spread more quickly through areas containing higher proportions of natural habitat. These multiple lines of evidence demonstrate how urbanisation can change patterns of contact-dependent pathogen transmission and provide insights into how continued urban development may influence the incidence and management of wildlife disease.

Agricultural land use shapes dispersal in white-tailed deer (Odocoileus virginianus)

BACKGROUND

Dispersal is a fundamental process to animal population dynamics and gene flow. In white-tailed deer (WTD; Odocoileus virginianus), dispersal also presents an increasingly relevant risk for the spread of infectious diseases. Across their wide range, WTD dispersal is believed to be driven by a suite of landscape and host behavioral factors, but these can vary by region, season, and sex. Our objectives were to (1) identify dispersal events in Wisconsin WTD and determine drivers of dispersal rates and distances, and (2) determine how landscape features (e.g., rivers, roads) structure deer dispersal paths.

METHODS

We developed an algorithmic approach to detect dispersal events from GPS collar data for 590 juvenile, yearling, and adult WTD. We used statistical models to identify host and landscape drivers of dispersal rates and distances, including the role of agricultural land use, the traversability of the landscape, and potential interactions between deer. We then performed a step selection analysis to determine how landscape features such as agricultural land use, elevation, rivers, and roads affected deer dispersal paths.

RESULTS

Dispersal predominantly occurred in juvenile males, of which 64.2% dispersed, with dispersal events uncommon in other sex and age classes. Juvenile male dispersal probability was positively associated with the proportion of the natal range that was classified as agricultural land use, but only during the spring. Dispersal distances were typically short (median 5.77 km, range: 1.3-68.3 km), especially in the fall. Further, dispersal distances were positively associated with agricultural land use in potential dispersal paths but negatively associated with the number of proximate deer in the natal range. Lastly, we found that, during dispersal, juvenile males typically avoided agricultural land use but selected for areas near rivers and streams.

CONCLUSION

Land use-particularly agricultural-was a key driver of dispersal rates, distances, and paths in Wisconsin WTD. In addition, our results support the importance of deer social environments in shaping dispersal behavior. Our findings reinforce knowledge of dispersal ecology in WTD and how landscape factors-including major rivers, roads, and land-use patterns-structure host gene flow and potential pathogen transmission.

Parasites as conservation tools

Parasite success typically depends on a close relationship with one or more hosts; therefore, attributes of parasitic infection have the potential to provide indirect details of host natural history and are biologically relevant to animal conservation. Characterization of parasite infections has been useful in delineating host populations and has served as a proxy for assessment of environmental quality. In other cases, the utility of parasites is just being explored, for example, as indicators of host connectivity. Innovative studies of parasite biology can provide information to manage major conservation threats by using parasite assemblage, prevalence, or genetic data to provide insights into the host. Overexploitation, habitat loss and fragmentation, invasive species, and climate change are major threats to animal conservation, and all of these can be informed by parasites.

Host, microbiome, and complex space: applying population and landscape genetic approaches to gut microbiome research in wild populations

In recent years, emerging sequencing technologies and computational tools have driven a tidal wave of research on host-associated microbiomes, particularly the gut microbiome. These studies demonstrate numerous connections between the gut microbiome and vital host functions, primarily in humans, model organisms, and domestic animals. As the adaptive importance of the gut microbiome becomes clearer, interest in studying the gut microbiomes of wild populations has increased, in part due to the potential for discovering conservation applications. The study of wildlife gut microbiomes holds many new challenges and opportunities due to the complex genetic, spatial, and environmental structure of wild host populations, and the potential for these factors to interact with the microbiome. The emerging picture of adaptive coevolution in host-microbiome relationships highlights the importance of understanding microbiome variation in the context of host population genetics and landscape heterogeneity across a wide range of host populations. We propose a conceptual framework for understanding wildlife gut microbiomes in relation to landscape variables and host population genetics, including the potential of approaches derived from landscape genetics. We use this framework to review current research, synthesize important trends, highlight implications for conservation, and recommend future directions for research. Specifically, we focus on how spatial structure and environmental variation interact with host population genetics and microbiome variation in natural populations, and what we can learn from how these patterns of covariation differ depending on host ecological and evolutionary traits.

Viral Sequences Recovered From Puma Tooth DNA Reconstruct Statewide Viral Phylogenies

Monitoring pathogens in wildlife populations is imperative for effective management, and for identifying locations for pathogen spillover among wildlife, domestic species and humans. Wildlife pathogen surveillance is challenging, however, as sampling often requires the capture of a significant proportion of the population to understand host pathogen dynamics. To address this challenge, we assessed the ability to use hunter-collected teeth from puma across Colorado to recover genetic data of two feline retroviruses, feline foamy virus (FFV) and feline immunodeficiency virus (FIVpco) and show they can be utilized for this purpose. Comparative phylogenetic analyses of FIVpco and FFV from tooth and blood samples to previous analyses conducted with blood samples collected over a nine-year period from two distinct areas was undertaken highlighting the value of tooth derived samples. We found less FIVpco phylogeographic structuring than observed from sampling only two regions and that FFV data confirmed previous findings of endemic infection, minimal geographic structuring, and supported frequent cross-species transmission from domestic cats to pumas. Viral analysis conducted using intentionally collected blood samples required extensive financial, capture and sampling efforts. This analysis illustrates that viral genomic data can be cost effectively obtained using tooth samples incidentally-collected from hunter harvested pumas, taking advantage of samples collected for morphological age identification. This technique should be considered as an opportunistic method to provide broad geographic sampling to define viral dynamics more accurately in wildlife.

MrIML: Multi-response interpretable machine learning to model genomic landscapes

We introduce a new R package "MrIML" ("Mister iml"; Multi-response Interpretable Machine Learning). MrIML provides a powerful and interpretable framework that enables users to harness recent advances in machine learning to quantify multilocus genomic relationships, to identify loci of interest for future landscape genetics studies, and to gain new insights into adaptation across environmental gradients. Relationships between genetic variation and environment are often nonlinear and interactive; these characteristics have been challenging to address using traditional landscape genetic approaches. Our package helps capture this complexity and offers functions that fit and interpret a wide range of highly flexible models that are routinely used for single-locus landscape genetics studies but are rarely extended to estimate response functions for multiple loci. To demonstrate the package's broad functionality, we test its ability to recover landscape relationships from simulated genomic data. We also apply the package to two empirical case studies. In the first, we model genetic variation of North American balsam poplar (Populus balsamifera, Salicaceae) populations across environmental gradients. In the second case study, we recover the landscape and host drivers of feline immunodeficiency virus genetic variation in bobcats (Lynx rufus). The ability to model thousands of loci collectively and compare models from linear regression to extreme gradient boosting, within the same analytical framework, has the potential to be transformative. The MrIML framework is also extendable and not limited to modelling genetic variation; for example, it can quantify the environmental drivers of microbiomes and coinfection dynamics.

Host relatedness and landscape connectivity shape pathogen spread in the puma, a large secretive carnivore

Urban expansion can fundamentally alter wildlife movement and gene flow, but how urbanization alters pathogen spread is poorly understood. Here, we combine high resolution host and viral genomic data with landscape variables to examine the context of viral spread in puma (Puma concolor) from two contrasting regions: one bounded by the wildland urban interface (WUI) and one unbounded with minimal anthropogenic development (UB). We found landscape variables and host gene flow explained significant amounts of variation of feline immunodeficiency virus (FIV) spread in the WUI, but not in the unbounded region. The most important predictors of viral spread also differed; host spatial proximity, host relatedness, and mountain ranges played a role in FIV spread in the WUI, whereas roads might have facilitated viral spread in the unbounded region. Our research demonstrates how anthropogenic landscapes can alter pathogen spread, providing a more nuanced understanding of host-pathogen relationships to inform disease ecology in free-ranging species.

Does the virus cross the road? Viral phylogeographic patterns among bobcat populations reflect a history of urban development

Urban development has major impacts on connectivity among wildlife populations and is thus likely an important factor shaping pathogen transmission in wildlife. However, most investigations of wildlife diseases in urban areas focus on prevalence and infection risk rather than potential effects of urbanization on transmission itself. Feline immunodeficiency virus (FIV) is a directly transmitted retrovirus that infects many felid species and can be used as a model for studying pathogen transmission at landscape scales. We investigated phylogenetic relationships among FIV isolates sampled from five bobcat (Lynx rufus) populations in coastal southern California that appear isolated due to major highways and dense urban development. Divergence dates among FIV phylogenetic lineages in several cases reflected historical urban growth and construction of major highways. We found strong FIV phylogeographic structure among three host populations north-west of Los Angeles, largely coincident with host genetic structure. In contrast, relatively little FIV phylogeographic structure existed among two genetically distinct host populations south-east of Los Angeles. Rates of FIV transfer among host populations did not vary significantly, with the lack of phylogenetic structure south-east of Los Angeles unlikely to reflect frequent contemporary transmission among populations. Our results indicate that major barriers to host gene flow can also act as barriers to pathogen spread, suggesting potentially reduced susceptibility of fragmented populations to novel directly transmitted pathogens. Infrequent exchange of FIV among host populations suggests that populations would best be managed as distinct units in the event of a severe disease outbreak. Phylogeographic inference of pathogen transmission is useful for estimating the ability of geographic barriers to constrain disease spread and can provide insights into contemporary and historical drivers of host population connectivity.

Comparative landscape genetics reveals differential effects of environment on host and pathogen genetic structure in Tasmanian devils (Sarcophilus harrisii) and their transmissible tumour

Genetic structure in host species is often used to predict disease spread. However, host and pathogen genetic variation may be incongruent. Understanding landscape factors that have either concordant or divergent influence on host and pathogen genetic structure is crucial for wildlife disease management. Devil facial tumour disease (DFTD) was first observed in 1996 and has spread throughout almost the entire Tasmanian devil geographic range, causing dramatic population declines. Whereas DFTD is predominantly spread via biting among adults, devils typically disperse as juveniles, which experience low DFTD prevalence. Thus, we predicted little association between devil and tumour population structure and that environmental factors influencing gene flow differ between devils and tumours. We employed a comparative landscape genetics framework to test the influence of environmental factors on patterns of isolation by resistance (IBR) and isolation by environment (IBE) in devils and DFTD. Although we found evidence for broad-scale costructuring between devils and tumours, we found no relationship between host and tumour individual genetic distances. Further, the factors driving the spatial distribution of genetic variation differed for each. Devils exhibited a strong IBR pattern driven by major roads, with no evidence of IBE. By contrast, tumours showed little evidence for IBR and a weak IBE pattern with respect to elevation in one of two tumour clusters we identify herein. Our results warrant caution when inferring pathogen spread using host population genetic structure and suggest that reliance on environmental barriers to host connectivity may be ineffective for managing the spread of wildlife diseases. Our findings demonstrate the utility of comparative landscape genetics for identifying differential factors driving host dispersal and pathogen transmission.

Population genetic structure and gene flow of rare and endangered Tetraena mongolica Maxim. revealed by reduced representation sequencing

BACKGROUND

Studying population genetic structure and gene flow of plant populations and their influencing factors is of particular significance in the field of conservation biology, especially important for species such as rare and endangered plants. Tetraena mongolica Maxim. (TM), belongs to Zygophyllaceae family, a rare and endangered plant with narrow distribution. However, for the last decade, due to excessive logging, urban expansion, industrial and tourism development, habitat fragmentation and loss of natural habitats have become major threats to the population of endangered plants.

RESULTS

In this study, genetic diversity, population genetic structure and gene flow of TM populations were evaluated by reduced representation sequencing technology, and a total of more than 133.45 GB high-quality clean reads and 38,097 high-quality SNPs were generated. Analysis based on multiple methods, we found that the existing TM populations have moderate levels of genetic diversity, and very low genetic differentiation as well as high levels of gene flow between populations. Population structure and principal coordinates analysis showed that 8 TM populations can be divided into two groups. The Mantel test detected no significant correlation between geographical distances and genetic distance for the whole sampling. Moreover, the migration model indicated that the gene flow is more of a north to south migration pattern in history.

CONCLUSIONS

This study demonstrates that the present genetic structure is mainly due to habitat fragmentation caused by urban sprawl, industrial development and coal mining. Our recommendation with respect to conservation management is that, all 8 populations should be preserved as a whole population, rather than just those in the core area of TM nature reserve. In particular, the populations near the edge of TM distribution in cities and industrial areas deserve our special protection.

Quantitative characterization of animal social organization: Applications for epidemiological modelling

Social organization is a key aspect of animal ecology, closely interlinked with all aspects of animal behaviour. The structure of animal assemblages is highly diverse, both within and between species. The complexity and variety of social systems and the dynamic nature of interactions and dependencies between members of social groups have long been major obstacles for developing operational characterizations of social organization. Here, social network analysis, a set of statistical tools rooted in graph theory, suggests itself as a potential solution for this problem, by offering quantitative measures for various aspects of social relationships. In this review I will first introduce network analysis as a tool to characterize the social organization of animal groups and population and, then, focus on the application of this method for epidemiological modelling, specifically the prediction of spreading patterns of pathogens in livestock and its potential for informing targeted surveillance and planning of intervention measures.

The Expectations and Challenges of Wildlife Disease Research in the Era of Genomics: Forecasting with a Horizon Scan-like Exercise

The outbreak and transmission of disease-causing pathogens are contributing to the unprecedented rate of biodiversity decline. Recent advances in genomics have coalesced into powerful tools to monitor, detect, and reconstruct the role of pathogens impacting wildlife populations. Wildlife researchers are thus uniquely positioned to merge ecological and evolutionary studies with genomic technologies to exploit unprecedented "Big Data" tools in disease research; however, many researchers lack the training and expertise required to use these computationally intensive methodologies. To address this disparity, the inaugural "Genomics of Disease in Wildlife" workshop assembled early to mid-career professionals with expertise across scientific disciplines (e.g., genomics, wildlife biology, veterinary sciences, and conservation management) for training in the application of genomic tools to wildlife disease research. A horizon scanning-like exercise, an activity to identify forthcoming trends and challenges, performed by the workshop participants identified and discussed 5 themes considered to be the most pressing to the application of genomics in wildlife disease research: 1) "Improving communication," 2) "Methodological and analytical advancements," 3) "Translation into practice," 4) "Integrating landscape ecology and genomics," and 5) "Emerging new questions." Wide-ranging solutions from the horizon scan were international in scope, itemized both deficiencies and strengths in wildlife genomic initiatives, promoted the use of genomic technologies to unite wildlife and human disease research, and advocated best practices for optimal use of genomic tools in wildlife disease projects. The results offer a glimpse of the potential revolution in human and wildlife disease research possible through multi-disciplinary collaborations at local, regional, and global scales.

When are pathogen dynamics likely to reflect host population genetic structure?

Does the structure and connectivity of host populations influence the dynamics and evolution of their pathogens? This topical question is the essence of research investigating the ecology of a Pteropus fruit bat and its zoonotic Nipah virus (NiV) published by Olival et al. in this issue of Molecular Ecology. Questioned less overtly, but nonetheless implicit to the study, is "what are the mechanisms underpinning intraspecific host-pathogen congruence (IHPC) of genetic structure?". Olival et al. investigated the phylogeographical structure of Pteropus medius and NiV isolates across Bangladesh, from areas inside and outside of the Nipah belt-an area where most human spillover events occur. A high degree of host panmixia was discovered, with some population differentiation east of the Nipah belt. NiV genetic structure was congruent with the host. The authors attributed the panmixia and structuring, respectively, to (a) the highly vagile nature of P. medius, and (b) possible differences between bioregions within and outside the Nipah belt. Other potential explanatory mechanisms were acknowledged, including hybridization and transmission mode. This study makes a valuable contribution to a growing body of literature examining IHPC. This has implications not only for pathogen spillover to humans and domestic animals, but more generally for thinking about the mechanisms that underlie patterns of host and pathogen genetic associations.

Climate connectivity of the bobcat in the Great Lakes region

The Great Lakes and the St. Lawrence River are imposing barriers for wildlife, and the additive effect of urban and agricultural development that dominates the lower Great Lakes region likely further reduces functional connectivity for many terrestrial species. As the climate warms, species will need to track climate across these barriers. It is important therefore to investigate land cover and bioclimatic hypotheses that may explain the northward expansion of species through the Great Lakes. We investigated the functional connectivity of a vagile generalist, the bobcat, as a representative generalist forest species common to the region. We genotyped tissue samples collected across the region at 14 microsatellite loci and compared different landscape hypotheses that might explain the observed gene flow or functional connectivity. We found that the Great Lakes and the additive influence of forest stands with either low or high canopy cover and deep lake-effect snow have disrupted gene flow, whereas intermediate forest cover has facilitated gene flow. Functional connectivity in southern Ontario is relatively low and was limited in part by the low amount of forest cover. Pathways across the Great Lakes were through the Niagara region and through the Lower Peninsula of Michigan over the Straits of Mackinac and the St. Marys River. These pathways are important routes for bobcat range expansion north of the Great Lakes and are also likely pathways that many other mobile habitat generalists must navigate to track the changing climate. The extent to which species can navigate these routes will be important for determining the future biodiversity of areas north of the Great Lakes.

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.

Population genetics of fruit bat reservoir informs the dynamics, distribution and diversity of Nipah virus

The structure and connectivity of wildlife host populations may influence zoonotic disease dynamics, evolution and therefore spillover risk to people. Fruit bats in the genus Pteropus, or flying foxes, are the primary natural reservoir for henipaviruses-a group of emerging paramyxoviruses that threaten livestock and public health. In Bangladesh, Pteropus medius is the reservoir for Nipah virus-and viral spillover has led to human fatalities nearly every year since 2001. Here, we use mitochondrial DNA and nuclear microsatellite markers to measure the population structure, demographic history and phylogeography of P. medius in Bangladesh. We combine this with a phylogeographic analysis of all known Nipah virus sequences and strains currently available to better inform the dynamics, distribution and evolutionary history of Nipah virus. We show that P. medius is primarily panmictic, but combined analysis of microsatellite and morphological data shows evidence for differentiation of two populations in eastern Bangladesh, corresponding to a divergent strain of Nipah virus also found in bats from eastern Bangladesh. Our demographic analyses indicate that a large, expanding population of flying foxes has existed in Bangladesh since the Late Pleistocene, coinciding with human population expansion in South Asia, suggesting repeated historical spillover of Nipah virus likely occurred. We present the first evidence of mitochondrial introgression, or hybridization, between P. medius and flying fox species found in South-East Asia (P. vampyrus and P. hypomelanus), which may help to explain the distribution of Nipah virus strains across the region.

The mystery of the phylogeographic structural pattern in rodent-borne hantaviruses

Hantaviruses (order Bunyavirales, family Hantaviridae) are important zoonotic pathogens. Because of the great diversity of their reservoir hosts, hantaviruses are excellent models to evaluate the dynamics of virus-host co-evolution. To understand the mechanisms behind the evolutionary history of hantaviruses through virus-reservoir interactions, it is important to know how the radiation and diversity of hantaviruses occurred. In this paper, we evaluate the pattern of hantavirus diversification based on a complete S segment representing major groups of hantaviruses found in the Americas. Phylogenetic analyses revealed a high degree of phylogeographic structure and a surprising pattern of geographical distribution of New World hantaviruses. The available data suggest that hantaviruses related to the Arvicolinae rodent subfamily in North America probably emerged and initially adapted from a shared common ancestor of the Tula virus. The first clade of hantaviruses associated with Neotominae occupied a stem lineage, especially those that emerged in Central America or Mexico. Hantaviruses from Central America and Mexico found in Neotominae rodents spread northward and probably gave rise to the first phylogroup of hantaviruses associated with Sigmodontinae in North America. Two preferential host-switching transmissions in hantaviruses apparently gave rise to two different paraphyletic group in Neotominae and Sigmodontinae. Our study supports a probable epicenter of diversification in Central America and/or Mexico for hantaviruses related to both the Neotominae and Sigmodontinae subfamilies.

Incorporating genomic methods into contact networks to reveal new insights into animal behavior and infectious disease dynamics

Utilization of contact networks has provided opportunities for assessing the dynamic interplay between pathogen transmission and host behavior. Genomic techniques have, in their own right, provided new insight into complex questions in disease ecology, and the increasing accessibility of genomic approaches means more researchers may seek out these tools. The integration of network and genomic approaches provides opportunities to examine the interaction between behavior and pathogen transmission in new ways and with greater resolution. While a number of studies have begun to incorporate both contact network and genomic approaches, a great deal of work has yet to be done to better integrate these techniques. In this review, we give a broad overview of how network and genomic approaches have each been used to address questions regarding the interaction of social behavior and infectious disease, and then discuss current work and future horizons for the merging of these techniques.

Pathogens in space: Advancing understanding of pathogen dynamics and disease ecology through landscape genetics

Landscape genetics has provided many insights into how heterogeneous landscape features drive processes influencing spatial genetic variation in free-living organisms. This rapidly developing field has focused heavily on vertebrates, and expansion of this scope to the study of infectious diseases holds great potential for landscape geneticists and disease ecologists alike. The potential application of landscape genetics to infectious agents has garnered attention at formative stages in the development of landscape genetics, but systematic examination is lacking. We comprehensively review how landscape genetics is being used to better understand pathogen dynamics. We characterize the field and evaluate the types of questions addressed, approaches used and systems studied. We also review the now established landscape genetic methods and their realized and potential applications to disease ecology. Lastly, we identify emerging frontiers in the landscape genetic study of infectious agents, including recent phylogeographic approaches and frameworks for studying complex multihost and host-vector systems. Our review emphasizes the expanding utility of landscape genetic methods available for elucidating key pathogen dynamics (particularly transmission and spread) and also how landscape genetic studies of pathogens can provide insight into host population dynamics. Through this review, we convey how increasing awareness of the complementarity of landscape genetics and disease ecology among practitioners of each field promises to drive important cross-disciplinary advances.

Viruses as indicators of contemporary host dispersal and phylogeography: an example of feline immunodeficiency virus (FIVP le ) in free-ranging African lion (Panthera leo)

Measuring contemporary dispersal in highly mobile terrestrial species is challenging, especially when species are characterized by low levels of population differentiation. Directly transmitted viruses can be used as a surrogate for traditional methods of tracking host movement. Feline immunodeficiency virus (FIV) is a species-specific lentivirus, which has an exceptionally high mutation rate and circulates naturally in wild felids. Using samples derived from 35 lion (Panthera leo) prides, we tested the prediction that FIV in lions (FIV_{P le} ) can be used to track the dispersal of individuals between prides. As FIV_{P le} subtypes are geographically structured throughout Africa, we predicted that this marker could be used to detect phylogeographic structure of lions at smaller spatial scales. Phylogenetic analyses of FIV_{P le} pol-RT sequences showed that core pride members (females and subadults) shared evolutionary close viral lineages which differed from neighbouring core prides, whereas sequences from sexually mature males associated with the same pride were always the most divergent. In six instances, natal pride associations of divergent male lions could be inferred, on the assumption that FIV_{P le} infections are acquired during early life stages. Congruence between the genetic pattern of FIV and pride structure suggests that vertical transmission plays an important role in lion FIV dynamics. At a fine spatial scale, significant viral geographic structuring was also detected between lions occurring north of the Olifants River within the Kruger National Park (KNP) and those occupying the southern and central regions. This pattern was further supported by phylogenetic analyses and the confinement of FIV_{P le} subtype E to the northern region of KNP. The study provides new insights into the use of retroviral sequences to predict host dispersal and fine-scale contemporary geographic structure in a social felid species.

High prevalence of Lynx rufus gammaherpesvirus 1 in wild Vermont bobcats

Gammaherpesviruses (GHVs) are host specific DNA viruses that infect a large range of mammalian species. These viruses preferentially target host lymphocyte cell populations and infection may lead to morbidity or mortality in immunocompromised, co-infected, or non-adapted hosts. In this study, we tested for the presence of Lynx rufus gammaherpesvirus 1 (LruGHV1) in a northeastern United States population of wild bobcats (L. rufus). We estimated prevalence of infection and viral load in infected individuals using quantitative real-time PCR analysis of spleen DNA from 64 Vermont bobcats. We observed an overall prevalence of 64% using this methodology. Bobcat age was significantly positively associated with GHV infection status, and we noted a trend for higher viral loads in young animals, but prevalence and viral load were similar in male and female bobcats. A single LruGHV1 variant was identified from the sequencing of the viral glycoprotein B gene of Vermont bobcats. This gene sequence was 100% similar to that reported in Florida bobcats and slightly variant from other isolates identified in the Western USA. Our work suggests broad geographic distribution and high prevalence of LruGHV1 in bobcat populations across the United States with infection attributes that suggest horizontal transmission of the agent. Geographic differences in viral genotype may reflect historical migration and expansion events among bobcat populations.

Evolution of life in urban environments

Our planet is an increasingly urbanized landscape, with over half of the human population residing in cities. Despite advances in urban ecology, we do not adequately understand how urbanization affects the evolution of organisms, nor how this evolution may affect ecosystems and human health. Here, we review evidence for the effects of urbanization on the evolution of microbes, plants, and animals that inhabit cities. Urbanization affects adaptive and nonadaptive evolutionary processes that shape the genetic diversity within and between populations. Rapid adaptation has facilitated the success of some native species in urban areas, but it has also allowed human pests and disease to spread more rapidly. The nascent field of urban evolution brings together efforts to understand evolution in response to environmental change while developing new hypotheses concerning adaptation to urban infrastructure and human socioeconomic activity. The next generation of research on urban evolution will provide critical insight into the importance of evolution for sustainable interactions between humans and our city environments.

Functional connectivity in replicated urban landscapes in the land snail (Cornu aspersum)

Urban areas are highly fragmented and thereby exert strong constraints on individual dispersal. Despite this, some species manage to persist in urban areas, such as the garden snail, Cornu aspersum, which is common in cityscapes despite its low mobility. Using landscape genetic approaches, we combined study area replication and multiscale analysis to determine how landscape composition, configuration and connectivity influence snail dispersal across urban areas. At the overall landscape scale, areas with a high percentage of roads decreased genetic differentiation between populations. At the population scale, genetic differentiation was positively linked with building surface, the proportion of borders where wooded patches and roads appeared side by side and the proportion of borders combining wooded patches and other impervious areas. Analyses based on pairwise genetic distances validated the isolation-by-distance and isolation-by-resistance models for this land snail, with an equal fit to least-cost paths and circuit-theory-based models. Each of the 12 landscapes analysed separately yielded specific relations to environmental features, whereas analyses integrating all replicates highlighted general common effects. Our results suggest that urban transport infrastructures facilitate passive snail dispersal. At a local scale, corresponding to active dispersal, unfavourable habitats (wooded and impervious areas) isolate populations. This work upholds the use of replicated landscapes to increase the generalizability of landscape genetics results and shows how multiscale analyses provide insight into scale-dependent processes.

Towards an eco-phylogenetic framework for infectious disease ecology

Identifying patterns and drivers of infectious disease dynamics across multiple scales is a fundamental challenge for modern science. There is growing awareness that it is necessary to incorporate multi-host and/or multi-parasite interactions to understand and predict current and future disease threats better, and new tools are needed to help address this task. Eco-phylogenetics (phylogenetic community ecology) provides one avenue for exploring multi-host multi-parasite systems, yet the incorporation of eco-phylogenetic concepts and methods into studies of host pathogen dynamics has lagged behind. Eco-phylogenetics is a transformative approach that uses evolutionary history to infer present-day dynamics. Here, we present an eco-phylogenetic framework to reveal insights into parasite communities and infectious disease dynamics across spatial and temporal scales. We illustrate how eco-phylogenetic methods can help untangle the mechanisms of host-parasite dynamics from individual (e.g. co-infection) to landscape scales (e.g. parasite/host community structure). An improved ecological understanding of multi-host and multi-pathogen dynamics across scales will increase our ability to predict disease threats.

Urban landscapes can change virus gene flow and evolution in a fragmentation-sensitive carnivore

Urban expansion has widespread impacts on wildlife species globally, including the transmission and emergence of infectious diseases. However, there is almost no information about how urban landscapes shape transmission dynamics in wildlife. Using an innovative phylodynamic approach combining host and pathogen molecular data with landscape characteristics and host traits, we untangle the complex factors that drive transmission networks of feline immunodeficiency virus (FIV) in bobcats (Lynx rufus). We found that the urban landscape played a significant role in shaping FIV transmission. Even though bobcats were often trapped within the urban matrix, FIV transmission events were more likely to occur in areas with more natural habitat elements. Urban fragmentation also resulted in lower rates of pathogen evolution, possibly owing to a narrower range of host genotypes in the fragmented area. Combined, our findings show that urban landscapes can have impacts on a pathogen and its evolution in a carnivore living in one of the most fragmented and urban systems in North America. The analytical approach used here can be broadly applied to other host-pathogen systems, including humans.

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.

Applying landscape genetics to the microbial world

Landscape genetics, which explicitly quantifies landscape effects on gene flow and adaptation, has largely focused on macroorganisms, with little attention given to microorganisms. This is despite overwhelming evidence that microorganisms exhibit spatial genetic structuring in relation to environmental variables. The increasing accessibility of genomic data has opened up the opportunity for landscape genetics to embrace the world of microorganisms, which may be thought of as 'the invisible regulators' of the macroecological world. Recent developments in bioinformatics and increased data accessibility have accelerated our ability to identify microbial taxa and characterize their genetic diversity. However, the influence of the landscape matrix and dynamic environmental factors on microorganism genetic dispersal and adaptation has been little explored. Also, because many microorganisms coinhabit or codisperse with macroorganisms, landscape genomic approaches may improve insights into how micro- and macroorganisms reciprocally interact to create spatial genetic structure. Conducting landscape genetic analyses on microorganisms requires that we accommodate shifts in spatial and temporal scales, presenting new conceptual and methodological challenges not yet explored in 'macro'-landscape genetics. We argue that there is much value to be gained for microbial ecologists from embracing landscape genetic approaches. We provide a case for integrating landscape genetic methods into microecological studies and discuss specific considerations associated with the novel challenges this brings. We anticipate that microorganism landscape genetic studies will provide new insights into both micro- and macroecological processes and expand our knowledge of species' distributions, adaptive mechanisms and species' interactions in changing environments.

Dental Pathology of the California Bobcat (Lynx rufus californicus)

Skulls from 277 California bobcats (Lynx rufus californicus) were examined macroscopically and by radiography. The majority of the skulls were from adult animals (79.8%). The skulls were from 128 male (46.2%) and 114 female (41.2%) animals and gender was unknown for the remainder. The majority (95.6%) of teeth were present for examination. Only 16 teeth were identified as absent congenitally and 15 of these were incisor teeth. Teeth with abnormal morphology were rare (0.5%). The two most common abnormalities were unusually large crowns of the maxillary first molar teeth and bigemination of the mandibular incisor teeth. Teeth with an abnormal number of roots were uncommon (n = 68). Sixty-three teeth had abnormal roots, mostly the presence of two roots instead of one for the maxillary first molar tooth. The most prevalent dental lesions found in the California bobcat were attrition/abrasion (85.2%), periodontitis (56.0%) and tooth fractures (50.9%). Less common dental lesions were endodontal disease (n = 114 teeth) and tooth resorption (n = 73 teeth).

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.

Population and genetic outcomes 20 years after reintroducing bobcats (Lynx rufus) to Cumberland Island, Georgia USA

In 1988-1989, 32 bobcats Lynx rufus were reintroduced to Cumberland Island (CUIS), Georgia, USA, from which they had previously been extirpated. They were monitored intensively for 3 years immediately post-reintroduction, but no estimation of the size or genetic diversity of the population had been conducted in over 20 years since reintroduction. We returned to CUIS in 2012 to estimate abundance and effective population size of the present-day population, as well as to quantify genetic diversity and inbreeding. We amplified 12 nuclear microsatellite loci from DNA isolated from scats to establish genetic profiles to identify individuals. We used spatially explicit capture-recapture population estimation to estimate abundance. From nine unique genetic profiles, we estimate a population size of 14.4 (SE = 3.052) bobcats, with an effective population size (N e) of 5-8 breeding individuals. This is consistent with predictions of a population viability analysis conducted at the time of reintroduction, which estimated the population would average 12-13 bobcats after 10 years. We identified several pairs of related bobcats (parent-offspring and full siblings), but ~75% of the pairwise comparisons were typical of unrelated individuals, and only one individual appeared inbred. Despite the small population size and other indications that it has likely experienced a genetic bottleneck, levels of genetic diversity in the CUIS bobcat population remain high compared to other mammalian carnivores. The reintroduction of bobcats to CUIS provides an opportunity to study changes in genetic diversity in an insular population without risk to this common species. Opportunities for natural immigration to the island are limited; therefore, continued monitoring and supplemental bobcat reintroductions could be used to evaluate the effect of different management strategies to maintain genetic diversity and population viability. The successful reintroduction and maintenance of a bobcat population on CUIS illustrates the suitability of translocation as a management tool for re-establishing felid populations.

Origin and genetic structure of a recovering bobcat (Lynx rufus) population

Genetic analyses can provide important insights into the demographic processes that underlie recovering populations of mammals of conservation concern such as felid species. To better understand the recent and rapid recovery of bobcats (Lynx rufus (Schreber, 1777)) in Ohio, we analyzed samples from four states in the lower Great Lakes Region using 12 microsatellite DNA loci and a portion of the mtDNA control region. Our results showed that a newly established population of bobcats in the eastern part of Ohio was genetically distinct from a multistate population distributed across Kentucky, southern Ohio, West Virginia, and western Pennsylvania. There was no direct genetic evidence of a bottleneck or inbreeding in this population. A lack of private alleles and only slightly lower levels of allelic richness and heterozygosity compared with its neighbors suggest that the eastern Ohio population likely originated from the migration of relatively large numbers of individuals from a source population rather than re-emerging from an undetected residual population. We recommend that a management plan should define the areas occupied by the two populations in Ohio as separate management units at least for the near future.

Infectious disease transmission and contact networks in wildlife and livestock

The use of social and contact networks to answer basic and applied questions about infectious disease transmission in wildlife and livestock is receiving increased attention. Through social network analysis, we understand that wild animal and livestock populations, including farmed fish and poultry, often have a heterogeneous contact structure owing to social structure or trade networks. Network modelling is a flexible tool used to capture the heterogeneous contacts of a population in order to test hypotheses about the mechanisms of disease transmission, simulate and predict disease spread, and test disease control strategies. This review highlights how to use animal contact data, including social networks, for network modelling, and emphasizes that researchers should have a pathogen of interest in mind before collecting or using contact data. This paper describes the rising popularity of network approaches for understanding transmission dynamics in wild animal and livestock populations; discusses the common mismatch between contact networks as measured in animal behaviour and relevant parasites to match those networks; and highlights knowledge gaps in how to collect and analyse contact data. Opportunities for the future include increased attention to experiments, pathogen genetic markers and novel computational tools.

Disease and freeways drive genetic change in urban bobcat populations

Urbanization profoundly impacts animal populations by causing isolation, increased susceptibility to disease, and exposure to toxicants. Genetic effects include reduced effective population size, increased population substructure, and decreased adaptive potential. We investigated the influence that urbanization and a disease epizootic had on the population genetics of bobcats (Lynx rufus) distributed across a highly fragmented urban landscape. We genotyped more than 300 bobcats, sampled from 1996 to 2012, for variation at nine neutral and seven immune gene-linked microsatellite loci. We found that two freeways are significant barriers to gene flow. Further, a 3-year disease epizootic, associated with secondary anticoagulant rodenticide exposure, caused a population bottleneck that led to significant genetic differentiation between pre- and post-disease populations that was greater than that between populations separated by major freeways for >60 years. However, balancing selection acted on immune-linked loci during the epizootic, maintaining variation at functional regions. Conservation assessments need to assay loci that are potentially under selection to better preserve the adaptive potential of populations at the urban–wildland interface. Further, interconnected regions that contain appropriate habitat for wildlife will be critical to the long-term viability of animal populations in urban landscapes.

Applying evolutionary concepts to wildlife disease ecology and management

Existing and emerging infectious diseases are among the most pressing global threats to biodiversity, food safety and human health. The complex interplay between host, pathogen and environment creates a challenge for conserving species, communities and ecosystem functions, while mediating the many known ecological and socio-economic negative effects of disease. Despite the clear ecological and evolutionary contexts of host-pathogen dynamics, approaches to managing wildlife disease remain predominantly reactionary, focusing on surveillance and some attempts at eradication. A few exceptional studies have heeded recent calls for better integration of ecological concepts in the study and management of wildlife disease; however, evolutionary concepts remain underused. Applied evolution consists of four principles: evolutionary history, genetic and phenotypic variation, selection and eco-evolutionary dynamics. In this article, we first update a classical framework for understanding wildlife disease to integrate better these principles. Within this framework, we explore the evolutionary implications of environment-disease interactions. Subsequently, we synthesize areas where applied evolution can be employed in wildlife disease management. Finally, we discuss some future directions and challenges. Here, we underscore that despite some evolutionary principles currently playing an important role in our understanding of disease in wild animals, considerable opportunities remain for fostering the practice of evolutionarily enlightened wildlife disease management.

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.