High adaptive variability and virus-driven selection on major histocompatibility complex (MHC) genes in invasive wild rabbits in Australia

Characterisation of Lagovirus europaeus GI-RHDVs (Rabbit Haemorrhagic Disease Viruses) in Terms of Their Pathogenicity and Immunogenicity

Rabbit haemorrhagic disease viruses (RHDV) belong to the family Caliciviridae, genus Lagovirus europaeus, genogroup GI, comprising four genotypes GI.1-GI.4, of which the genotypes GI.1 and GI.2 are pathogenic RHD viruses, while the genotypes GI.3 and GI.4 are non-pathogenic RCV (Rabbit calicivirus) viruses. Among the pathogenic genotypes GI.1 and GI.2 of RHD viruses, an antigenic variant of RHDV, named RHDVa-now GI.1a-RHDVa, was distinguished in 1996; and in 2010, a variant of RHDV-named RHDVb, later RHDV2 and now GI.2-RHDV2/b-was described; and recombinants of these viruses were registered. Pathogenic viruses of the genotype GI.1 were the cause of a disease described in 1984 in China in domestic (Oryctolagus (O.) cuniculus domesticus) and wild (O. cuniculus) rabbits, characterised by a very rapid course and a mortality rate of 90-100%, which spread in countries all over the world and which has been defined since 1989 as rabbit haemorrhagic disease. It is now accepted that GI.1-RHDV, including GI.1a-RHDVa, cause the predetermined primary haemorrhagic disease in domestic and wild rabbits, while GI.2-RHDV2/b cause it not only in rabbits, including domestic rabbits' young up to 4 weeks and rabbits immunised with rabbit haemorrhagic disease vaccine, but also in five various species of wild rabbits and seven different species of hares, as well as wild ruminants: mountain muskoxen and European badger. Among these viruses, haemagglutination-positive, doubtful and harmful viruses have been recorded and described and have been shown to form phylogenogroups, immunotypes, haematotypes and pathotypes, which, together with traits that alter and expand their infectious spectrum (rabbit, hare, wild ruminant, badger and various rabbit and hare species), are the determinants of their pathogenicity (infectivity) and immunogenicity and thus shape their virulence. These relationships are the aim of our consideration in this article.

Invasive Bullfrogs Maintain MHC Polymorphism Including Alleles Associated with Chytrid Fungal Infection

Maintenance of genetic diversity at adaptive loci may facilitate invasions by non-native species by allowing populations to adapt to novel environments, despite the loss of diversity at neutral loci that typically occurs during founder events. To evaluate this prediction, we compared genetic diversity at major histocompatibility complex (MHC) and cytochrome b (cytb) loci from 20 populations of the American bullfrog (Rana catesbeiana) across theinvasive and native ranges in North America and quantified the presence of the pathogen Batrachochytrium dendrobatidis (Bd). Compared to native populations, invasive populations had significantly higher Bd prevalence and intensity, significantly higher pairwise MHC and cytb FST, and significantly lower cytb diversity, but maintained similar levels of MHC diversity. The two most common MHC alleles (LiCA_B and Rapi_33) were associated with a significant decreased risk of Bd infection, and we detected positive selection acting on four peptide binding residues. Phylogenetic analysis suggested invasive populations likely arose from a single founding population in the American Midwest with a possible subsequent invasion in the northwest. Overall, our study suggests that the maintenance of diversity at adaptive loci may contribute to invasion success and highlights the importance of quantifying diversity at functional loci to assess the evolutionary potential of invasive populations.

Towards a more healthy conservation paradigm: integrating disease and molecular ecology to aid biological conservation†

Parasites, and the diseases they cause, are important from an ecological and evolutionary perspective because they can negatively affect host fitness and can regulate host populations. Consequently, conservation biology has long recognized the vital role that parasites can play in the process of species endangerment and recovery. However, we are only beginning to understand how deeply parasites are embedded in ecological systems, and there is a growing recognition of the important ways in which parasites affect ecosystem structure and function. Thus, there is an urgent need to revisit how parasites are viewed from a conservation perspective and broaden the role that disease ecology plays in conservation-related research and outcomes. This review broadly focusses on the role that disease ecology can play in biological conservation. Our review specifically emphasizes on how the integration of tools and analytical approaches associated with both disease and molecular ecology can be leveraged to aid conservation biology. Our review first concentrates on disease-mediated extinctions and wildlife epidemics. We then focus on elucidating how host–parasite interactions has improved our understanding of the eco-evolutionary dynamics affecting hosts at the individual, population, community and ecosystem scales. We believe that the role of parasites as drivers and indicators of ecosystem health is especially an exciting area of research that has the potential to fundamentally alter our view of parasites and their role in biological conservation. The review concludes with a broad overview of the current and potential applications of modern genomic tools in disease ecology to aid biological conservation.

Genomic insights into a population of introduced European rabbits Oryctolagus cuniculus in Australia and the development of genetic resistance to rabbit hemorrhagic disease virus

The European rabbit (Oryctolagus cuniculus) is one of the most devastating invasive species in Australia. Since the 1950s, myxoma virus (MYXV) and rabbit haemorrhagic disease virus (RHDV) have been used to manage overabundant rabbit populations. Resistance to MYXV was observed within a few years of the release. More recently, resistance to lethal RHDV infection has also been reported, undermining the efficiency of landscape-scale rabbit control. Previous studies suggest that genetic resistance to lethal RHDV infection may differ locally between populations, yet the mechanisms of genetic resistance remain poorly understood. Here, we used genotyping by sequencing (GBS) data representing a reduced representation of the genome, to investigate Australian rabbit populations. Our aims were to understand the relationship between populations and identify possible genomic signatures of selection for RHDV resistance. One population we investigated had previously been reported to show levels of resistance to lethal RHDV infection. This population was compared to three other populations with lower or no previously reported RHDV resistance. We identified a set of novel candidate genes that could be involved in host-pathogen interactions such as virus binding and infection processes. These genes did not overlap with previous studies on RHDV resistance carried out in different rabbit populations, suggesting that multiple mechanisms are feasible. These findings provide useful insights into the different potential mechanisms of genetic resistance to RHDV virus which will inform future functional studies in this area.

Pathogen-associated selection on innate immunity genes (TLR4, TLR7) in a neotropical rodent in landscapes differing in anthropogenic disturbance

Toll-like receptors (TLRs) form part of the innate immune system and can recognize structurally conserved pathogen-associated molecular pattern (PAMP) molecules. Their functional importance in the resistance to pathogens has been documented in laboratory experimental settings and in humans. TLR diversity, however, has been rarely investigated in wildlife species. How the genetic diversity of TLRs is associated with various pathogens and how it is shaped by habitat disturbance are understudied. Therefore, we investigated the role of genetic diversity in the functionally important parts of TLR4 and TLR7 genes in resistance towards gastrointestinal nematodes and Hepacivirus infection. We chose a generalist study species, the rodent Proechimys semispinosus, because it is highly abundant in three Panamanian landscapes that differ in their degree of anthropogenic modification. We detected only two TLR7 haplotypes that differed by one synonymous single-nucleotide polymorphism (SNP) position. The TLR4 variability was higher, and we detected four TLR4 haplotypes that differed at one synonymous SNP and at three amino acid positions within the leucine-rich repeat region. Only TLR4 haplotypes had different frequencies in each landscape. Using generalized linear models, we found evidence that nematode loads and virus prevalence were influenced by both specific TLR4 haplotypes and landscape. Here, the variable “landscape” served as a surrogate for the important influential ecological factors distinguishing landscapes in our study, i.e. species diversity and host population density. Individuals carrying the common TLR4_Ht1 haplotype were less intensely infected by the most abundant strongyle nematode. Individuals carrying the rare TLR4_Ht3 haplotype were all Hepacivirus-positive, where those carrying the rare haplotype TLR4_Ht4 were less often infected by Hepacivirus than individuals with other haplotypes. Our study highlights the role of TLR diversity in pathogen resistance and the importance of considering immune genetic as well as ecological factors in order to understand the effects of anthropogenic changes on wildlife health.

Adaptive changes in the genomes of wild rabbits after 16 years of viral epidemics

Since its introduction to control overabundant invasive European rabbits (Oryctolagus cuniculus), the highly virulent rabbit haemorrhagic disease virus (RHDV) has caused regular annual disease outbreaks in Australian rabbit populations. Although initially reducing rabbit abundance by 60%, continent-wide, experimental evidence has since indicated increased genetic resistance in wild rabbits that have experienced RHDV-driven selection. To identify genetic adaptations, which explain the increased resistance to this biocontrol virus, we investigated genome-wide SNP (single nucleotide polymorphism) allele frequency changes in a South Australian rabbit population that was sampled in 1996 (pre-RHD genomes) and after 16 years of RHDV outbreaks. We identified several SNPs with changed allele frequencies within or close to genes potentially important for increased RHD resistance. The identified genes are known to be involved in virus infections and immune reactions or had previously been identified as being differentially expressed in healthy versus acutely RHDV-infected rabbits. Furthermore, we show in a simulation study that the allele/genotype frequency changes cannot be explained by drift alone and that several candidate genes had also been identified as being associated with surviving RHD in a different Australian rabbit population. Our unique data set allowed us to identify candidate genes for RHDV resistance that have evolved under natural conditions, and over a time span that would not have been feasible in an experimental setting. Moreover, it provides a rare example of host genetic adaptations to virus-driven selection in response to a suddenly emerging infectious disease.

Specific MHC class I supertype associated with parasite infection and color morph in a wild lizard population

The major histocompatibility complex (MHC) is a large gene family that plays a central role in the immune system of all jawed vertebrates. Nonavian reptiles are underrepresented within the MHC literature and little is understood regarding the mechanisms maintaining MHC diversity in this vertebrate group. Here, we examined the relative roles of parasite-mediated selection and sexual selection in maintaining MHC class I diversity of a color polymorphic lizard. We discovered evidence for parasite-mediated selection acting via rare-allele advantage or fluctuating selection as ectoparasite load was significantly lower in the presence of a specific MHC supertype (functional clustering of alleles): supertype four. Based on comparisons between ectoparasite prevalence and load, and assessment of the impact of ectoparasite load on host fitness, we suggest that supertype four confers quantitative resistance to ticks or an intracellular tickborne parasite. We found no evidence for MHC-associated mating in terms of pair genetic distance, number of alleles, or specific supertypes. An association was uncovered between supertype four and male throat color morph. However, it is unlikely that male throat coloration acts as a signal of MHC genotype to conspecifics because we found no evidence to suggest that male throat coloration predicts male mating status. Overall, our results suggest that parasite-mediated selection plays a role in maintaining MHC diversity in this population via rare-allele advantage and/or fluctuating selection. Further work is required to determine whether sexual selection also plays a role in maintaining MHC diversity in agamid lizards.

Robust Innate Immunity of Young Rabbits Mediates Resistance to Rabbit Hemorrhagic Disease Caused by Lagovirus Europaeus GI.1 But Not GI.2

The rabbit caliciviruses Lagovirus europaeus GI.1 and GI.2 both cause acute necrotizing hepatitis in European rabbits (Oryctolagus cuniculus). Whilst GI.2 is highly virulent in both young and adult rabbits, rabbits younger than eight weeks of age are highly resistant to disease caused by GI.1, although they are still permissive to infection and viral replication. To investigate the underlying mechanism(s) of this age related resistance to GI.1, we compared liver transcriptomes of young rabbits infected with GI.1 to those of adult rabbits infected with GI.1 and young rabbits infected with GI.2. Our data suggest that kittens have constitutively heightened innate immune responses compared to adult rabbits, particularly associated with increased expression of major histocompatibility class II molecules and activity of natural killer cells, macrophages, and cholangiocytes. This enables them to respond more rapidly to GI.1 infection than adult rabbits and thus limit virus-induced pathology. In contrast, these responses were not fully developed during GI.2 infection. We speculate that the observed downregulation of multiple genes associated with innate immunity in kittens during GI.2 infection may be due to virally-mediated immunomodulation, permitting fatal disease to develop. Our study provides insight into the fundamental host⁻pathogen interactions responsible for the differences in age-related susceptibility, which likely plays a critical role in defining the success of GI.2 in outcompeting GI.1 in the field.

Pathogen richness and abundance predict patterns of adaptive major histocompatibility complex variation in insular amphibians

The identification of the factors responsible for genetic variation and differentiation at adaptive loci can provide important insights into the evolutionary process and is crucial for the effective management of threatened species. We studied the impact of environmental viral richness and abundance on functional diversity and differentiation of the MHC class Ia locus in populations of the black-spotted pond frog (Pelophylax nigromaculatus), an IUCN-listed species, on 24 land-bridge islands of the Zhoushan Archipelago and three nearby mainland sites. We found a high proportion of private MHC alleles in mainland and insular populations, corresponding to 32 distinct functional supertypes, and strong positive selection on MHC antigen-binding sites in all populations. Viral pathogen diversity and abundance were reduced at island sites relative to the mainland, and islands housed distinctive viral communities. Standardized MHC diversity at island sites exceeded that found at neutral microsatellites, and the representation of key functional supertypes was positively correlated with the abundance of specific viruses in the environment (Frog virus 3 and Ambystoma tigrinum virus). These results indicate that pathogen-driven diversifying selection can play an important role in maintaining functionally important MHC variation following island isolation, highlighting the importance of considering functionally important genetic variation and host-pathogen associations in conservation planning and management.

Resistance to RHD virus in wild Australian rabbits: Comparison of susceptible and resistant individuals using a genomewide approach

Deciphering the genes involved in disease resistance is essential if we are to understand host-pathogen coevolutionary processes. The rabbit haemorrhagic disease virus (RHDV) was imported into Australia in 1995 as a biocontrol agent to manage one of the most successful and devastating invasive species, the European rabbit (Oryctolagus cuniculus). During the first outbreaks of the disease, RHDV caused mortality rates of up to 97%. Recently, however, increased genetic resistance to RHDV has been reported. Here, we have aimed to identify genomic differences between rabbits that survived a natural infection with RHDV and those that died in the field using a genomewide next-generation sequencing (NGS) approach. We detected 72 SNPs corresponding to 133 genes associated with survival of a RHD infection. Most of the identified genes have known functions in virus infections and replication, immune responses or apoptosis, or have previously been found to be regulated during RHD. Some of the genes identified in experimental studies, however, did not seem to play a role under natural selection regimes, highlighting the importance of field studies to complement the genomic background of wildlife diseases. Our study provides a set of candidate markers as a tool for the future scanning of wild rabbits for their resistance to RHDV. This is important both for wild rabbit populations in southern Europe where RHD is regarded as a serious problem decimating the prey of endangered predator species and for assessing the success of currently planned RHDV variant biocontrol releases in Australia.