'All I do is fight fires': Qualitative exploration of UK veterinarians' attitudes towards and experiences of pre-purchase consultations regarding brachycephalic dogs.

BACKGROUND: Demand for brachycephalic dogs has dramatically increased over the past decade, despite growing evidence of strong associations between brachycephaly and severe, chronic disease. Positive post-purchase attitudes of owners towards brachycephalic dogs are highly intractable, facilitated by normalisation of breed-related disease and strong dog-owner bonds. Veterinarian-led pre-purchase consultations (PPCs) have been promoted as an opportunity to influence acquisition decisions regarding brachycephalic dogs. METHODS: Using reflexive thematic analysis, this qualitative study investigated veterinarians' attitudes towards, and experiences of, engaging in PPCs for brachycephalic dogs. RESULTS: Veterinarians described significant structural barriers, such as lack of time and staff, and perceptual barriers, such as perceived futility, negative outcomes of PPCs and public distrust in the profession, that combine to inhibit many veterinarians from offering or delivering brachycephalic PPCs. Some of these barriers are intractable at an individual-veterinarian level and lead to moral distress due to conflicts between veterinarians' responsibilities to animal welfare versus clients' and business needs, with resultant censorship of authentic views on brachycephalic ownership compromising professional integrity and autonomy. LIMITATIONS: Wider views and/or experiences may exist outside of the sampled population. CONCLUSION: Stronger veterinary leadership, including from the RCVS, BVA and veterinary corporates, is now urgently required to build profession-wide solutions that overcome the barriers identified to create public-facing consensus and innovative solutions to the 'brachycephalic crisis'.

A phylogenetic approach identifies patterns of beta diversity and floristic subregions of the Qinghai-Tibet Plateau.

Patterns of taxonomic and phylogenetic beta diversity and their relationships with environmental correlates can help reveal the origin and evolutionary history of regional biota. The Qinghai-Tibet Plateau (QTP) harbors an exceptionally diverse flora, however, a phylogenetic perspective has rarely been used to investigate its beta diversity and floristic regions. In this study, we used a phylogenetic approach to identify patterns of beta diversity and quantitatively delimit floristic regions on the Qinghai-Tibet Plateau. We also examined the relationships between multifaceted beta diversity, geographical distance, and climatic difference, and evaluated the relative importance of various factors (i.e., climate, topography and history) in shaping patterns of beta diversity. Sørensen dissimilarity indices indicated that patterns of species turnover among sites dominated the QTP. We also found that patterns of both taxonomic and phylogenetic beta diversity were significantly related to geographical distance and climatic difference. The environmental factors that contributed most to these patterns of beta diversity include annual precipitation, mean annual temperature, climatic gradients and climatic instability. Hierarchical dendrograms of dissimilarity and non-metric multidimensional scaling ordination based on phylogenetic beta diversity data identified ten floristic subregions in the QTP. Our results suggest that the contemporary environment and historical climate changes have filtered species composition among sites and eventually determined beta diversity patterns of plants in the QTP.

Biotic homogenisation and differentiation as directional change in beta diversity: synthesising driver-response relationships to develop conceptual models across ecosystems.

Biotic homogenisation is defined as decreasing dissimilarity among ecological assemblages sampled within a given spatial area over time. Biotic differentiation, in turn, is defined as increasing dissimilarity over time. Overall, changes in the spatial dissimilarities among assemblages (termed 'beta diversity') is an increasingly recognised feature of broader biodiversity change in the Anthropocene. Empirical evidence of biotic homogenisation and biotic differentiation remains scattered across different ecosystems. Most meta-analyses quantify the prevalence and direction of change in beta diversity, rather than attempting to identify underlying ecological drivers of such changes. By conceptualising the mechanisms that contribute to decreasing or increasing dissimilarity in the composition of ecological assemblages across space, environmental managers and conservation practitioners can make informed decisions about what interventions may be required to sustain biodiversity and can predict potential biodiversity outcomes of future disturbances. We systematically reviewed and synthesised published empirical evidence for ecological drivers of biotic homogenisation and differentiation across terrestrial, marine, and freshwater realms to derive conceptual models that explain changes in spatial beta diversity. We pursued five key themes in our review: (i) temporal environmental change; (ii) disturbance regime; (iii) connectivity alteration and species redistribution; (iv) habitat change; and (v) biotic and trophic interactions. Our first conceptual model highlights how biotic homogenisation and differentiation can occur as a function of changes in local (alpha) diversity or regional (gamma) diversity, independently of species invasions and losses due to changes in species occurrence among assemblages. Second, the direction and magnitude of change in beta diversity depends on the interaction between spatial variation (patchiness) and temporal variation (synchronicity) of disturbance events. Third, in the context of connectivity and species redistribution, divergent beta diversity outcomes occur as different species have different dispersal characteristics, and the magnitude of beta diversity change associated with species invasions also depends strongly on alpha and gamma diversity prior to species invasion. Fourth, beta diversity is positively linked with spatial environmental variability, such that biotic homogenisation and differentiation occur when environmental heterogeneity decreases or increases, respectively. Fifth, species interactions can influence beta diversity via habitat modification, disease, consumption (trophic dynamics), competition, and by altering ecosystem productivity. Our synthesis highlights the multitude of mechanisms that cause assemblages to be more or less spatially similar in composition (taxonomically, functionally, phylogenetically) through time. We consider that future studies should aim to enhance our collective understanding of ecological systems by clarifying the underlying mechanisms driving homogenisation or differentiation, rather than focusing only on reporting the prevalence and direction of change in beta diversity, per se.

