Next Article in Journal
Ebola Virus Causes Intestinal Tract Architectural Disruption and Bacterial Invasion in Non-Human Primates
Previous Article in Journal
Baculovirus as a Tool for Gene Delivery and Gene Therapy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
CSIRO Australian Animal Health Laboratory, Geelong, VIC 3220, Australia
2
CSIRO Health and Biosecurity, Canberra, ACT 2601, Australia
3
CSIRO Agriculture and Food, Canberra, ACT 2601, Australia
*
Author to whom correspondence should be addressed.
Submission received: 16 August 2018 / Revised: 12 September 2018 / Accepted: 15 September 2018 / Published: 19 September 2018
(This article belongs to the Section Animal Viruses)

Abstract

:
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.

1. Introduction

Rabbit hemorrhagic disease (RHD) is a peracute and often lethal hepatitis caused by the rabbit caliciviruses Lagovirus europaeus GI.1 (previously called RHDV) and Lagovirus europaeus GI.2 (previously called RHDV2 or RHDVb) [1,2,3]. Both viruses are positive-sense RNA viruses belonging to the genus Lagovirus, family Caliciviridae. GI.1 was first reported in China in 1984, and spread rapidly through European rabbit (Oryctolagus cuniculus) populations in Asia and Europe, with outbreaks also reported from North Africa and the Americas in the 1980s [1]. GI.1 was released in Australia in the mid-1990s to control wild rabbit populations, and spread to New Zealand in 1997 [4]. It is now considered endemic in wild rabbit populations in Asia, Europe, Australia, New Zealand, and parts of Africa [1]. Subsequently, a novel lagovirus, GI.2, emerged in France in 2010 [2]. GI.2 has since been detected throughout Europe and nearby islands [2,5,6,7,8,9,10,11,12], and in Australia in 2015 [13]. GI.2 is now endemic in Europe and Australia, and appears to be replacing GI.1 strains in these regions [5,10,14,15].
GI.2 is antigenically and genetically distinct from GI.1, and is able to cause disease in GI.1-vaccinated rabbits and in various hare species, including Italian hares (Lepus corsicanus), Sardinian cape hares (Lepus capensis mediterraneus), and European brown hares (Lepus europaeus) [2,5,16,17,18,19,20]. Both viruses have an identical genome organisation comprising two open reading frames (ORFs). ORF1 encodes seven non-structural proteins and the major capsid protein VP60, while ORF2 encodes a minor structural protein, VP10 [1]. Interestingly, multiple GI.2 recombinants have been reported, containing non-structural gene sequences of either GI.2, GI.1, or GI.4 (a group largely comprised of non-pathogenic lagoviruses) with a recombination breakpoint located at the junction between the non-structural genes and VP60 [15,21]. All of these recombinants have been reported to be virulent in rabbits as well as hares [22]. This suggests that determinants of pathogenicity, as well as the broader host range of GI.2 viruses, lie within the structural gene sequences.
In natural infections, the virus usually enters the host via mucosal membranes after ingestion [23]. However, experimentally, animals can also be infected through parenteral inoculation. Orally-infected rabbits showed an increased survival time compared to animals infected by intramuscular or intradermal injection; however, it is not known whether this is due to a prolonged disease duration or to delayed virus entry in orally-infected animals [24]. Typically, over 90% of adult rabbits infected with GI.1 develop pyrexia within 24 to 36 h post infection (hpi) and die 12 to 36 h later [1]. In contrast, disease associated with early (2010 to 2011) strains of GI.2 was frequently subacute, with mortalities occurring three to nine days post-infection, and highly variable case fatality rates between experiments [16]. Pathogenicity studies on more recent GI.2 strains (2014 to 2015) showed case fatality rates comparable with GI.1, although low numbers of animals were tested [25]. Typical necropsy lesions associated with both GI.1 and GI.2 viruses include extensive hepatic necrosis, initially with a predominantly periportal distribution, splenomegaly, and hemorrhages throughout most solid organs resulting from disseminated intravascular coagulation (DIC; reviewed in Reference [1]).
Studies on lagovirus replication have thus far been impeded by the lack of a reliable cell culture system, and in vivo studies on pathogenicity mechanisms have largely been limited to GI.1. The pathogenesis of GI.2 is presumed to follow a similar pattern to that of GI.1, although the chronology of infection may differ. Cellular receptors have not yet been identified for either virus, although histo-blood group antigens (HBGAs), complex carbohydrates linked to glycoproteins or glycolipids that are expressed on epithelial cells of the trachea and duodenum in rabbits, have been identified as co-receptors/attachment factors for lagoviruses [26,27,28,29,30,31]. Viral replication is presumed to occur predominantly in hepatocytes, and GI.1 antigen can be detected in these cells as early as 12 hpi [32,33]. Viral antigens have also been detected in cells of the monocyte/macrophage lineage, including Kupffer cells and intravascular macrophages, most commonly in the lungs and spleen [33,34,35,36,37,38,39]. Infected hepatocytes initially show hydropic change and ultimately become apoptotic and necrotic [34,37,40,41,42]. Necrotic foci are associated with infiltration of inflammatory cells, predominantly heterophils [18,37,41,43]. Necrotic foci coalesce as infection proceeds and this widespread hepatic damage leads to typical biochemical changes associated with hepatic disease, including marked elevations in serum transaminases and bilirubin [43,44]. Subsequently, infected rabbits develop DIC, characterized by thrombocytopenia and prolongations in blood clotting times [44,45,46].
In contrast to the severe pathology induced in adult rabbits, young rabbits are innately resistant to disease caused by GI.1, although they can be infected and seroconvert, even after the inoculation of low infectious doses [47,48,49]. This resistance to disease is independent of maternal antibody status and is gradually lost over time, with case-fatality rates increasing in a non-linear fashion until approximately nine weeks of age, when case-fatality rates are comparable to those for adult rabbits [23,50]. The viral antigen can only be detected in a limited number (0.01 to 0.2%) of hepatocytes in young rabbits, compared to over 60% of hepatocytes in adult animals [33,41,50,51,52]. Furthermore, the levels of virus RNA in the liver of young rabbits sacrificed at five dpi are approximately 1000-fold lower than those in moribund adult rabbits [49]. Microscopically, foci of infected hepatocytes show a similar appearance to that reported in adults, and are surrounded by small infiltrates of inflammatory cells [33,41,50]. Viral antigen has also been detected in macrophages from infected kittens and adults [33,50]. In accordance with the limited liver pathology, infected kittens show only mild increases in serum transaminases compared to those observed in adult rabbits, and other hematological and biochemical parameters, such as bilirubin, blood glucose concentration, and measures of coagulation, remain unchanged [41,51]. The mechanism underlying this resistance of young rabbits is not understood, although various suggestions have been proposed. These include (i) a developmental antigen that functions as the cellular receptor for the virus [30,53]; (ii) structural or functional changes to the liver, possibly due to dietary changes at weaning, that induce susceptibility [33,50]; or (iii) differences in innate immune responses between young and adult rabbits that drive pathogenesis [51,52].
To further investigate the mechanisms underlying the innate resistance of kittens to GI.1-induced disease, we conducted a genome-wide analysis of liver transcripts from young and adult rabbits infected with GI.1, and kittens infected with GI.2. This allowed us to assess both potential virus-related factors (i.e., difference in responses between GI.1 and GI.2 infection in kittens) as well as potential age-related host factors (i.e., difference in responses between adults and kittens during GI.1 infection) contributing to the resistance observed in kittens to GI.1 RHD. Here, we discuss those genes showing strong differential regulation between the experimental groups and the corresponding pathways involved in host responses to GI.1 and GI.2 infections, and compare transcriptional responses between young and adult rabbits.

2. Materials and Methods

2.1. Experimental Treatments and Virus Inoculation

