James C. Bull

A Virulent Strain of Deformed Wing Virus (DWV) of Honeybees (Apis mellifera) Prevails after Varroa destructor-Mediated, or In Vitro, Transmission

The globally distributed ectoparasite Varroa destructor is a vector for viral pathogens of the Western honeybee (Apis mellifera), in particular the Iflavirus Deformed Wing Virus (DWV). In the absence of Varroa low levels DWV occur, generally causing asymptomatic infections. Conversely, Varroa-infested colonies show markedly elevated virus levels, increased overwintering colony losses, with impairment of pupal development and symptomatic workers. To determine whether changes in the virus population were due Varroa amplifying and introducing virulent virus strains and/or suppressing the host immune responses, we exposed Varroa-naïve larvae to oral and Varroa-transmitted DWV. We monitored virus levels and diversity in developing pupae and associated Varroa, the resulting RNAi response and transcriptome changes in the host. Exposed pupae were stratified by Varroa association (presence/absence) and virus levels (low/high) into three groups. Varroa-free pupae all exhibited low levels of a highly diverse DWV population, with those exposed per os (group NV) exhibiting changes in the population composition. Varroa-associated pupae exhibited either low levels of a diverse DWV population (group VL) or high levels of a near-clonal virulent variant of DWV (group VH). These groups and unexposed controls (C) could be also discriminated by principal component analysis of the transcriptome changes observed, which included several genes involved in development and the immune response. All Varroa tested contained a diverse replicating DWV population implying the virulent variant present in group VH, and predominating in RNA-seq analysis of temporally and geographically separate Varroa-infested colonies, was selected upon transmission from Varroa, a conclusion supported by direct injection of pupae in vitro with mixed virus populations. Identification of a virulent variant of DWV, the role of Varroa in its transmission and the resulting host transcriptome changes furthers our understanding of this important viral pathogen of honeybees.

Persistence of an occlusion-negative recombinant nucleopolyhedrovirus in Trichoplusia ni indicates high multiplicity of cellular infection

ABSTRACT We use data from the serial passage of co-occluded recombinantAutographa californica nuclear polyhedrosis virus (AcMNPV) to estimate the viral multiplicity of infection of cells within infected insects. Co-occlusion, the incorporation of wild-type and mutant virus genomes in the same occlusion body, has been proposed as a strategy to deliver genetically modified viruses as insecticides in a way that contains their spread in the environment. It may also serve as a means whereby naturally occurring mutant forms of NPVs can be maintained in a stable polymorphism. Here, a recombinant strain of AcMNPV was constructed with a deletion of itspolyhedrin gene, rendering it incapable of producing occlusion bodies (i.e., occlusion negative). This was co-occluded with wild-type AcMNPV and used to infect fifth-instarTrichoplusia ni larvae. The fate of both genotypes was monitored over several rounds of insect infection. Levels of the occlusion-negative virus genome declined slowly over successive rounds of infection. We applied these data to a model of NPV population genetics to derive an estimate of 4.3 ± 0.3 viral genomes per occlusion body-producing cell.

A Strong Immune Response in Young Adult Honeybees Masks Their Increased Susceptibility to Infection Compared to Older Bees

Honeybees, Apis mellifera, show age-related division of labor in which young adults perform maintenance (“housekeeping”) tasks inside the colony before switching to outside foraging at approximately 23 days old. Disease resistance is an important feature of honeybee biology, but little is known about the interaction of pathogens and age-related division of labor. We tested a hypothesis that older forager bees and younger “house” bees differ in susceptibility to infection. We coupled an infection bioassay with a functional analysis of gene expression in individual bees using a whole genome microarray. Forager bees treated with the entomopathogenic fungus Metarhizium anisopliae s.l. survived for significantly longer than house bees. This was concomitant with substantial differences in gene expression including genes associated with immune function. In house bees, infection was associated with differential expression of 35 candidate immune genes contrasted with differential expression of only two candidate immune genes in forager bees. For control bees (i.e. not treated with M. anisopliae) the development from the house to the forager stage was associated with differential expression of 49 candidate immune genes, including up-regulation of the antimicrobial peptide gene abaecin, plus major components of the Toll pathway, serine proteases, and serpins. We infer that reduced pathogen susceptibility in forager bees was associated with age-related activation of specific immune system pathways. Our findings contrast with the view that the immunocompetence in social insects declines with the onset of foraging as a result of a trade-off in the allocation of resources for foraging. The up-regulation of immune-related genes in young adult bees in response to M. anisopliae infection was an indicator of disease susceptibility; this also challenges previous research in social insects, in which an elevated immune status has been used as a marker of increased disease resistance and fitness without considering the effects of age-related development.