Two dominant forms of multisite similarity decline - Their origins and interpretation.

The number of species shared by two or more sites is a fundamental measure of spatial variation in species composition. As more sites are included in the comparison of species composition, the average number of species shared across them declines, with a rate increasingly dependent on only the most widespread species. In over 80% of empirical communities, models of decline in shared species across multiple sites (multisite similarity decline) follow one of two distinct forms. An exponential form is assumed to reflect stochastic assembly and a power law form niche-based sorting, yet these explanations are largely untested, and little is known of how the two forms arise in nature. Using simulations, we first show that the distribution of the most widespread species largely differentiates the two forms, with the power law increasingly favored where such species occupy more than ~75% of sites. We reasoned the less cosmopolitan distribution of widespread species within exponential communities would manifest as differences in community biodiversity properties, specifically more aggregated within-species distributions, less even relative abundance distributions, and weaker between-species spatial associations. We tested and largely confirmed these relationships using 80 empirical datasets, suggesting that the form of multisite similarity decline offers a basis to predict how landscape-scale loss or gain of widespread species is reflected in different local-scale community structures. Such understanding could, for example, be used to predict changes in local-scale competitive interactions following shifts in widespread species' distributions. We propose multiple explanations for the origin of exponential decline, including high among-site abiotic variation, sampling of highly specialized (narrow niche width) taxa, and strong dispersal limitation. We recommend these are evaluated as alternative hypotheses to stochastic assembly.

Species accumulation in small-large vs large-small order: more species but not all species?

Although groups of small habitat patches often support more species than large patches of equal total area, their biodiversity value remains controversial. An important line of evidence in this debate compares species accumulation curves, where patches are ordered from small-large and large-small (aka 'SLOSS analysis'). However, this method counts species equally and is unable to distinguish patch size dependence in species' occupancies. Moreover, because of the species-area relationship, richness differences typically only contribute to accumulation in small-large order, maximizing the probability of adding species in this direction. Using a null model to control for this, I tested 202 published datasets from archipelagos, habitat islands and fragments for patch size dependence in species accumulation and compared conclusions regarding relative species accumulation with SLOSS analysis. Relative to null model expectations, species accumulation was on average 2.7% higher in large-small than small-large order. The effect was strongest in archipelagos (5%), intermediate for fragments (1.5%) and smallest for habitat islands (1.1%). There was no difference in effect size among taxonomic groups, but each shared this same trend. Results suggest most meta-communities include species that either prefer, or depend upon, larger habitat patches. Relative to SLOSS analysis, null models found lower frequency of greater small-patch importance for species representation (e.g., for fragments: 69 vs 16% respectively) and increased frequency for large patches (fragments: 3 vs 25%). I suggest SLOSS analysis provides unreliable inference on species accumulation and the outcome largely depends on island species-area relationships, not the relative diversity value of small vs large patches.

The Effects of Species Abundance, Spatial Distribution, and Phylogeny on a Plant-Ectomycorrhizal Fungal Network.

