Next Article in Journal
Symmetry Breaking and Emergence of Directional Flows in Minimal Actomyosin Cortices
Next Article in Special Issue
Vasculopathy and Coagulopathy Associated with SARS-CoV-2 Infection
Previous Article in Journal
High Throughput strategies Aimed at Closing the GAP in Our Knowledge of Rho GTPase Signaling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Host Gene Expression of Macrophages in Response to Feline Coronavirus Infection

by
Yvonne Drechsler
1,*,†,
Elton J. R. Vasconcelos
2,*,†,
Lisa M. Griggs
1,
Pedro P. P. V. Diniz
1 and
Ellen Collisson
1
1
College of Veterinary Medicine, Western University of Health Sciences, Pomona, CA 91766, USA
2
Leeds Omics, University of Leeds, Leeds LS2 9JT, UK
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Cells 2020, 9(6), 1431; https://doi.org/10.3390/cells9061431
Submission received: 3 May 2020 / Revised: 5 June 2020 / Accepted: 6 June 2020 / Published: 9 June 2020
(This article belongs to the Special Issue The Cell Biology of Coronavirus Infection)

Abstract

:
Feline coronavirus is a highly contagious virus potentially resulting in feline infectious peritonitis (FIP), while the pathogenesis of FIP remains not well understood, particularly in the events leading to the disease. A predominant theory is that the pathogenic FIPV arises from a mutation, so that it could replicate not only in enterocytes of the intestines but also in monocytes, subsequently systemically transporting the virus. The immune status and genetics of affected cats certainly play an important role in the pathogenesis. Considering the importance of genetics and host immune responses in viral infections, the goal of this study was to elucidate host gene expression in macrophages using RNA sequencing. Macrophages from healthy male cats infected with FIPV 79-1146 ex vivo displayed a differential host gene expression. Despite the virus uptake, aligned viral reads did not increase from 2 to 17 h. The overlap of host gene expression among macrophages from different cats was limited, even though viral transcripts were detected in the cells. Interestingly, some of the downregulated genes in all macrophages were involved in immune signaling, while some upregulated genes common for all cats were found to be inhibiting immune activation. Our results highlight individual host responses playing an important role, consistent with the fact that few cats develop feline infectious peritonitis despite a common presence of enteric FCoV.

1. Introduction

Feline coronavirus (FCoV) is a highly contagious virus that is distributed worldwide and is ubiquitous in virtually all cat populations, particularly in multi-cat environments, such as shelters and rescues [1,2,3,4]. FCoV exists as two pathotypes, feline enteric coronavirus (FECV) and feline infectious peritonitis virus (FIPV). The enteric virus, FECV, commonly causes an asymptomatic infection presenting with mild gastrointestinal signs, and can persist in certain individuals [5]. FECV is readily transmitted via the fecal-oral route; therefore, the prevalence of FECV infection is generally associated with the number and density of cats housed together [1,4,6].
Sporadically and unpredictably, the infection can turn pathogenic, in which case the FIPV infection is associated with the highly fatal, systemic immune-mediated disease, feline infectious peritonitis (FIP) [5]. In this form, the virus develops a 100% fatal syndrome with two possible presentations: A pyogranulomatous disease of the visceral serosa and omentum, with consequent cavitary effusions, termed the “wet” or “effusive” form; and the dry form with granulomatous inflammation of parenchymatous organs such as kidneys, mesenteric lymph nodes, liver, pancreas, and central nervous system [7].
The pathogenesis of FIPV strains is still not fully understood, and several different theories are being discussed. A predominant theory has been that the pathogenic FIPV was caused by an FECV-like virus that had mutated such that it could acquire tropism for monocytes/macrophages. The implication was that FECV was present only in the intestines, while FIPV would now replicate mainly in monocytes, which subsequently transport the virus systemically, causing the disease. A study investigating strains FIPV 79-1146 and FECV 79-1183 comparing an ex vivo replication in cultured monocytes/macrophages of cats indicated a differential replication between the two viruses [8]. The authors speculated that a mutation potentially enables the FIPV to replicate more efficiently in monocytes and consequently, aids in spreading systemically in the cat. However, FECV has been shown to be present systemically in monocytes [9,10]. Pedersen et al. (2012) further showed that FIPV was present in feces of experimentally infected cats, but the virus was not infectious when tested via the oral-nasal route [11]. Pedersen inferred that intestinal and monocyte-macrophage tropisms of FECVs and FIPVs are related but not identical and that the mutational transformation of an FECV to FIPV might rather be happening in the monocytes/macrophages. These results show that analyzing responses of macrophages to the feline coronavirus is crucial for a better understanding of virus-host interactions.
A strong argument can also be made for the role of the immune status and genetics of affected cats in the pathogenesis. The latter theory is supported by the higher incidence of infection in very young or geriatric cats, purebreds, specifically Siamese and Burmese, and immune-compromised animals, such as those previously infected with the feline leukemia virus or feline immunodeficiency virus [4,12,13,14]. Several studies also have shown an important role for overall immune status and functioning cell-mediated immunity [15,16], antibody-mediated enhanced uptake of the virus [17,18], and numerous other factors such as destruction of lymphocytes in infected cats [19].
Another indication of the significance of the immune response in the pathogenesis is highlighted by the fact that inflammatory cytokine patterns play a role in the development of the dry or wet form of the disease [20,21,22]. However, little is known about the immunological pathways involved in early infection and the cells involved. The FIPV infection leads to a T-cell depletion in cats, although the virus is not directly infecting CD4+ or CD8+ T-cells [23,24]. The apoptosis of T-cells is probably caused by signaling mediators from other cells, such as infected macrophages or even epithelial cells [19,24].
Next-generation sequencing of RNA has increased our understanding of gene expression in response to pathogens in other viral infections and therefore, is a valuable tool to investigate early events in macrophage responses to the feline coronavirus. While more recent studies have investigated host responses to FIPV via transcriptome studies in Crandell-Rees feline kidney (CRFK) cells early after infection [25,26], peripheral blood mononuclear cells [26] or peritoneal cells including monocytes/macrophages [27] from infected cats, there are no reported transcriptome studies on primary monocytes/macrophages infected in the cell culture (ex vivo). In the current study, the transcriptome analysis of macrophages infected with FIPV 79-1146 was performed and analyzed for differential gene expression of ex vivo infected macrophages from six healthy male cats.

2. Materials and Methods

2.1. Cell Culture

CRFK cells (ATCC, Manassas, Virginia, cat#: CCL-94) were grown as monolayers in Dulbecco’s modified essential medium (DMEM) containing 10% of fetal bovine serum (FBS) and 1% of penicillin/streptomycin at 37 °C and 5% CO2.

2.2. Animal Procedures

All animal procedures were conducted and approved under the guidelines of the IACUC of Western University of Health Sciences, protocol approval number R10/IACUC/017. Peripheral blood for transcriptome studies was taken from six healthy male, specific pathogen-free (SPF) cats residing in an existing colony at the University of California, Davis. 30–40 mL of blood, equivalent to 1% of body weight or less was collected in heparinized tubes. The ages of cats at the time of blood draw were five months up to two years (five cats), and four years (one cat).

2.3. Monocyte Isolation

Monocytes from peripheral blood were isolated as previously described for canine monocytes [28] with some modifications. Briefly, the gradient centrifugation steps occurred at 450× g without break at deceleration and the subsequent washes to remove platelets were performed at 200× g for a total of three washes. PBMCs were counted, resuspended in RPMI 1640, containing 10% of FBS, 1% of penicillin and streptomycin, and 1× of non-essential amino acids, and plated at 5 × 106 in 6-well plates. Non-adherent cells were removed after 24 h by vigorously washing with a culture medium and cells were infected the following day.

2.4. Viral Infection for Transcriptome

For host transcriptome studies, macrophages were infected with FIPV 79-1146 (ATCC VR2128). The viruses were incubated at a multiplicity of infection (MOI) of 2 in a serum-free OptiMEM (Gibco, Thermo Fisher Scientific, Waltham, MA, USA) for 1 h for virus attachment, washed with OptiMEM, and incubated with fresh supplemented RPMI1640 for an additional 2 or 17 h. Technical replicates for the control, 2 and 17 h for macrophages from each cat were plated and incubated with PBS or the virus, respectively. CRFK cells (including technical replicates) were also infected as the control at an MOI of 1 in OptiMEM, followed by incubation in a supplemented DMEM. Uninfected controls underwent the same process with PBS without the virus. After incubation, the cell culture medium was completely removed and 600 μL of TRIzol (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) was added to each well, followed by RNA extraction with the ZymoResearch RNA kit (ZymoResearch, Irvine, CA) according to the manufacturer’s instructions. The RNA quality was evaluated via the Bioanalyzer (Agilent, Santa Clara, CA) and sent (1 μg per sample) for mRNA sequencing to Novogene, Inc. (Sacramento, CA, USA).

2.5. Quality Control of RNA Sequence Data

RNA paired-end sequencing quality control was assessed through FastQC (www.bioinformatics.babraham.ac.uk/projects/fastqc). An average of 35 million paired reads were sequenced per sample. Both adapters and low-quality bases (QV < 20) were trimmed from the reads’ extremities with Trimmomatic [29].

2.6. Alignment Against Reference Transcriptomes

Kallisto [30] was the algorithm of choice for performing the alignment of all paired reads against the whole Felis catus reference transcriptome (F. catus NCBI-RefSeq-9.0). An average of 87.5% of the total reads from each sample was mapped onto the cat’s annotated transcriptome. Alternatively, we also attempted to retrieve viral reads for both macrophages and CRFK from the sequenced libraries using Kallisto to align reads against the 11 protein-coding genes from the feline coronavirus (FCoV). An average of 0.03% of the total reads per sample was aligned against the FCoV annotated transcripts. Viral read counts were normalized as fragments per kilobase per million (FPKM): [read_counts / (gene_length_in_kb × total_reads_in_sample)] × 1,000,000.

2.7. Differential Expression

Tables generated by Kallisto were used as input for differential expression (DE) analyses. Due to the unique host responses in the macrophage dataset, the NOISeq version 2.14.1 [31] (Ctrl, 2, and 17 h) was employed to assess differentially expressed genes (DEGs) from each cat from which the ex vivo infected macrophages were derived. NOISeq output tables contained DEGs for each comparison (2 h versus control and 17 h vs. control) per cat. vennCounts and vennDiagram functions from the limma R package [32] were used for combining DEGs from each cat. Datasets were submitted to a multidimensional scaling (MDS) analysis, with the plotMDS function from EdgeR, to identify distinct samples clustered in a two dimensions-reduction landscape prior to the start of DE analyses. All tools described above for differential expression were run within the R environment version 3.5.2.

2.8. Gene Enrichment Analyses Using Both GO Terms and KEGG Pathways

After generating a list of differentially expressed genes (DEGs), we used ClueGO [33] under the Cytoscape version 3.7.1 [34,35] for a gene enrichment analysis relying on the Felis catus annotation from both gene ontology (http://geneontology.org/) and KEGG pathways (https://www.genome.jp/kegg/pathway.html) consortia. Both enrichment analyses adopted the Hypergeometric test along with the Benjamini and Hochberg p-value adjustment method. A 0.05 threshold was set for the latter.

3. Results

3.1. Host Responses Are Unique to Individual Cats with Several Clusters Present

To analyze host responses to FIPV, mRNA from feline macrophages (n = 6), infected ex vivo, was isolated and processed for next-generation sequencing. A unique bioinformatic analysis of differentially expressed genes (NOISeq) was necessary to compare gene expression. Individual cats displayed very different patterns of gene activation, with macrophages from some cats showing a clear separation of uninfected control versus infected samples (#1 and 2), while macrophages from other cats (#5 and 6) did cluster closely together independent of infection status or infection time (Figure 1).

3.2. Analysis of the Presence of Viral RNA in Infected Macrophages and CRFK

Cultured macrophages were infected with the FIPV 79-1146 for 2 and 17 h before RNA was collected for the expression analysis of host genes and viral reads. The viral presence was very low in the samples as indicated by the aligned reads for most macrophage samples at 2 h. However, the presence in most samples did indicate that viral particles were taken up by the cells (Table 1). No viral reads were detected in macrophages from cat #3.
Comparing 2 to 17 h samples, there was no indication of viral amplification despite viral uptake, as no increase of viral reads was observed in macrophages at the later time point. The viral reads at 17 h were either similar or lower than the number of reads at 2 h. Macrophages from one cat (cat #3) did not show any viral presence at 2 h. In contrast to the macrophages exposed to the virus, RNA sequencing showed several log-fold increases of viral isolates in the CRFK cells used as replication controls from 2 to 17 h (Table 1). This indicated that very few viral particles entering CRFK cells result in a significant amplification of the virus. This is also associated with pronounced cytopathic effects, which was observed in CRFK, but not macrophages.

3.3. Host Transcriptome Analysis Shows Cat Specific Responses to FIPV

Macrophages from each of the individual cats did differentially express several hundred genes at 2 and 17 h, compared to uninfected macrophages, highlighting the unique host responses of the cells from individual cats. Cat #3 was excluded from the bioinformatic analysis due to the lack of viral RNA in the macrophages. Since no control sample was available for cat #4, no NOISeq individual analysis of 2 or 17 h samples of cat #4 was possible in comparison to its own uninfected control. At 2 h after infection, macrophages from cats #1, 2, 5, and 6 expressed 1787, 499, 455, and 608 differentially downregulated genes, while 1620, 673, 416, and 639 genes were differentially upregulated (Figure 2).
It is of interest that, while the macrophages clearly show a differential regulation of genes, there are relatively few common DEGs expressed in all of the macrophages from the four cats that were compared. Only 10 genes were significantly upregulated in macrophages from all cats, while one gene was downregulated (Figure 2). The gene significantly downregulated in all macrophages was ATPase phospholipid transporting 8B4 (ATP8B4), which is involved in the cation transport and biosynthesis of ATP. Among the 10 upregulated genes, several are involved in immune signaling of the viral infection, such as SMAD4, ATP binding cassette, HAUS8, FCH domain, ubiquitin peptidase 20, and heteronuclear Ribonuclear Protein C (hnRNP C). The remaining genes in this group are involved in biosynthesis and other non-immune cellular functions (Table 2).
Monocyte isolations yielded different amounts of cells for each cat resulting in sufficient macrophages from only four cats to conduct the 17 h infection time point, with cat #3 not showing viral reads and thus not included in the analysis. Therefore, macrophage gene expression from cats #1, 2, and 5 was compared. At 17 h, macrophages from cats #1, 2, and 5 differentially downregulated 573, 1203, and 1452 genes, respectively, with 120 of the same genes downregulated in macrophages from all three cats (Figure 2). The number for upregulation of genes was similar, with 582, 1407, and 1611 genes for macrophages from cats #1, 2, and 5, respectively. One hundred and thirty-two genes were differentially upregulated in macrophages from all three cats (Figure 2). Downregulated immune genes at 17 h tumor necrosis factor (TNF) receptor superfamily 4, interleukin 10, transforming growth factor (TGF)-beta, signal transducer and transcription activator 3 (STAT3), and transcription factors interferon regulatory factor (IRF) 4, and SMAD family member 7 (Table A1). One hundred and thirty-four upregulated genes common for all cats included several other proteins involved in cell cycle and metabolism. Among these were several in the centromere or centrosomal category, nuclear body protein SP140, or 7-dehydrocholesterol reductase. There were only two classical cytokine genes among the upregulated, namely leukemia inhibitory factor (LIF), from the interleukin 6 cytokine family, and the transcription factor IRF 3 (Table A2).

3.4. Gene Enrichment for Pathway Analysis and Ontology

As indicated in Figure 2, at 2 h after infection, there was no significant overlap of DEGs common to all infected macrophages in order to allow gene enrichment analyses (both GO terms and KEGG pathways). There was more commonality of DEGs expressed in response to FIPV at 17 h when macrophage responses of individual cats were compared, therefore gene enrichment was performed. Overall, 120 genes were differentially downregulated in all cats in response to FIPV with one KEGG pathway enriched, “Arginine and proline metabolism” (Table 3).
Gene ontology for enriched genes downregulated at 17 h showed the terms “extracellular matrix binding” and “collagen binding” for the molecular function, and “regulation of endothelial cell proliferation”, “regulation of epithelial cell proliferation”, “positive regulation of vasculature development”, “positive regulation of angiogenesis”, and “regulation of peptidase activity” for biological processes (Table 4).
KEGG pathways for enriched genes upregulated at 17 h were “Steroid biosynthesis”, “Valine, leucine and isoleucine degradation”, “Butanoate metabolism”, “Terpenoid backbone biosynthesis”, “Cell cycle”, “p53 signaling pathway”, and “Progesterone-mediated oocyte maturation” (Table 3). Gene ontology for this same set of DEGs resulted in terms related to several nuclear and mitotic cellular processes, including cytoskeletal and spindle organization (Table 4).
The Cytoscape network analysis (by clueGO) shows how several of these pathways interact with each other, yielding two networks (Figure 3).

4. Discussion

While isolates of both FCoV serotypes I and II have been shown to be pathogenic, only serotype II has been shown to efficiently replicate in the cell culture, using CRFK or felis catus whole fetus (FCWF) cell lines. Therefore, serotype II strains have been the subject of expanded in vitro and ex vivo investigations, including those focusing on viral entry and replication in monocytes/macrophages. It remains to be satisfactorily explained how both of these serotypes can arise from FECV in individual cats, with a yet unidentified mutation potentially resulting in a differential replication in monocytes and subsequently, leading to de novo pathogenesis in the affected cats. It is, however, likely that host immune responses play a significant, but not sufficiently elucidated, role in the pathogenesis, and the question remains: What are the responses by the macrophages when infected with FCoV?
The presence of the virus has been evaluated directly in studies in macrophages [8,36,37] using PCR or immunofluorescence. Most studies showing an increased replication of FIPV in vitro, particularly in comparison to FECV, use indirect methods such as TCID50 in CRFK cells incubated with virus-infected macrophage extracts or supernatants [8,37]. However, to our knowledge, no studies so far have verified replication in macrophages by more state-of-the-art sensitive techniques, such as RNA sequencing.
In this study, RNA sequencing and analysis of viral reads indicated that at 2 h, the virus is taken up by the macrophages, as evidenced by the presence of viral RNA. However, there was no significant increase in viral RNA at 17 h, a time close to peak replication of the virus, which is usually at 24 h. Instead of a significant increase in viral reads, which would be expected with viral amplification, a similar amount, or even decrease of the virus was observed in the samples. In contrast to this, just a limited uptake of the virus into CRFK cells at 2 h led to several log-fold of replication, indicated by several thousands to hundreds of thousands of viral transcripts at 17 h. While the samples were not enriched for the pathogen, but rather the host RNA, there is a striking difference between replication of the virus in the macrophages compared to CRFK cells. Since subtle differences of replication in macrophages might not be recognized and could be below the threshold of detection, further investigations are warranted.
Macrophages of one cat #3 did not show the presence of viral RNA at all. This might be a technical issue, such as a higher MOI needed for successful infection. It is also possible that this particular host was more resistant to infection, which has been shown before both in vivo and in vitro [38,39]. Macrophages are competent immune cells and it has been shown that in vitro or ex vivo infection of macrophages with other viruses is challenging even if viremia exists in vivo. For example, there is evidence that Marek’s disease virus, a herpesvirus, is phagocytized by macrophages in vivo and then disseminated to infect T and B lymphocytes. In contrast, in vitro/ex vivo infection of macrophages has been shown to be difficult [40]. The investigators were successful when virally infected target cells were incubated with the macrophages, which then phagocytized the infected cells and therefore, took up the virus much more effectively. Similarly, it might be worth exploring the incubation of feline macrophages with infected epithelial target cells for studies requiring a higher ratio of macrophages positive for the viral antigen. In addition, methods such as RNA sequencing could be explored to better quantify differences between FIPV and FECV uptake, viral presence, and amplification in macrophages.
Gene expression of monocytes/macrophages following exposure to the feline coronavirus so far has not been investigated with next-generation sequencing. The PCR analysis of mRNA expression of individual cytokines and other immune-related proteins has been done in vitro and in vivo [19,20,21,41], but next-generation sequencing is poised to give a more complete picture of gene expression after infection. In particular, differences in responses to FIPV and FECV can be better defined.
Transcriptome analyses of FIPV-infected feline macrophages were done in this study at 2 h (early phase of infection) and 17 h (closer to the known peak of amplification of FCoV), and samples were positive for a limited presence of viral transcripts. However, no enrichment of pathways or gene ontology common for all cats was possible in the early phase after infection. The only gene downregulated in all analyzed samples at 2 h was ATPase phospholipid transporting 8B4 (ATP8B4). This protein is involved in phospholipid transport, but there are no studies to its involvement in viral infection and replication. From the 10 upregulated genes, several are involved in immune-related host responses (Table A1). SMAD4 is a transcription factor involved in immune signaling, particularly in TGF-beta signaling. TGF-beta is an anti-inflammatory or regulatory cytokine that has been demonstrated to act as a proviral factor in epithelial cells during an influenza infection [42]. Haus8 is involved in maintaining cellular spindle integrity, but appears to have a role in the RIG-1 like antiviral signaling pathway in a Sendai virus infection [43]. ATP binding cassette subfamily A member 1, alternatively named RNase L inhibitor, blocks ribonuclease L. The interferon-regulated 2-5A/RNase L pathway plays a major role in the antiviral innate immune response. Consequently, several viruses are known to inhibit this pathway, including HIV which also induces the RNase L inhibitor [44]. The FCH domain is a phosphoprotein, associated with a viral infection. It is linked to endocytosis in an avian influenza virus replication [45]. Heterogeneous nuclear ribonucleoprotein C associated with pre-mRNA processing, RNA metabolism, and transport is a host factor important in the replication of positive-strand RNA viruses [46,47]. Finally, ubiquitin-specific peptidase 20 is an enzyme that has been shown to negatively regulate NFkB in an HTLV infection [48]. Other genes in this group are involved in cellular processes not linked to viral replication. Taken together, we identified several genes that are downregulated in all macrophages infected with FIPV and are linked to antiviral signaling. However, further studies are needed to elucidate the role of these genes in viral host interactions.
Interestingly, several downregulated genes that did overlap among cats were involved in immune signaling, which might confirm that the immune responses of macrophages are negatively altered by FIPV. While a few of these genes are easily identifiable for involvement in antiviral or inflammatory responses, such as interferon regulatory factor 4 (IRF 4) or TGF beta, a further investigation shows that several of the transcripts do, in fact, play a role in responses to infection. For example, heterogeneous nuclear ribonucleoprotein A2/B1 (hNRP A2/B1), which is involved in RNA binding and trafficking, has previously been shown to bind to the NP of Avian Influenza [46]. Other proteins downregulated are purinergic receptors (P2Y), that are involved in antiviral responses, affecting cytokine responses and T-cell activation [49]. These receptors also contribute to the direct elimination of the virus by inhibiting their intracellular replication [50]. Intersectin 2 plays important role in the regulation of the adaptive immune response in viral infection [51]. Ubiquitin-like proteins have been shown to modify proteins, thus conferring functions related to programmed cell death, autophagy and regulation of the immune system [52], and IRF 4, which is critical to T-cell effector function [53].
In regards to the viral infection of macrophages, it would be expected to see a strong inflammatory response of the cell in response to the virus, such as upregulation of Toll-like receptor pathways, interferon type I signaling, etc. In our study, the absence of typical inflammatory signals is notable. While immune molecules IRF3 and LIF were upregulated, no classical virally activated signaling is discernable or enriched in KEGG pathways. LIF has not been well investigated in regards to viral infections, but it has been implicated in suppressing the replication of HIV [54]. IRF3 is necessary for IFN beta induction and SARS-CoV has been shown to block a step between the nuclear transport of IRF3 and its phosphorylation, which is necessary for an IFN induction, but not with the induction of IRF3 itself [55]. Our data is consistent with FCoV employing a similar strategy, however, further studies are needed.
Some differentially upregulated genes common for all cats, after a macrophage infection, were found to be within unexpected cellular functions, such as centromere proteins and TPX2 that are involved in specific phases of the cell cycle. Recent studies have indicated that centrosomes, as well as spindle organization, are actually involved in responses to a viral infection [56,57,58]. On the other hand, some upregulated transcripts included the protein tyrosine phosphatase, which is involved in attenuating T-cell activation [59] or nuclear body protein SP140. The latter has been shown to act as important repressor of genes involved in the regulation of cytokine production, inflammatory response, and cell-cell adhesion [60].
The picture emerging from these differentially expressed genes is that genes involved in antiviral responses and immune activation are depressed by FIPV uptake, while other genes involved in cell cycle and proteins repressing immune responses are upregulated. Even without an active viral replication, these changes most likely influence the pathogenesis and might explain how the monocytes carrying the virus potentially affect other cells, specifically lymphocytes.

5. Conclusions

The FIPV exposure and uptake leads to a limited differential gene expression in feline macrophages that might affect cell function. It will be of importance to further investigate cellular responses to FCoV isolates, in order to understand virus interactions with the host macrophages. A comparison between infection with FIPV and FECV strains focusing on host responses might yield insights into the pathogenesis of the virus. On the other hand, individual cat responses may be found to be of significant relevance to the pathogenesis and the number of samples required for analyses might be very high. The problem remains that very few cells take up the virus, and thus it will be difficult to analyze the transcriptome of a low number of infected cells. However, with high sample numbers, it could still be possible to identify differential gene expression that can provide valuable insight into the pathogenesis. An alternative is to use bone marrow-derived macrophages for viral infection since these macrophages have been shown to support replication at a higher rate, especially at a high MOI. An attempt to increase the number of virus-positive cells by incubating them with infected target cells may also identify more infected macrophages. Another possible source of infected macrophages is to isolate and analyze cells directly from viremic cats for transcriptome studies, however, immunologically those cells are representative of an ongoing disease process, and gene expression probably would be very different than early in the pathogenesis. In any case, further analyzing the host responses with next-generation molecular techniques will increase our understanding of FIPV pathogenesis.

Author Contributions

Conceptualization, Y.D., E.C., and P.P.P.V.D.; methodology, Y.D. and L.M.G.; software, E.J.R.V.; validation, Y.D. and L.G.; formal analysis, E.J.R.V. and Y.D.; investigation, L.M.G. and Y.D.; data curation, E.J.R.V.; writing—original draft preparation, Y.D.; writing—review and editing, Y.D., E.J.R.V., E.C., L.M.G., and P.P.P.V.D.; visualization, E.J.R.V. and Y.D.; supervision, Y.D.; project administration, Y.D.; funding acquisition, Y.D. and P.P.P.V.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by funding from the Morris Animal Foundation (D10FE-303) and Winn Feline Foundation (W13-019).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Genes differentially downregulated in all macrophages 17 h after infection with FIPV.
Table A1. Genes differentially downregulated in all macrophages 17 h after infection with FIPV.
Gene IdentifierGene Name
NM_001009200.1TNF receptor superfamily member 4
NM_001009209.1interleukin 10
NM_001009253.1CD80 molecule
XM_003980782.4transforming growth factor beta induced
XM_003981485.5CD74 molecule, transcript variant X1
XM_003981616.5growth arrest and DNA damage inducible beta
XM_003984354.4limb bud and heart development, transcript variant X1
XM_003986416.5PR/SET domain 1, transcript variant X1
XM_003987780.4serine and arginine rich splicing factor 5, transcript variant X2
XM_003987833.5NPC intracellular cholesterol transporter 2
XM_003990035.4peroxiredoxin 1
XM_003990688.5glycophorin C (Gerbich blood group)
XM_003994088.5ribosomal protein S24, transcript variant X6
XM_003995273.5carnosine dipeptidase 2, transcript variant X1
XM_003996781.5post-GPI attachment to proteins 3, transcript variant X1
XM_003998055.5metallothionein-1
XM_003999680.5T-cell surface glycoprotein CD1b
XM_006928964.4bromodomain and PHD finger containing 1, transcript variant X9
XM_006929396.3transmembrane protein 140, transcript variant X2
XM_006929637.3ribonucleic acid export 1, transcript variant X1
XM_006929792.2AAR2 splicing factor homolog, transcript variant X1
XM_006930307.4BCL2-like 11, transcript variant X1
XM_006930765.3cathepsin B, transcript variant X2
XM_006931414.4nucleoporin 153, transcript variant X3
XM_006931511.2interferon regulatory factor 4, transcript variant X3
XM_006932723.3mitochondrial ribosomal protein L52, transcript variant X5
XM_006932762.3ring finger protein 31, transcript variant X2
XM_006933673.2microspherule protein 1, transcript variant X4
XM_006934848.4Rab geranylgeranyltransferase beta subunit, transcript variant X2
XM_006934929.4coagulation factor III, tissue factor
XM_006935700.4diacylglycerol kinase delta, transcript variant X2
XM_006938342.3ankyrin repeat and LEM domain containing 2, transcript variant X2
XM_006938804.4elastin microfibril interfacer 2, transcript variant X1
XM_006940068.4zinc finger protein 207, transcript variant X5
XM_006941006.4CD37 molecule, transcript variant X5
XM_006943029.3Rho GTPase activating protein 30, transcript variant X2
XM_011282426.3ELOVL fatty acid elongase 5, transcript variant X3
XM_011282594.3BCL2 associated transcription factor 1, transcript variant X3
XM_011283077.3thrombospondin 1
XM_011283790.3C-type lectin domain family 4 member A, transcript variant X2
XM_011285121.3extracellular matrix protein 1, transcript variant X4
XM_011285878.3ATPase 13A3, transcript variant X1
XM_011286809.3CD44 molecule (Indian blood group), transcript variant X13
XM_011290903.3coagulation factor V
XM_011291023.2quinone oxidoreductase-like protein 2, transcript variant X1
XM_011291505.3poly(U) binding splicing factor 60, transcript variant X5
XM_011291709.2RNA binding motif protein 3, transcript variant X1
XM_019811387.2hypoxia upregulated 1, transcript variant X5
XM_019811602.2family with sequence similarity 76 member B, transcript variant X1
XM_019812261.2CD44 molecule (Indian blood group), transcript variant X2
XM_019812263.2CD44 molecule (Indian blood group), transcript variant X3
XM_019812267.2CD44 molecule (Indian blood group), transcript variant X7
XM_019812269.2CD44 molecule (Indian blood group), transcript variant X10
XM_019812271.2CD44 molecule (Indian blood group), transcript variant X12
XM_019815319.2cold inducible RNA binding protein, transcript variant X2
XM_019816708.2TSC complex subunit 1, transcript variant X3
XM_019817540.2HEAT repeat containing 6, transcript variant X4
XM_019817905.2signal transducer and activator of transcription 3, transcript variant X2
XM_019819573.1heterogeneous nuclear ribonucleoprotein L, transcript variant X3
XM_019819600.2ubiquitin-like modifier activating enzyme 2, transcript variant X2
XM_019820314.1Sad1 and UNC84 domain containing 1, transcript variant X5
XM_019821314.2signaling lymphocytic activation molecule family member 1, transcript variant X2
XM_019821551.2cellular repressor of E1A stimulated genes 1, transcript variant X2
XM_019822490.2ubiquitin associated protein 2-like, transcript variant X12
XM_019822721.2YTH N6-methyladenosine RNA binding protein 3, transcript variant X5
XM_019823162.2Src-like adaptor, transcript variant X1
XM_019823774.1RNA binding motif protein 3, transcript variant X2
XM_019823775.1RNA binding motif protein 3, transcript variant X4
XM_019824901.1ubiquitin specific peptidase-like 1, transcript variant X5
XM_019825689.2septin 7, transcript variant X7
XM_019829713.2tubulin gamma complex associated protein 3, transcript variant X3
XM_019832464.2EP300 interacting inhibitor of differentiation 1
XM_019836585.2myocyte enhancer factor 2C, transcript variant X18
XM_019836980.2Fas associated factor 1, transcript variant X1
XM_019837265.1zinc finger ZZ-type containing 3, transcript variant X2
XM_023238977.1pleckstrin homology-like domain family A member 2
XM_023239126.1cell adhesion molecule 1, transcript variant X3
XM_023240047.1calpain 1, transcript variant X1
XM_023241058.1ornithine aminotransferase, transcript variant X2
XM_023241314.1dynein axonemal heavy chain 10
XM_023241586.1neurofibromin 2, transcript variant X2
XM_023241792.1elastin microfibril interfacer 2, transcript variant X2
XM_023241843.1SS18, nBAF chromatin remodeling complex subunit, transcript variant X3
XM_023241980.1SMAD family member 7, transcript variant X2
XM_023242133.1nuclear factor of activated T-cells 1, transcript variant X6
XM_023242558.1thioredoxin
XM_023242850.1cold inducible RNA binding protein, transcript variant X3
XM_023243143.1transcription factor 3, transcript variant X8
XM_023243448.1lysine demethylase 6B, transcript variant X7
XM_023244057.1methyltransferase-like 23, transcript variant X2
XM_023244808.1TATA-box binding protein associated factor 15, transcript variant X3
XM_023245178.1POU class 2 homeobox 2, transcript variant X8
XM_023246400.1mitogen-activated protein kinase 8 interacting protein 3, transcript variant X16
XM_023247978.1lamin A/C
XM_023248580.1family with sequence similarity 49 member B, transcript variant X6
XM_023249097.1P2Y receptor family member 8, transcript variant X2
XM_023249251.1dedicator of cytokinesis 11, transcript variant X1
XM_023251817.1thymosin beta 10, transcript variant X2
XM_023252054.1intersectin 2, transcript variant X7
XM_023252112.1ornithine decarboxylase 1, transcript variant X2
XM_023253670.1CD83 molecule, transcript variant X1
XM_023253671.1CD83 molecule, transcript variant X2
XM_023253685.1human immunodeficiency virus type I enhancer binding protein 1
XM_023254138.1CD109 molecule, transcript variant X1
XM_023254166.1synaptotagmin binding cytoplasmic RNA interacting protein, transcript variant X4
XM_023255110.1transducing-like enhancer of split 3, transcript variant X7
XM_023255350.1mitogen-activated protein kinase 6, transcript variant X2
XM_023255486.1regulator of microtubule dynamics 3, transcript variant X3
XM_023255629.1HECT domain E3 ubiquitin protein ligase 1, transcript variant X5
XM_023256852.1cancer susceptibility 1, transcript variant X7
XM_023256920.1dynamin-1-like, transcript variant X6
XM_023257772.1adaptor-related protein complex 3 beta 1 subunit, transcript variant X1
XM_023258773.1ubiquitin specific peptidase 1, transcript variant X3
XM_023258859.1G protein subunit gamma 5
XM_023260612.1nuclear receptor subfamily 1 group D member 2
XM_023260684.1trafficking kinesin protein 1, transcript variant X8
XR_002736856.1tetraspanin 14, transcript variant X2
XR_002737850.1uncharacterized LOC111557626, transcript variant X1
XR_002740181.1Eukaryotic 18S ribosomal RNA
XR_890054.3heterogeneous nuclear ribonucleoprotein A2/B1, transcript variant X6
Table A2. Genes differentially upregulated in all macrophages 17 h after infection with FIPV.
Table A2. Genes differentially upregulated in all macrophages 17 h after infection with FIPV.
Gene IdentifierGene Name
XM_003981024.5splicing regulatory glutamic acid and lysine rich protein 1, transcript variant X1
XM_003981036.5cyclin B1
XM_003981430.53-hydroxy-3-methylglutaryl-CoA synthase 1, transcript variant X2
XM_003981592.4centrosomal protein 72, transcript variant X2
XM_003982742.5lanosterol 14-alpha demethylase, transcript variant X1
XM_003982984.3aminoadipate-semialdehyde synthase, transcript variant X1
XM_003984244.4cytoskeleton associated protein-2-like, transcript variant X1
XM_003984844.5farnesyl-diphosphate farnesyltransferase 1, transcript variant X2
XM_003984882.5methylsterol monooxygenase 1, transcript variant X1
XM_003985082.5ELOVL fatty acid elongase 6, transcript variant X1
XM_003986073.4FYVE, RhoGEF and PH domain containing 2, transcript variant X1
XM_003987096.4cyclin B2, transcript variant X2
XM_003988112.5optineurin, transcript variant X7
XM_003988333.5cell division cycle associated 3, transcript variant X1
XM_003988595.5Rap guanine nucleotide exchange factor 3, transcript variant X1
XM_003988664.3trophinin associated protein
XM_003990222.5zinc finger RANBP2-type containing 2, transcript variant X1
XM_003992405.5myelin protein zero like-3, transcript variant X1
XM_003993771.37-dehydrocholesterol reductase, transcript variant X1
XM_003995047.4NDC80, kinetochore complex component
XM_003995401.5transmembrane protein 2, transcript variant X2
XM_003995673.4Rho guanine nucleotide exchange factor 39
XM_003996487.4sperm associated antigen 5, transcript variant X2
XM_003996771.5SH3 and cysteine rich domain 2, transcript variant X1
XM_003998762.5Polo-like kinase 1
XM_003999323.5NIMA related kinase 2
XM_003999347.5abnormal spindle microtubule assembly, transcript variant X2
XM_004000594.4kinesin family member 4A
XM_004000604.5mediator complex subunit 12
XM_006927499.4fms-related tyrosine kinase 4, transcript variant X2
XM_006928617.3phosphodiesterase 4C, transcript variant X3
XM_006929641.4aurora kinase A
XM_006932356.4Nei-like DNA glycosylase 1, transcript variant X1
XM_006934172.3G2 and S-phase expressed 1, transcript variant X1
XM_006934388.4zinc finger and BTB domain containing 40, transcript variant X2
XM_006935662.3nuclear body protein SP140, transcript variant X7
XM_006935975.4ribosomal oxygenase 2, transcript variant X1
XM_006936407.2DNA topoisomerase II binding protein 1, transcript variant X3
XM_006937682.4Ras association domain family member 7, transcript variant X2
XM_006938581.4LIF, interleukin 6 family cytokine, transcript variant X2
XM_006939451.4NIMA related kinase 6, transcript variant X2
XM_006939718.4KIAA0753 ortholog, transcript variant X7
XM_006943087.3farnesyl diphosphate synthase, transcript variant X3
XM_011281000.3VPS54, GARP complex subunit, transcript variant X1
XM_011281524.3fibroblast growth factor receptor 1, transcript variant X6
XM_011283442.3tyrosyl-DNA phosphodiesterase 1, transcript variant X1
XM_011284205.3cytoskeleton associated protein 2
XM_011284366.2DENN domain containing 6B, transcript variant X2
XM_011287426.3centrosomal protein 55, transcript variant X3
XM_011287724.3checkpoint kinase 2, transcript variant X2
XM_011288029.3regulation of nuclear pre-mRNA domain containing 1A, transcript variant X1
XM_011288421.2syntaxin binding protein 1, transcript variant X2
XM_011289261.3sclerostin
XM_011289687.3carnitine palmitoyltransferase 1C, transcript variant X1
XM_011290536.3F-box and leucine rich repeat protein 19, transcript variant X1
XM_011290681.3brain-specific serine protease 4, transcript variant X4
XM_011291526.3kinesin family member C2, transcript variant X2
XM_019812837.27-dehydrocholesterol reductase, transcript variant X3
XM_019813506.2poly(ADP-ribose) glycohydrolase, transcript variant X3
XM_019816169.1talin 1, transcript variant X2
XM_019817080.2discs large MAGUK scaffold protein 4, transcript variant X9
XM_019817127.2aurora kinase B, transcript variant X3
XM_019817247.2angiopoietin-like 6, transcript variant X2
XM_019817546.2synergin gamma, transcript variant X8
XM_019817554.1acetyl-CoA carboxylase alpha, transcript variant X5
XM_019817596.2proline rich 11, transcript variant X3
XM_019817950.2ETS variant 4, transcript variant X5
XM_019819207.2interferon regulatory factor 3, transcript variant X3
XM_019820028.2fucokinase, transcript variant X1
XM_019821899.1complement C3d receptor 2, transcript variant X3
XM_019821902.1complement C3d receptor 2, transcript variant X6
XM_019822306.2upstream transcription factor 1, transcript variant X4
XM_019824070.1centromere protein I
XM_019826681.2TPX2, microtubule nucleation factor, transcript variant X7
XM_019829073.2alpha kinase 1, transcript variant X4
XM_019830601.2peroxisomal biogenesis factor 6, transcript variant X3
XM_019830890.2TTK protein kinase, transcript variant X6
XM_019831446.2acetyl-CoA acetyltransferase 2
XM_019834535.2apolipoprotein L domain containing 1, transcript variant X1
XM_019834543.2apolipoprotein L domain containing 1, transcript variant X8
XM_019835370.2apoptotic peptidase activating factor 1, transcript variant X1
XM_019836058.2nephrocystin 4, transcript variant X2
XM_019839257.2lanosterol synthase, transcript variant X3
XM_023239847.1telomerase reverse transcriptase, transcript variant X11
XM_023239925.1fatty acid desaturase 2, transcript variant X1
XM_023240698.1chondroitin sulfate N-acetylgalactosaminyltransferase 2, transcript variant X8
XM_023241280.1acetoacetyl-CoA synthetase, transcript variant X2
XM_023241432.1HECT domain E3 ubiquitin protein ligase 4, transcript variant X2
XM_023242979.1vav guanine nucleotide exchange factor 2, transcript variant X7
XM_023243604.1zinc finger and BTB domain containing 4, transcript variant X3
XM_023243877.1myosin XVIIIA, transcript variant X20
XM_023244268.1ring finger protein 213, transcript variant X5
XM_023244594.1centrobin, centriole duplication and spindle assembly protein, transcript variant X3
XM_023244664.1acetyl-CoA carboxylase alpha, transcript variant X9
XM_023244901.1mevalonate diphosphate decarboxylase, transcript variant X1
XM_023245382.1dynamin 2, transcript variant X5
XM_023245498.1dedicator of cytokinesis 6
XM_023245844.1lipase E, hormone sensitive type, transcript variant X4
XM_023245914.1DM1 locus, WD repeat containing, transcript variant X1
XM_023246136.1LSM14A, mRNA processing body assembly factor, transcript variant X1
XM_023246494.1cyclin F
XM_023247155.1chromosome E3 C16orf62 homolog, transcript variant X4
XM_023248638.1protein tyrosine kinase 2, transcript variant X19
XM_023249981.1tRNA nucleotidyl transferase 1, transcript variant X3
XM_023250648.1EPH receptor B6, transcript variant X1
XM_023250847.1protein tyrosine phosphatase, receptor type S, transcript variant X10
XM_023250896.1insulin induced gene 1
XM_023252016.1elastin microfibril interfacer 1
XM_023252392.1pericentriolar material 1, transcript variant X15
XM_023252891.1centromere protein E, transcript variant X7
XM_023252992.1heterogeneous nuclear ribonucleoprotein D, transcript variant X6
XM_023253738.1discoidin domain receptor tyrosine kinase 1, transcript variant X1
XM_023253977.1ubiquitin protein ligase E3 component n-recognin 2, transcript variant X3
XM_023254254.1cyclin dependent kinase 19, transcript variant X4
XM_023254265.1FYN proto-oncogene, Src family tyrosine kinase, transcript variant X1
XM_023255117.1kinesin family member 23, transcript variant X2
XM_023255227.1HECT and RLD domain containing E3 ubiquitin protein ligase family member 1, transcript variant X3
XM_023255500.1BUB1 mitotic checkpoint serine/threonine kinase B
XM_023255706.1DLG associated protein 5, transcript variant X1
XM_023256274.1isopentenyl-diphosphate delta-isomerase 1, transcript variant X2
XM_023256524.1cullin 2, transcript variant X1
XM_023257105.1integrin subunit beta 7, transcript variant X5
XM_023257126.1diacylglycerol kinase alpha, transcript variant X2
XM_023257348.1thymopoietin, transcript variant X4
XM_023257700.1casein kinase I, transcript variant X4
XM_023257872.1Holliday junction recognition protein
XM_023259342.1plakophilin 4, transcript variant X6
XM_023259713.1Obscurin-like 1, transcript variant X6
XM_023260160.1cytoplasmic FMR1 interacting protein 2, transcript variant X2
XR_002150997.2uncharacterized LOC109495548
XR_002160083.2MFNG O-fucosylpeptide 3-beta-N-acetylglucosaminyltransferase, transcript variant X2
XR_002737671.1phosphatidylinositol-4-phosphate 5-kinase-like 1, transcript variant X2

References

  1. Addie, D.D.; Jarrett, O. A study of naturally occurring feline coronavirus infections in kittens. Vet. Rec. 1992, 130, 133–137. [Google Scholar] [CrossRef]
  2. Pedersen, N.C. Serologic studies of naturally occurring feline infectious peritonitis. Am. J. Vet. Res. 1976, 37, 1449–1453. [Google Scholar]
  3. Cave, T.A.; Golder, M.C.; Simpson, J.; Addie, D.D. Risk factors for feline coronavirus seropositivity in cats relinquished to a UK rescue charity. J. Feline Med. Surg. 2004, 6, 53–58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Drechsler, Y.; Alcaraz, A.; Bossong, F.J.; Collisson, E.W.; Diniz, P.P. Feline coronavirus in multicat environments. Vet. Clin. North. Am. Small Anim. Pract. 2011, 41, 1133–1169. [Google Scholar] [CrossRef] [PubMed]
  5. Pedersen, N.C.; Boyle, J.F.; Floyd, K.; Fudge, A.; Barker, J. An enteric coronavirus infection of cats and its relationship to feline infectious peritonitis. Am. J. Vet. Res. 1981, 42, 368–377. [Google Scholar]
  6. Foley, J.E.; Poland, A.; Carlson, J.; Pedersen, N.C. Risk factors for feline infectious peritonitis among cats in multiple-cat environments with endemic feline enteric coronavirus. J. Am. Vet. Med. Assoc. 1997, 210, 1313–1318. [Google Scholar] [PubMed]
  7. Pedersen, N.C. Feline Infectious Peritonitis. In Comparative Pathobiology of Viral Diseases; Olsen, R.G., Krakowka, S., Blakeslee, J.R., Jr., Eds.; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar]
  8. Dewerchin, H.L.; Cornelissen, E.; Nauwynck, H.J. Replication of feline coronaviruses in peripheral blood monocytes. Arch. Virol. 2005, 150, 2483–2500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Kipar, A.; Meli, M.L.; Baptiste, K.E.; Bowker, L.J.; Lutz, H. Sites of feline coronavirus persistence in healthy cats. J. Gen. Virol. 2010, 91, 1698–1707. [Google Scholar] [CrossRef] [Green Version]
  10. Vogel, L.; Van der Lubben, M.; te Lintelo, E.G.; Bekker, C.P.; Geerts, T.; Schuijff, L.S.; Grinwis, G.C.; Egberink, H.F.; Rottier, P.J. Pathogenic characteristics of persistent feline enteric coronavirus infection in cats. Vet. Res. 2010, 41, 71. [Google Scholar] [CrossRef] [Green Version]
  11. Pedersen, N.C.; Liu, H.; Scarlett, J.; Leutenegger, C.M.; Golovko, L.; Kennedy, H.; Kamal, F.M. Feline infectious peritonitis: Role of the feline coronavirus 3c gene in intestinal tropism and pathogenicity based upon isolates from resident and adopted shelter cats. Virus Res. 2012, 165, 17–28. [Google Scholar] [CrossRef]
  12. Pesteanu-Somogyi, L.D.; Radzai, C.; Pressler, B.M. Prevalence of feline infectious peritonitis in specific cat breeds. J. Feline Med. Surg. 2006, 8, 1–5. [Google Scholar] [CrossRef] [PubMed]
  13. Norris, J.M.; Bosward, K.L.; White, J.D.; Baral, R.M.; Catt, M.J.; Malik, R. Clinicopathological findings associated with feline infectious peritonitis in Sydney, Australia: 42 cases (1990–2002). Aust. Vet. J. 2005, 83, 666–673. [Google Scholar] [CrossRef] [PubMed]
  14. Tsai, H.Y.; Chueh, L.L.; Lin, C.N.; Su, B.L. Clinicopathological findings and disease staging of feline infectious peritonitis: 51 cases from 2003 to 2009 in Taiwan. J. Feline Med. Surg 2011, 13, 74–80. [Google Scholar] [CrossRef] [PubMed]
  15. Kiss, I.; Kecskemeti, S.; Tanyi, J.; Klingeborn, B.; Belak, S. Preliminary studies on feline coronavirus distribution in naturally and experimentally infected cats. Res. Vet. Sci 2000, 68, 237–242. [Google Scholar] [CrossRef] [PubMed]
  16. Pedersen, N.C. An overview of feline enteric coronavirus and infectious peritonitis virus infections. Feline Pract. 1995, 23, 7–20. [Google Scholar]
  17. Takano, T.; Nakaguchi, M.; Doki, T.; Hohdatsu, T. Antibody-dependent enhancement of serotype II feline enteric coronavirus infection in primary feline monocytes. Arch. Virol. 2017, 162, 3339–3345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Hohdatsu, T.; Yamada, M.; Tominaga, R.; Makino, K.; Kida, K.; Koyama, H. Antibody-dependent enhancement of feline infectious peritonitis virus infection in feline alveolar macrophages and human monocyte cell line U937 by serum of cats experimentally or naturally infected with feline coronavirus. J. Vet. Med. Sci. 1998, 60, 49–55. [Google Scholar] [CrossRef] [Green Version]
  19. Takano, T.; Hohdatsu, T.; Hashida, Y.; Kaneko, Y.; Tanabe, M.; Koyama, H. A “possible” involvement of TNF-alpha in apoptosis induction in peripheral blood lymphocytes of cats with feline infectious peritonitis. Vet. Microbiol. 2007, 119, 121–131. [Google Scholar] [CrossRef]
  20. Kipar, A.; Meli, M.L.; Failing, K.; Euler, T.; Gomes-Keller, M.A.; Schwartz, D.; Lutz, H.; Reinacher, M. Natural feline coronavirus infection: Differences in cytokine patterns in association with the outcome of infection. Vet. Immunol. Immunopathol. 2006, 112, 141–155. [Google Scholar] [CrossRef]
  21. Dean, G.A.; Olivry, T.; Stanton, C.; Pedersen, N.C. In vivo cytokine response to experimental feline infectious peritonitis virus infection. Vet. Microbiol. 2003, 97, 1–12. [Google Scholar] [CrossRef]
  22. Kiss, I.; Poland, A.M.; Pedersen, N.C. Disease outcome and cytokine responses in cats immunized with an avirulent feline infectious peritonitis virus (FIPV)-UCD1 and challenge-exposed with virulent FIPV-UCD8. J. Feline Med. Surg. 2004, 6, 89–97. [Google Scholar] [CrossRef] [PubMed]
  23. de Groot-Mijnes, J.D.; van Dun, J.M.; van der Most, R.G.; de Groot, R.J. Natural history of a recurrent feline coronavirus infection and the role of cellular immunity in survival and disease. J. Virol. 2005, 79, 1036–1044. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Haagmans, B.L.; Egberink, H.F.; Horzinek, M.C. Apoptosis and T-cell depletion during feline infectious peritonitis. J. Virol. 1996, 70, 8977–8983. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Mehrbod, P.; Harun, M.S.; Shuid, A.N.; Omar, A.R. Transcriptome analysis of feline infectious peritonitis virus infection. Methods Mol. Biol. 2015, 1282, 241–250. [Google Scholar] [CrossRef] [PubMed]
  26. Harun, M.S.; Kuan, C.O.; Selvarajah, G.T.; Wei, T.S.; Arshad, S.S.; Hair Bejo, M.; Omar, A.R. Transcriptional profiling of feline infectious peritonitis virus infection in CRFK cells and in PBMCs from FIP diagnosed cats. Virol. J. 2013, 10, 329. [Google Scholar] [CrossRef] [Green Version]
  27. Watanabe, R.; Eckstrand, C.; Liu, H.; Pedersen, N.C. Characterization of peritoneal cells from cats with experimentally-induced feline infectious peritonitis (FIP) using RNA-seq. Vet. Res. 2018, 49, 81. [Google Scholar] [CrossRef] [Green Version]
  28. Goto-Koshino, Y.; Ohno, K.; Nakajima, M.; Mochizuki, H.; Kanemoto, H.; Tsujimoto, H. A rapid and simple method to obtain canine peripheral blood-derived macrophages. J. Vet. Med. Sci 2011, 73, 773–778. [Google Scholar] [CrossRef] [Green Version]
  29. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  30. Bray, N.L.; Pimentel, H.; Melsted, P.; Pachter, L. Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol. 2016, 34, 525–527. [Google Scholar] [CrossRef]
  31. Tarazona, S.; Furio-Tari, P.; Turra, D.; Pietro, A.D.; Nueda, M.J.; Ferrer, A.; Conesa, A. Data quality aware analysis of differential expression in RNA-seq with NOISeq R/Bioc package. Nucleic Acids Res. 2015, 43, e140. [Google Scholar] [CrossRef] [Green Version]
  32. 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, e47. [Google Scholar] [CrossRef] [PubMed]
  33. Bindea, G.; Mlecnik, B.; Hackl, H.; Charoentong, P.; Tosolini, M.; Kirilovsky, A.; Fridman, W.H.; Pages, F.; Trajanoski, Z.; Galon, J. ClueGO: A Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 2009, 25, 1091–1093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  35. Smoot, M.E.; Ono, K.; Ruscheinski, J.; Wang, P.L.; Ideker, T. Cytoscape 2.8: New features for data integration and network visualization. Bioinformatics 2011, 27, 431–432. [Google Scholar] [CrossRef] [Green Version]
  36. Shirato, K.; Chang, H.W.; Rottier, P.J.M. Differential susceptibility of macrophages to serotype II feline coronaviruses correlates with differences in the viral spike protein. Virus Res. 2018, 255, 14–23. [Google Scholar] [CrossRef]
  37. Stoddart, C.A.; Scott, F.W. Intrinsic resistance of feline peritoneal macrophages to coronavirus infection correlates with in vivo virulence. J. Virol. 1989, 63, 436–440. [Google Scholar] [CrossRef] [Green Version]
  38. Pedersen, N.C. A review of feline infectious peritonitis virus infection: 1963-2008. J. Feline Med. Surg. 2009, 11, 225–258. [Google Scholar] [CrossRef]
  39. Poland, A.M.; Vennema, H.; Foley, J.E.; Pedersen, N.C. Two related strains of feline infectious peritonitis virus isolated from immunocompromised cats infected with a feline enteric coronavirus. J. Clin. Microbiol. 1996, 34, 3180–3184. [Google Scholar] [CrossRef]
  40. Chakraborty, P.; Vervelde, L.; Dalziel, R.G.; Wasson, P.S.; Nair, V.; Dutia, B.M.; Kaiser, P. Marek’s disease virus infection of phagocytes: A de novo in vitro infection model. J. Gen. Virol. 2017, 98, 1080–1088. [Google Scholar] [CrossRef]
  41. Takano, T.; Ohyama, T.; Kokumoto, A.; Satoh, R.; Hohdatsu, T. Vascular endothelial growth factor (VEGF), produced by feline infectious peritonitis (FIP) virus-infected monocytes and macrophages, induces vascular permeability and effusion in cats with FIP. Virus Res. 2011, 158, 161–168. [Google Scholar] [CrossRef]
  42. Denney, L.; Branchett, W.; Gregory, L.G.; Oliver, R.A.; Lloyd, C.M. Epithelial-derived TGF-beta1 acts as a pro-viral factor in the lung during influenza A infection. Mucosal Immunol. 2018, 11, 523–535. [Google Scholar] [CrossRef] [Green Version]
  43. He, T.S.; Chen, T.; Wang, D.D.; Xu, L.G. HAUS8 regulates RLRVISA antiviral signaling positively by targeting VISA. Mol. Med. Rep. 2018, 18, 2458–2466. [Google Scholar] [CrossRef]
  44. Martinand, C.; Montavon, C.; Salehzada, T.; Silhol, M.; Lebleu, B.; Bisbal, C. RNase L inhibitor is induced during human immunodeficiency virus type 1 infection and down regulates the 2-5A/RNase L pathway in human T cells. J. Virol. 1999, 73, 290–296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Soderholm, S.; Kainov, D.E.; Ohman, T.; Denisova, O.V.; Schepens, B.; Kulesskiy, E.; Imanishi, S.Y.; Corthals, G.; Hintsanen, P.; Aittokallio, T.; et al. Phosphoproteomics to Characterize Host Response During Influenza A Virus Infection of Human Macrophages. Mol. Cell Proteomics 2016, 15, 3203–3219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Chang, C.K.; Chen, C.J.; Wu, C.C.; Chen, S.W.; Shih, S.R.; Kuo, R.L. Cellular hnRNP A2/B1 interacts with the NP of influenza A virus and impacts viral replication. PLoS ONE 2017, 12, e0188214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Dave, P.; George, B.; Balakrishnan, S.; Sharma, D.K.; Raheja, H.; Dixit, N.M.; Das, S. Strand-specific affinity of host factor hnRNP C1/C2 guides positive to negative-strand ratio in Coxsackievirus B3 infection. RNA Biol. 2019, 16, 1286–1299. [Google Scholar] [CrossRef]
  48. Yasunaga, J.; Lin, F.C.; Lu, X.; Jeang, K.T. Ubiquitin-specific peptidase 20 targets TRAF6 and human T cell leukemia virus type 1 tax to negatively regulate NF-kappaB signaling. J. Virol. 2011, 85, 6212–6219. [Google Scholar] [CrossRef] [Green Version]
  49. Junger, W.G. Immune cell regulation by autocrine purinergic signalling. Nat. Rev. Immunol. 2011, 11, 201–212. [Google Scholar] [CrossRef] [Green Version]
  50. Paoletti, A.; Raza, S.Q.; Voisin, L.; Law, F.; Pipoli da Fonseca, J.; Caillet, M.; Kroemer, G.; Perfettini, J.L. Multifaceted roles of purinergic receptors in viral infection. Microbes Infect. 2012, 14, 1278–1283. [Google Scholar] [CrossRef]
  51. Burbage, M.; Gasparrini, F.; Aggarwal, S.; Gaya, M.; Arnold, J.; Nair, U.; Way, M.; Bruckbauer, A.; Batista, F.D. Tuning of in vivo cognate B-T cell interactions by Intersectin 2 is required for effective anti-viral B cell immunity. Elife 2018, 7. [Google Scholar] [CrossRef] [Green Version]
  52. Cajee, U.F.; Hull, R.; Ntwasa, M. Modification by ubiquitin-like proteins: Significance in apoptosis and autophagy pathways. Int. J. Mol. Sci 2012, 13, 11804–11831. [Google Scholar] [CrossRef] [PubMed]
  53. Mahnke, J.; Schumacher, V.; Ahrens, S.; Kading, N.; Feldhoff, L.M.; Huber, M.; Rupp, J.; Raczkowski, F.; Mittrucker, H.W. Interferon Regulatory Factor 4 controls TH1 cell effector function and metabolism. Sci. Rep. 2016, 6, 35521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Tjernlund, A.; Walther-Jallow, L.; Behbahani, H.; Screpanti, V.; Nowak, P.; Grandien, A.; Andersson, J.; Patterson, B.K. Leukemia inhibitor factor (LIF) inhibits HIV-1 replication via restriction of stat 3 activation. AIDS Res. Hum. Retroviruses 2007, 23, 398–406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Spiegel, M.; Pichlmair, A.; Martinez-Sobrido, L.; Cros, J.; Garcia-Sastre, A.; Haller, O.; Weber, F. Inhibition of Beta interferon induction by severe acute respiratory syndrome coronavirus suggests a two-step model for activation of interferon regulatory factor 3. J. Virol. 2005, 79, 2079–2086. [Google Scholar] [CrossRef] [Green Version]
  56. Afonso, P.V.; Zamborlini, A.; Saib, A.; Mahieux, R. Centrosome and retroviruses: The dangerous liaisons. Retrovirology 2007, 4, 27. [Google Scholar] [CrossRef] [Green Version]
  57. Full, F.; van Gent, M.; Sparrer, K.M.J.; Chiang, C.; Zurenski, M.A.; Scherer, M.; Brockmeyer, N.H.; Heinzerling, L.; Stürzl, M.; Korn, K.; et al. Centrosomal protein TRIM43 restricts herpesvirus infection by regulating nuclear lamina integrity. Nat. Microbiol. 2019, 4, 164–176. [Google Scholar] [CrossRef]
  58. Naghavi, M.H.; Walsh, D. Microtubule Regulation and Function during Virus Infection. J. Virol. 2017, 91. [Google Scholar] [CrossRef] [Green Version]
  59. Wiede, F.; Shields, B.J.; Chew, S.H.; Kyparissoudis, K.; van Vliet, C.; Galic, S.; Tremblay, M.L.; Russell, S.M.; Godfrey, D.I.; Tiganis, T. T cell protein tyrosine phosphatase attenuates T cell signaling to maintain tolerance in mice. J. Clin. Investig. 2011, 121, 4758–4774. [Google Scholar] [CrossRef] [Green Version]
  60. Karaky, M.; Fedetz, M.; Potenciano, V.; Andres-Leon, E.; Codina, A.E.; Barrionuevo, C.; Alcina, A.; Matesanz, F. SP140 regulates the expression of immune-related genes associated with multiple sclerosis and other autoimmune diseases by NF-kappaB inhibition. Hum. Mol. Genet. 2018, 27, 4012–4023. [Google Scholar] [CrossRef]
Figure 1. Multidimensional scaling (MDS) analysis on 27 RNA-seq samples from feline macrophages infected with feline infectious peritonitis virus (FIPV) at two different time points (2 and 17 h). Colors indicate infection status: Non-infected (green), 2 h of infection (blue), and 17 h of infection (red). Macrophages were derived from six different cats that are represented by distinct symbol shapes on the chart: Cat #1 (circle), cat #2 (square), cat #3 (asterisk), cat #4 (triangle), cat #5 (diamond), and cat #6 (cross symbol).
Figure 1. Multidimensional scaling (MDS) analysis on 27 RNA-seq samples from feline macrophages infected with feline infectious peritonitis virus (FIPV) at two different time points (2 and 17 h). Colors indicate infection status: Non-infected (green), 2 h of infection (blue), and 17 h of infection (red). Macrophages were derived from six different cats that are represented by distinct symbol shapes on the chart: Cat #1 (circle), cat #2 (square), cat #3 (asterisk), cat #4 (triangle), cat #5 (diamond), and cat #6 (cross symbol).
Cells 09 01431 g001
Figure 2. Venn diagram of differentially expressed genes in macrophages after infection with the feline coronavirus. (a) Number of genes downregulated at 2 h. (b) Number of genes upregulated at 2 h. (c) Number of genes downregulated at 17 h. (d) Number of genes upregulated at 17 h.
Figure 2. Venn diagram of differentially expressed genes in macrophages after infection with the feline coronavirus. (a) Number of genes downregulated at 2 h. (b) Number of genes upregulated at 2 h. (c) Number of genes downregulated at 17 h. (d) Number of genes upregulated at 17 h.
Cells 09 01431 g002
Figure 3. Cytoscape network output from a KEGG-based gene enrichment analysis performed by clueGO on the 132 differentially upregulated genes (main intersection from Figure 2d) in macrophages infected with the feline coronavirus (FCoV) for 17 h.
Figure 3. Cytoscape network output from a KEGG-based gene enrichment analysis performed by clueGO on the 132 differentially upregulated genes (main intersection from Figure 2d) in macrophages infected with the feline coronavirus (FCoV) for 17 h.
Cells 09 01431 g003
Table 1. Viral RNA load per sample in both raw read count summation and normalized fragments per kilobase million (FPKM) average from all virus’ genes combined, obtained in cellular extracts of either macrophages or CRFK cells infected with FIPV, at 2 and 17 h post-infection. Cats #1–6 were healthy cats from the UC Davis cat colony, ages five months (5M) to four years (4Y). In macrophages, there was a limited viral presence/uptake of FIPV at 2 h, or not detected (nd), and no evidence of significant replication at 17 h (cat #6 n/a (no sample available)). In CRFK cells, viral RNA increased several log-fold from 2 to 17 h.
Table 1. Viral RNA load per sample in both raw read count summation and normalized fragments per kilobase million (FPKM) average from all virus’ genes combined, obtained in cellular extracts of either macrophages or CRFK cells infected with FIPV, at 2 and 17 h post-infection. Cats #1–6 were healthy cats from the UC Davis cat colony, ages five months (5M) to four years (4Y). In macrophages, there was a limited viral presence/uptake of FIPV at 2 h, or not detected (nd), and no evidence of significant replication at 17 h (cat #6 n/a (no sample available)). In CRFK cells, viral RNA increased several log-fold from 2 to 17 h.
MØ OriginRead CountsRead CountsFPKM AverageFPKM Average
(age of cat)2 h17 h2 h17 h
cat #1 (2Y)72724.46 × 10−21.38 × 10−4
cat #2 (2Y)375nd1.82 × 10−2nd
cat #3 (4Y)ndndndnd
cat #4 (5M)179815559.93 × 10−29.05 × 10−2
cat #5 (2Y)6732023.62 × 10−28.69 × 10−3
cat #6 (2Y)245n/a1.4 × 10−1n/a
CRFK 13715459,8991.38 × 10−117.4
CRFK 2348712.19 × 10−41.54 × 10−1
CRFK 380845,2084.71 × 10−21.70
CRFK 4157671,7469.04 × 10−22.91
Table 2. Genes #1–10 were upregulated at 2 h in all samples, gene #11 was downregulated at 2 h in all samples.
Table 2. Genes #1–10 were upregulated at 2 h in all samples, gene #11 was downregulated at 2 h in all samples.
#Gene IdentifierGene Name
1XM_019815278.2SMAD family member 4, transcript variant X3
2XM_019816257.2ATP binding cassette subfamily A member 1, transcript variant X2
3XM_019819160.2HAUS augmin-like complex subunit 8, transcript variant X1
4XM_019819460.2FCH domain only 1, transcript variant X4
5XM_019832772.2heterogeneous nuclear ribonucleoprotein C (C1/C2), transcript variant X1
6XM_023240709.1pantothenate kinase 1, transcript variant X2
7XM_023240796.1ligand-dependent corepressor, transcript variant X4
8XM_023242837.1ubiquitin specific peptidase 20, transcript variant X4
9XM_023248013.1RAB interacting factor, transcript variant X2
10XM_023255193.1transcription factor 7, transcript variant X5
11XM_019832447.2ATPase phospholipid transporting 8B4 (putative), transcript variant X2
Table 3. Enriched KEGG pathways of differentially expressed genes (DEGs) in macrophages 17 h after infection with FIPV. Term P-Value corrected: Corrected with Benjamini and Hochberg.
Table 3. Enriched KEGG pathways of differentially expressed genes (DEGs) in macrophages 17 h after infection with FIPV. Term P-Value corrected: Corrected with Benjamini and Hochberg.
KEGG Pathways 17 h Downregulated
KEGG IDKEGG PathwayTerm p-ValueTerm p-Value CorrectedFold EnrichmentAssociated Genes Found
KO 0000330Arginine and proline metabolism6.73 × 10−46.73 × 10−45.88[CNDP2, OAT, ODC1]
KEGG Pathways 17 h Upregulated
KEGG IDKEGG PathwayTerm p-ValueTerm p-Value CorrectedFold EnrichmentAssociated Genes Found
KO 0000100Steroid biosynthesis4.13 × 10−62.89 × 10−518.2[DHCR7, FDFT1, LOC101083499, LSS]
KO 0000280Valine, leucine, and isoleucine degradation2.35 × 10−32.35 × 10−35.77[AACS, ACAT2, HMGCS1]
KO 0000650Butanoate metabolism3.40 × 10−45.94 × 10−411.1[AACS, ACAT2, HMGCS1]
KO 0000900Terpenoid backbone biosynthesis1.58 × 10−43.68 × 10−414.3[ACAT2, HMGCS1, MVD]
KO 0004110Cell cycle3.89 × 10−45.45 × 10−44.07[CCNB1, CCNB2, CHEK2, PLK1, TTK]
KO 0004115p53 signaling pathway3.27 × 10−51.14 × 10−46.85[APAF1, CCNB1, CCNB2, CHEK2, GTSE1]
KO 0004914Progesterone-mediated oocyte maturation1.08 × 10−31.25 × 10−34.49[AURKA, CCNB1, CCNB2, PLK1]
Table 4. Gene ontology enrichment for differentially expressed genes in macrophages 17 h after infection with FIPV.
Table 4. Gene ontology enrichment for differentially expressed genes in macrophages 17 h after infection with FIPV.
Gene Ontology for Genes Differentially Downregulated at 17 h
GO Molecular Function CompleteFold EnrichmentRaw p-ValueFDR
extracellular matrix binding (GO:0050840)90.546.35 × 10−62.37 × 10−2
collagen binding (GO:0005518)61.831.85 × 10−53.45 × 10−2
GO Biological Process CompleteFold EnrichmentRaw p-ValueFDR
regulation of endothelial cell proliferation (GO:0001936)56.333.29 × 10−84.49 × 10−4
positive regulation of angiogenesis (GO:0045766)34.146.07 × 10−62.76 × 10−2
positive regulation of vasculature development (GO:1904018)31.018.77 × 10−62.99 × 10−2
regulation of epithelial cell proliferation (GO:0050678)19.844.81 × 10−63.28 × 10−2
regulation of peptidase activity (GO:0052547)16.311.22 × 10−53.33 × 10−2
Gene Ontology for Genes Differentially Upregulated at 17 h
GO Molecular Function CompleteFold EnrichmentRaw p-ValueFDR
condensed nuclear chromosome outer kinetochore (GO:0000942)>1001.99 × 10−53.32 × 10−2
GO Molecular Function CompleteFold EnrichmentRaw p-ValueFDR
regulation of spindle organization (GO:0090224)>1004.63 × 10−61.26 × 10−2
establishment of spindle orientation (GO:0051294)>1004.63 × 10−61.05 × 10−2
establishment of spindle localization (GO:0051293)83.298.53 × 10−61.16 × 10−2
spindle localization (GO:0051653)69.411.41 × 10−51.75 × 10−2
mitotic spindle organization (GO:0007052)48.433.88 × 10−53.79 × 10−2
mitotic sister chromatid segregation (GO:0000070)45.522.18 × 10−61.49 × 10−2
microtubule cytoskeleton organization involved in mitosis (GO:1902850)45.522.18 × 10−69.94 × 10−3
mitotic nuclear division (GO:0140014)41.321.61 × 10−72.20 × 10−3
sister chromatid segregation (GO:0000819)35.65.55 × 10−69.47 × 10−3
spindle organization (GO:0007051)35.145.82 × 10−68.83 × 10−3
nuclear chromosome segregation (GO:0098813)22.043.44 × 10−53.91 × 10−2
nuclear division (GO:0000280)21.693.46 × 10−61.18 × 10−2
organelle fission (GO:0048285)19.835.29 × 10−61.03 × 10−2
mitotic cell cycle process (GO:1903047)13.33.49 × 10−53.66 × 10−2
microtubule cytoskeleton organization (GO:0000226)12.484.69 × 10−54.27 × 10−2

Share and Cite

MDPI and ACS Style

Drechsler, Y.; Vasconcelos, E.J.R.; Griggs, L.M.; Diniz, P.P.P.V.; Collisson, E. Host Gene Expression of Macrophages in Response to Feline Coronavirus Infection. Cells 2020, 9, 1431. https://doi.org/10.3390/cells9061431

AMA Style

Drechsler Y, Vasconcelos EJR, Griggs LM, Diniz PPPV, Collisson E. Host Gene Expression of Macrophages in Response to Feline Coronavirus Infection. Cells. 2020; 9(6):1431. https://doi.org/10.3390/cells9061431

Chicago/Turabian Style

Drechsler, Yvonne, Elton J. R. Vasconcelos, Lisa M. Griggs, Pedro P. P. V. Diniz, and Ellen Collisson. 2020. "Host Gene Expression of Macrophages in Response to Feline Coronavirus Infection" Cells 9, no. 6: 1431. https://doi.org/10.3390/cells9061431

APA Style

Drechsler, Y., Vasconcelos, E. J. R., Griggs, L. M., Diniz, P. P. P. V., & Collisson, E. (2020). Host Gene Expression of Macrophages in Response to Feline Coronavirus Infection. Cells, 9(6), 1431. https://doi.org/10.3390/cells9061431

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