A few-polyhedra mutant and wild-type nucleopolyhedrovirus remain as a stable polymorphism during serial coinfection in Trichoplusia ni.

ABSTRACT Few-polyhedra (FP) mutants of nucleopolyhedroviruses (NPVs) are a well-known phenomenon during serial passage of virus in cell culture. Under these circumstances such mutants produce low yields of occlusion bodies (OBs) and poorly occlude virions, but they are selected for through advantageous rates of budded virus replication. Spontaneous insertion of transposable elements originating from host cell DNA into the viral fp25 gene has been shown to be a common cause of the phenotype. A model of NPV population genetics predicts that mutants with these characteristics might persist within stable polymorphisms in viral populations during serial passage of virus in vivo. However, this hypothesis was previously untested, and FP mutants have not been recovered from field isolates of NPVs. We isolated and characterized an FP mutant that arose during routine passage of Autographa californica multinucleocapsid NPV (AcMNPV) in cell culture and identified a transposable element within the fp25 gene. We tracked the fates of coinfecting wild-type and FP mutant AcMNPV strains through serial passage in fifth-instar Trichoplusia ni larvae. The levels of both strains remained stable during successive rounds of infection. We applied the data obtained to a model of NPV population genetics in order to derive the frequency distribution of the multiplicity of cell infection in infected insects and estimated that 4.3 baculovirus genomes per OB-producing cell would account for this equilibrium.

Indirect effects and spatial scaling affect the persistence of multispecies metapopulations

Quantifying the role of space and spatial scale on the population dynamics of ecological assemblages is a contemporary challenge in ecology. Here, we evaluate the role of metapopulation dynamics on the persistence and dynamics of a multispecies predator–prey assemblage where two prey species shared a common natural enemy (apparent competition). By partitioning the effects of increased resource availability from the effects of metapopulation structure on regional population persistence we show that space has a marked impact on the dynamics of apparent competition in multispecies predator–prey assemblages. Further, the role of habitat size and stochasticity are also shown to influence the dynamics and persistence of this multispecies interaction. The broader consequences of these processes are discussed.

Unity in defence: honeybee workers exhibit conserved molecular responses to diverse pathogens

BackgroundOrganisms typically face infection by diverse pathogens, and hosts are thought to have developed specific responses to each type of pathogen they encounter. The advent of transcriptomics now makes it possible to test this hypothesis and compare host gene expression responses to multiple pathogens at a genome-wide scale. Here, we performed a meta-analysis of multiple published and new transcriptomes using a newly developed bioinformatics approach that filters genes based on their expression profile across datasets. Thereby, we identified common and unique molecular responses of a model host species, the honey bee (Apis mellifera), to its major pathogens and parasites: the Microsporidia Nosema apis and Nosema ceranae, RNA viruses, and the ectoparasitic mite Varroa destructor, which transmits viruses.ResultsWe identified a common suite of genes and conserved molecular pathways that respond to all investigated pathogens, a result that suggests a commonality in response mechanisms to diverse pathogens. We found that genes differentially expressed after infection exhibit a higher evolutionary rate than non-differentially expressed genes. Using our new bioinformatics approach, we unveiled additional pathogen-specific responses of honey bees; we found that apoptosis appeared to be an important response following microsporidian infection, while genes from the immune signalling pathways, Toll and Imd, were differentially expressed after Varroa/virus infection. Finally, we applied our bioinformatics approach and generated a gene co-expression network to identify highly connected (hub) genes that may represent important mediators and regulators of anti-pathogen responses.ConclusionsOur meta-analysis generated a comprehensive overview of the host metabolic and other biological processes that mediate interactions between insects and their pathogens. We identified key host genes and pathways that respond to phylogenetically diverse pathogens, representing an important source for future functional studies as well as offering new routes to identify or generate pathogen resilient honey bee stocks. The statistical and bioinformatics approaches that were developed for this study are broadly applicable to synthesize information across transcriptomic datasets. These approaches will likely have utility in addressing a variety of biological questions.