Animal experiments were conducted at the Commonwealth Scientific and Industrial Research Organisation (CSIRO) Black Mountain Laboratories following the Australian Code for the Care and Use of Animals for Scientific Purposes (2013) and approved by the CSIRO Ecosystem Sciences Animal Ethics Committee (permit identifiers: CESAEC DOMRAB, SEAEC 10-12, ESAEC 13-10). New Zealand white rabbits (36 days old or >10 weeks old) bred from GI.4 and GI.1 antibody-free does were housed either individually or in litter groups. All rabbits were confirmed negative for antibodies to the non-pathogenic calicivirus GI.4 [54] and were healthy prior to inclusion in the study. Five experimental groups were delineated: (i) six control kittens sham inoculated with phosphate-buffered saline (PBS) and with tissues harvested at 12 (n = 3) and 24 (n = 3) hpi; (ii) eight kittens infected with GI.1 and tissues harvested at 12 (n = 4) and 24 (n = 4) hpi; (iii) seven kittens infected with GI.2 and tissues harvested at 12 (n = 3) and 24 (n = 4) hpi; (iv) seven adult rabbits infected with GI.1 and tissues harvested at 12 (n = 3) and 24 (n = 4) hpi; and (v) five control adult rabbits, four of which were sham inoculated with PBS with tissues harvested at 12 (n = 2) and 24 (n = 2) hpi, as well as a retired breeding doe that was used as an additional control (Table 1). Animals were infected orally with clarified liver homogenate produced from infected rabbits diluted to either 3 × 108 capsid gene copies of GI.2 or 1500 LD50 of GI.1 (equivalent to approximately 1.5 × 108 capsid gene copies). These doses are both considered to be high infectious doses. Control animals received 1 mL of PBS. The GI.2 and GI.1 strains used for infection were the Australian field isolate BlMt-1 (GenBank #KT280060), and a commercially available titrated preparation of strain Czech-351 (GenBank #KF594475, Elizabeth Macarthur Agriculture Institute, Menangle, Australia), respectively. The Czech-351 isolate has been reported to cause a high case-fatality rate in adult domestic rabbits and a low case-fatality rate in 5 week old kittens [1,49]. A recent study reported a very high case-fatality rate of BlMt-1 in 5 week old domestic rabbits [55]. Tissues were harvested at 12 or 24 hpi, after rabbits had been anesthetised by intramuscular injection of either Zoletil 100 (Virbac, Peakhurst, NSW, Australia) or a combination of 20 mg/mL xylazine hydrochloride (Troy Laboratories, Smithfield, NSW, Australia) and 100 mg/mL ketamine hydrochloride (Mavlab, Logan, QLD, Australia), and then euthanized by intravenous or intracardiac injection of 325 mg/mL sodium pentobarbitone (Virbac). Infection times of the 12 h and 24 h groups were staggered such that all animals were sacrificed at similar times of the day to avoid any differences in gene expression caused by diurnal fluctuations in rabbit metabolism. Tissues were collected at necropsy and stored in RNAlater (Qiagen, Chadstone Centre, VIC, Australia) at −80 °C until processing.

2.2. RNA Extraction

Twenty to 30 mg of liver from each rabbit was homogenized with 1-mm glass beads (Daintree Scientific, St Helens, TAS, Australia) using a Precellys 24-dual tissue homogenizer (Bertin Technologies, Montigny-le-Bretonneux, France). RNA was then extracted using either the Qiagen RNeasy mini kit (Qiagen, Chadstone Centre, VIC, Australia) or the Maxwell 16 LEV simplyRNA tissue kit (Promega, Sydney, NSW, Australia) as per the manufacturers’ instructions.

2.3. Virus Quantification

Viral load was quantified in terms of the “capsid gene copy number” using the quantitative reverse transcription polymerase chain reaction (RT-qPCR) method that detects both GI.1 and GI.2, as described previously [22]. Briefly, reactions were performed in duplicate using the SensiFAST SYBR No-ROX One-Step kit (Bioline, Alexandra, NSW, Australia) on a BioRad CFX96/C1000 thermal cycler platform using the primers RHDV all 3 Fw 5’-TTGACRTACGCCCTGTGGGACC-3′ and RHDV all 3b Rv 5’-TCAGACATAAGARAAGCCATTRGYTG-3′ [22]. Absolute quantification was performed using a standard curve generated from in vitro generated full length GI.1 RNA transcripts of known concentration [22].

2.4. RNA Sequencing and Analysis

Messenger RNA (mRNA) was enriched in the total RNA samples using polyA selection and sequenced using two 75 bp single-end NextSeq lanes (Illumina, Scoresby, VIC, Australia) at the Australian Cancer Research Foundation (ACRF) Biomolecular Resource Facility (The John Curtin School of Medical Research, Australian National University, Canberra, Australia). RNA-Seq reads were cleaned by removing Illumina TruSeq adapters, trimming sequences when the quality score dropped below 20, and discarding reads shorter than 50 bp using Trimmomatic v.0.32 [56]. The cleaned reads were then mapped to the rabbit genome (OryCun2.0; http://www.ensembl.org/Oryctolagus_cuniculus) and rabbit calicivirus genomes (GI.1, GenBank #KF594473; GI.2, GenBank #KT280060) using TopHat v.2.1.0 [57]. Counts for each gene in the rabbit genome were calculated using HTSeq v.0.6.0 [58] with the “union” mode for processing overlapping and multi-mapping reads. The counts were then imported into edgeR v.3.12.1 [59] for testing of significant expression differences across the treatments. The control rabbits within each age category had very similar expression profiles at both 12 and 24 h time points (Figure 1) and, therefore, these time points were combined to generate a single control pool for each age group for differential expression testing. The differentially expressed genes were tested for functional pathway enrichment using the Gene Ontology vocabulary [60] with Fisher’s exact tests in the R package topGO v.2.22.0 [61] and the KEGG framework [62] using the kegga function within limma v.3.32.1 [63]. Reads per kilobase per million (rpkm) values were calculated using the fpkm function in edgeR v.2.12.1 [59] and plotted with ggplot v.2.1.0 [64]. Raw coverage of the viral genomes was obtained by mapping the cleaned reads to the GI genomes using Bowtie v.2.2.9 [65], extracting the coverage profile with Samtools v.1.3.1, and plotting with Gviz v.1.14.2 [66].
Gene co-expression networks were constructed using the R package “WGCNA” (Weighted Gene Co-expression Network Analysis), following the authors’ recommendations [67]. Briefly, significantly differentially expressed genes were filtered by removing genes with low counts and log transformed. The function “pickSoftThreshold” was used to plot the data, and a soft thresholding power of 10 was chosen based on examination of the plot. The genes were clustered using the function “hclust” and the tree was cut into modules using “cutreeDynamic” with a minimum cluster size of 10. The modules were further refined using “mergeCloseModules” and manual examination of dendrograms using “plotDendroAndColors”. Co-expression modules were then tested for significant associations with the rabbit treatment groups using the R base function “cor” and plotted using “labeledHeatmap” [67].

3. Results

3.1. Genome Mapping and Significance Testing

The transcriptional response of rabbits to lagovirus infection was investigated 12 and 24 hpi in juveniles and adults, and compared to control animals for each age group. Messenger and viral RNA were enriched from the liver of each rabbit and RNA-seq was conducted, producing several million high-quality reads per sample (Table 1). More than 90% of these reads were successfully aligned to the rabbit genome, providing good coverage for the calculation of expression values (Table 1). Viral reads were only detected in rabbits at 24 hpi, indicating that viral titres at 12 hpi were too low to be detected using the sequencing depths obtained here (Table 1). Although there was considerable variation between individual animals, adults 24 h post-GI.1 infection, and kittens 24 h post-GI.2 infection, contained the highest viral loads (Table 1), which is consistent with the higher levels of virulence described for GI.1 in adult compared to young rabbits, and the recently described high case-fatality rate reported for 5 week old rabbits infected with the Australian GI.2 isolate [55]. These broad patterns were validated using RT-qPCR assays (Table 1). Viral genome copy numbers varied between individuals in each treatment group, likely representing individual variation in outbred rabbits very early during infection, although we cannot be certain that all GI.1 infected adults and all GI.2 infected kittens would have succumbed to infection. In all of the samples, viral genome coverage was consistent for approximately 5.3 kb, before a large increase in coverage at the 3’ end (Figure S1). This was likely due to the sequencing of multiple 2.1 kb subgenomic viral RNAs (co-linear with the 3’ end of the viral genome) that are produced by lagoviruses in large quantities early in infection [68]. Despite the individual variation between animals in each group, replicate rabbits within each of the treatment groups had broadly similar expression profiles, which provided high power for significance testing across the groups (Figure 1). The control samples were also very similar to each other, regardless of whether they were collected at 12 or 24 hpi (Figure 1). For this reason, the control samples at each time-point were combined into a single treatment group for the remainder of the analysis.
The number of differentially expressed genes in GI-infected animals correlated with expected disease severity (Table 2; Figure S2). For example, kittens that were resistant to disease caused by GI.1 infection (no disease/low severity), only differentially expressed 62 genes at 24 hpi. The overall transcriptional profile of this treatment group was similar to that in control animals (Figure 1). On the other hand, adult rabbits, which were susceptible to GI.1 (high severity), differentially expressed 322 genes at 24 hpi, and the transcriptional profile of infected animals at 24 hpi was generally distinct from the controls (Figure 1). In contrast to GI.1, GI.2 reached high titres in kittens (high severity; Table 1). In this case, changes in the transcriptomic profiles could be seen after only 12 hpi, and, by 24 hpi, there were 4211 dysregulated genes in the kittens, the most of any treatment group (Table 2).

3.2. Differential Gene Expression and Pathway Enrichment in Adults and Kittens

Several hundred genes were differentially expressed between uninfected kittens and uninfected adults (Table 2). Unsurprisingly, many of the genes upregulated in kittens compared to adults were related to growth and maturation, such as genes involved in the cell cycle, e.g., HN1 and CDK1, and DNA replication, e.g., NCAPH (Table 3). Interestingly, the polymeric immunoglobulin receptor PIGR was significantly upregulated in adults compared to kittens (Table 3).
During GI.1 infection, kittens upregulated several immune related genes, including interferon induced proteins (MX1, IFI44, and IFIT5) and EPCAM (Table 4). Similarly, many of the most differentially expressed genes in GI.1-infected adults were also involved in immune processes, including MX2, CD80 (B7-1), IFIH1, and CXCL10 (Table 4). Functionally, both adults and kittens differentially regulated pathways such as “defense response to virus,” “negative regulation of viral genome replication,” and “immune response” (Table 5), likely reflecting the alteration of broad generic stress response pathways. At 12 hpi, an insufficient number of genes were differentially expressed to conduct pathway analyses.

3.3. Kittens Upregulate Important Components of Innate Immunity Compared to Adults, Which May Limit GI.1-Induced Pathology

To examine differences in the response of adults and kittens to GI infection, we used gene co-expression networks to find patterns in the data. One particular co-expression cluster (the “pink” module) contained genes that were more highly expressed by uninfected kittens compared to uninfected adults, further upregulated in GI.1-infected kittens, but markedly downregulated in GI.2-infected kittens (Figure 2; Figure S3). This pattern may suggest that the co-expressed genes were involved in the resistance of kittens to disease caused by GI.1. Functionally, the cluster was enriched for signal transduction, natural killer (NK) cell regulation, and positive GTPase regulation (Figure 2), including genes such as HLA-DPB1, TNFRSF11B, DOK3, and OSGIN1 (Figure S3). Several other genes in the cluster encoded Rho-GTPases, which have a key regulatory role in immune responses [69].
Of particular interest was the upregulation of major histocompatibility (MHC) class II genes (e.g., HLA-DPB1), which are an important component of the immune system. In fact, when genes of the MHC complex were analyzed specifically, MHC class II genes tended to be expressed at a higher rate in uninfected kittens compared to adults, and their expression levels further increased in GI.1-infected kittens, beginning as early as 12 hpi (Figure 3). In contrast, the expression of many MHC II genes decreased in GI.2-infected kittens (Figure 3).
NK cells can express MHC class II genes, and NK cell regulation was enriched in the “pink” co-expression cluster [70]. Consequently, we specifically examined the expression of rabbit genes involved in NK cell regulation and found that several key genes, including PTPN22, VAV1, and ARRB2, tended to be upregulated in control kittens compared to control adults (as seen at time-point 0), and were upregulated or maintained their expression levels in GI.1-infected kittens, but were downregulated in GI.2-infected kittens (Figure 4), reflecting the co-expression pattern observed for MHC class II genes (Figure 3). Moreover, natural cytotoxicity triggering receptor 3 (NCR3) was also upregulated in GI.1-infected kittens compared to GI.2-infected kittens (Figure 4). Macrophages also contribute to a higher expression of MHC class II genes [71]. Thus, we examined the expression of rabbit genes likely to be associated with macrophage function, and found that a number of these genes, including CSF1R, ZBTB46, and CD68, tended to be upregulated by 24 hpi (if not earlier) in GI.1-infected kittens compared with GI.1-infected adults and GI.2-infected kittens (Figure 4). Another type of cell that may have contributed to MHC class II expression in GI.1-infected kittens are cholangiocytes (biliary epithelial cells). These cells not only express MHC class II genes, but also many other genes that were upregulated in GI.1-infected kittens, such as genes encoding MX proteins, interferons and adhesion molecules [72] (Table 4 and Table 5). In fact, EPCAM, which is a biomarker of cholangiocytes and their progenitor cells [72,73], was the most up-regulated gene at 12 hpi in GI.1-infected kittens (Table 4; Figure 4). In addition, the gene encoding FUT1, an enzyme belonging to the family of alpha 1,2 fucosyltransferases, was found to be upregulated in control kittens compared to adults and further upregulated in GI.1-infected kittens, but downregulated in GI.2-infected kittens (Figure 4).
In the context of the KEGG pathway “Antigen Processing and Presentation”, the up-regulation of MHC class I and II genes in GI.1-infected kittens coincided with increased expression of markers associated with NK cells, such as KLRD1 and NKG2-A/NKG2-B type II integral membrane protein, and CD4 T-cells, including the T-cell surface glycoprotein, CD4 (Figure 5). Adults infected with GI.1 also tended to slightly up-regulate MHC genes, and these were associated with minor increases in NK cell markers, however, the magnitude was much less than for kittens (Figure 3 and Figure 5).

3.4. In Contrast to GI.1, GI.2 Infection Restricts the Activation of Several Innate Immune Pathways

GI.2-infected kittens differentially expressed far more genes than any other treatment group (Table 2; Figure S2). Similar to other experimental groups, however, GI.2-infected kittens upregulated generic viral response pathways involved in the negative regulation of viral genome replication, viral defense responses, and more general immune and inflammatory responses (Tables S1 and S2). A number of genes involved in immune responses were amongst the most up-regulated in GI.2-infected kittens, including FLVCR1, IRAK2, and the cell adhesion molecules VCAM1 and ICAM1 (Table S1). However, many of the innate immunity genes that were upregulated in GI.1-infected kittens were downregulated in GI.2-infected kittens (Figure 3, Figure 4 and Figure 5).
The downregulation of genes encoding MHC complex molecules in GI.2-infected kittens coincided with decreased expression of markers associated with CD8 T cells, NK cells, macrophages, and CD4 T cells (Figure 4 and Figure 5).

4. Discussion

We undertook a genome-wide transcriptome study to investigate differences in the immune responses of young and adult rabbits during very early stages of infection with two caliciviruses, GI.1 and GI.2. Interestingly, while adults are susceptible to GI.1, young rabbits are resistant [1]. In contrast, both age groups are susceptible to fulminant disease caused by GI.2 [2,5].
Several components of the innate immune system were identified as potentially important in the resistance of kittens to GI.1, including MHC class II genes. These genes were constitutively expressed at a higher rate in uninfected kittens compared to uninfected adults, suggesting that kittens have a “primed” innate immune system. Upon infection of kittens with GI.1, a further upregulation of MHC class II alleles was seen from as early as 12 hpi, indicating a rapid and coordinated innate immune response. The precise nature of this immune response, however, depends on the type of cell expressing the MHC II molecules. For example, MHC class II molecules are cell surface proteins generally expressed by antigen presenting cells (APCs), such as macrophages, natural killer (NK) cells and B cells [70,74]. In this study, macrophages and NK cells were likely to be the primary MHC II-producing cells due to the early sampling times, which precluded the development of B cell populations specific to viral antigens. We also found several markers suggesting that cholangiocytes, which are also MHC II-producing cells and potentially APCs [72,75], may play an important role in the innate resistance of kittens to GI-induced disease. The upregulation of MHC II genes indicates that either individual APCs expressed more of these molecules, or alternatively, that a greater number of APCs were present. Although high throughput sequencing cannot discriminate between these possibilities, a combination of the two processes is likely. In either case, our findings suggest that the increased activity of NK cells, macrophages, and/or cholangiocytes, have a critical role in limiting GI.1-induced pathology in young rabbits.
NK cells provide front-line innate immunity and are well known for controlling viral spread. Indeed, individuals with lower NK cell numbers or reduced NK cell-mediated cytotoxicity are far more susceptible to certain viral infections [76,77]. Here, in addition to higher MHC II expression, we found several other indications that NK cells play an important role in the resistance of kittens to GI.1. For example, a co-expressed cluster of genes containing several involved in NK cell activity, were highly expressed in GI.1-infected kittens compared to GI.1-infected adults and GI.2-infected kittens. Genes in this cluster included natural cytotoxicity triggering receptor 3 (NCR3 or NKp30), which is a specialized receptor on the surface of NK cells that plays a key role in the recognition and destruction of virus-infected cells [78]. Several other potentially important genes in this cluster may be expressed as part of the same process, including PTPN22, VAV1, and EPCAM. For example, PTPN22 encodes lymphoid protein tyrosine phosphatase, a protein expressed only in immune cells, and most highly expressed by NK cells [79,80]. PTPN22 also induces cytotoxic activity of NK cells through the activation of an array of proteins including the guanine nucleotide-exchange factor, VAV1 [81]. VAV1 facilitates the reorganization of actin cytoskeletal proteins and augments the expression of cell adhesion molecules on cell surfaces [81,82]. In line with this, many of the most differentially expressed genes in infected rabbits were cell adhesion molecules, including EPCAM (epithelial cell adhesion molecule). In fact, EPCAM was the most differentially expressed gene in GI.1-infected kittens at 12 hpi, correlating with increases in VAV1 expression. Moreover, a number of different GTPases were part of the co-expressed gene cluster, which may be used to boost NK cell cytotoxicity [82]. Overall, the coordinated expression of these genes suggests an orchestrated effort to increase the activity and cytotoxicity of NK cells in kittens and points to an important role for NK cells in GI.1-resistant kittens.
Macrophages may also play a role in the resistance of kittens to disease caused by GI.1. Indeed, a greater abundance of macrophages in the liver of control kittens compared to control adults, with further increases in kittens after GI.1 infection, could account for the MHC II expression patterns that we observed. Moreover, the expression of CSF1R (macrophage colony-stimulating factor 1 receptor) was higher in uninfected kittens compared to uninfected adults, and steadily increased over the GI.1 infection time-points, but remained stable in GI.2-infected kittens. CSF1R is a transmembrane protein that is almost ubiquitously expressed on mononuclear phagocytes, and can be used to differentiate macrophages from most dendritic-cell types [83,84]. We also found that the “classical” macrophage marker, CD68 [85], and a more recently recognized marker, ZBTB46 [86], tended to be most highly expressed in GI.1-infected kittens at 24 hpi, further suggesting a role for macrophages in clearing GI.1 virus from kittens.
In the present study, there were several indications that cholangiocytes might also be important in the resistance of kittens to GI.1-induced disease. Cholangiocytes can express MHC II molecules and comprise a relatively large proportion of cell mass in the liver (4–5%) [72]. Traditionally, the immunological contribution of these cells was thought to be restricted to immunoglobulin secretion. More recently, however, cholangiocytes have been implicated in many innate and adaptive immune processes, such as the recognition of pathogen-associated molecular patterns (PAMPs), secretion of cytokines and antimicrobial peptides, and interactions with other immune cells through the expression of cell-surface adhesion molecules (for a review, see Reference [72]). EPCAM, which was highly expressed in GI.1-infected kittens at 12 hpi, is frequently used as a biomarker of cholangiocytes [72,73], and EPCAM receptors are present on immune modulatory cells such as lymphocytes, monocytes, dendritic cells, and NK cells [87]. Moreover, cholangiocytes express MX proteins [72], and MX1 was one of the most upregulated genes in GI.1-infected kittens at 24 hpi, along with several interferons. MX proteins, which are induced by interferons, have broad antiviral activity against RNA viruses through the recognition of viral nucleocapsid proteins and subsequent inhibition of viral replication [88]. Although antiviral activity of MX proteins has never been demonstrated for GI viruses, it is tempting to speculate that the upregulation of MX-1 constitutes a possible mechanism contributing to the observed reduction of GI.1 replication in kittens. Cholangiocytes may also be associated with other immune responses, such as the observed upregulation of macrophage markers. For instance, cholangiocytes can secrete molecules that attract monocytes and macrophages during liver disease [72,89].
Of note was the differential expression of FUT1, which encodes an enzyme belonging to the family of alpha 1,2 fucosyltransferases. These enzymes are required for the synthesis of the histo-blood group antigen (HGBA) H type 2, which has been described as a co-receptor for lagovirus binding [28]. HBGA H type 2 is present on mucosal surfaces, such as epithelial cells of the trachea and duodenum, but is not expressed in liver parenchyma [29,30]. FUT1 expression levels were elevated in kittens compared to adults, upregulated in GI.1-infected kittens at 24 hpi, but downregulated in GI.2-infected kittens as early as 12 hpi. This could suggest that differential fucosylation patterns of immune cell ligands and receptors may contribute to the attenuated GI-1 infections in kittens.
We propose that young kittens have a primed innate immune system that perhaps compensates for the lack of a fully developed adaptive immune system. Possible causes for elevated innate immune responses in rabbit kittens are unclear but warrant further investigations. Exposure to a suite of different antigens while the rabbit’s complex gut microbiome is establishing after birth is a possibility. It is also of note that the waning of age-related resistance to lethal GI.1 infection coincides with weaning and a change from a high-fat diet to a vegetarian diet. Links between high-fat diets and inflammation have been demonstrated in adult rabbits used as a model for human diabetes [90], although no such links have yet been reported in very young rabbits.
The consequences of this primed innate immunity includes coordinated increases in the expression of genes associated with NK cell, macrophage, and cholangiocyte activity, which may allow kittens to rapidly respond to and limit infection with GI.1. Indeed, previous studies support this assertion. First, treatment of adult rabbits with poly(I:C) to artificially “prime” the innate immune system by stimulating type 1 interferon responses can protect against GI.1 challenge for up to six hours post-treatment [91]. Second, adult rabbits treated with the pro-inflammatory cytokine cardiotrophin-1, a member of the IL-6 cytokine family, prior to infection with GI.1 had lower case fatality rates, improved biochemical parameters of liver function, and less severe hepatic necrosis based on histopathology, although virus loads in the liver reached similar titres to untreated rabbits [92]. Thirdly, GI.1 infection of kittens was associated with increased levels of circulating pro-inflammatory cytokines as early as 6 hpi, along with a rapid influx of macrophages and lymphocytes into the liver, although these results were not compared to adult rabbits [4]. Perhaps most convincingly, immunosuppression of young rabbits with corticosteroid treatment negated their innate resistance to GI.1-induced disease, with kittens developing widespread hepatic necrosis indistinguishable from that of adult animals and succumbing to disease in 24 to 72 hpi [52]. It should be noted, however, that corticosteroids also have wide-ranging systemic effects on multiple metabolic pathways, and treatment may lead to functional changes in the liver that could increase susceptibility to infection. This innate resistance of young rabbits to GI.1-induced disease is not observed with GI.2; indeed, mortalities have been observed in rabbit kittens as young as 11 days old [93].
Compared to kittens infected with GI.1, GI.2-infected kittens had decreased expression of several MHC class I genes, including those encoding the MHC class I structural protein, B2M (beta-2-microglobulin). B2M forms part of the alpha chain of the MHC class I complex, and is required for the assembly and cell surface expression of functional MHC class I molecules [74]. The limited expression of these pathways during GI.2 infection suggests that GI.2 may suppress host innate immune responses, allowing the virus to replicate rapidly. Indeed, there were more viral genome copies at 24 hpi in GI.2-infected kittens compared to any other treatment group (Table 1). Some viruses intentionally downregulate host MHC class I molecules to avoid recognition by CD8 T cells [94]. The consequence of this downregulation, however, is that the virus-infected cells become susceptible to NK cell attack [95]. In turn, certain viruses have evolved the ability to preferentially express, or to mimic the expression of, ligands that inhibit NK cell receptors [96]. While these examples involve large double-stranded DNA viruses, small RNA viruses have also evolved mechanisms to avoid host immune responses. For example, murine norovirus (MNV), which belongs to the same family as lagoviruses (Caliciviridae), expresses a protein that delays the upregulation of genes involved in host innate immune responses [97]. Thus, although host immune modulation by GI.2 has not been demonstrated, caliciviruses as a family are known to subvert host responses, and the functions of several of the non-structural proteins of the lagoviruses are as yet unknown [98]. Moreover, variability in MHC class I genes has been implicated in the genetic resistance of rabbits to GI.1 viruses [99]. It is feasible that when genetic resistance against GI.2 begins to develop in wild rabbit populations, the mechanism may be different. These examples provide some suggestion that GI.2 may have the ability to modulate the expression of MHC class I molecules, and enable the virus to replicate and spread before the development of an immune response. This scenario should be tested in future studies.

5. Conclusions

Overall, we found that kittens infected with GI.1 increased the expression of multiple genes encoding components of the innate immune response compared with adult rabbits, particularly those associated with MHC II genes, such as natural killer cells, macrophages, and cholangiocytes. In contrast, these genes were downregulated in kittens during GI.2 infection, suggesting that these genes, including MHC class II alleles and those involved in NK cell regulation, play a role in the pathogenicity differences of GI.1 and GI.2. We propose a model where young kittens have a primed innate immune system, perhaps due to the increased load of environmental antigens they encounter at this age, or due to a high-fat milk diet stimulating innate immune responses. This primed innate immunity allows them to respond rapidly to the incoming viral pathogen. These innate responses may limit GI.1 replication to a point where only individual hepatocytes become infected, precluding progression to the coalescing hepatitis seen in adult animals. In contrast, GI.2 may be able to suppress these innate responses sufficiently to permit extensive viral replication, leading to the fulminant hepatitis observed in GI.2-infected kittens, although the mechanism by which GI.2 induces this suppression requires further elucidation. This model is supported by our data and earlier work done by others demonstrating the importance of a robust immune system in the natural resistance of young rabbits to GI.1-induced disease [52]. However, these findings need to be validated experimentally in the future, for example, by demonstrating changes in expression levels of key proteins implicated in this study. The fundamental host and viral mechanisms underlying the resistance of young rabbits to disease induced by GI.1, but not GI.2, also has important epidemiological implications. The ability of GI.2 to infect new cohorts of rabbits at a much younger age likely represents a key competitive advantage over GI.1. Our findings therefore represent an important step in increasing our understanding of how an emerging pathogen (i.e., GI.2) can displace a highly successful existing pathogen (GI.1) in a widely dispersed host population. This may ultimately lead to improved management of wild rabbit populations, either to better protect domestic and wild rabbits in their native range, or to improve the management of overabundant rabbit populations in parts of the world where these are considered an invasive pest.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/10/9/512/s1, Figure S1: Read coverage of the GI.1 and GI.2 genomes, Figure S2: Significantly differentially expressed genes in infected rabbits compared to controls. Figure S3: Expression of genes in the “pink” co-expression module and change in selected genes over time. Table S1: The top 10 most differentially expressed genes in GI.2-infected rabbits. Table S2: Significantly enriched Gene Ontology (GO) terms in the category “Biological Process” for up regulated genes in GI.2-infected rabbits.

Author Contributions

Conceptualization, M.J.N., P.K., J.T., and T.S.; Data curation, M.J.N.; Formal analysis, M.J.N.; Investigation, N.H., P.K., M.H., and T.S.; Methodology, M.J.N., R.N.H., K.A.M., P.K., J.T., and T.S.; Supervision, R.N.H. and T.S.; Validation, M.J.N.; Visualization, M.J.N.; Writing—original draft, M.J.N. and R.N.H.; Writing—review and editing, K.A.M., P.K., and T.S.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Data Availability

Raw RNA-seq reads have been deposited in the NCBI Sequence Read Archive under BioProject accession PRJNA434149. A javascript app is also publically available for the exploration and visualisation of this dataset (https://neavemj.github.io/posts/D3_rab_app).

References

  1. Abrantes, J.; van der Loo, W.; Le Pendu, J.; Esteves, P.J. Rabbit haemorrhagic disease (RHD) and rabbit haemorrhagic disease virus (RHDV): A review. Vet. Res. 2012, 43, 12. [Google Scholar] [CrossRef] [PubMed]
  2. Le Gall-Recule, G.; Zwingelstein, F.; Boucher, S.; Le Normand, B.; Plassiart, G.; Portejoie, Y.; Decors, A.; Bertagnoli, S.; Guerin, J.L.; Marchandeau, S. Detection of a new variant of rabbit haemorrhagic disease virus in France. Vet. Rec. 2011, 168, 137–138. [Google Scholar] [CrossRef] [PubMed]
  3. Le Pendu, J.; Abrantes, J.; Bertagnoli, S.; Guitton, J.-S.; Le Gall-Reculé, G.; Lopes, A.M.; Marchandeau, S.; Alda, F.; Almeida, T.; Célio, A.P.; et al. Proposal for a unified classification system and nomenclature of lagoviruses. J. Gen. Virol. 2017, 98, 1658–1666. [Google Scholar] [CrossRef] [PubMed]
  4. Cooke, B.D.; Fenner, F. Rabbit haemorrhagic disease and the biological control of wild rabbits, Oryctolagus cuniculus, in Australia and New Zealand. Wildl. Res. 2002, 29, 689–706. [Google Scholar] [CrossRef]
  5. Dalton, K.P.; Nicieza, I.; Balseiro, A.; Muguerza, M.A.; Rosell, J.M.; Casais, R.; Alvarez, A.L.; Parra, F. Variant rabbit hemorrhagic disease virus in young rabbits, Spain. Emerg. Infect. Dis. 2012, 18, 2009–2012. [Google Scholar] [CrossRef] [PubMed]
  6. Abrantes, J.; Lopes, A.M.; Dalton, K.P.; Melo, P.; Correia, J.J.; Ramada, M.; Alves, P.C.; Parra, F.; Esteves, P.J. New variant of rabbit hemorrhagic disease virus, Portugal, 2012–2013. Emerg. Infect. Dis. 2013, 19, 1900–1902. [Google Scholar] [CrossRef] [PubMed]
  7. Westcott, D.G.; Frossard, J.P.; Everest, D.; Dastjerdi, A.; Duff, J.P.; Steinbach, F.; Choudhury, B. Incursion of RHDV2-like variant in Great Britain. Vet. Rec. 2014, 174, 333. [Google Scholar] [CrossRef] [PubMed]
  8. Baily, J.L.; Dagleish, M.P.; Graham, M.; Maley, M.; Rocchi, M.S. RHDV variant 2 presence detected in Scotland. Vet. Rec. 2014, 174, 411. [Google Scholar] [CrossRef] [PubMed]
  9. Duarte, M.; Carvalho, C.; Bernardo, S.; Barros, S.V.; Benevides, S.; Flor, L.; Monteiro, M.; Marques, I.; Henriques, M.; Barros, S.C.; et al. Rabbit haemorrhagic disease virus 2 (RHDV2) outbreak in Azores: Disclosure of common genetic markers and phylogenetic segregation within the European strains. Infect. Genet. Evol. 2015, 35, 163–171. [Google Scholar] [CrossRef] [PubMed]
  10. Martin-Alonso, A.; Martin-Carrillo, N.; Garcia-Livia, K.; Valladares, B.; Foronda, P. Emerging rabbit haemorrhagic disease virus 2 (RHDV2) at the gates of the African continent. Infect. Genet. Evol. 2016, 44, 46–50. [Google Scholar] [CrossRef] [PubMed]
  11. Rocchi, M.; Maley, M.; Dagleish, M.; Vick, C.; Ryan, D.; Lee, A.; Jahns, H. RHDV-2 on the Isle of Man and in the Republic of Ireland. Vet. Rec. 2016, 179, 389–390. [Google Scholar] [CrossRef] [PubMed]
  12. Neimanis, A.S.; Ahola, H.; Zohari, S.; Larsson Pettersson, U.; Brojer, C.; Capucci, L.; Gavier-Widen, D. Arrival of rabbit haemorrhagic disease virus 2 to northern Europe: Emergence and outbreaks in wild and domestic rabbits (Oryctolagus cuniculus) in Sweden. Transbound. Emerg. Dis. 2017, 65, 213–220. [Google Scholar] [CrossRef] [PubMed]
  13. Hall, R.N.; Mahar, J.E.; Haboury, S.; Stevens, V.; Holmes, E.C.; Strive, T. Emerging rabbit hemorrhagic disease virus 2 (RHDVb), Australia. Emerg. Infect. Dis. 2015, 21, 2276–2278. [Google Scholar] [CrossRef] [PubMed]
  14. Lopes, A.M.; Correia, J.; Abrantes, J.; Melo, P.; Ramada, M.; Magalhaes, M.J.; Alves, P.C.; Esteves, P.J. Is the new variant RHDV replacing genogroup 1 in Portuguese wild rabbit populations? Viruses 2014, 7, 27–36. [Google Scholar] [CrossRef] [PubMed]
  15. Mahar, J.E.; Hall, R.N.; Peacock, D.; Kovaliski, J.; Piper, M.; Mourant, R.; Huang, N.; Campbell, S.; Gu, X.; Read, A.; et al. Rabbit haemorrhagic disease virus 2 (GI.2) is replacing endemic strains of RHDV in the Australian landscape within 18 months of its arrival. J. Virol. 2017, 92. [Google Scholar] [CrossRef] [PubMed]
  16. Le Gall-Recule, G.; Lavazza, A.; Marchandeau, S.; Bertagnoli, S.; Zwingelstein, F.; Cavadini, P.; Martinelli, N.; Lombardi, G.; Guerin, J.L.; Lemaitre, E.; et al. Emergence of a new lagovirus related to rabbit haemorrhagic disease virus. Vet. Res. 2013, 44, 81. [Google Scholar] [CrossRef] [PubMed]
  17. Puggioni, G.; Cavadini, P.; Maestrale, C.; Scivoli, R.; Botti, G.; Ligios, C.; Le Gall-Recule, G.; Lavazza, A.; Capucci, L. The new French 2010 rabbit hemorrhagic disease virus causes an RHD-like disease in the Sardinian Cape hare (Lepus capensis mediterraneus). Vet. Res. 2013, 44, 96. [Google Scholar] [CrossRef] [PubMed]
  18. Camarda, A.; Pugliese, N.; Cavadini, P.; Circella, E.; Capucci, L.; Caroli, A.; Legretto, M.; Mallia, E.; Lavazza, A. Detection of the new emerging rabbit haemorrhagic disease type 2 virus (RHDV2) in Sicily from rabbit (Oryctolagus cuniculus) and Italian hare (Lepus corsicanus). Res. Vet. Sci. 2014, 97, 642–645. [Google Scholar] [CrossRef] [PubMed]
  19. Velarde, R.; Cavadini, P.; Neimanis, A.; Cabezon, O.; Chiari, M.; Gaffuri, A.; Lavin, S.; Grilli, G.; Gavier-Widen, D.; Lavazza, A.; et al. Spillover events of infection of brown hares (Lepus europaeus) with rabbit haemorrhagic disease type 2 virus (RHDV2) caused sporadic cases of an European brown hare syndrome-like disease in Italy and Spain. Transbound. Emerg. Dis. 2016, 64, 1750–1761. [Google Scholar] [CrossRef] [PubMed]
  20. Hall, R.N.; Peacock, D.E.; Kovaliski, J.; Mahar, J.E.; Mourant, R.; Piper, M.; Strive, T. Detection of RHDV2 in European brown hares (Lepus europaeus) in Australia. Vet. Rec. 2016, 180, 121. [Google Scholar] [CrossRef] [PubMed]
  21. Lopes, A.M.; Dalton, K.P.; Magalhaes, M.J.; Parra, F.; Esteves, P.J.; Holmes, E.C.; Abrantes, J. Full genomic analysis of new variant rabbit hemorrhagic disease virus revealed multiple recombination events. J. Gen. Virol. 2015, 96, 1309–1319. [Google Scholar] [CrossRef] [PubMed]
  22. Hall, R.N.; Mahar, J.E.; Read, A.J.; Mourant, R.; Piper, M.; Huang, N.; Strive, T. A strain-specific multiplex RT-PCR for Australian rabbit haemorrhagic disease viruses uncovers a new recombinant virus variant in rabbits and hares. Transbound. Emerg. Dis. 2017, 65, e444–e456. [Google Scholar] [CrossRef] [PubMed]
  23. Morisse, J.P.; Le Gall, G.; Boilletot, E. Hepatitis of viral origin in Leporidae: Introduction and aetiological hypotheses. Rev. Sci. Tech. 1991, 10, 269–310. [Google Scholar] [CrossRef] [PubMed]
  24. Cooke, B.D.; Berman, D. Effect of inoculation route and ambient temperature on the survival time of rabbits, Oryctolagus cuniculus, infected with Rabbit Haemorrhagic Disease Virus. Wildl. Res. 2000, 27, 137–142. [Google Scholar] [CrossRef]
  25. Capucci, L.; Cavadini, P.; Schiavitto, M.; Lombardi, G.; Lavazza, A. Increased pathogenicity in rabbit haemorrhagic disease virus type 2 (RHDV2). Vet. Rec. 2017, 180, 426. [Google Scholar] [CrossRef] [PubMed]
  26. Song, Y.; Wang, F.; Fan, Z.; Hu, B.; Liu, X.; Wei, H.; Xue, J.; Xu, W.; Qiu, R. Identification of novel rabbit hemorrhagic disease virus B-cell epitopes and their interaction with host histo-blood group antigens. J. Gen. Virol. 2016, 97, 356–365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Leuthold, M.M.; Dalton, K.P.; Hansman, G.S. Structural analysis of a rabbit hemorrhagic disease virus binding to histo-blood group antigens. J. Virol. 2015, 89, 2378–2387. [Google Scholar] [CrossRef] [PubMed]
  28. Nystrom, K.; Le Gall-Recule, G.; Grassi, P.; Abrantes, J.; Ruvoen-Clouet, N.; Le Moullac-Vaidye, B.; Lopes, A.M.; Esteves, P.J.; Strive, T.; Marchandeau, S.; et al. Histo-blood group antigens act as attachment factors of rabbit hemorrhagic disease virus infection in a virus strain-dependent manner. PLoS Pathog. 2011, 7, e1002188. [Google Scholar] [CrossRef] [PubMed]
  29. Guillon, P.; Ruvoen-Clouet, N.; Le Moullac-Vaidye, B.; Marchandeau, S.; Le Pendu, J. Association between expression of the H histo-blood group antigen, alpha1,2fucosyltransferases polymorphism of wild rabbits, and sensitivity to rabbit hemorrhagic disease virus. Glycobiology 2009, 19, 21–28. [Google Scholar] [CrossRef] [PubMed]
  30. Ruvoen-Clouet, N.; Ganiere, J.P.; Andre-Fontaine, G.; Blanchard, D.; Le Pendu, J. Binding of rabbit hemorrhagic disease virus to antigens of the ABH histo-blood group family. J. Virol. 2000, 74, 11950–11954. [Google Scholar] [CrossRef] [PubMed]
  31. Lopes, A.M.; Breiman, A.; Lora, M.; Le Moullac-Vaidye, B.; Galanina, O.; Nyström, K.; Marchandeau, S.; Le Gall-Reculé, G.; Strive, T.; Neimanis, A.; et al. Host specific glycans are correlated with susceptibility to infection by lagoviruses, but not with their virulence. J. Virol. 2017, 92. [Google Scholar] [CrossRef] [PubMed]
  32. Jung, J.Y.; Lee, B.J.; Tai, J.H.; Park, J.H.; Lee, Y.S. Apoptosis in rabbit haemorrhagic disease. J. Comp. Pathol. 2000, 123, 135–140. [Google Scholar] [CrossRef] [PubMed]
  33. Prieto, J.M.; Fernandez, F.; Alvarez, V.; Espi, A.; Garcia Marin, J.F.; Alvarez, M.; Martin, J.M.; Parra, F. Immunohistochemical localisation of rabbit haemorrhagic disease virus VP-60 antigen in early infection of young and adult rabbits. Res. Vet. Sci. 2000, 68, 181–187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Marcato, P.S.; Benazzi, C.; Vecchi, G.; Galeotti, M.; Della Salda, L.; Sarli, G.; Lucidi, P. Clinical and pathological features of viral haemorrhagic disease of rabbits and the European brown hare syndrome. Rev. Sci. Tech. 1991, 10, 371–392. [Google Scholar] [CrossRef] [PubMed]
  35. Alexandrov, M.; Peshev, R.; Yanchev, I.; Bozhkov, S.; Doumanova, L.; Dimitrov, T.; Zacharieva, S. Immunohistochemical localization of the rabbit haemorrhagic disease viral antigen. Arch. Virol. 1992, 127, 355–363. [Google Scholar] [CrossRef] [PubMed]
  36. Stoerckle-Berger, N.; Keller-Berger, B.; Ackermann, M.; Ehrensperger, F. Immunohistological diagnosis of rabbit haemorrhagic disease (RHD). Zbl. Veterinarmed. B 1992, 39, 237–245. [Google Scholar] [CrossRef]
  37. Park, J.H.; Itakura, C. Detection of rabbit haemorrhagic disease virus antigen in tissues by immunohistochemistry. Res. Vet. Sci. 1992, 52, 299–306. [Google Scholar] [CrossRef]
  38. Ramiro-Ibanez, F.; Martin-Alonso, J.M.; Garcia Palencia, P.; Parra, F.; Alonso, C. Macrophage tropism of rabbit hemorrhagic disease virus is associated with vascular pathology. Virus Res. 1999, 60, 21–28. [Google Scholar] [CrossRef]
  39. Gelmetti, D.; Grieco, V.; Rossi, C.; Capucci, L.; Lavazza, A. Detection of rabbit haemorrhagic disease virus (RHDV) by in situ hybridisation with a digoxigenin labelled RNA probe. J. Virol. Methods 1998, 72, 219–226. [Google Scholar] [CrossRef]
  40. Fuchs, A.; Weissenbock, H. Comparative histopathological study of rabbit haemorrhagic disease (RHD) and European brown hare syndrome (EBHS). J. Comp. Pathol. 1992, 107, 103–113. [Google Scholar] [CrossRef]
  41. Ferreira, P.G.; Costa-e-Silva, A.; Monteiro, E.; Oliveira, M.J.; Aguas, A.P. Transient decrease in blood heterophils and sustained liver damage caused by calicivirus infection of young rabbits that are naturally resistant to rabbit haemorrhagic disease. Res. Vet. Sci. 2004, 76, 83–94. [Google Scholar] [CrossRef] [PubMed]
  42. Ferreira, P.G.; Costa, E.S.A.; Monteiro, E.; Oliveira, M.J.; Aguas, A.P. Liver enzymes and ultrastructure in rabbit haemorrhagic disease (RHD). Vet. Res. Commun. 2006, 30, 393–401. [Google Scholar] [CrossRef] [PubMed]
  43. Tunon, M.J.; Sanchez-Campos, S.; Garcia-Ferreras, J.; Alvarez, M.; Jorquera, F.; Gonzalez-Gallego, J. Rabbit hemorrhagic viral disease: Characterization of a new animal model of fulminant liver failure. J. Lab. Clin. Med. 2003, 141, 272–278. [Google Scholar] [CrossRef] [PubMed]
  44. Ferreira, P.G.; Costa-e-Silva, A.; Oliveira, M.J.; Monteiro, E.; Cunha, E.M.; Aguas, A.P. Severe leukopenia and liver biochemistry changes in adult rabbits after calicivirus infection. Res. Vet. Sci. 2006, 80, 218–225. [Google Scholar] [CrossRef] [PubMed]
  45. Plassiart, G.; Guelfi, J.F.; Ganiere, J.P.; Wang, B.; Andre-Fontaine, G.; Wyers, M. Hematological parameters and visceral lesions relationships in rabbit viral hemorrhagic disease. Zbl. Veterinarmed. B 1992, 39, 443–453. [Google Scholar] [CrossRef]
  46. Ueda, K.; Park, J.H.; Ochiai, K.; Itakura, C. Disseminated intravascular coagulation (DIC) in rabbit haemorrhagic disease. Jpn. J. Vet. Res. 1992, 40, 133–141. [Google Scholar] [PubMed]
  47. Robinson, A.J.; So, P.T.M.; Müller, W.J.; Cooke, B.D.; Capucci, L. Statistical models for the effect of age and maternal antibodies on the development of rabbit haemorrhagic disease in Australian wild rabbits. Wildl. Res. 2002, 29, 663–671. [Google Scholar] [CrossRef]
  48. Ferreira, P.G.; Dinis, M.; Costa, E.S.A.; Aguas, A.P. Adult rabbits acquire resistance to lethal calicivirus infection by adoptive transfer of sera from infected young rabbits. Vet. Immunol. Immunopathol. 2008, 121, 364–369. [Google Scholar] [CrossRef] [PubMed]
  49. Matthaei, M.; Kerr, P.J.; Read, A.J.; Hick, P.; Haboury, S.; Wright, J.D.; Strive, T. Comparative quantitative monitoring of rabbit haemorrhagic disease viruses in rabbit kittens. Virol. J. 2014, 11, 109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Mikami, O.; Park, J.H.; Kimura, T.; Ochiai, K.; Itakura, C. Hepatic lesions in young rabbits experimentally infected with rabbit haemorrhagic disease virus. Res. Vet. Sci. 1999, 66, 237–242. [Google Scholar] [CrossRef] [PubMed]
  51. Marques, R.M.; Costa, E.S.A.; Aguas, A.P.; Teixeira, L.; Ferreira, P.G. Early inflammatory response of young rabbits attending natural resistance to calicivirus (RHDV) infection. Vet. Immunol. Immunopathol. 2012, 150, 181–188. [Google Scholar] [CrossRef] [PubMed]
  52. Marques, R.M.; Teixeira, L.; Aguas, A.P.; Ribeiro, J.C.; Costa-e-Silva, A.; Ferreira, P.G. Immunosuppression abrogates resistance of young rabbits to rabbit haemorrhagic disease (RHD). Vet. Res. 2014, 45, 14. [Google Scholar] [CrossRef] [PubMed]
  53. Ruvoen-Clouet, N.; Blanchard, D.; Andre-Fontaine, G.; Ganiere, J.P. Partial characterization of the human erythrocyte receptor for rabbit haemorrhagic disease virus. Res. Virol. 1995, 146, 33–41. [Google Scholar] [CrossRef]
  54. Liu, J.; Kerr, P.J.; Wright, J.D.; Strive, T. Serological assays to discriminate rabbit haemorrhagic disease virus from Australian non-pathogenic rabbit calicivirus. Vet. Microbiol. 2012, 157, 345–354. [Google Scholar] [CrossRef] [PubMed]
  55. Neimanis, A.; Pettersson, U.L.; Huang, N.; Gavier-Widen, D.; Strive, T. Elucidation of the pathology and tissue distribution of Lagovirus europaeus GI.2/RHDV2 (rabbit haemorrhagic disease virus 2) in young and adult rabbits (Oryctolagus cuniculus). Vet. Res. 2018, 49, 46. [Google Scholar] [CrossRef] [PubMed]
  56. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  57. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  59. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed]
  60. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene Ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [PubMed]
  61. Alexa, A.; Rahnenfuhrer, J. topGO: Enrichment Analysis for Gene Ontology. Bioconductor. Available online: http://www.bioconductor.org/packages//2.10/bioc/html/topGO.html (accessed on 1 October 2010).
  62. Kanehisa, M.; Furumichi, M.; Tanabe, M.; Sato, Y.; Morishima, K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017, 45, D353–D361. [Google Scholar] [CrossRef] [PubMed]
  63. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43. [Google Scholar] [CrossRef] [PubMed]
  64. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer Publishing Company: New York, NY, USA, 2009; p. 216. [Google Scholar]
  65. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Hahne, F.; Ivanek, R. Visualizing genomic data using Gviz and Bioconductor. Methods Mol. Biol. 2016, 1418, 335–351. [Google Scholar] [PubMed]
  67. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef] [PubMed]
  68. Meyers, G.; Wirblich, C.; Thiel, H.-J. Genomic and subgenomic RNAs of rabbit hemorrhagic disease virus are both protein-linked and packaged into particles. Virology 1991, 184, 677–686. [Google Scholar] [CrossRef]
  69. Van den Broeke, C.; Jacob, T.; Favoreel, H.W. Rho’ing in and out of cells: Viral interactions with Rho GTPase signaling. Small GTPases 2014, 5, e28318. [Google Scholar] [CrossRef] [PubMed]
  70. Nakayama, M.; Takeda, K.; Kawano, M.; Takai, T.; Ishii, N.; Ogasawara, K. Natural killer (NK)–dendritic cell interactions generate MHC class II-dressed NK cells that regulate CD4(+) T cells. Proc. Natl. Acad. Sci. USA 2011, 108, 18360–18365. [Google Scholar] [CrossRef] [PubMed]
  71. Roche, P.A.; Furuta, K. The ins and outs of MHC class II-mediated antigen processing and presentation. Nat. Rev. Immunol. 2015, 15, 203–216. [Google Scholar] [CrossRef] [PubMed]
  72. Syal, G.; Fausther, M.; Dranoff, J.A. Advances in cholangiocyte immunobiology. Am. J. Physiol. Gastrointest. Liver Physiol. 2012, 303, G1077–G1086. [Google Scholar] [CrossRef] [PubMed]
  73. Hattoum, A.; Rubin, E.; Orr, A.; Michalopoulos, G.K. Expression of hepatocyte epidermal growth factor receptor, FAS and glypican 3 in EpCAM-positive regenerative clusters of hepatocytes, cholangiocytes, and progenitor cells in human liver failure. Hum. Pathol. 2013, 44, 743–749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Neefjes, J.; Jongsma, M.L.; Paul, P.; Bakke, O. Towards a systems understanding of MHC class I and MHC class II antigen presentation. Nat. Rev. Immunol. 2011, 11, 823–836. [Google Scholar] [CrossRef] [PubMed]
  75. Cruickshank, S.M.; Southgate, J.; Selby, P.J.; Trejdosiewicz, L.K. Expression and cytokine regulation of immune recognition elements by normal human biliary epithelial and established liver cell lines in vitro. J. Hepatol. 1998, 29, 550–558. [Google Scholar] [CrossRef]
  76. Biron, C.A.; Byron, K.S.; Sullivan, J.L. Severe herpesvirus infections in an adolescent without natural killer cells. N. Engl. J. Med. 1989, 320, 1731–1735. [Google Scholar] [CrossRef] [PubMed]
  77. Shaw, R.K.; Issekutz, A.C.; Fraser, R.; Schmit, P.; Morash, B.; Monaco-Shawver, L.; Orange, J.S.; Fernandez, C.V. Bilateral adrenal EBV-associated smooth muscle tumors in a child with a natural killer cell deficiency. Blood 2012, 119, 4009–4012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Horton, N.C.; Mathew, P.A. NKp44 and natural cytotoxicity receptors as damage-associated molecular pattern recognition receptors. Front. Immunol. 2015, 6, 31. [Google Scholar] [CrossRef] [PubMed]
  79. Burn, G.L.; Svensson, L.; Sanchez-Blanco, C.; Saini, M.; Cope, A.P. Why is PTPN22 a good candidate susceptibility gene for autoimmune disease? FEBS Lett. 2011, 585, 3689–3698. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Begovich, A.B.; Carlton, V.E.H.; Honigberg, L.A.; Schrodi, S.J.; Chokkalingam, A.P.; Alexander, H.C.; Ardlie, K.G.; Huang, Q.; Smith, A.M.; Spoerke, J.M.; et al. A missense single-nucleotide polymorphism in a gene encoding a protein tyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis. Am. J. Hum. Genet. 2004, 75, 330–337. [Google Scholar] [CrossRef] [PubMed]
  81. Pisegna, S.; Zingoni, A.; Pirozzi, G.; Cinque, B.; Cifone, M.G.; Morrone, S.; Piccoli, M.; Frati, L.; Palmieri, G.; Santoni, A. Src-dependent Syk activation controls CD69-mediated signaling and function on human NK cells. J. Immunol. 2002, 169, 68–74. [Google Scholar] [CrossRef] [PubMed]
  82. Malorni, W.; Quaranta, M.G.; Straface, E.; Falzano, L.; Fabbri, A.; Viora, M.; Fiorentini, C. The Rac-activating toxin cytotoxic necrotizing factor 1 oversees NK cell-mediated activity by regulating the actin/microtubule interplay. J. Immunol. 2003, 171, 4195–4202. [Google Scholar] [CrossRef] [PubMed]
  83. Wynn, T.A.; Chawla, A.; Pollard, J.W. Macrophage biology in development, homeostasis and disease. Nature 2013, 496, 445. [Google Scholar] [CrossRef] [PubMed]
  84. Satpathy, A.T.; Wu, X.; Albring, J.C.; Murphy, K.M. Re(de)fining the dendritic cell lineage. Nat. Immunol. 2012, 13, 1145–1154. [Google Scholar] [CrossRef] [PubMed]
  85. Holness, C.L.; Simmons, D.L. Molecular cloning of CD68, a human macrophage marker related to lysosomal glycoproteins. Blood 1993, 81, 1607. [Google Scholar] [PubMed]
  86. Satpathy, A.T.; Kc, W.; Albring, J.C.; Edelson, B.T.; Kretzer, N.M.; Bhattacharya, D.; Murphy, T.L.; Murphy, K.M. Zbtb46 expression distinguishes classical dendritic cells and their committed progenitors from other immune lineages. J. Exp. Med. 2012, 209, 1135–1152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Armstrong, A.; Eck, S.L. EpCAM: A new therapeutic target for an old cancer antigen. Cancer Biol. Ther. 2003, 2, 320–325. [Google Scholar] [CrossRef] [PubMed]
  88. Haller, O.; Kochs, G. Interferon-induced Mx proteins: Dynamin-like GTPases with antiviral activity. Traffic 2002, 3, 710–717. [Google Scholar] [CrossRef] [PubMed]
  89. Marra, F.; DeFranco, R.; Grappone, C.; Milani, S.; Pastacaldi, S.; Pinzani, M.; Romanelli, R.G.; Laffi, G.; Gentilini, P. Increased expression of monocyte chemotactic protein-1 during active hepatic fibrogenesis: Correlation with monocyte infiltration. Am. J. Pathol. 1998, 152, 423–430. [Google Scholar] [PubMed]
  90. Otogawa, K.; Kinoshita, K.; Fujii, H.; Sakabe, M.; Shiga, R.; Nakatani, K.; Ikeda, K.; Nakajima, Y.; Ikura, Y.; Ueda, M.; et al. Erythrophagocytosis by liver macrophages (Kupffer cells) promotes oxidative stress, inflammation, and fibrosis in a rabbit model of steatohepatitis. Am. J. Pathol. 2007, 170, 967–980. [Google Scholar] [CrossRef] [PubMed]
  91. Huang, H.B. Vaccination against and immune response to viral haemorrhagic disease of rabbits: A review of research in the People’s Republic of China. Rev. Sci. Tech. 1991, 10, 481–498. [Google Scholar] [PubMed]
  92. Tunon, M.J.; San Miguel, B.; Crespo, I.; Riezu-Boj, J.I.; Larrea, E.; Alvarez, M.; Gonzalez, I.; Bustos, M.; Gonzalez-Gallego, J.; Prieto, J. Cardiotrophin-1 promotes a high survival rate in rabbits with lethal fulminant hepatitis of viral origin. J. Virol. 2011, 85, 13124–13132. [Google Scholar] [CrossRef] [PubMed]
  93. Dalton, K.P.; Nicieza, I.; Abrantes, J.; Esteves, P.J.; Parra, F. Spread of new variant RHDV in domestic rabbits on the Iberian Peninsula. Vet. Microbiol. 2014, 169, 67–73. [Google Scholar] [CrossRef] [PubMed]
  94. Alcami, A.; Koszinowski, U.H. Viral mechanisms of immune evasion. Immunol. Today 2000, 21, 447–455. [Google Scholar] [CrossRef]
  95. Karre, K.; Ljunggren, H.G.; Piontek, G.; Kiessling, R. Selective rejection of H-2-deficient lymphoma variants suggests alternative immune defence strategy. Nature 1986, 319, 675–678. [Google Scholar] [CrossRef] [PubMed]
  96. Gumá, M.; Angulo, A.; López-Botet, M. NK Cell receptors involved in the response to human cytomegalovirus infection. In Immunobiology of Natural Killer Cell Receptors; Compans, R.W., Cooper, M.D., Honjo, T., Koprowski, H., Melchers, F., Oldstone, M.B.A., Olsnes, S., Potter, M., Vogt, P.K., Wagner, H., et al., Eds.; Springer: Berlin/Heidelberg, Germany, 2006; pp. 207–223. [Google Scholar]
  97. McFadden, N.; Bailey, D.; Carrara, G.; Benson, A.; Chaudhry, Y.; Shortland, A.; Heeney, J.; Yarovinsky, F.; Simmonds, P.; Macdonald, A.; et al. Norovirus regulation of the innate immune response and apoptosis occurs via the product of the alternative open reading frame 4. PLoS Pathog. 2011, 7, e1002413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  98. Urakova, N.; Frese, M.; Hall, R.N.; Liu, J.; Matthaei, M.; Strive, T. Expression and partial characterisation of rabbit haemorrhagic disease virus non-structural proteins. Virology 2015, 484, 69–79. [Google Scholar] [CrossRef] [PubMed]
  99. Schwensow, N.; Mazzoni, C.J.; Marmesat, E.; Fickel, J.; Peacock, D.; Kovaliski, J.; Sinclair, R.; Cassey, P.; Cooke, B.; Sommer, S. High adaptive variability and virus-driven selection on major histocompatibility complex (MHC) genes in invasive wild rabbits in Australia. Biol. Invasions 2017, 19, 1255–1271. [Google Scholar] [CrossRef]
Figure 1. Similarity of gene expression profiles across rabbit age cohorts during infection with different GI viruses. GI.1-infected kittens had the smallest transcriptomic response after infection, while large changes were observed from 12 hpi in GI.2-infected kittens, and in GI.1-infected adults at 24 hpi. Notably, expression profiles varied considerably between adult and kitten control rabbits. Expression values from the top 500 most variable genes across the samples were normalized using regularized log transformation and plotted using principal component analysis.
Figure 1. Similarity of gene expression profiles across rabbit age cohorts during infection with different GI viruses. GI.1-infected kittens had the smallest transcriptomic response after infection, while large changes were observed from 12 hpi in GI.2-infected kittens, and in GI.1-infected adults at 24 hpi. Notably, expression profiles varied considerably between adult and kitten control rabbits. Expression values from the top 500 most variable genes across the samples were normalized using regularized log transformation and plotted using principal component analysis.
Viruses 10 00512 g001
Figure 2. Correlation between gene co-expression modules and the rabbit treatments (a) and functional enrichment of select modules (b). The pink co-expression module contained genes that tended to be downregulated in control adults, upregulated in control kittens, further upregulated in GI.1-infected kittens, but markedly downregulated in GI.2-infected kittens. Many of the other modules were significantly correlated with the large expression changes in GI.2-infected kittens. In (a), the correlation coefficient is given, and the p-value is in brackets. In (b), significantly enriched gene ontology terms are given with the p-value in brackets.
Figure 2. Correlation between gene co-expression modules and the rabbit treatments (a) and functional enrichment of select modules (b). The pink co-expression module contained genes that tended to be downregulated in control adults, upregulated in control kittens, further upregulated in GI.1-infected kittens, but markedly downregulated in GI.2-infected kittens. Many of the other modules were significantly correlated with the large expression changes in GI.2-infected kittens. In (a), the correlation coefficient is given, and the p-value is in brackets. In (b), significantly enriched gene ontology terms are given with the p-value in brackets.
Viruses 10 00512 g002
Figure 3. Regulation of genes in the major histocompatibility complex (MHC) during GI.1 and GI.2 infections. Several MHC genes tended to be expressed at a higher level in uninfected kittens compared to uninfected adults. Frequently, these were further upregulated in GI.1-infected kittens, but downregulated in GI.2-infected kittens (e.g., RLA-DRB1, SLA-DQD1, HLA-DQB2, CD74, HLA-DRA, HLA-DMA, and RLA-DMB). The header in each plot gives the Ensembl gene accession followed by the gene symbol.
Figure 3. Regulation of genes in the major histocompatibility complex (MHC) during GI.1 and GI.2 infections. Several MHC genes tended to be expressed at a higher level in uninfected kittens compared to uninfected adults. Frequently, these were further upregulated in GI.1-infected kittens, but downregulated in GI.2-infected kittens (e.g., RLA-DRB1, SLA-DQD1, HLA-DQB2, CD74, HLA-DRA, HLA-DMA, and RLA-DMB). The header in each plot gives the Ensembl gene accession followed by the gene symbol.
Viruses 10 00512 g003
Figure 4. Regulation of key markers for natural killer cell (PTPN22, VAV1, ARRB2, and NCR3), macrophage (CSF1R, ZBTB46, and CD68), and cholangiocyte development (EPCAM) during GI.1 and GI.2 infections. The header in each plot gives the Ensembl gene accession followed by the gene symbol.
Figure 4. Regulation of key markers for natural killer cell (PTPN22, VAV1, ARRB2, and NCR3), macrophage (CSF1R, ZBTB46, and CD68), and cholangiocyte development (EPCAM) during GI.1 and GI.2 infections. The header in each plot gives the Ensembl gene accession followed by the gene symbol.
Viruses 10 00512 g004
Figure 5. Regulation of antigen processing and presentation genes during GI.1 and GI.2 infections. GI.1-infected kittens upregulated pathways involved in natural killer cell and CD4 T-cell activation, such as MHC class II genes, while the same pathways were downregulated in GI.2-infected kittens. GI.1-infected adults also upregulated MHC class II genes, but to a lesser degree.
Figure 5. Regulation of antigen processing and presentation genes during GI.1 and GI.2 infections. GI.1-infected kittens upregulated pathways involved in natural killer cell and CD4 T-cell activation, such as MHC class II genes, while the same pathways were downregulated in GI.2-infected kittens. GI.1-infected adults also upregulated MHC class II genes, but to a lesser degree.
Viruses 10 00512 g005
Table 1. Experimental design and RNA-Seq mapping results.
Table 1. Experimental design and RNA-Seq mapping results.
IDVirushpiAgeTotal ReadsCleaned ReadsMapped to Rabbit Genome *Mapped to GI Genome ^RT-qPCR (Copies/mg Tissue)
A130contN/A **A36,608,81133,019,05730,372,424 (92.0%)0 (0.00%)55
A440cont12A31,544,53929,690,68827,851,164 (93.8%)0 (0.00%)0
A445cont12A33,024,74230,839,50428,557,666 (92.6%)0 (0.00%)0
K11cont12K30,376,75028,436,07426,697,303 (93.9%)0 (0.00%)0
K2cont12K26,466,33624,621,22123,121,670 (93.9%)0 (0.00%)0
K9cont12K38,879,76035,838,44933,677,252 (94.0%)0 (0.00%)0
A446cont24A26,875,05625,190,31423,565,537 (93.5%)0 (0.00%)0
A447cont24A30,801,40328,647,13026,648,792 (93.0%)0 (0.00%)0
K12cont24K21,522,68219,510,40418,264,437 (93.6%)0 (0.00%)0
K14cont24K32,273,37530,146,42228,183,027 (93.5%)0 (0.00%)0
K3cont24K34,270,68031,856,52229,771,558 (93.5%)0 (0.00%)0
A444GI.112A41,554,57438,521,84035,819,416 (93.0%)0 (0.00%)94
A451GI.112A34,269,02931,413,21729,237,190 (93.1%)0 (0.00%)28
A452GI.112A30,744,54428,741,20226,469,612 (92.1%)0 (0.00%)0
K1GI.112K35,648,93233,188,89531,053,518 (93.6%)0 (0.00%)0
K13GI.112K28,834,63026,852,29925,283,903 (94.2%)0 (0.00%)0
K5GI.112K32,247,35130,065,55928,089,931 (93.4%)0 (0.00%)0
K7GI.112K34,983,90532,474,26530,104,239 (92.7%)0 (0.00%)85
A443GI.124A32,025,31729,266,84327,437,173 (93.7%)4321 (0.01%)1,359,870
A448GI.124A41,956,00238,817,34935,919,297 (92.5%)817 (0.00%)213,792
A449GI.124A32,829,76629,594,92727,422,586 (92.7%)352 (0.00%)83,234
A450GI.124A39,743,58837,093,31034,584,907 (93.2%)5 (0.00%)1205
K10GI.124K26,707,07424,485,69822,779,842 (93.0%)2 (0.00%)307
K4GI.124K35,102,44032,387,79830,179,957 (93.2%)49 (0.00%)18,929
K6GI.124K30,314,74827,870,81625,865,930 (92.8%)19 (0.00%)6405
K8GI.124K30,187,24128,311,53926,206,586 (92.6%)5 (0.00%)1490
K379GI.212K32,179,57130,150,29728,485,814 (94.5%)0 (0.00%)81
K380GI.212K30,937,89429,113,19627,110,731 (93.1%)0 (0.00%)66
K381GI.212K24,208,54122,521,03020,774,527 (92.2%)0 (0.00%)0
K375GI.224K31,299,46529,544,13227,913,303 (94.5%)26,558 (0.09%)20,966,939
K376GI.224K30,626,98328,656,60226,667,507 (93.1%)1410 (0.00%)3,634,306
K377GI.224K35,624,01733,340,44731,470,468 (94.4%)1947 (0.00%)323,906
K378GI.224K26,510,09124,505,53422,803,925 (93.1%)22,528 (0.09%)4,857,896
* rabbit genome is OryCun2.0 (http://www.ensembl.org/Oryctolagus_cuniculus). ^ GI.1 is GenBank#KF594473, and GI.2 is GenBank #KT280060. ** A130 was an additional negative control animal that was not part of the infection experiment. A = adult, k = kitten, cont = noninfected control, hpi = hours post infection.
Table 2. Number of significantly differentially expressed genes (false discovery rate < 0.05).
Table 2. Number of significantly differentially expressed genes (false discovery rate < 0.05).
ComparisonUpregulatedDownregulated
Control adults vs. control kittens556 (up in adults)817 (up in kittens)
GI.1-infected kittens (12 hpi) * 21
GI.1-infected kittens (24 hpi) *602
GI.1-infected adults (12 hpi) *20
GI.1-infected adults (24 hpi) *29626
GI.2-infected kittens (12 hpi) *347293
GI.2-infected kittens (24 hpi) *21522059
* Compared to control animals.
Table 3. The top 10 most differentially expressed genes between control adults and control kittens.
Table 3. The top 10 most differentially expressed genes between control adults and control kittens.
Ensembl Gene #AnnotationLog Fold Change ^FDR *
ENSOCUG00000005528HAL (histidine ammonia-lyase)2.93.69 × 10−19
ENSOCUG00000014137NAALAD2 (N-acetylated alpha-linked acidic dipeptidase 2)−3.84.71 × 10−16
ENSOCUG00000002950HN1 (hematological and neurological expressed 1)−1.91.99 × 10−15
ENSOCUG00000003663ACTG2 (actin, gamma 2, smooth muscle, enteric)−2.37.09 × 10−13
ENSOCUG00000014725N/A6.98.92 × 10−13
ENSOCUG00000014317PIGR (polymeric immunoglobulin receptor)2.31.70 × 10−11
ENSOCUG00000023455N/A−3.01.70 × 10−11
ENSOCUG00000002858CDK1 (cyclin dependent kinase 1)−4.72.35 × 10−11
ENSOCUG00000009882NCAPH (condensing complex subunit 2)−2.92.35 × 10−11
ENSOCUG00000012902N/A−4.02.35 × 10−11
* FDR is the false discovery rate. ^ Positive fold change was upregulated in adults. N/A indicates unannotated gene.
Table 4. The top differentially expressed genes in GI.1-infected rabbits.
Table 4. The top differentially expressed genes in GI.1-infected rabbits.
Ensembl Gene #AnnotationLog Fold ChangeFDR *
Adults after 12 h
ENSOCUG00000023455N/A2.40.0003
ENSOCUG00000025530N/A2.30.0350
Adults after 24 h
ENSOCUG00000026233HRASLS2 (HRAS like suppressor 2)7.48.32 × 10−13
ENSOCUG00000029154N/A5.43.30 × 10−12
ENSOCUG00000002863IFIH1 (interferon induced with helicase C domain 1)3.41.31 × 10−11
ENSOCUG00000009504CD80 (T-lymphocyte activation antigen CD80)3.51.31 × 10−11
ENSOCUG00000020931SAMD9 (sterile alpha motif domain containing 9)4.65.12 × 10−11
ENSOCUG00000021037MX2 (MX dynamin like GTPase 2)5.41.81 × 10−10
ENSOCUG00000010311NXPE4 (neurexophilin and PC-esterase domain family member 4)3.73.05 × 10−10
ENSOCUG00000009811ZNFX1 (zinc finger NFX1-type containing 1)2.63.61 × 10−10
ENSOCUG00000000716EPSTI1 (epithelial stromal interaction 1)4.75.53 × 10−10
ENSOCUG00000016280CXCL10 (C-X-C motif chemokine ligand 10)6.85.53 × 10−10
Kittens after 12 h
ENSOCUG00000006223EPCAM (epithelial cell adhesion molecule)2.32.84 × 10−5
ENSOCUG00000009559KIAA18411.30.0472
ENSOCUG00000011988N/A−3.80.0472
Kittens after 24 h
ENSOCUG00000004501N/A4.14.62 × 10−6
ENSOCUG00000006595UBA7 (ubiquitin like modifier activating enzyme 7)2.91.32 × 10−5
ENSOCUG00000021924MX1 (MX dynamin like GTPase 1)4.21.84 × 10−5
ENSOCUG00000006482IFI44 (interferon induced protein 44)2.90.0001
ENSOCUG00000015823OAS2 (2’-5’-oligoadenylate synthetase 2)3.30.0001
ENSOCUG00000013278DHX58 (DExH-box helicase 58)2.50.0002
ENSOCUG00000027981ISG15 (ubiquitin-like modifier)5.40.0002
ENSOCUG00000015872HERC5 (HECT and RLD domain containing E3 ubiquitin protein ligase 5)2.60.0002
ENSOCUG00000024570IFIT5 (interferon induced protein with tetratricopeptide repeats 5)2.60.0002
ENSOCUG00000024734N/A2.70.0002
A maximum of ten differentially expressed genes is shown per experimental group. * FDR is the false discovery rate. N/A indicates unannotated genes.
Table 5. Significantly enriched Gene Ontology (GO) terms in the category “Biological Process” for upregulated genes in GI.1-infected rabbits.
Table 5. Significantly enriched Gene Ontology (GO) terms in the category “Biological Process” for upregulated genes in GI.1-infected rabbits.
GO #TermAnnotatedSignificantExpectedp-Value
Adults after 24 h
GO:0051607defense response to virus106291.722.40 × 10−19
GO:0045071negative regulation of viral genome replication27120.443.50 × 10−15
GO:0006955immune response6727010.925.70 × 10−13
GO:0071360cellular response to exogenous dsRNA1150.184.60 × 10−7
GO:0035455response to interferon-alpha1250.197.80 × 10−7
GO:0002474antigen processing and presentation of peptide antigen via MHC class 11250.197.80 × 10−7
GO:0060333interferon-gamma-mediated signaling pathway1250.197.80 × 10−7
GO:0034123positive regulation of toll-like receptor signaling pathway1350.211.30 × 10−6
GO:0050688regulation of defense response to virus3890.622.10 × 10−6
Kittens after 24 h
GO:0051607defense response to virus106170.374.60 × 10−17
GO:0045071negative regulation of viral genome replication2790.091.30 × 10−16
GO:0035455response to interferon-alpha1260.041.10 × 10−12
GO:0032727positive regulation of interferon-alpha production1140.044.00 × 10−8
GO:0032728positive regulation of interferon-beta production2240.088.60 × 10−7
GO:0009615response to virus155210.532.00 × 10−6
GO:0039529RIG-I signaling pathway1030.034.50 × 10−6
GO:0006955immune response672212.312.60 × 10−5
GO:0035456response to interferon-beta1430.056.60 × 10−5
GO:0019941modification-dependent protein catabolic process35441.220.00011

Share and Cite

MDPI and ACS Style

Neave, M.J.; Hall, R.N.; Huang, N.; McColl, K.A.; Kerr, P.; Hoehn, M.; Taylor, J.; Strive, T. Robust Innate Immunity of Young Rabbits Mediates Resistance to Rabbit Hemorrhagic Disease Caused by Lagovirus Europaeus GI.1 But Not GI.2. Viruses 2018, 10, 512. https://0-doi-org.brum.beds.ac.uk/10.3390/v10090512

AMA Style

Neave MJ, Hall RN, Huang N, McColl KA, Kerr P, Hoehn M, Taylor J, Strive T. Robust Innate Immunity of Young Rabbits Mediates Resistance to Rabbit Hemorrhagic Disease Caused by Lagovirus Europaeus GI.1 But Not GI.2. Viruses. 2018; 10(9):512. https://0-doi-org.brum.beds.ac.uk/10.3390/v10090512

Chicago/Turabian Style

Neave, Matthew J., Robyn N. Hall, Nina Huang, Kenneth A. McColl, Peter Kerr, Marion Hoehn, Jennifer Taylor, and Tanja Strive. 2018. "Robust Innate Immunity of Young Rabbits Mediates Resistance to Rabbit Hemorrhagic Disease Caused by Lagovirus Europaeus GI.1 But Not GI.2" Viruses 10, no. 9: 512. https://0-doi-org.brum.beds.ac.uk/10.3390/v10090512

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop