Coevolution of host and virus: the pathogenesis of virulent and attenuated strains of myxoma virus in resistant and susceptible European rabbits

The potential for evolutionary rescue in an Arctic seashore plant threatened by climate change

The impacts of climate change may be particularly severe for geographically isolated populations, which must adjust through plastic responses or evolve. Here, we study an endangered Arctic plant, Primula nutans ssp. finmarchica, confined to Fennoscandian seashores and showing indications of maladaptation to warming climate. We evaluate the potential of these populations to evolve to facilitate survival in the rapidly warming Arctic (i.e. evolutionary rescue) by utilizing manual crossing experiments in a nested half-sibling breeding design. We estimate G-matrices, evolvability and genetic constraints in traits with potentially conflicting selection pressures. To explicitly evaluate the potential for climate change adaptation, we infer the expected time to evolve from a northern to a southern phenotype under different selection scenarios, using demographic and climatic data to relate expected evolutionary rates to projected rates of climate change. Our results indicate that, given the nearly 10-fold greater evolvability of vegetative than of floral traits, adaptation in these traits may take place nearly in concert with changing climate, given effective climate mitigation. However, the comparatively slow expected evolutionary modification of floral traits may hamper the evolution of floral traits to track climate-induced changes in pollination environment, compromising sexual reproduction and thus reducing the likelihood of evolutionary rescue.

Quantifying resistance to myxomatosis in wild rabbits produces novel evolutionary insights

Wild rabbits in Australia developed genetic resistance to the myxoma virus, which was introduced as a biological control agent. However, little is known about the rate at which this evolutionary change occurred. We collated data from challenge trials that estimated rabbit resistance to myxomatosis in Australia and expressed resistance on a continuous scale, enabling trends in its development to be assessed over 45 years up to 1995. Resistance initially increased rapidly, followed by a plateau lasting ten years, before a second rapid increase occurred associated with the introduction of European rabbit fleas as myxoma virus vectors. By contrast, in the United Kingdom, where rabbit flea vectors were already present when the myxoma virus initially spread, resistance developed more slowly. No estimates of rabbit resistance to myxomatosis have been made for almost 30 years, despite other highly lethal rabbit pathogens becoming established worldwide. Continued testing of wild-caught rabbits in Australia to determine current levels of resistance to myxomatosis is recommended to assess its current effectiveness for managing pest rabbits. Given the economic and environmental significance of invasive rabbits, it would be remiss to manage such biological resources and ecosystem services poorly.

Host-pathogen coevolution promotes the evolution of general, broad-spectrum resistance and reduces foreign pathogen spillover risk

Genetic variation for disease resistance within host populations can strongly impact the spread of endemic pathogens. In plants, recent work has shown that within-population variation in resistance can also affect the transmission of foreign spillover pathogens if that resistance is general. However, most hosts also possess specific resistance mechanisms that provide strong defenses against coevolved endemic pathogens. Here we use a modeling approach to ask how antagonistic coevolution between hosts and their endemic pathogen at the specific resistance locus can affect the frequency of general resistance, and therefore a host's vulnerability to foreign pathogens. We develop a two-locus model with variable recombination that incorporates both general (resistance to all pathogens) and specific (resistance to endemic pathogens only). We find that introducing coevolution into our model greatly expands the regions where general resistance can evolve, decreasing the risk of foreign pathogen invasion. Furthermore, coevolution greatly expands which conditions maintain polymorphisms at both resistance loci, thereby driving greater genetic diversity within host populations. This genetic diversity often leads to positive correlations between host resistance to foreign and endemic pathogens, similar to those observed in natural populations. However, if resistance loci become linked, the resistance correlations can shift to negative. If we include a third, linkage modifying locus into our model, we find that selection often favors complete linkage. Our model demonstrates how coevolutionary dynamics with an endemic pathogen can mold the resistance structure of host populations in ways that affect its susceptibility to foreign pathogen spillovers, and that the nature of these outcomes depends on resistance costs, as well as the degree of linkage between resistance genes.

Polygenic plague resistance in the great gerbil uncovered by population sequencing

Pathogens can elicit high selective pressure on hosts, potentially altering genetic diversity over short evolutionary timescales. Intraspecific variation in immune response is observable as variable survivability from specific infections. The great gerbil (Rhombomys opimus) is a rodent plague host with a heterogenic but highly resistant phenotype. Here, we investigate the genomic basis for plague-resistant phenotypes by exposing wild-caught great gerbils to plague (Yersinia pestis). Whole genome sequencing of 10 survivors and 10 moribund individuals revealed a subset of genomic regions showing elevated differentiation. Gene ontology analysis of candidate genes in these regions demonstrated enrichment of genes directly involved in immune functions, cellular metabolism and the regulation of apoptosis as well as pathways involved in transcription, translation, and gene regulation. Transcriptomic analysis revealed that the early activated great gerbil immune response to plague consisted of classical components of the innate immune system. Our approach combining challenge experiments with transcriptomics and population level sequencing, provides new insight into the genetic background of plague-resistance and confirms its complex nature, most likely involving multiple genes and pathways of both the immune system and regulation of basic cellular functions.

Spatiotemporal monitoring of myxomatosis in European wild rabbit (Oryctolagus cuniculus) in Spanish Mediterranean ecosystems

A long-term active epidemiological surveillance programme was conducted to determine seroprevalence to myxoma virus (MYXV), infection prevalence and spatiotemporal patterns and factors associated with MYXV circulation in wild rabbits (Oryctolagus cuniculus) in Spanish Mediterranean ecosystems. A total of 2376 animals were sampled over four study periods: 2009-2012 (P1), 2012-2015 (P2), 2015-2018 (P3) and 2018-2021 (P4). Antibodies against MYXV were detected by a commercial indirect ELISA in 59.9% (1424/2376; 95% CI: 58.0-61.9) of wild rabbits. At least one seropositive animal was detected on 131 (96.3%) of 136 game estates sampled. MYXV infection was confirmed by PCR in 94 of 1063 (8.8%; 95% CI: 7.3-10.7) wild rabbits. Circulation of the novel recombinant MYXV (ha-MYXV) was not found in wild rabbits analysed during P4. Five statistically significant spatiotemporal clusters of high MYXV seroprevalence were identified using a Bernoulli model: one in P2 and four in P3. A generalized linear mixed model (GLMM) analysis identified sampling season (autumn), age (adult and juvenile), outbreaks of myxomatosis in the month prior to sampling, mean annual temperature, humidity and seropositivity to rabbit haemorrhagic disease virus as factors potentially linked with MYXV seropositivity. GLMM analysis identified outbreaks of myxomatosis in the month prior to sampling, MYXV seropositivity and presence of lesions compatible with myxomatosis as factors associated with MYXV infection. The results indicate high exposure, widespread but non-homogeneous distribution, and endemic circulation of MYXV in wild rabbit populations in southern Spain during the last decade. Prevalence of antibodies against MYXV showed fluctuations both within the year and over the study periods, revealing variations in the immunity of wild rabbit populations in Mediterranean ecosystems that could increase the risk of MYXV re-emergence in immunologically naïve populations. The present study highlights the importance of long-term surveillance to better understand the epidemiology of MYXV in wild lagomorphs.

Divergent Evolutionary Pathways of Myxoma Virus in Australia: Virulence Phenotypes in Susceptible and Partially Resistant Rabbits Indicate Possible Selection for Transmissibility

To characterize the ongoing evolution of myxoma virus in Australian rabbits, we used experimental infections of laboratory rabbits to determine the virulence and disease phenotypes of recent virus isolates. The viruses, collected between 2012 and 2015, fell into three lineages, one of which, lineage c, experienced a punctuated increase in evolutionary rate. All viruses were capable of causing acute death with aspects of neutropenic septicemia, characterized by minimal signs of myxomatosis, the occurrence of pulmonary edema and bacteria invasions throughout internal organs, but with no inflammatory response. For the viruses of highest virulence all rabbits usually died at this point. In more attenuated viruses, some rabbits died acutely, while others developed an amyxomatous phenotype. Rabbits that survived for longer periods developed greatly swollen cutaneous tissues with very high virus titers. This was particularly true of lineage c viruses. Unexpectedly, we identified a line of laboratory rabbits with some innate resistance to myxomatosis and used these in direct comparisons with the fully susceptible rabbit line. Importantly, the same disease phenotype occurred in both susceptible and resistant rabbits, although virulence was shifted toward more attenuated grades in resistant animals. We propose that selection against inflammation at cutaneous sites prolongs virus replication and enhances transmission, leading to the amyxomatous phenotype. In some virus backgrounds this creates an immunosuppressive state that predisposes to high virulence and acute death. The alterations in disease pathogenesis, particularly the overwhelming bacterial invasions that characterize the modern viruses, suggest that their virulence grades are not directly comparable with earlier studies.

IMPORTANCE The evolution of the myxoma virus (MYXV) following its release as a biological control for European rabbits in Australia is the textbook example of the coevolution of virus virulence and host resistance. However, most of our knowledge of MYXV evolution only covers the first few decades of its spread in Australia and often with little direct connection between how changes in virus phenotype relate to those in the underlying virus genotype. By conducting detailed experimental infections of recent isolates of MYXV in different lines of laboratory rabbits, we examined the ongoing evolution of MYXV disease phenotypes. Our results reveal a wide range of phenotypes, including an amyxomatous type, as well as the impact of invasive bacteria, that in part depended on the level of rabbit host resistance. These results provide a unique insight into the complex virus and host factors that combine to shape disease phenotype and viral evolution.

Social interactions of juvenile rabbits (Oryctolagus cuniculus) and their potential role in lagovirus transmission

Rabbit Haemorrhagic Disease Virus (RHDV), which is a calicivirus, is used as a biocontrol agent to suppress European wild rabbit populations in Australia. The transmission of RHDV can be influenced by social interactions of rabbits; however, there is a paucity of this knowledge about juvenile rabbits and the roles they may play in the transmission of RHDV. We aimed to quantify the social interactions of juvenile (< 900 g) and adult (> 1200 g) rabbits in a locally abundant population in the Central Tablelands of New South Wales, Australia. Twenty-six juvenile and 16 adult rabbits were fitted with VHF proximity loggers to monitor intra- and inter-group pairings. Use of multiple warrens by these rabbits was investigated using VHF base stations at nine warrens and on foot with a hand-held Yagi antenna. Juvenile rabbits were strongly interconnected with both juveniles and adults within and outside their warren of capture, and almost all juveniles were well-connected to other individuals within their own social group. Inter-group pairings were infrequent and fleeting between adults. Both juvenile and adult rabbits used multiple warrens. However, visits to warrens outside their warren of capture, particularly those within 50 m, were more common and longer in duration in juveniles than in adults. The high connectivity of juveniles within and between warrens in close proximity increases potential pathogen exchange between warrens. Therefore, juvenile rabbits could be of greater importance in lagovirus transmission than adult rabbits. The strength of juvenile rabbit inter- and intra-group pairings, and their tendency to use multiple warrens, highlight their potential to act as ‘superspreaders’ of both infection and immunity for lagoviruses and other pathogens with similar lifecycles. Confirmation of this potential is required through examination of disease progress and rabbit age-related immune responses during outbreaks.

Assessing the risk of human-to-wildlife pathogen transmission for conservation and public health

The SARS-CoV-2 pandemic has led to increased concern over transmission of pathogens from humans to animals, and its potential to threaten conservation and public health. To assess this threat, we reviewed published evidence of human-to-wildlife transmission events, with a focus on how such events could threaten animal and human health. We identified 97 verified examples, involving a wide range of pathogens; however, reported hosts were mostly non-human primates or large, long-lived captive animals. Relatively few documented examples resulted in morbidity and mortality, and very few led to maintenance of a human pathogen in a new reservoir or subsequent "secondary spillover" back into humans. We discuss limitations in the literature surrounding these phenomena, including strong evidence of sampling bias towards non-human primates and human-proximate mammals and the possibility of systematic bias against reporting human parasites in wildlife, both of which limit our ability to assess the risk of human-to-wildlife pathogen transmission. We outline how researchers can collect experimental and observational evidence that will expand our capacity for risk assessment for human-to-wildlife pathogen transmission.

The PB1 gene from H9N2 avian influenza virus showed high compatibility and increased mutation rate after reassorting with a human H1N1 influenza virus

BACKGROUND

Reassortment between human and avian influenza viruses (AIV) may result in novel viruses with new characteristics that may threaten human health when causing the next flu pandemic. A particular risk may be posed by avian influenza viruses of subtype H9N2 that are currently massively circulating in domestic poultry in Asia and have been shown to infect humans. In this study, we investigate the characteristics and compatibility of a human H1N1 virus with avian H9N2 derived genes.

METHODS

The polymerase activity of the viral ribonucleoprotein (RNP) complex as combinations of polymerase-related gene segments derived from different reassortment events was tested in luciferase reporter assays. Reassortant viruses were generated by reverse genetics. Gene segments of the human WSN-H1N1 virus (A/WSN/1933) were replaced by gene segments of the avian A2093-H9N2 virus (A/chicken/Jiangsu/A2093/2011), which were both the Hemagglutinin (HA) and Neuraminidase (NA) gene segments in combination with one of the genes involved in the RNP complex (either PB2, PB1, PA or NP). The growth kinetics and virulence of reassortant viruses were tested on cell lines and mice. The reassortant viruses were then passaged for five generations in MDCK cells and mice lungs. The HA gene of progeny viruses from different passaging paths was analyzed using Next-Generation Sequencing (NGS).

RESULTS

We discovered that the avian PB1 gene of H9N2 increased the polymerase activity of the RNP complex in backbone of H1N1. Reassortant viruses were able to replicate in MDCK and DF1 cells and mice. Analysis of the NGS data showed a higher substitution rate for the PB1-reassortant virus. In particular, for the PB1-reassortant virus, increased virulence for mice was measured by increased body weight loss after infection in mice.

CONCLUSIONS

The higher polymerase activity and increased mutation frequency measured for the PB1-reassortant virus suggests that the avian PB1 gene of H9N2 may drive the evolution and adaptation of reassortant viruses to the human host. This study provides novel insights in the characteristics of viruses that may arise by reassortment of human and avian influenza viruses. Surveillance for infections with H9N2 viruses and the emergence of the reassortant viruses in humans is important for pandemic preparedness.

Controlling long-term SARS-CoV-2 infections can slow viral evolution and reduce the risk of treatment failure

The rapid emergence and expansion of novel SARS-CoV-2 variants threatens our ability to achieve herd immunity for COVID-19. These novel SARS-CoV-2 variants often harbor multiple point mutations, conferring one or more evolutionarily advantageous traits, such as increased transmissibility, immune evasion and longer infection duration. In a number of cases, variant emergence has been linked to long-term infections in individuals who were either immunocompromised or treated with convalescent plasma. In this paper, we used a stochastic evolutionary modeling framework to explore the emergence of fitter variants of SARS-CoV-2 during long-term infections. We found that increased viral load and infection duration favor emergence of such variants. While the overall probability of emergence and subsequent transmission from any given infection is low, on a population level these events occur fairly frequently. Targeting these low-probability stochastic events that lead to the establishment of novel advantageous viral variants might allow us to slow the rate at which they emerge in the patient population, and prevent them from spreading deterministically due to natural selection. Our work thus suggests practical ways to achieve control of long-term SARS-CoV-2 infections, which will be critical for slowing the rate of viral evolution.

An ancient viral epidemic involving host coronavirus interacting genes more than 20,000 years ago in East Asia

The current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has emphasized the vulnerability of human populations to novel viral pressures, despite the vast array of epidemiological and biomedical tools now available. Notably, modern human genomes contain evolutionary information tracing back tens of thousands of years, which may help identify the viruses that have impacted our ancestors-pointing to which viruses have future pandemic potential. Here, we apply evolutionary analyses to human genomic datasets to recover selection events involving tens of human genes that interact with coronaviruses, including SARS-CoV-2, that likely started more than 20,000 years ago. These adaptive events were limited to the population ancestral to East Asian populations. Multiple lines of functional evidence support an ancient viral selective pressure, and East Asia is the geographical origin of several modern coronavirus epidemics. An arms race with an ancient coronavirus, or with a different virus that happened to use similar interactions as coronaviruses with human hosts, may thus have taken place in ancestral East Asian populations. By learning more about our ancient viral foes, our study highlights the promise of evolutionary information to better predict the pandemics of the future. Importantly, adaptation to ancient viral epidemics in specific human populations does not necessarily imply any difference in genetic susceptibility between different human populations, and the current evidence points toward an overwhelming impact of socioeconomic factors in the case of coronavirus disease 2019 (COVID-19).

Resistance Correlations Influence Infection by Foreign Pathogens

AbstractReciprocal selection promotes the specificity of host-pathogen associations and resistance polymorphisms in response to disease. However, plants and animals also vary in response to pathogen species not previously encountered in nature, with potential effects on new disease emergence. Using anther smut disease, we show that resistance (measured as infection rates) to foreign pathogens can be correlated with standing variation in resistance to an endemic pathogen. In Silene vulgaris, genetic variation in resistance to its endemic anther smut pathogen correlated positively with resistance variation to an anther smut pathogen from another host, but the relationship was negative between anther smut and a necrotrophic pathogen. We present models describing the genetic basis for assessing resistance relationships between endemic and foreign pathogens and for quantifying infection probabilities on foreign pathogen introduction. We show that even when the foreign pathogen has a lower average infection ability than the endemic pathogen, infection outcomes are determined by the sign and strength of the regression of the host's genetic variation in infection rates by a foreign pathogen on variation in infection rates by an endemic pathogen as well as by resistance allele frequencies. Given that preinvasion equilibria of resistance are determined by factors including resistance costs, we show that protection against foreign pathogens afforded by positively correlated resistances can be lessened or even result in elevated infection risk at the population level, depending on local dynamics. Therefore, a pathogen's emergence potential could be influenced not only by its average infection rate but also by resistance variation resulting from prior selection imposed by endemic diseases.

A globally ubiquitous symbiont can drive experimental host evolution

Organisms harbour myriad microbes which can be parasitic or protective against harm. The costs and benefits resulting from these symbiotic relationships can be context-dependent, but the evolutionary consequences to hosts of these transitions remain unclear. Here, we mapped the Leucobacter genus across 13,715 microbiome samples (163 studies) to reveal a global distribution as a free-living microbe or a symbiont of animals and plants. We showed that across geographically distant locations (South Africa, France, Cape Verde), Leucobacter isolates vary substantially in their virulence to an associated animal host, Caenorhabditis nematodes. We further found that multiple Leucobacter sequence variants co-occur in wild Caenorhabditis spp. which combined with natural variation in virulence provides real-world potential for Leucobacter community composition to influence host fitness. We examined this by competing C. elegans genotypes that differed in susceptibility to different Leucobacter species in an evolution experiment. One Leucobacter species was found to be host-protective against another, virulent parasitic species. We tested the impact of host genetic background and Leucobacter community composition on patterns of host-based defence evolution. We found host genotypes conferring defence against the parasitic species were maintained during infection. However, when hosts were protected during coinfection, host-based defences were nearly lost from the population. Overall, our results provide insight into the role of community context in shaping host evolution during symbioses.

Evolutionary rescue in host-pathogen systems

Natural populations encounter a variety of threats that can increase their risk of extinction. Populations can avoid extinction through evolutionary rescue (ER), which occurs when an adaptive, genetic response to selection allows a population to recover from an environmental change that would otherwise cause extinction. While the traditional framework for ER was developed with abiotic risk factors in mind, ER may also occur in response to a biotic source of demographic change, such as the introduction of a novel pathogen. We first describe how ER in response to a pathogen differs from the traditional ER framework; density-dependent transmission, pathogen evolution, and pathogen extinction can change the strength of selection imposed by a pathogen and make host population persistence more likely. We also discuss several variables that affect traditional ER (abundance, genetic diversity, population connectivity, and community composition) that also directly affect disease risk resulting in diverse outcomes for ER in host-pathogen systems. Thus, generalizations developed in studies of traditional ER may not be relevant for ER in response to the introduction of a pathogen. Incorporating pathogens into the framework of ER will lead to a better understanding of how and when populations can avoid extinction in response to novel pathogens.

Virus shedding kinetics and unconventional virulence tradeoffs

Tradeoff theory, which postulates that virulence provides both transmission costs and benefits for pathogens, has become widely adopted by the scientific community. Although theoretical literature exploring virulence-tradeoffs is vast, empirical studies validating various assumptions still remain sparse. In particular, truncation of transmission duration as a cost of virulence has been difficult to quantify with robust controlled in vivo studies. We sought to fill this knowledge gap by investigating how transmission rate and duration were associated with virulence for infectious hematopoietic necrosis virus (IHNV) in rainbow trout (Oncorhynchus mykiss). Using host mortality to quantify virulence and viral shedding to quantify transmission, we found that IHNV did not conform to classical tradeoff theory. More virulent genotypes of the virus were found to have longer transmission durations due to lower recovery rates of infected hosts, but the relationship was not saturating as assumed by tradeoff theory. Furthermore, the impact of host mortality on limiting transmission duration was minimal and greatly outweighed by recovery. Transmission rate differences between high and low virulence genotypes were also small and inconsistent. Ultimately, more virulent genotypes were found to have the overall fitness advantage, and there was no apparent constraint on the evolution of increased virulence for IHNV. However, using a mathematical model parameterized with experimental data, it was found that host culling resurrected the virulence tradeoff and provided low virulence genotypes with the advantage. Human-induced or natural culling, as well as host population fragmentation, may be some of the mechanisms by which virulence diversity is maintained in nature. This work highlights the importance of considering non-classical virulence tradeoffs.

Evaluating alternative hypotheses to explain the downward trend in the indices of the COVID-19 pandemic death rate

BACKGROUND

In the ongoing Covid-19 pandemic, in the global data on the case fatality ratio (CFR) and other indices reflecting death rate, there is a consistent downward trend from mid-April to mid-November. The downward trend can be an illusion caused by biases and limitations of data or it could faithfully reflect a declining death rate. A variety of explanations for this trend are possible, but a systematic analysis of the testable predictions of the alternative hypotheses has not yet been attempted.

METHODOLOGY

We state six testable alternative hypotheses, analyze their testable predictions using public domain data and evaluate their relative contributions to the downward trend.

RESULTS

We show that a decline in the death rate is real; changing age structure of the infected population and evolution of the virus towards reduced virulence are the most supported hypotheses and together contribute to major part of the trend. The testable predictions from other explanations including altered testing efficiency, time lag, improved treatment protocols and herd immunity are not consistently supported, or do not appear to make a major contribution to this trend although they may influence some other patterns of the epidemic.

CONCLUSION

The fatality of the infection showed a robust declining time trend between mid April to mid November. Changing age class of the infected and decreasing virulence of the pathogen were found to be the strongest contributors to the trend.

Viruses for Landscape-Scale Therapy: Biological Control of Rabbits in Australia

Viral diseases, whether of animals or humans, are normally considered as problems to be managed. However, in Australia, two viruses have been used as landscape-scale therapeutics to control European rabbits (Oryctolagus cuniculus), the preeminent invasive vertebrate pest species. Rabbits have caused major environmental and agricultural losses and contributed to extinction of native species. It was not until the introduction of Myxoma virus that effective control of this pest was obtained at a continental scale. Subsequent coevolution of rabbit and virus saw a gradual reduction in the effectiveness of biological control that was partially ameliorated by the introduction of the European rabbit flea to act as an additional vector for the virus. In 1995, a completely different virus, Rabbit hemorrhagic disease virus (RHDV), escaped from testing and spread through the Australian rabbit population and again significantly reduced rabbit numbers and environmental impacts. The evolutionary pressures on this virus appear to be producing quite different outcomes to those that occurred with myxoma virus and the emergence and invasion of a novel genotype of RHDV in 2014 have further augmented control. Molecular studies on myxoma virus have demonstrated multiple proteins that manipulate the host innate and adaptive immune response; however the molecular basis of virus attenuation and reversion to virulence are not yet understood.

A Mouse Model of Acute Liver Injury by Warm, Partial Ischemia-Reperfusion for Testing the Efficacy of Virus-Derived Therapeutics

Ischemia-reperfusion injury (IRI) drives early and long-term damage to organs as well as compounding damage from acute transplant rejection and surgical trauma. IRI initiates an aggressive and prolonged inflammation leading to tissue injury, organ failure, and death. However, there are few effective therapeutic interventions for IRI. The destructive inflammatory cell activity in IRI is part of an aberrant innate immune response that triggers multiple pathways. Hence, immune-modulating treatments to control pathways triggered by IRI hold great therapeutic potential. Viruses, especially large DNA viruses, have evolved highly effective immune-modulating proteins for the purpose of immune evasion and to protect the virus from the host immune defenses. A number of these immune-modulating proteins have proven therapeutically effective in preclinical models, many with function targeting pathways known to be involved in IRI. The use of virus-derived immune-modulating proteins thus represents a promising source for new treatments to target ischemia-reperfusion injury. Laboratory small animal models of IRI are well established and are able to reproduce many aspects of ischemia-reperfusion injury seen in humans. This chapter will discuss the methods used to perform the IRI procedure in mice, as well as clinically relevant diagnostic tests to evaluate liver injury and approaches for assessing histological damage while testing novel immune modulating protein treatments.

COVID-19, varying genetic resistance to viral disease and immune tolerance checkpoints

Coronavirus disease 2019 (COVID‐19) is a zoonosis like most of the great plagues sculpting human history, from smallpox to pandemic influenza and human immunodeficiency virus. When viruses jump into a new species the outcome of infection ranges from asymptomatic to lethal, historically ascribed to “genetic resistance to viral disease.” People have exploited these differences for good and bad, for developing vaccines from cowpox and horsepox virus, controlling rabbit plagues with myxoma virus and introducing smallpox during colonization of America and Australia. Differences in resistance to viral disease are at the core of the severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) crisis, yet our understanding of the mechanisms in any interspecies leap falls short of the mark. Here I review how the two key parameters of viral disease are countered by fundamentally different genetic mechanisms for resistance: (1) virus transmission, countered primarily by activation of innate and adaptive immune responses; and (2) pathology, countered primarily by tolerance checkpoints to limit innate and adaptive immune responses. I discuss tolerance thresholds and the role of CD8 T cells to limit pathological immune responses, the problems posed by tolerant superspreaders and the signature coronavirus evasion strategy of eliciting only short‐lived neutralizing antibody responses. Pinpointing and targeting the mechanisms responsible for varying pathology and short‐lived antibody were beyond reach in previous zoonoses, but this time we are armed with genomic technologies and more knowledge of immune checkpoint genes. These known unknowns must now be tackled to solve the current COVID‐19 crisis and the inevitable zoonoses to follow.

Recombinant Myxoma Virus-Derived Immune Modulator M-T7 Accelerates Cutaneous Wound Healing and Improves Tissue Remodeling

Complex dermal wounds represent major medical and financial burdens, especially in the context of comorbidities such as diabetes, infection and advanced age. New approaches to accelerate and improve, or "fine tune" the healing process, so as to improve the quality of cutaneous wound healing and management, are the focus of intense investigation. Here, we investigate the topical application of a recombinant immune modulating protein which inhibits the interactions of chemokines with glycosaminoglycans, reducing damaging or excess inflammation responses in a splinted full-thickness excisional wound model in mice. M-T7 is a 37 kDa-secreted, virus-derived glycoprotein that has demonstrated therapeutic efficacy in numerous animal models of inflammatory immunopathology. Topical treatment with recombinant M-T7 significantly accelerated wound healing when compared to saline treatment alone. Healed wounds exhibited properties of improved tissue remodeling, as determined by collagen maturation. M-T7 treatment accelerated the rate of peri-wound angiogenesis in the healing wounds with increased levels of TNF, VEGF and CD31. The immune cell response after M-T7 treatment was associated with a retention of CCL2 levels, and increased abundances of arginase-1-expressing M2 macrophages and CD4 T cells. Thus, topical treatment with recombinant M-T7 promotes a pro-resolution environment in healing wounds, and has potential as a novel treatment approach for cutaneous tissue repair.

Mechanisms underlying host persistence following amphibian disease emergence determine appropriate management strategies

Emerging infectious diseases have caused many species declines, changes in communities and even extinctions. There are also many species that persist following devastating declines due to disease. The broad mechanisms that enable host persistence following declines include evolution of resistance or tolerance, changes in immunity and behaviour, compensatory recruitment, pathogen attenuation, environmental refugia, density-dependent transmission and changes in community composition. Here we examine the case of chytridiomycosis, the most important wildlife disease of the past century. We review the full breadth of mechanisms allowing host persistence, and synthesise research on host, pathogen, environmental and community factors driving persistence following chytridiomycosis-related declines and overview the current evidence and the information required to support each mechanism. We found that for most species the mechanisms facilitating persistence have not been identified. We illustrate how the mechanisms that drive long-term host population dynamics determine the most effective conservation management strategies. Therefore, understanding mechanisms of host persistence is important because many species continue to be threatened by disease, some of which will require intervention. The conceptual framework we describe is broadly applicable to other novel disease systems.

Detection of recombinant Hare Myxoma Virus in wild rabbits (Oryctolagus cuniculus algirus)

In late 2018, an epidemic myxomatosis outbreak emerged on the Iberian Peninsula leading to high mortality in Iberian hare populations. A recombinant Myxoma virus (strains MYXV-Tol and ha-MYXV) was rapidly identified, harbouring a 2.8 kbp insertion containing evolved duplicates of M060L, M061L, M064L, and M065L genes from myxoma virus (MYXV) or other Poxviruses. Since 2017, 1616 rabbits and 125 hares were tested by a qPCR directed to M000.5L/R gene, conserved in MYXV and MYXV-Tol/ha-MYXV strains. A subset of the positive samples (20%) from both species was tested for the insert with MYXV being detected in rabbits and the recombinant MYXV in hares. Recently, three wild rabbits were found dead South of mainland Portugal, showing skin oedema and pulmonary lesions that tested positive for the 2.8 kbp insert. Sequencing analysis showed 100% similarity with the insert sequences described in Iberian hares from Spain. Viral particles were observed in the lungs and eyelids of rabbits by electron microscopy, and isolation in RK13 cells attested virus infectivity. Despite that the analysis of complete genomes may predict the recombinant MYXV strains’ ability to infect rabbit, routine analyses showed species segregation for the circulation of MYXV and recombinant MYXV in wild rabbit and in Iberian hares, respectively. This study demonstrates, however, that recombinant MYXV can effectively infect and cause myxomatosis in wild rabbits and domestic rabbits, raising serious concerns for the future of the Iberian wild leporids while emphasises the need for the continuous monitoring of MYXV and recombinant MYXV in both species.

Seasonality and risk factors for myxomatosis in pet rabbits in Great Britain

Myxomatosis is a highly contagious, frequently fatal viral disease affecting both wild and domesticated European rabbits across many areas of the world. Here we used electronic health records (EHRs) collected from pet rabbits attending a sentinel voluntary network of 191 veterinary practices across Great Britain (GB) between March 2014 and June 2019 to identify new features of this disease's epidemiology. From a total of 89,408 rabbit consultations, text mining verified by domain experts identified 207 (0.23 %) cases where myxomatosis was the only differential diagnosis recorded by the attending practitioner. Cases occurred in all months but February and were distributed across the country. Consistent with studies in wild rabbits, the majority of cases occurred between August and November. However, there was also evidence for considerable variation between years. A nested case control study identified important risk factors for myxomatosis within this pet animal population including season, sex, age, vaccination status and distance to likely wild rabbit habitats. Female entire rabbits were twice as likely to be a case (odds ratio (OR) 1.98, 95 % confidence interval (CI) 1.26-3.13, p = 0.003), suggesting a novel role for behaviour in driving transmission from wild to domesticated rabbits. Vaccination had the largest protective effect with vaccinated rabbits being 8.3 times less likely to be a case than unvaccinated rabbits (OR = 0.12, 95 % CI 0.06-0.21, p = <0.001).

First outbreak of myxomatosis in Iberian hares (Lepus granatensis)

Myxomatosis is an infectious disease caused by myxoma virus (MYXV; genus Leporipoxvirus), which affects the European wild rabbit (Oryctolagus cuniculus) and sporadically brown hares (Lepus europaeus). Here, we describe the first outbreak of myxomatosis in Iberian hares (Lepus granatensis). Between mid-July and the end of September 2018, around 530 dead animals were detected in Iberian hare populations in southern Spain. The apparent mean mortality rate was 56.7%, and the estimated mean case fatality rate was 69.2%. Histopathological and molecular results confirmed MYXV infections in all hares analysed. To the authors' knowledge, this is the first myxomatosis outbreak causing a high mortality in hares and the first detailed characterization of a myxomatosis outbreak in the Iberian hare. The absence of cases in sympatric wild rabbits suggests differences in the susceptibility between both lagomorph species to the virus strain implicated in the outbreak. After the first case, the number of affected areas increased sharply affecting most of the Iberian Peninsula where the Iberian hare is present. Further studies are required to elucidate the origin of the implicated MYXV strain as well as to assess the impact of this outbreak on the Iberian hare populations.

Serp-2, a virus-derived apoptosis and inflammasome inhibitor, attenuates liver ischemia-reperfusion injury in mice

BACKGROUND

Ischemia-reperfusion injury (IRI) is an antigen-independent, innate immune response to arterial occlusion and ischemia with subsequent paradoxical exacerbation after reperfusion. IRI remains a critical problem after vessel occlusion and infarction or during harvest and surgery in transplants. After transplant, liver IRI (LIRI) contributes to increased acute and chronic rejection and graft loss. Tissue loss during LIRI has been attributed to local macrophage activation and invasion with excessive inflammation together with hepatocyte apoptosis and necrosis. Inflammatory and apoptotic signaling are key targets for reducing post-ischemic liver injury.Myxomavirus is a rabbit-specific leporipoxvirus that encodes a suite of immune suppressing proteins, often with extensive function in other mammalian species. Serp-2 is a cross-class serine protease inhibitor (serpin) which inhibits the inflammasome effector protease caspase-1 as well as the apoptotic proteases granzyme B and caspases 8 and 10. In prior work, Serp-2 reduced inflammatory cell invasion after angioplasty injury and after aortic transplantation in rodents. In this report, we explore the potential for therapeutic treatment with Serp-2 in a mouse model of LIRI.

METHODS

Wildtype (C57BL/6 J) mice were subjected to warm, partial (70%) hepatic ischemia for 90 min followed by treatment with saline or Serp-2 or M-T7, 100 ng/g/day given by intraperitoneal injection on alternate days for 5 days. M-T7 is a Myxomavirus-derived inhibitor of chemokine-GAG interactions and was used in this study for comparative analysis of an unrelated viral protein with an alternative immunomodulating mechanism of action. Survival, serum ALT levels and histopathology were assessed 24 h and 10 days post-LIRI.

RESULTS

Serp-2 treatment significantly improved survival to 85.7% percent versus saline-treated wildtype mice (p = 0.0135), while M-T7 treatment did not significantly improve survival (p = 0.2584). Liver viability was preserved by Serp-2 treatment with a significant reduction in serum ALT levels (p = 0.0343) and infarct scar thickness (p = 0.0016), but with no significant improvement with M-T7 treatment. Suzuki scoring by pathologists blinded with respect to treatment group indicated that Serp-2 significantly reduced hepatocyte necrosis (p = 0.0057) and improved overall pathology score (p = 0.0046) compared to saline. Immunohistochemistry revealed that Serp-2 treatment reduced macrophage infiltration into the infarcted liver tissue (p = 0.0197).

CONCLUSIONS

Treatment with Serp-2, a virus-derived inflammasome and apoptotic pathway inhibitor, improves survival after liver ischemia-reperfusion injury in mouse models. Treatment with a cross-class immune modulator provides a promising new approach designed to reduce ischemia-reperfusion injury, improving survival and reducing chronic transplant damage.

Eco-evolutionary rescue promotes host-pathogen coexistence

Emerging infectious pathogens are responsible for some of the most severe host mass mortality events in wild populations. Yet, effective pathogen control strategies are notoriously difficult to identify, in part because quantifying and forecasting pathogen spread and disease dynamics is challenging. Following an outbreak, hosts must cope with the presence of the pathogen, leading to host-pathogen coexistence or extirpation. Despite decades of research, little is known about host-pathogen coexistence post-outbreak when low host abundances and cryptic species make these interactions difficult to study. Using a novel disease-structured N-mixture model, we evaluate empirical support for three host-pathogen coexistence hypotheses (source-sink, eco-evolutionary rescue, and spatial variation in pathogen transmission) in a Neotropical amphibian community decimated by Batrachochytrium dendrobatidis (Bd) in 2004. During 2010-2014, we surveyed amphibians in Parque Nacional G. D. Omar Torríjos Herrera, Coclé Province, El Copé, Panama. We found that the primary driver of host-pathogen coexistence was eco-evolutionary rescue, as evidenced by similar amphibian survival and recruitment rates between infected and uninfected hosts. Average apparent monthly survival rates of uninfected and infected hosts were both close to 96%, and the expected number of uninfected and infected hosts recruited (via immigration/reproduction) was less than one host per disease state per 20-m site. The secondary driver of host-pathogen coexistence was spatial variation in pathogen transmission as we found that transmission was highest in areas of low abundance but there was no support for the source-sink hypothesis. Our results indicate that changes in the host community (i.e., through genetic or species composition) can reduce the impacts of emerging infectious disease post-outbreak. Our disease-structured N-mixture model represents a valuable advancement for conservation managers trying to understand underlying host-pathogen interactions and provides new opportunities to study disease dynamics in remnant host populations decimated by virulent pathogens.

Parallel evolution of gene classes, but not genes: Evidence from Hawai'ian honeycreeper populations exposed to avian malaria

Adaptation in nature is ubiquitous, yet characterizing its genomic basis is difficult because population demographics cause correlations with nonadaptive loci. Introduction events provide opportunities to observe adaptation over known spatial and temporal scales, facilitating the identification of genes involved in adaptation. The pathogen causing avian malaria, Plasmodium relictum, was introduced to Hawai'i in the 1930s and elicited extinctions and precipitous population declines in native honeycreepers. After a sharp initial population decline, the Hawai'i 'amakihi (Chlorodrepanis virens) has evolved tolerance to the parasite at low elevations where P. relictum exists, and can sustain infection without major fitness consequences. High-elevation, unexposed populations of 'amakihi display little to no tolerance. To explore the genomic basis of adaptation to P. relictum in low-elevation 'amakihi, we genotyped 125 'amakihi from the island of Hawai'i via hybridization capture to 40,000 oligonucleotide baits containing SNPs and used the reference 'amakihi genome to identify genes potentially under selection from malaria. We tested for outlier loci between low- and high-elevation population pairs and identified loci with signatures of selection within low-elevation populations. In some cases, genes commonly involved in the immune response (e.g., major histocompatibility complex) were associated with malaria presence in the population. We also detected several novel candidate loci that may be implicated in surviving malaria infection (e.g., beta-defensin, glycoproteins and interleukin-related genes). Our results suggest that rapid adaptation to pathogens may occur through changes in different immune genes, but in the same classes of genes, across populations.

Quantitative host resistance drives the evolution of increased virulence in an emerging pathogen

Emergent infectious diseases can have a devastating impact on host populations. The high selective pressures on both the hosts and the pathogens frequently lead to rapid adaptations not only in pathogen virulence but also host resistance following an initial outbreak. However, it is often unclear whether hosts will evolve to avoid infection-associated fitness costs by preventing the establishment of infection (here referred to as qualitative resistance) or by limiting its deleterious effects through immune functioning (here referred to as quantitative resistance). Equally, the evolutionary repercussions these different resistance mechanisms have for the pathogen are often unknown. Here, we investigate the co-evolutionary dynamics of pathogen virulence and host resistance following the epizootic outbreak of the highly pathogenic bacterium Mycoplasma gallisepticum in North American house finches (Haemorhous mexicanus). Using an evolutionary modelling approach and with a specific emphasis on the evolved resistance trait, we demonstrate that the rapid increase in the frequency of resistant birds following the outbreak is indicative of strong selection pressure to reduce infection-associated mortality. This, in turn, created the ecological conditions that selected for increased bacterial virulence. Our results thus suggest that quantitative host resistance was the key factor underlying the evolutionary interactions in this natural host-pathogen system.

Humoral immune response to different routes of myxomatosis vaccine application

The aim of our study was to monitor the dynamics of the serological response to different application routes of live attenuated myxomatosis vaccine. The study included 42 Californian breed rabbits, aged 3 mo, of both sexes. They were separated into 7 groups: 6 experimental and 1 control. All experimental groups were vaccinated on day 0 with a single dose of myxomatosis vaccine (min 10^3.3 tissue culture infective dose 50 [TCID₅₀], max 10^5.8 TCID₅₀). Three of the groups were injected with monovalent attenuated myxomatosis vaccine using different types of application: intradermal (i.d.), intramuscular (i.m.) and subcutaneous (s.c.). The other 3 groups were injected with bivalent attenuated vaccine against myxomatosis and rabbit haemorrhagic disease; again the routes of administration were i.d., i.m. and s.c.. There were no clinical signs or serious side effects after vaccination. The serological response was evaluated on days 7, 15 and 30 with a monoclonal antibody based-competition enzyme-linked immunosorbent assay (cELISA). More rapid and potent humoral response was detected in groups with i.d. inoculation in comparison to i.m. and s.c. routes. Vaccination with monovalent vaccine against myxomatosis induced higher antibody titre in comparison to bivalent vaccine. Our study showed that the vaccine application route and the type of vaccine used influence the speed and intensity of antibody response.

Reverse Engineering Field Isolates of Myxoma Virus Demonstrates that Some Gene Disruptions or Losses of Function Do Not Explain Virulence Changes Observed in the Field

The coevolution of myxoma virus (MYXV) and wild European rabbits in Australia and Europe is a paradigm for the evolution of a pathogen in a new host species. Genomic analyses have identified the mutations that have characterized this evolutionary process, but defining causal mutations in the pathways from virulence to attenuation and back to virulence has not been possible. Using reverse genetics, we examined the roles of six selected mutations found in Australian field isolates of MYXV that fall in known or potential virulence genes. Several of these mutations occurred in genes previously identified as virulence genes in whole-gene knockout studies. Strikingly, no single or double mutation among the mutations tested had an appreciable impact on virulence. This suggests either that virulence evolution was defined by amino acid changes other than those analyzed here or that combinations of multiple mutations, possibly involving epistatic interactions or noncoding sequences, have been critical in the ongoing evolution of MYXV virulence. In sum, our results show that single-gene knockout studies of a progenitor virus can have little power to predict the impact of individual mutations seen in the field. The genetic determinants responsible for this canonical case of virulence evolution remain to be determined.IMPORTANCE The species jump of myxoma virus (MYXV) from the South American tapeti to the European rabbit populations of Australia and Europe is a canonical example of host-pathogen coevolution. Detailed molecular studies have identified multiple genes in MYXV that are critical for virulence, and genome sequencing has revealed the evolutionary history of MYXV in Australia and Europe. However, it has not been possible to categorically identify the key mutations responsible for the attenuation of or reversion to virulence during this evolutionary process. Here we use reverse genetics to examine the role of mutations in viruses isolated early and late in the Australian radiation of MYXV. Surprisingly, none of the candidate mutations that we identified as likely having roles in attenuation proved to be important for virulence. This indicates that considerable caution is warranted when interpreting the possible role of individual mutations during virulence evolution.

Next step in the ongoing arms race between myxoma virus and wild rabbits in Australia is a novel disease phenotype

In host-pathogen arms races, increases in host resistance prompt counteradaptation by pathogens, but the nature of that counteradaptation is seldom directly observed outside of laboratory models. The best-documented field example is the coevolution of myxoma virus (MYXV) in European rabbits. To understand how MYXV in Australia has continued to evolve in wild rabbits under intense selection for genetic resistance to myxomatosis, we compared the phenotypes of the progenitor MYXV and viral isolates from the 1950s and the 1990s in laboratory rabbits with no resistance. Strikingly, and unlike their 1950s counterparts, most virus isolates from the 1990s induced a highly lethal immune collapse syndrome similar to septic shock. Thus, the next step in this canonical case of coevolution after a species jump has been further escalation by the virus in the face of widespread host resistance.

Genomic and phenotypic characterization of myxoma virus from Great Britain reveals multiple evolutionary pathways distinct from those in Australia

The co-evolution of myxoma virus (MYXV) and the European rabbit occurred independently in Australia and Europe from different progenitor viruses. Although this is the canonical study of the evolution of virulence, whether the genomic and phenotypic outcomes of MYXV evolution in Europe mirror those observed in Australia is unknown. We addressed this question using viruses isolated in the United Kingdom early in the MYXV epizootic (1954-1955) and between 2008-2013. The later UK viruses fell into three distinct lineages indicative of a long period of separation and independent evolution. Although rates of evolutionary change were almost identical to those previously described for MYXV in Australia and strongly clock-like, genome evolution in the UK and Australia showed little convergence. The phenotypes of eight UK viruses from three lineages were characterized in laboratory rabbits and compared to the progenitor (release) Lausanne strain. Inferred virulence ranged from highly virulent (grade 1) to highly attenuated (grade 5). Two broad disease types were seen: cutaneous nodular myxomatosis characterized by multiple raised secondary cutaneous lesions, or an amyxomatous phenotype with few or no secondary lesions. A novel clinical outcome was acute death with pulmonary oedema and haemorrhage, often associated with bacteria in many tissues but an absence of inflammatory cells. Notably, reading frame disruptions in genes defined as essential for virulence in the progenitor Lausanne strain were compatible with the acquisition of high virulence. Combined, these data support a model of ongoing host-pathogen co-evolution in which multiple genetic pathways can produce successful outcomes in the field that involve both different virulence grades and disease phenotypes, with alterations in tissue tropism and disease mechanisms.

Overlooking the smallest matter: viruses impact biological invasions

Parasites and pathogens have recently received considerable attention for their ability to affect biological invasions, however, researchers have largely overlooked the distinct role of viruses afforded by their unique ability to rapidly mutate and adapt to new hosts. With high mutation and genomic substitution rates, RNA and single-stranded DNA (ssDNA) viruses may be important constituents of invaded ecosystems, and could potentially behave quite differently from other pathogens. We review evidence suggesting that rapidly evolving viruses impact invasion dynamics in three key ways: (1) Rapidly evolving viruses may prevent exotic species from establishing self-sustaining populations. (2) Viruses can cause population collapses of exotic species in the introduced range. (3) Viruses can alter the consequences of biological invasions by causing population collapses and extinctions of native species. The ubiquity and frequent host shifting of viruses make their ability to influence invasion events likely. Eludicating the viral ecology of biological invasions will lead to an improved understanding of the causes and consequences of invasions, particularly as regards establishment success and changes to community structure that cannot be explained by direct interspecific interactions among native and exotic species.

Comparative study of host response to chytridiomycosis in a susceptible and a resistant toad species

In the past century, recently emerged infectious diseases have become major drivers of species decline and extinction. The fungal disease chytridiomycosis has devastated many amphibian populations and exacerbated the amphibian conservation crisis. Biologists are beginning to understand what host traits contribute to disease susceptibility, but more work is needed to determine why some species succumb to chytridiomycosis while others do not. We conducted an integrative laboratory experiment to examine how two toad species respond to infection with the pathogen Batrachochytrium dendrobatidis in a controlled environment. We selected two toad species thought to differ in susceptibility - Bufo marinus (an invasive and putatively resistant species) and Bufo boreas (an endangered and putatively susceptible species). We measured infection intensity, body weight, histological changes and genomewide gene expression using a custom assay developed from transcriptome sequencing. Our results confirmed that the two species differ in susceptibility with the more susceptible species, B. boreas, showing higher infection intensities, loss in body weight, more dramatic histological changes and larger perturbations in gene expression. We found key differences in skin expression responses in multiple pathways including upregulation of skin integrity-related genes in the resistant B. marinus. Together, our results show intrinsic differences in host response between related species, which are likely to be important in explaining variation in response to a deadly emerging pathogen in wild populations. Our study also underscores the importance of understanding differences among host species to better predict disease outcomes and reveal generalities about host response to emerging infectious diseases of wildlife.

Tissue-specific transcription profile of cytokine and chemokine genes associated with flavivirus control and non-lethal neuropathogenesis in rabbits

We previously showed that New Zealand White (NZWRs) and cottontail rabbits (CTRs) are a suitable model for studying immune mechanisms behind virus control and non-lethal neuropathogenesis associated with West Nile virus (WNV) and Murray Valley encephalitis virus (MVEV) infections. In the current study, we observed that MVEV infection induced high IFNα, TNFα, IL6, and CXCL10 transcript levels in the brains of weanling NZWRs, unlike infection with the less virulent WNVNSW2011. These transcript levels also correlated with encephalitis severity. Widespread STAT1 protein expression in brain with moderate neuropathology suggests that IFN-I signaling is crucial for limiting neural infection and mediating non-lethal neuropathogenesis. Unlike NZWRs, CTRs limit neuroinvasion without upregulation of many cytokine/chemokine transcripts, suggesting a species-dependent virus control mechanism. However, the common IFNγ, TNFα and IL6 transcript upregulation in specific lymphoid organs suggest some conserved elements in the response against flaviviruses, unique to all rabbits.

Targeted gene flow for conservation

Anthropogenic threats often impose strong selection on affected populations, causing rapid evolutionary responses. Unfortunately, these adaptive responses are rarely harnessed for conservation. We suggest that conservation managers pay close attention to adaptive processes and geographic variation, with an eye to using them for conservation goals. Translocating pre-adapted individuals into recipient populations is currently considered a potentially important management tool in the face of climate change. Targeted gene flow, which involves moving individuals with favorable traits to areas where these traits would have a conservation benefit, could have a much broader application in conservation. Across a species' range there may be long-standing geographic variation in traits or variation may have rapidly developed in response to a threatening process. Targeted gene flow could be used to promote natural resistance to threats to increase species resilience. We suggest that targeted gene flow is a currently underappreciated strategy in conservation that has applications ranging from the management of invasive species and their impacts to controlling the impact and virulence of pathogens.

An overview of the lagomorph immune system and its genetic diversity