Plant and root fungal interactions are among the most important belowground ecological interactions, however, the mechanisms underlying pairwise interactions and network patterns of rhizosphere fungi and host plants remain unknown. We tested whether neutral process or spatial constraints individually or jointly best explained quantitative plant-ectomycorrhizal fungal network assembly in a subtropical forest in southern China. Results showed that the observed plant-ectomycorrhizal fungal network had low connectivity, high interaction evenness, and an intermediate level of specialization, with nestedness and modularity both greater than random expectation. Incorporating information on the relative abundance and spatial overlap of plants and fungi well predicted network nestedness and connectance, but not necessarily explained other network metrics such as specificity. Spatial overlap better predicted pairwise species interactions of plants and ectomycorrhizal fungi than species abundance or a combination of species abundance and spatial overlap. There was a significant phylogenetic signal on species degree and interaction strength for ectomycorrhizal fungal but not for plant species. Our study suggests that neutral processes (species abundance matching) and niche/dispersal-related processes (implied by spatial overlap and phylogeny) jointly drive the shaping of a plant-ectomycorrhizal fungal network.

Rare, common, alien and native species follow different rules in an understory plant community.

Biological invasions are a leading threat to biodiversity globally. Increasingly, ecosystems experience multiple introductions, which can have significant effects on patterns of diversity. The way these communities assemble will depend partly on whether rare and common alien species respond to environmental predictors in the same manner as rare and common native species, but this is not well understood. To examine this question across four national parks in south-eastern Australia, we sampled the understory plant community of eucalypt-dominated dry forest subject to multiple plant introductions. The drivers of diversity and turnover in alien and native species of contrasting frequency of occurrence (low, intermediate, and high) were each tested individually. We found alien species diversity and turnover were both strongly associated with abiotic conditions (e.g., soil pH), while distance had little influence because of the greater extent of occurrence and more homogeneous composition of common aliens. In contrast, native species diversity was not associated with abiotic conditions and their turnover was as strongly influenced by distance as by abiotic conditions. In both alien and native species, however, the most important predictors of turnover changed with frequency of occurrence. Although local coexistence appears to be facilitated by life history trade-offs, species richness of aliens and natives was negatively correlated and native species might face greater competition in areas with more neutral soils (e.g., pH > ~5.5) where alien richness and relative frequency were both highest. We conclude that diversity and turnover in the generally more widespread alien species are mainly driven by species sorting along an environmental gradient associated with pH and nutrient availability, whereas turnover of native species is driven by more neutral processes associated with dispersal limitation. We show alien and native plant species respond to different environmental factors, as do rare and common species within each component.

Chronic Bilateral Symmetric Anterior Shoulder Dislocation Secondary to Seizures in Chikungunya Encephalitis.

Bilateral shoulder dislocations are a rare occurrence and can be categorized as either symmetric (both humeral heads dislocate in the same direction) or asymmetric (wherein the humeral heads dislocate in different directions). Shoulder dislocations may be overlooked if they are the result of systemic injury; if diagnosed >21 days after occurring, they are considered chronic dislocations. We describe the case of a 31-year-old male who presented with an eight-week history of bilateral shoulder pain. His onset of pain coincided with a seizure secondary to Chikungunya encephalitis. Clinical and radiological examination demonstrated bilateral symmetric anterior shoulder dislocation with associated greater tuberosity fractures and extensive callus formation bilaterally. Open surgical management was performed first on the left shoulder via the deltopectoral approach. The callus was removed, the greater tuberosity fragment lifted off, reattached to the original position, and held in place with sutures and proximal humeral locking plates. The right shoulder was reduced six weeks after the left shoulder due to patient preference; the reduction utilized the same approach as with the left shoulder. Post-operatively the patient was immobilized, and physiotherapy commenced. He achieved a satisfactory range of motion four months post-operation. Physicians should be cognizant that shoulder pain after a convulsive seizure may signify shoulder dislocation. Thorough clinical and radiological examinations are warranted in such an instance. There exists no consensus on the treatment of chronic shoulder dislocations, but it is recommended that closed reduction only be attempted up to six weeks post-dislocation due to the high risk of iatrogenic fractures and neurovascular damage beyond this time.

Direct effects of selection on aboveground biomass contrast with indirect structure-mediated effects of complementarity in a subtropical forest.

Understanding the multiple biotic and abiotic controls of aboveground biomass (AGB) is important for projecting the consequences of global change and to effectively manage carbon storage. Although large-scale studies have identified the major environmental and biological controls of AGB, drivers of local-scale variation are less well known. Additionally, involvement of multiple causal paths and scale dependence in effect sizes potentially confounds comparisons among studies differing in methodology and sampling grain. We tested for scale dependence in evidence supporting selection, complementarity and environmental factors as the main determinants of AGB variation over a 50 ha study extent in subtropical China, modelling this at four sampling grains (0.01, 0.04, 0.25 and 1 ha). At each grain, we used piecewise structural equation models to quantify the direct and indirect effects of environmental (topographic and edaphic properties) and forest attributes (structure, diversity and functional traits) on AGB, while controlling for spatial autocorrelation. Direct scale-invariant effects on AGB were evident for structure and community-mean traits, supporting dominance of selection effects. However, diversity had strong indirect effects on AGB via forest structure, particularly at larger sampling grains (≥ 0.25 ha), while direct effects only emerged at the smallest grain size (0.01 ha). The direct and indirect effects of edaphic and topographic factors were also important for explaining both forest attributes and AGB across all scales. Although selection effects appeared to be more influential on ecosystem function, ignoring indirect causal pathways for diversity via structural attributes risks overlooking the importance of complementarity on ecosystem functioning, particularly as sampling grain increases.

Effects of host phylogeny, habitat and spatial proximity on host specificity and diversity of pathogenic and mycorrhizal fungi in a subtropical forest.

Soil plant-pathogenic (PF) and mycorrhizal fungi (MF) are both important in maintaining plant diversity, for example via host-specialized effects. However, empirical knowledge on the degree of host specificity and possible factors affecting the fungal assemblages is lacking. We identified PF and MF in fine roots of 519 individuals across 45 subtropical tree species in southern China in order to quantify the importance of host phylogeny (including via its effects on functional traits), habitat and space in determining fungal communities. We also compared host specificity in PF and MF at different host-phylogenetic scales. In both PF and MF, host phylogeny independently accounted for > 19% of the variation in fungal richness and composition, whereas environmental and spatial factors each explained no more than 4% of the variation. Over 77% of the variation explained by phylogeny was attributable to covariation in plant functional traits. Host specificity was phylogenetically scale-dependent, being stronger in PF than in MF at low host-phylogenetic scales (e.g. within genus) but similar at larger scales. Our study suggests that host-phylogenetic effects dominate the assembly of both PF and MF communities, resulting from phylogenetically clustered plant traits. The scale-dependent host specificity implies that PF were specialized at lower-level and MF at higher-level host taxa.

Loss of only the smallest patches will reduce species diversity in most discrete habitat networks.

Under many global-change scenarios, small habitat patches are the most vulnerable to destruction. For example, smaller ponds are at greater risk in a drying climate and their loss would remove any obligate aquatic individuals present. We asked what proportional loss of species diversity from metacommunities comprised of discrete habitat patches should be expected from attrition (complete loss) of only the smallest patches under such a premise. We analyzed 175 published datasets for different taxonomic groups (vertebrates, invertebrates, and plants) and habitat types (islands, habitat islands, and fragments). We simulated the destruction of only the smallest patches to an approximate 20% of total area (range: 15.2%-24.2%) and analyzed species loss. Mean [± 95% CI] species loss was 12.7% [10.8, 14.6], although 18.3% of datasets lost no species. Four broad patterns of species loss were evident, reflecting underlying differences in minimum area requirements and the degree of species turnover among patches. Regression modeling showed species loss increased with greater species turnover among patches (ßSIM ) and decreased with greater area scaling of diversity (i.e., larger power-law island species-area relationship exponents). Losses also increased with greater numbers of single-patch endemics and with increasing proportions of patches destroyed. After accounting for these predictors, neither taxonomic group nor habitat type increased explained variation in species loss. Attrition of the smallest patches removed species in >80% of metacommunities, despite all larger patches and >75% of total area remaining intact. At both 10% and 20% area reduction, median species loss across all datasets was around 50% higher than predicted from methods based on the species-area relationship. We conclude that any mechanism of global change that selectively destroys small habitat patches will lead to imminent extinctions in most discrete metacommunities.

Hydrological-niche models predict water plant functional group distributions in diverse wetland types.

Human use of water resources threatens environmental water supplies. If resource managers are to develop policies that avoid unacceptable ecological impacts, some means to predict ecosystem response to changes in water availability is necessary. This is difficult to achieve at spatial scales relevant for water resource management because of the high natural variability in ecosystem hydrology and ecology. Water plant functional groups classify species with similar hydrological niche preferences together, allowing a qualitative means to generalize community responses to changes in hydrology. We tested the potential for functional groups in making quantitative prediction of water plant functional group distributions across diverse wetland types over a large geographical extent. We sampled wetlands covering a broad range of hydrogeomorphic and salinity conditions in South Australia, collecting both hydrological and floristic data from 687 quadrats across 28 wetland hydrological gradients. We built hydrological-niche models for eight water plant functional groups using a range of candidate models combining different surface inundation metrics. We then tested the predictive performance of top-ranked individual and averaged models for each functional group. Cross validation showed that models achieved acceptable predictive performance, with correct classification rates in the range 0.68-0.95. Model predictions can be made at any spatial scale that hydrological data are available and could be implemented in a geographical information system. We show the response of water plant functional groups to inundation is consistent enough across diverse wetland types to quantify the probability of hydrological impacts over regional spatial scales.

Alcelaphine herpesvirus 1 glycoprotein B: recombinant expression and antibody recognition.

The gammaherpesvirus alcelaphine herpesvirus 1 (AlHV-1) causes fatal malignant catarrhal fever (MCF) in susceptible species including cattle, but infects its reservoir host, wildebeest, without causing disease. Pathology in cattle may be influenced by virus-host cell interactions mediated by the virus glycoproteins. Cloning and expression of a haemagglutinin-tagged version of the AlHV-1 glycoprotein B (gB) was used to demonstrate that the AlHV-1-specific monoclonal antibody 12B5 recognised gB and that gB was the main component of the gp115 complex of AlHV-1, a glycoprotein complex of five components identified on the surface of AlHV-1 by immunoprecipitation and radiolabelling. Analysis of AlHV-1 virus particles showed that the native form of gB was detected by mAb 12B5 as a band of about 70 kDa, whilst recombinant gB expressed by transfected HEK293T cells appeared to be subject to additional cleavage and incomplete post-translational processing. Antibody 12B5 recognised an epitope on the N-terminal furin-cleaved fragment of gB on AlHV-1 virus particles. It could be used to detect recombinant and virus-expressed gB on western blots and on the surface of infected cells by flow cytometry, whilst recombinant gB was detected on the surface of transfected cells by immunofluorescence. Recombinant gB has potential as an antigen for ELISA detection of MCF virus infection and as a candidate vaccine antigen.

Assessment of the rabbit as a wildlife reservoir of bovine viral diarrhea virus: serological analysis and generation of trans-placentally infected offspring.

Eradication of bovine viral diarrhea virus (BVDV) is ongoing in many European countries and is based on removal of persistently infected (PI) cattle. In this context, low-level risks, including alternative reservoirs of infection, may become more important as the number of BVDV-free herds increases. Alternative reservoirs include livestock, such as sheep and goats, as well as wildlife, including deer and rabbits. Due to the extensive nature of the beef industry in Scotland, where an eradication program started in 2010, contact between cattle and alternative reservoir hosts is common. Seroprevalence to BVDV in rabbit populations can be high. In addition, rabbits can be infected with BVDV by natural routes, indicating that they could be a wildlife reservoir of infection. We analyzed the potential risk to livestock from rabbit populations in the UK by two approaches. First, ∼260 serum samples from free-ranging wild rabbits in Scotland and northern England were tested for BVDV-specific antibodies by ELISA. Only three samples exhibited low level BVDV-specific reactivity, suggesting that BVDV infection of rabbits was not frequent. Second, rabbits were challenged with BVDV at day 7 or 12 of pregnancy. This did not lead to any clinical signs in the infected animals or obvious increases in abortion or stillbirth in the infected dams. Samples from the dams, placental material and ∼130 offspring were tested by BVDV-specific RT-PCR and antibody ELISA. Positive PCR results in the placentas and in the tissues and body fluids of rabbits up to 10 days old showed that trans-placental infection of rabbits with BVDV had occurred. Many of the offspring had BVDV-specific antibodies. These data support the view that a wildlife reservoir of BVDV in rabbit poses a small but non-zero risk of re-infection for BVDV-free cattle herds. Rabbits are susceptible to infection with BVDV but only a small proportion of free-living rabbits in the UK appear to have been infected.

Identification of immuno-reactive capsid proteins of malignant catarrhal fever viruses.

Malignant catarrhal fever (MCF) is a fatal disease of cattle and other ungulates caused by certain gamma-herpesviruses including alcelaphine herpesvirus-1 (AlHV-1) and ovine herpesvirus-2 (OvHV-2). An attenuated virus vaccine based on AlHV-1 has been shown to induce virus-neutralising antibodies in plasma and nasal secretions of protected cattle but the targets of virus-specific antibodies are unknown. Proteomic analysis and western blotting of virus extracts allowed the identification of eight candidate AlHV-1 virion antigens. Recombinant expression of selected candidates and their OvHV-2 orthologues confirmed that two polypeptides, the products of the ORF17.5 and ORF65 genes, were antigens recognised by antibodies from natural MCF cases or from AlHV-1 vaccinated cattle. These proteins have potential as diagnostic and/or vaccine antigens.

The effect of the TLR9 ligand CpG-oligodeoxynucleotide on the protective immune response to alcelaphine herpesvirus-1-mediated malignant catarrhal fever in cattle.

We wished to determine the effect of of CpG ODN adjuvant on the magnitude and duration of protective immunity against alcelaphine herpesvirus-1 (AlHV-1) malignant catarrhal fever (MCF), a fatal lymphoproliferative disease of cattle. Immunity was associated with a mucosal barrier of virus-neutralising antibody. The results showed that CpG ODN included either with emulsigen adjuvant and attenuated AlHV-1 (atAlHV-1) or alone with atAlHV-1 did not affect the overall protection from clinical disease or duration of immunity achieved using emulsigen and atAlHV-1. This is in contrast to other similar studies in cattle with BoHV-1 or cattle and pigs with various other immunogens. In addition to this, several other novel observations were made, not reported previously. Firstly, we were able to statistically verify that vaccine protection against MCF was associated with virus-neutralising antibodies (nAbs) in nasal secretions but was not associated with antibodies in blood plasma, nor with total virus-specific antibody (tAb) titres in either nasal secretions or blood plasma. Furthermore, CpG ODN alone as adjuvant did not support the generation of virus-neutralising antibodies. Secondly, there was a significant boost in tAb in animals with MCF comparing titres before and after challenge. This was not seen with protected animals. Finally, there was a strong IFN-γ response in animals with emulsigen and atAlHV-1 immunisation, as measured by IFN-γ secreting PBMC in culture (and a lack of IL-4) that was not affected by the inclusion of CpG ODN. This suggests that nAbs at the oro-nasal-pharyngeal region are important in protection against AlHV-1 MCF.

Preliminary characterisation of Pentlands paramyxovirus-1, -2 and -3, three new paramyxoviruses of rodents.

A paramyxovirus was discovered by chance during the primary culture of grey squirrel (Sciurus carolinensis) kidney cells from the UK. Amplification, sequencing and phylogenetic analysis of part of the genome encoding a region of the RNA polymerase (L gene) confirmed that the virus was a member of the Paramyxovirinae subfamily, but that it did not partition with any of the currently recognised paramyxovirus genera and instead segregated with the unclassified rodent viruses, J-virus, Beilong virus and Tailam virus as well as paramyxoviruses recently detected in rodents in Africa. A subsequent examination of kidney samples from red squirrels (Sciurus vulgaris) revealed that they too harboured a paramyxovirus, but sequence analysis of the corresponding region of the L gene revealed that it was approximately 67% identical to the grey squirrel virus, suggesting the presence of a second species of virus. In addition, one of the red squirrels examined harboured a second virus with approximately 69% identity to the grey squirrel virus, but only approximately 63% identity to the other red squirrel viruses, signifying the presence of a third species of paramyxovirus. In a sample of 22 red and grey squirrels 68% of those examined were found to harbour virus suggesting that paramyxovirus infection in squirrels may be common within the UK.

Effects of parasitism and morphology on squirrelpox virus seroprevalence in grey squirrels (Sciurus carolinensis).

Invasive species have been cited as major causes of population extinctions in several animal and plant classes worldwide. The North American grey squirrel (Sciurus carolinensis) has a major detrimental effect on native red squirrel (Sciurus vulgaris) populations across Britain and Ireland, in part because it can be a reservoir host for the deadly squirrelpox virus (SQPV). Whilst various researchers have investigated the epizootiology of SQPV disease in grey squirrels and have modelled the consequent effects on red squirrel populations, less work has examined morphological and physiological characteristics that might make individual grey squirrels more susceptible to contracting SQPV. The current study investigated the putative relationships between morphology, parasitism, and SQPV exposure in grey squirrels. We found geographical, sex, and morphological differences in SQPV seroprevalence. In particular, larger animals, those with wide zygomatic arch widths (ZAW), males with large testes, and individuals with concurrent nematode and/or coccidial infections had an increased seroprevalence of SQPV. In addition, males with larger spleens, particularly those with narrow ZAW, were more likely to be exposed to SQPV. Overall these results show that there is variation in SQPV seroprevalence in grey squirrels and that, consequently, certain individual, or populations of, grey squirrels might be more responsible for transmitting SQPV to native red squirrel populations.