Wasting disease regulates long-term population dynamics in a threatened seagrass

The role of disease in the long-term dynamics of threatened species is poorly quantified, as well as being under-represented in ecology and conservation management. To understand persistent host–pathogen interaction operating in a vulnerable habitat, we quantified dynamics driving patterns of seagrass density using a longitudinal study in a relatively pristine site (Isles of Scilly, UK). Replicated samples of eelgrass (Zostera marina) density and wasting disease prevalence, presumably caused by Labyrinthula zosterae, were taken from five meadows at the height of the growing season, over the years 1997–2010. Data were used to parameterise a population dynamic model, incorporating density-dependent factors and sea temperature records. We found that direct density and disease-mediated feedback operate within a network of local populations. Furthermore, our results indicate that the strength of limitation to seagrass growth by disease was increased at higher temperatures. This modification of the coupled host–pathogen dynamics forms a novel hypothesis to account for dramatic die-backs of Z. marina widely reported elsewhere. Our findings highlight the importance of disease in structuring distributions of vulnerable species, as well as the application of population modelling in order to reveal ecological processes and prioritise future mechanistic investigation.

Metapopulation extinction risk is increased by environmental stochasticity and assemblage complexity

Extinction risk is a key area of investigation for contemporary ecologists and conservation biologists. Practical conservation efforts for vulnerable species can be considerably enhanced by thoroughly understanding the ecological processes that interact to determine species persistence or extinction. Theory has highlighted the importance of both extrinsic environmental factors and intrinsic demographic processes. In laboratory microcosms, single-species single-habitat patch experimental designs have been widely used to validate the theoretical prediction that environmental heterogeneity can increase extinction risk. Here, we develop on this theme by testing the effects of fluctuating resource levels in experimental multispecies metapopulations. We compare a three-species host–parasitoid assemblage that exhibits apparent competition to the individual pairwise, host–parasitoid interactions. Existing theory is broadly supported for two-species assemblages: environmental stochasticity reduces trophic interaction persistence time, while metapopulation structure increases persistence time. However, with increasing assemblage complexity, the effects of trophic interactions mask environmental impacts and persistence time is further reduced, regardless of resource renewal regime. We relate our findings to recent theory, highlighting the importance of taking into account both intrinsic and extrinsic factors, over a range of spatial scales, in order to understand resource–consumer dynamics.

The relationship between polychlorinated biphenyls in blubber and levels of nematode infestations in harbour porpoises, Phocoena phocoena

Post-mortem examinations of harbour porpoises, Phocoena phocoena, regularly reveal heavy parasitic worm burdens. These same post-mortem records show varying levels of polychlorinated biphenyls (PCBs) accumulating in the blubber of porpoises. Although a number of papers have documented geospatial and temporal changes of PCBs and their detrimental effects on marine mammal health, as yet none have examined their role in determining nematode burdens in wild marine mammal populations. Using a data set consisting of harbour porpoises stranded in the UK between 1989 and 2002, we found a significant, positive association between PCB levels and nematode burdens, although the nature of the relationship was confounded with porpoise sex, age and cause of death. It was also apparent that individuals with the heaviest infestations of nematodes did not have the highest PCB level: while PCBs are important, they are clearly not the sole determinants of nematode burdens in wild populations of the harbour porpoise around the UK.

Modelling the spread of American foulbrood in honeybees

We investigate the spread of American foulbrood (AFB), a disease caused by the bacterium Paenibacillus larvae, that affects bees and can be extremely damaging to beehives. Our dataset comes from an inspection period carried out during an AFB epidemic of honeybee colonies on the island of Jersey during the summer of 2010. The data include the number of hives of honeybees, location and owner of honeybee apiaries across the island. We use a spatial SIR model with an underlying owner network to simulate the epidemic and characterize the epidemic using a Markov chain Monte Carlo (MCMC) scheme to determine model parameters and infection times (including undetected ‘occult’ infections). Likely methods of infection spread can be inferred from the analysis, with both distance- and owner-based transmissions being found to contribute to the spread of AFB. The results of the MCMC are corroborated by simulating the epidemic using a stochastic SIR model, resulting in aggregate levels of infection that are comparable to the data. We use this stochastic SIR model to simulate the impact of different control strategies on controlling the epidemic. It is found that earlier inspections result in smaller epidemics and a higher likelihood of AFB extinction.