Our knowledge of the lagomorph immune system remains largely based upon studies of the European rabbit (Oryctolagus cuniculus), a major model for studies of immunology. Two important and devastating viral diseases, rabbit hemorrhagic disease and myxomatosis, are affecting European rabbit populations. In this context, we discuss the genetic diversity of the European rabbit immune system and extend to available information about other lagomorphs. Regarding innate immunity, we review the most recent advances in identifying interleukins, chemokines and chemokine receptors, Toll-like receptors, antiviral proteins (RIG-I and Trim5), and the genes encoding fucosyltransferases that are utilized by rabbit hemorrhagic disease virus as a portal for invading host respiratory and gut epithelial cells. Evolutionary studies showed that several genes of innate immunity are evolving by strong natural selection. Studies of the leporid CCR5 gene revealed a very dramatic change unique in mammals at the second extracellular loop of CCR5 resulting from a gene conversion event with the paralogous CCR2. For the adaptive immune system, we review genetic diversity at the loci encoding antibody variable and constant regions, the major histocompatibility complex (RLA) and T cells. Studies of IGHV and IGKC genes expressed in leporids are two of the few examples of trans-species polymorphism observed outside of the major histocompatibility complex. In addition, we review some endogenous viruses of lagomorph genomes, the importance of the European rabbit as a model for human disease studies, and the anticipated role of next-generation sequencing in extending knowledge of lagomorph immune systems and their evolution.

Myxoma virus and the Leporipoxviruses: an evolutionary paradigm

Myxoma virus (MYXV) is the type species of the Leporipoxviruses, a genus of Chordopoxvirinae, double stranded DNA viruses, whose members infect leporids and squirrels, inducing cutaneous fibromas from which virus is mechanically transmitted by biting arthropods. However, in the European rabbit (Oryctolagus cuniculus), MYXV causes the lethal disease myxomatosis. The release of MYXV as a biological control for the wild European rabbit population in Australia, initiated one of the great experiments in evolution. The subsequent coevolution of MYXV and rabbits is a classic example of natural selection acting on virulence as a pathogen adapts to a novel host species. Slightly attenuated mutants of the progenitor virus were more readily transmitted by the mosquito vector because the infected rabbit survived longer, while highly attenuated viruses could be controlled by the rabbit immune response. As a consequence, moderately attenuated viruses came to dominate. This evolution of the virus was accompanied by selection for genetic resistance in the wild rabbit population, which may have created an ongoing co-evolutionary dynamic between resistance and virulence for efficient transmission. This natural experiment was repeated on a continental scale with the release of a separate strain of MYXV in France and its subsequent spread throughout Europe. The selection of attenuated strains of virus and resistant rabbits mirrored the experience in Australia in a very different environment, albeit with somewhat different rates. Genome sequencing of the progenitor virus and the early radiation, as well as those from the 1990s in Australia and Europe, has shown that although MYXV evolved at high rates there was no conserved route to attenuation or back to virulence. In contrast, it seems that these relatively large viral genomes have the flexibility for multiple pathways that converge on a similar phenotype.

The α2,3-sialyltransferase encoded by myxoma virus is a virulence factor that contributes to immunosuppression

Myxoma virus (MYXV) induces a lethal disease called Myxomatosis in European rabbits. MYXV is one of the rare viruses that encodes an α2,3-sialyltransferase through its M138L gene. In this study, we showed that although the absence of the enzyme was not associated with any in vitro deficit, the M138L deficient strains are highly attenuated in vivo. Indeed, while all rabbits infected with the parental and the revertant strains died within 9 days post-infection from severe myxomatosis, all but one rabbit inoculated with the M138L deficient strains survived the infection. In primary lesions, this resistance to the infection was associated with an increased ability of innate immune cells, mostly neutrophils, to migrate to the site of virus replication at 4 days post-infection. This was followed by the development of a better specific immune response against MYXV. Indeed, at day 9 post-infection, we observed an important proliferation of lymphocytes and an intense congestion of blood vessels in lymph nodes after M138L knockouts infection. Accordingly, in these rabbits, we observed an intense mononuclear cell infiltration throughout the dermis in primary lesions and higher titers of neutralizing antibodies. Finally, this adaptive immune response provided protection to these surviving rabbits against a challenge with the MYXV WT strain. Altogether, these results show that expression of the M138L gene contributes directly or indirectly to immune evasion by MYXV. In the future, these results could help us to better understand the pathogenesis of myxomatosis but also the importance of glycans in regulation of immune responses.

Viral biocontrol: grand experiments in disease emergence and evolution

Although viral emergence is commonly associated with cross-species transmission, the processes and determinants of viral evolution in a novel host environment are poorly understood. We address key questions in virus emergence and evolution using data generated from two unique natural experiments: the deliberate release of myxoma virus (MYXV) and rabbit hemorrhagic disease virus (RHDV) as biological control (biocontrol) agents against the European rabbit in Australia, and which have been of enormous benefit to Australia's ecosystem and agricultural industries. Notably, although virulence evolution in MYXV and RHDV followed different trajectories, a strongly parallel evolutionary process was observed in Australia and Europe. These biocontrol agents were also characterized by a lack of transmission to nontarget host species, suggesting that there are major barriers to successful emergence.

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.

Malaria infections reinforce competitive asymmetry between two Ficedula flycatchers in a recent contact zone

Parasites may influence the outcome of interspecific competition between closely related host species through lower parasite virulence in the host with which they share the longer evolutionary history. We tested this idea by comparing the prevalence of avian malaria (Haemosporidia) lineages and their association with survival in pied and collared flycatchers (Ficedula hypoleuca and F. albicollis) breeding in a recent contact zone on the Swedish island of Öland. A nested PCR protocol amplifying haemosporidian fragments of mtDNA was used to screen the presence of malaria lineages in 1048 blood samples collected during 6 years. Competitively inferior pied flycatchers had a higher prevalence of blood parasites, including the lineages that were shared between the two flycatcher species. Multistate mark-recapture models revealed a lower survival of infected versus uninfected female pied flycatchers, while no such effects were detected in male pied flycatchers or in collared flycatchers of either sex. Our results show that a comparatively new host, the collared flycatcher, appears to be less susceptible to a local northern European malarial lineage where the collared flycatchers have recently expanded their distribution. Pied flycatchers experience strong reproductive interference from collared flycatchers, and the additional impact of species-specific blood parasite effects adds to this competitive exclusion. These results support the idea that parasites can strongly influence the outcome of interspecific competition between closely related host species, but that the invading species need not necessarily be more susceptible to local parasites.

How does human-induced environmental change influence host-parasite interactions?

Host-parasite interactions are an integral part of ecosystems that influence both ecological and evolutionary processes. Humans are currently altering environments the world over, often with drastic consequences for host-parasite interactions and the prevalence of parasites. The mechanisms behind the changes are, however, poorly known. Here, we explain how host-parasite interactions depend on two crucial steps – encounter rate and host-parasite compatibility – and how human activities are altering them and thereby host-parasite interactions. By drawing on examples from the literature, we show that changes in the two steps depend on the influence of human activities on a range of factors, such as the density and diversity of hosts and parasites, the search strategy of the parasite, and the avoidance strategy of the host. Thus, to unravel the mechanisms behind human-induced changes in host-parasite interactions, we have to consider the characteristics of all three parts of the interaction: the host, the parasite and the environment. More attention should now be directed to unfold these mechanisms, focusing on effects of environmental change on the factors that determine encounter rate and compatibility. We end with identifying several areas in urgent need of more investigations.