Skip to main content

Transcriptional profiling of feline infectious peritonitis virus infection in CRFK cells and in PBMCs from FIP diagnosed cats

Abstract

Background

Feline Infectious Peritonitis (FIP) is a lethal systemic disease, caused by the FIP Virus (FIPV); a virulent mutant of Feline Enteric Coronavirus (FECV). Currently, the viruses virulence determinants and host gene expressions during FIPV infection are not fully understood.

Methods

RNA sequencing of Crandell Rees Feline Kidney (CRFK) cells, infected with FIPV strain 79–1146 at 3 hours post infection (h.p.i), were sequenced using the Illumina next generation sequencing approach. Bioinformatic’s analysis, based on Felis catus 2X annotated shotgun reference genome, using CLC bio Genome Workbench mapped both control and infected cell reads to 18899 genes out of 19046 annotated genes. Kal’s Z test statistical analysis was used to analyse the differentially expressed genes from the infected CRFK cells. Real time RT-qPCR was developed for further transcriptional profiling of three genes (PD-1, PD-L1 and A3H) in infected CRFK cells and Peripheral Blood Mononuclear Cells (PBMCs) from healthy and FIP-diseased cats.

Results

Based on Kal’s Z-test, with False Discovery Rate (FDR) <0.05 and >1.99 fold change on gene expressions, a total of 61 genes were differentially expressed by both samples, where 44 genes were up-regulated and the remainder were down-regulated. Most genes were closely clustered together, suggesting a homogeneous expression. The majority of the genes that were significantly regulated, were those associated with monocytes-macrophage and Th1 cell functions, and the regulation of apoptosis. Real time RT-qPCR developed focusing on 2 up-regulated genes (PD-L1 and A3H) together with an apoptosis associated gene PD-1 expressions in FIPV infected CRFK cells and in PBMCs from healthy and FIP diagnosed cats produced concordant results with transcriptome data.

Conclusion

The possible roles of these genes, and their importance in feline coronaviruses infection, are discussed.

Background

Feline coronaviruses are enveloped, positive sense RNA viruses that can be classified into two biotypes, namely low virulent Feline Enteric Coronavirus (FECV) and highly virulent Feline Infectious Peritonitis Virus (FIPV). FECV is very common in the cat population worldwide, and has been shown to have infected 20-60% pet cats and shed by 75-100% cats in multi-cat environments[1, 2]. Of those shedding the virus, 1-5% will develop Feline Infectious Peritonitis (FIP) disease[3]. Within the biotypes, the viruses are differentiated further into serotype I and serotype II, based on virus neutralizing antibodies. Type I FECV and FIPV strains are more ubiquitous worldwide and are more likely to cause clinical FIP, while type II strains are less common, but more adaptable to cell culture[2].

It has been suggested that FIPV, the causative agent for FIP, is a mutant form of FECV[4, 5]; where several possible nature of mutation responsible for the increase in virulence has been characterized. Studies have shown that several mutations throughout the FIPV genome were detected, but mutations at 3c membrane protein and 7b secretory glycoprotein genes were suggested to be responsible for transforming FECV to FIPV[4, 5]. A recent study revealed that mutation of the S1/S2 locus and modulation of a furin recognition site normally present in the S gene of FECVs is a critical contributing factor for development of FIP[6]. Furthermore, it was found that FIPV infection is associated with T cell depletion by apoptosis; although the virus cannot infect CD4+ and CD8+ T cells[7, 8]. Therefore, apoptosis of CD4+ and CD8+ T cells is probably caused by mediators from infected macrophages and/or intestinal epithelial cells[8, 9]. Hence, little is known about the interaction of the virus and host cells; especially the early cellular transcriptional responses towards virus infection, virus mechanism of inducing T cell apoptosis, and the absence of cell-mediated immunity (CMI) response in FIP infected cats.

The use of a next generation sequencing approach in RNA sequencing has facilitated understanding in defining the expression profiles of cellular responses during pathogen infections. This method has been proven to be helpful in explaining the pathogenesis of various viruses[10, 11], including Feline Immunodeficiency Virus (FIV) infection[12, 13]. Furthermore, the availability of complete 1.9X of cat genome, using the Whole Genome Shotgun (WGS) approach, provides valuable information for the bioinformatic’s analysis of feline host responses, following pathogen infection. Moreover, the cat genome contigs were aligned, mapped, and annotated to NCBI annotated genome sequence of six index mammalian genomes (human, chimpanzee, mouse, rat, dog, and cow) using MegaBLAST[14].

Previous study has shown that more than 70% of FIPV strain 79–1146 were internalized by CRFK cells at 3 hours post infection[15]. Hence, in this study, mRNA from CRFK cells infected with FIPV strain 79–1146 at 3 h.p.i were sequenced using Illumina next generation sequencing technology. The generated data was then analyzed using CLC bio Genomic Workbench, where the genes were compared to Felis catus 1.9X annotated shotgun reference genome. Kal’s Z-test on expression proportions[16] was used to determine significantly expressed genes. Genes expressed with a False Discovery Rate (FDR) <0.05 and >1.99 fold change were considered for further analysis.

Results

Early gene expression of FIPV infected CRFK cells

Overall, the trimmed sequence reads match to 25,689 annotated transcripts; where only 215 (0.8%) were statistically significant (Kal’s Z test, p < 0.05), and out of the significant matched, only 96 (44.7%) transcripts were expressed with fold absolute change of 2 or more. Of these 96 transcripts, 76 (79.2%) were up-regulated and the remainder were down-regulated. After BLAST analysis, 76 up-regulated transcripts were matched to 44 genes while 20 down-regulated transcripts were matched to 17 genes. Of the 44 up-regulated genes, there were 2 transcripts per gene for 32 genes, but only one transcript per gene for the remaining 12 genes. Meanwhile, of the 17 down-regulated genes, there was one transcript per gene for 14 genes and 2 transcripts per gene for 3 genes.

As shown in Figure 1, the RPKM of control samples was plotted against the RPKM of infected samples; and as such, genes with equal expression will line-up on the diagonal identity line while genes with different expression values will either be over or under the diagonal line. The further away the point is from the identity line, the larger is the difference between its expressions in one experiment compared with the other. Except for 4 genes (MX1, RSAD2, PLIN2, and SERPINB1), most genes from both samples were closely clustered together, thus suggesting a homogeneous expression. The plot excluded 3 genes (CCL8, RNF7, and RPL39) that had infinite fold change expressions.

Figure 1
figure 1

Scatter plot of control RNA-seq RPKM versus infected RNA-seq RPKM of significant transcripts (p < 0.05) with absolute fold change value of 2 or more. Most genes from both samples are closely clustered except for 4 genes label 1 to 6. 1 and 2 are MX1 transcripts 1 and 2, respectively, 3 and 4 = RSAD2 transcripts 1 and 2, respectively, 5 = PLIN2 and 6 = SERPINB1.

Table 1 shows the top 20 up-regulated genes (in decreasing order) and their functions. The majority of the genes were those associated with immune response, while the remainder were associated with apoptosis, cell cycle, cytokine, and ubiquitination activities. Furthermore, there were also 5 Interferon Stimulated Genes (ISGs) coded for 6 proteins (RSAD2, A3C, A3H, MDA5, IFI35, and MX1) that were involved in inhibiting viral entry, replication, and production. Interestingly, one gene (PD-L1), which negatively regulates immune response in viral infection, was also found to be highly up-regulated. Meanwhile, the majority of the down-regulated genes were involved in pro-inflammatory cytokine’s activation, CMI, and anti-apoptosis activities (Table 2). Two unique down-regulated genes (RNF7 and RPL39) were found to be expressed in control uninfected cells only, where the former had anti-apoptotic effect and the latter translated RNA to protein.

Table 1 List of top 20 up regulated genes and their functions
Table 2 List of all 17 down regulated genes and their functions

Real-time RT-qPCR analysis of CRFK infected cells

Three host genes (A3H, PD-1, and PD-L1) were selected for real-time RT-qPCR analysis, because they were highly up-regulated and may play important roles during FIPV induced disease; judging by their functions. A3H was involved in viral RNA and DNA editing, causing mutation, while PD-1 and its ligand (PD-L1) were involved in programmed cell death that was associated with negative regulation of immune response. Comparisons of fold change results of real-time PCR and transcriptome study for A3H and PD-L1 genes revealed almost similar levels of fold changes. The transcriptome resulted in 4.07 and 5.65 fold change for 2 transcripts of A3H while real-time PCR resulted in 5.23 ± 1.15 fold changes (Tables 1 and3). Meanwhile, transcriptome resulted in 3.29 and 4.24 fold changes for 2 transcripts of PD-L1 gene while real-time PCR resulted in 3.97 ± 0.29 fold changes. In the case of the PD-1 gene, RNA sequencing was unable to detect the gene expression, due to low coverage (i.e., data not shown), although real-time PCR was able to detect an up-regulation of the gene at 3 h.p.i.

Table 3 Fold changes of A3H, PD-1 and PD-L1 genes in CRFK cells at different time points following infection with FIPV

FIPV induced high and exceptionally high A3H expression at 3 and 48 h.p.i. respectively, but from 6 to 24 h.p.i., the gene was down-regulated. FIPV infected cells also showed high up-regulation of PD-1 expression at 3 h.p.i. and moderately up-regulated at 12 h.p.i but were being down-regulated at 6, 9, 24 and 48 h.p.i. Meanwhile, PD-L1 gene was consistently down-regulated from 3 hours to 48 h.p.i.

Fold change expression analysis of PBMCs of FIP diagnosed cats

Peripheral Blood Mononuclear Cells (PBMCs), obtained from cats with clinical signs associated with FIP (Table 4), were purified and analysed with real-time PCR. In general, all of the FIP diagnosed cats expressed the PD-1 and PD-L1 genes more than 2 folds, while only 2 cats expressed A3H gene more than 2 folds (Table 5). The highest expression fold for A3H, PD-1 and PD-L1 genes was detected from cat no. 3 at 3.44 ± 0.36, 68.13 ± 19.45, and 96.94 ± 21.54 fold changes, respectively. Meanwhile, cat no. 5 showed less fold changes of the genes compared to the other FIP diagnosed cats.

Table 4 Clinical and serological parameters of FIPV diagnosed and normal cats used in this study
Table 5 Fold changes of A3H, PD-1 and PD-L1 genes from FIP diagnosed cats

Discussion

The pathogenesis of feline coronavirus infection is unclear. The reference feline genome sequence assembly of transcriptome analysis of early infection (3 h.p.i.) of CRFK cells with FIPV 79–1146 showed that the expressions of 215 transcripts (0.8% of the trimmed annotated) were statistically significant, based on Kal’s Z test. Only 96 transcripts, which consisted of 44 up-regulated genes and 17 down-regulated genes, were expressed with fold absolute changes of 2 or more. Since only one sample per group was analysed, Kal’s Z test was used to determine the significant differences in the expression profiles. Study has shown that this test evaluates single sample against another single sample where each group in an experiment has only one sample[16]. This test is based on an approximation of the binomial distribution by the normal distribution considering proportions rather than raw counts so that it can be used reliably on libraries of different size.

The transcriptional profiles of selected genes in FIPV in vitro infected cells, as well as cats diagnosed with FIP, were explored. The expressions of A3H, which involved in viral RNA and DNA editing causing mutation during RNA virus infection[17] and PD-1 and its ligand (PD-L1), which are involved in programmed cell death and negative regulation of T cells immune response[18], were characterised. A3H has antiretroviral activity by generating lethal hypermutations in viral genomes and is associated with increased resistance to HIV-1 infection in certain populations[19]. In the case of felines, A3H (but not A3C) has been found to reduce the infectivity of feline leukemia virus[17]. It is interesting to detect that the expression of A3H is readily expressed at higher levels in PBMCs than in FIPV infected cells (Tables 3 and5), indicating the possible involvement of the gene in antiviral activity. Furthermore, the up-regulation of the gene is less, compared to PD-1 and PD-L1 in FIPV diagnosed cats. In addition, the expressions of A3H are significantly up-regulated in FIPV infected CRFK cells at 48 h.p.i. The actual reason for this high expression of A3H is unclear, but suggests that this gene is essential in restricting viral replication or forming part of the type 1 interferon-induced innate response, since recent study has indicated the involvement of A3H in restricting virus replication[17, 20].

Our results show that up-regulation of PD-1 and PD-L1 gene’s expression in PBMCs occurred in cats diagnosed with the FIP disease. In general, their expressions are correlated to each other. Similar patterns were also observed in chronic FIV infection[21] and in HIV infection in humans[22]; where they were associated with increasing immune dysfunction and T cell depletion. Previous studies have shown that although FIPV cannot infect CD4+ and CD8+ T-cells[8], cats infected with the virus showed T cell depletion by apoptosis resulting in an acute immunodeficiency[7]. Hence, we hypothesized that PD-L1 could be a mediator that mediates apoptosis of CD4+ and CD8+ T-cells, since its expression is found in a wide range of nonhematopoietic cells[23]. Furthermore, the blockade of PD-L1 expression was found to enhance T cell immunity and cytokine production[24]. Nevertheless, more studies are required to confirm our hypothesis.

It has been established that cats infected with FIPV undergo extensive tissue destruction due to inflammatory reaction[2]. Even though the study was performed on infected CRFK cells, a kidney epithelial cells, it is interesting to note that transcriptome analysis of cells proposed that the inflammation process was associated with proinflammatory and Th1-like cytokines production, due to the up-regulation of several chemokine genes, such as CCL8, CXCL10, and CCL17; and genes associated with innate immune responses, such as PHF11, ATF3, and IRF1. Further study on infection of FIPV in particularly type I FIPV on macrophage-like cells namely fcwf-4 cells and samples from FIP diagnosed cats will add more value to our findings. Furthermore, the down-regulation of CRIP1 (a T helper regulatory gene) (Table 2) also suggests that FIPV infection is associated with Th1 response, based on a study on CRIP1 gene in mice, transgenic mice, and in murine cell line[25]. In that study, they found that the down-regulation of CRIP was associated with the expression of IL-2, IFN-γ, and TNF-α. Previous studies have shown that in vitro and in vivo FIPV infections were associated with TNF-α and IFN-γ expressions[9, 26].

Based on previous studies on the growth of FIPV 79–1146 in CRFK cells, production of progeny virus start between 3 and 6 hours post inoculation and increased rapidly until 12 hours post inoculation[27]. In our study, we found that at 3 hours, a few complete virus genomes that completely aligned to FIPV 79–1146 reference genome sequence (data not shown) has already been assembled. MX1 expression was up-regulated in FIPV infection (Table 1), similar to other RNA virus infections, its role in FIP pathogenesis is still unclear and requires further investigation. Previous studies have shown that the expression of MX1 gene inhibits viral replication during various RNA virus infections[28, 29]. Early in vitro FIPV infection is also associated with a marked increase in the expression of RSAD2 (radical S-adenosyl methionine domain-containing protein 2) also known as viperin (Table 1). Previous studies have shown viperin involvement in inhibiting viral RNA and/or protein synthesis during different virus infections, such as the West Nile virus, Dengue virus, and hepatitis virus[30, 31]. Hence, further study is required to define the role of viperin in FIPV replication and infection. Serine Proteinase Inhibitor clade B member 1 (SERPINB1) is the only gene that was markedly down-regulated (Table 2). SERPINB1 functions as an inhibitor of the neutrophil serine proteases, found at inflammatory sites where the inhibition of the gene prevents tissue destruction by phagocytic cells during the virus clearance process by infiltrating neutrophils and monocytes[32, 33]. Thus, SERPINB1 down-regulation could possibly be part of innate immune response in recruiting phagocytic cells for the proteolytic destruction of infected cells. Besides genes that modulate T cell functions, several genes with pro-apoptotic and/or anti-apoptotic are also differentially regulated; thus highlighting their functions in regulating the apoptosis of virus infected cells (Additional file1: Table S1 and Additional file1: Table S2). However, the pro-apoptotic gene YAP-65 or YAP-1, which was found up-regulated in FIV infection in CRFK cells[34], was not up-regulated in this study.

Conclusion

In conclusion, this present study has described the transcriptional profiles of cellular genes in vitro system can be applied to an in vivo situation and the possible involvement of A3H, PD-1, and PD-L1 genes during FIPV infection. However, further studies are required to elucidate the roles of these genes, and their interactions with other genes, during FIPV replication and infection especially in in vivo model.

Materials and methods

Virus and cell lines

Monolayers of Crandell Rees Feline Kidney (CRFK) cells (ATCC® no. CCL-94™) were grown in a base media consisting of Minimum Essential Medium (MEM), 10% Fetal Bovine Serum (FBS), 2% penicillin and streptomycin, 2% amphotericin B, and non-essential amino acids at 37°C and 5% CO2. For transcriptome study, CRFK cells were infected with FIPV 79–1146 (ATCC® no. VR-2126™) at a Multiplicity Of Infection (MOI) of 2. The virus was incubated for one hour for virus attachment. After incubation, 1% FBS MEM was added and the cells were incubated for a further 3 hours. At the end of incubation, the cells were harvested using TrypLE™ and centrifuged twice in a PBS at 4°C for 5 minutes at 1000rpm. Cell pellets were stored at -80°C until RNA extraction. For the control sample, the same process was applied, with the exception that 2 ml of sterile PBS was used to replace the virus.

RNA extraction and sequencing

An RNeasy® Mini Kit (Qiagen®, USA) was used to extract and purify RNA samples as per the methods recommended by the manufacturer. The quality of the extracted RNA was determined by an Ultrospec 3000 Pro UV/Visible spectrophotometer (GE Healthcare, UK), where samples with an absorbance ratio value (A260/A280) of 1.8 to 2.0 were considered for further analysis with an Agilent® 2100 Bioanalyzer. Samples with RNA Integrity Numbers (RIN) 9 to 10, and concentrations higher than 500 ng/μl per sample, were sent for Illumina GAII sequence analysis.

Bioinformatics analysis

A total of 17.3 Gb of sequencing data, comprised of both control and infected samples, was imported into the CLC bio Genomic Workbench (GWB). The sequences were trimmed for adapter sequences and low quality base. The trimmed raw sequences were subjected to RNA-sequence analysis, by mapping them to an annotated feline genome reference sequence[35] accounting for a maximum of two gaps or mismatches in each sequence. Unpaired group comparisons, based on Reads Per Kilo base per Million (RPKM)[36], were chosen as expression values for comparison. Kal’s Z test statistical analyses, based on False Discovery Rate (FDR) <0.05 and fold change >1.99 were used to filter the expressed transcripts. The resulting list was then BLAST at NCBI servers (http://www.ncbi.nlm.nih.gov/) using GWB’s built-in BLAST (blastn, refseq_rna or nr databases, mammals only). Homologous sequences with the lowest e-value, highest score, and lowest percentage of gaps to the query sequence, was selected as the transcript identity.

Real-time RT-qPCR analysis of FIPV infected CRFK cells

In order to validate the transcriptome results, the expression profiles of 3 genes (A3H, PD-1, and PD-L1) were analysed using real-time PCR. Briefly, viral RNA from FIPV strain 79–1146 infected CRFK cells at 3, 6, 9, 12, 24, and 48 hours post infection (h.p.i.) were collected and processed as described previously. Control cells, inoculated with PBS only, were used as a control. Primers were designed using Primer-BLAST (http://www.ncbi.nlm.nih.gov/tools/primer-blast/) and synthesized by AITbiotech PTE LTD (Singapore) (Table 6).

Table 6 Primers sequences used in this study

The reactions were performed using SensiFAST™ SYBR No-ROX One Step kit (Bioline Ltd, UK) on Bio-Rad CFX 96™ Real-Time System, with C1000™ Thermal Cycler (Bio-Rad Laboratories, USA). Briefly, the reaction mixture of 20 μl contained 10 μl 2× SensiFAST SYBR No-ROX One-Step mix, 0.5 μl forward & reverse primers (5 nM for GADPH, PD-L1 and A3H, 3 nM for PD-1 and 10 nM for YWHAZ), 0.2 μl RT, 0.4 μl RiboSafe RNase inhibitor, 2.4 μl H2O, and 6 μl extracted RNA. The RT-qPCR reaction conditions were as follows; one cycle at 45°C for 10mins, one cycle at 95°C for 2 mins, and 35 cycles at 95°C for 5 secs; then 57°C (YWHAZ), 58°C (PD-L1), 59°C (GAPDH), 64°C (A3H), and 65°C (PD-1) for 20 secs; and finally, at 72°C for 5 secs. One cycle for the dissociation curve for all reactions was added and melting curve analysis was performed. Data generated from the technical triplicate experiment was analysed with 2–ΔΔCT method[37] using Bio-Rad CFX Manager version 2.0. GAPDH and/or YWHAZ genes were chosen as reference genes, based on previous studies[38, 39].

Real-time RT-qPCR analysis of peripheral blood mononuclear cells from FIP diagnosed cats

Besides FIPV infected cell cultures, Peripheral Blood Mononuclear Cells (PBMCs) were also used to analyse the transcriptome results. Six FIP diagnosed domestic short hair breed cats, with ages ranging from 8 months to adult, that were admitted to University Veterinary Hospital, UPM, were considered for this study (Table 4). The cats tested negative for FeLV and FIV antibodies, but positive for FCoV antibodies, and showed abdominal/pleural effusion. Meanwhile, 5 healthy cats, with negative results for FCoV, FeLV, and FIV antibodies, were selected as controls. The kits for FCoV, FeLV and FIV antibody tests originated from Biogal’s feline coronavirus (FCoV) [FIP] ImmunoComb® Antibody test kit (Biogal Galed Laboratories, USA) and IDEXX’s SNAP® Combo FeLV Ag/FIV Antibody test kit (IDEXX Laboratories, USA), respectively. The tests were performed as per the methods recommended in their respective manuals. The health assessment and blood collection of the cats were performed by a trained and certified veterinarian (GTS). The sampling were performed according to internationally recognized guidelines and recommended by the Animal Care and Use Committee at the Faculty of Veterinary Medicine, Universiti Putra Malaysia.

Two to 5 ml of cat blood was drawn and stored at 4°C in BD Vacutainer® (BD USA) EDTA-K2 tubes. Parts of the blood were used for the test kits, while the rest was processed for PBMCs extraction. PBMCs were isolated using the Ficoll-Paque™ Plus (GE Healthcare, USA) method, according to the manufacturers protocol. Total RNA from PBMCs was isolated using an RNeasy mini plus kit (Qiagen, Germany), as described by the manufacturer. RNA quantity and purity was measured and assessed using a Nanodrop Nanophotometer P-class (Implen GmbH, Germany). The isolated RNA samples were either kept at -80°C for further analysis, or immediately used for real-time RT-qPCR analysis.

Abbreviations

BLAST:

Basic Local Alignment Search Tool

CRFK:

Crandell Rees Feline Kidney

CMI:

Cell-mediated immunity

FCOV:

Feline coronavirus

FDR:

False discovery rate

FECV:

Feline enteric coronavirus

FIPV:

Feline infectious peritonitis virus

FIV:

Feline immunodeficiency virus

kb:

Kilobase

min:

Minute

ml:

Mililitre

MOI:

Multiplicity of infection

NGS:

Next generation sequencing

PBMCs:

Peripheral Blood Mononuclear Cells

PCR:

Polymerase chain reaction

RPKM:

Reads per kilobase of exon model per million mapped reads.

References

  1. Sharif S, Arshad SS, Hair-Bejo M, Omar AR, Zeenathul NA, Hafidz MA: Prevalence of feline coronavirus in two cat populations in Malaysia. J Feline Med Surg 2009, 12: 1031-1034.

    Article  Google Scholar 

  2. Pedersen NC: A review of feline infectious peritonitis virus infection: 1963–2008. J Feline Med Surg 2009, 11: 225-258. 10.1016/j.jfms.2008.09.008

    Article  PubMed  Google Scholar 

  3. Hartmann K: Feline infectious peritonitis. Vet Clin North Am Small Anim Pract 2005,35(1):39-79. 10.1016/j.cvsm.2004.10.011

    Article  PubMed  Google Scholar 

  4. Herrewegh AA, Vennema H, Horzinek MC, Rottier PJ, De Groot RJ: The molecular genetics of feline coronaviruses: comparative sequence analysis of the ORF7a/7b transcription unit of different biotypes. Virol 1995, 212: 622-631. 10.1006/viro.1995.1520

    Article  CAS  Google Scholar 

  5. Vennema H, Poland A, Foley J, Pedersen NC: Feline infectious peritonitis viruses arise by mutation from endemic feline enteric coronaviruses. Virol 1998, 243: 150-157. 10.1006/viro.1998.9045

    Article  CAS  Google Scholar 

  6. Licitra BN, Millet JK, Regan AD, Hamilton BS, Rinaldi VD, Duhamel GE, et al.: Mutation in spike protein cleavage site and pathogenesis of feline coronavirus. Emerg Infect Dis 2013. http://dx.doi.org/10.3201/eid1907.121094/

    Google Scholar 

  7. De Groot Mijnes JD, Van Dun JM, Van Der Most RG, De Groot RJ: Natural history of a recurrent feline coronavirus infection and the role of cellular immunity in survival and disease. J Virol 2005, 79: 1036-1044. 10.1128/JVI.79.2.1036-1044.2005

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  8. Haagmans BL, Egberink HF, Horzinek MC: Apoptosis and T-cell depletion during feline infectious peritonitis. J Virol 1996,70(12):8977-8983.

    PubMed  CAS  PubMed Central  Google Scholar 

  9. 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(2–4):121-131.

    Article  PubMed  CAS  Google Scholar 

  10. Assarsson E, Greenbaum JA, Sundström M, Schaffer L, Hammond JA, Pasquetto V, Oseroff C, Hendrickson RC, Lefkowitz EJ, Tscharke DC, Sidney J, Grey HM, Head SR, Peters B, Sette A: Kinetic analysis of a complete poxvirus transcriptome an immediate early class of gene. Proc Natl Acad Sci U S A 2008,105(6):2140-2145. 10.1073/pnas.0711573105

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Nanda S, Havert MB, Calderón GM, Thomson M, Jacobson C, et al.: Hepatic transcriptome analysis of hepatitis C virus infection in chimpanzees defines unique gene expression patterns associated with viral clearance. PLoS ONE 2008,3(10):e3442. doi:10.1371/journal.pone.0003442 10.1371/journal.pone.0003442

    Article  PubMed  PubMed Central  Google Scholar 

  12. Dowling RJ, Bienzle D: Gene-expression changes induced by Feline immunodeficiency virus infection differ in epithelial cells and lymphocytes. J Gen Virol 2005,86(8):2239-2248. 10.1099/vir.0.80735-0

    Article  PubMed  CAS  Google Scholar 

  13. Ertl R, Birzele F, Hildebrandt T, Klein D: Viral transcriptome analysis of feline immunodeficiency virus infected cells using second generation sequencing technology. Vet Immunol Immunopathol 2011,143(3–4):314-324.

    Article  PubMed  CAS  Google Scholar 

  14. Zhang Z, Schwartz S, Wagner L, Miller W: A greedy algorithm for aligning DNA sequences. J Comput Biol 2000, 7: 203-214. 10.1089/10665270050081478

    Article  PubMed  CAS  Google Scholar 

  15. van Hamme E, Dewerchin HL, Cornelissen E, Nauwynck HJ: Attachment and internalization of feline infectious peritonitis virus in feline blood monocytes and Crandell feline kidney cells. J Gen Virol 2007, 88: 2527-2532. 10.1099/vir.0.82991-0

    Article  PubMed  CAS  Google Scholar 

  16. Kal AJ, Van Zonneveld AJ, Benes V, Van Den Berg M, Koerkamp MG, Albermann K, Strack N, Ruijter JM, Richter A, Dujon B, Ansorge W, Tabak HF: Dynamics of gene expression revealed by comparison of serial analysis of gene expression transcript profiles from yeast grown on two different carbon sources. Mol Biol Cell 1999,10(6):1859-1872. 10.1091/mbc.10.6.1859

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Münk C, Beck T, Zielonka J, Hotz-Wagenblatt A, Chareza S, Battenberg M, Thielebein J, Cichutek K, Bravo IG, O’Brien SJ, Löchelt M, Yuhki N: Functions, structure, and read-through alternative splicing of feline APOBEC3 genes. Genome Biol 2008,9(3):R48. 10.1186/gb-2008-9-3-r48

    Article  PubMed  PubMed Central  Google Scholar 

  18. Shen T, Chen X, Chen Y, Xu Q, Lu F, Liu S: Increased PD-L1 expression and PD-L1/CD86 ratio on dendritic cells were associated with impaired dendritic cells function in HCV infection. J Med Virol 2010,82(7):1152-1159. 10.1002/jmv.21809

    Article  PubMed  CAS  Google Scholar 

  19. Li MM, Emerman M: Polymorphism in human APOBEC3H affects a phenotype dominant for subcellular localization and antiviral activity. J Virol 2011,85(16):8197-8207. 10.1128/JVI.00624-11

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  20. Zhen A, Du J, Zhou X, Xiong Y, Yu XF: Reduced APOBEC3H variant anti-viral activities are associated with altered RNA binding activities. PLoS ONE 2012,7(7):e38771. doi:10.1371/journal.pone.0038771 10.1371/journal.pone.0038771

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  21. Folkl A, Wen X, Kuczynski E, Clark ME, Bienzle D: Feline programmed death and its ligand: characterization and changes with feline immunodeficiency virus infection. Vet Immunol Immunopathol 2010,134(1–2):107-114.

    Article  PubMed  CAS  Google Scholar 

  22. Trautmann L, Janbazian L, Chomont N, Said EA, Gimmig S, Bessette B, Boulassel MR, Delwart E, Sepulveda H, Balderas RS, Routy JP, Haddad EK, Sekaly RP: Upregulation of PD-1 expression on HIV-specific CD8+ T cells leads to reversible immune dysfunction. Nat Med 2006, 12: 1198-1202. 10.1038/nm1482

    Article  PubMed  CAS  Google Scholar 

  23. Keir ME, Butte MJ, Freeman GJ, Sharpe AH: PD-1 and its ligands in tolerance and immunity. Annu Rev Immunol 2008,2008(26):677-704.

    Article  Google Scholar 

  24. Brown JA, Dorfman DM, Ma FR, Sullivan EL, Munoz O, Wood CR, Greenfield EA, Freeman GJ: Blockade of programmed death-1 ligands on dendritic cells enhances T cell activation and cytokine production. J Immunol 2003,170(3):1257-1266.

    Article  PubMed  CAS  Google Scholar 

  25. Cousins RJ, Foster LL: Regulation of cysteine-rich intestinal protein, a zinc finger protein, by mediators of the immune response. J Infect Dis 2000,182(1):S81-S84. doi:10.1086/315917

    Article  PubMed  CAS  Google Scholar 

  26. Berg AL, Ekmanb K, Belákb S, Bergb M: Cellular composition and interferon-γ expression of the local inflammatory response in feline infectious peritonitis (FIP). Vet Microbiol 2005,111(1–2):15-23.

    Article  PubMed  CAS  Google Scholar 

  27. Dewerchin HL, Cornelissen E, Nauwynck HJ: Replication of feline coronaviruses in peripheral blood monocytes. Arch Virol 2005, 150: 2483-2500. 10.1007/s00705-005-0598-6

    Article  PubMed  CAS  Google Scholar 

  28. Khaiboullina SF, Rizvanov AA, Deyde VM, St Jeor SC: Andes virus stimulates interferon-inducible MxA protein expression in endothelial cells. J Med Virol 2005,75(2):267-275. 10.1002/jmv.20266

    Article  PubMed  CAS  Google Scholar 

  29. Hoenen A, Liu W, Kochs G, Khromykh AA, Mackenzie JM: West Nile virus-induced cytoplasmic membrane structures provide partial protection against the interferon-induced antiviral MxA protein. J Gen Virol 2007,88(11):3013-3017. 10.1099/vir.0.83125-0

    Article  PubMed  CAS  Google Scholar 

  30. Jiang D, Weidner JM, Qing M, Pan XB, Guo H, Xu C, Zhang X, Birk A, Chang J, Shi PY, Block TM, Guo JT: Identification of five interferon-induced cellular proteins that inhibit west nile virus and dengue virus infections. J Virol 2010,84(16):8332-8341. 10.1128/JVI.02199-09

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  31. Helbig KJ, Lau DT, Semendric L, Harley HA, Beard MR: Analysis of ISG expression in chronic hepatitis C identifies viperin as a potential antiviral effector. Hepatology 2005,42(3):702-710. 10.1002/hep.20844

    Article  PubMed  CAS  Google Scholar 

  32. Gong D, Benarafa C, Hartshorn KL, Remold-O’Donnell E: The neutrophil serine protease inhibitor SerpinB1 protects against inflammatory lung injury and morbidity in influenza virus infection. J Immunol 2009,182(43):10.

    Google Scholar 

  33. Benarafa C, LeCuyer TE, Baumann M, Stolley JM, Cremona TP, Remold-O’Donnell E: SerpinB1 protects the mature neutrophil reserve in the bone marrow. J Leukoc Biol 2011,90(1):21-29. 10.1189/jlb.0810461

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Basu S, Totty NF, Irwin MS, Sudol M, Downward J: Akt phosphorylates the Yes-associated protein, YAP, to induce interaction with 14-3-3 and attenuation of p73-mediated apoptosis. Mol Cell 2003,11(1):11-23. 10.1016/S1097-2765(02)00776-1

    Article  PubMed  CAS  Google Scholar 

  35. The UCSC felCat4 data from the Dec. 2008 catChrV17e draft assembly. . html#cat_credits http://genome.ucsc.edu/goldenPath/credits

  36. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 2008,5(7):621-628. 10.1038/nmeth.1226

    Article  PubMed  CAS  Google Scholar 

  37. Livak K, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001,25(4):402-408. 10.1006/meth.2001.1262

    Article  PubMed  CAS  Google Scholar 

  38. Penning LC, Vrieling HE, Brinkhof B, Riemers FM, Rothuizen J, Rutteman GR, Hazewinkel HA: A validation of 10 feline reference genes for gene expression measurements in snap-frozen tissues. Vet Immunol Immunopathol 2007, 120: 212-222. 10.1016/j.vetimm.2007.08.006

    Article  PubMed  CAS  Google Scholar 

  39. Kessler Y, Helfer-Hungerbuehler AK, Cattori V, Meli ML, Zellweger B, Ossent P, Riond B, Reusch CE, Lutz H, Hofmann-Lehmann R: Quantitative TaqMan® real-time PCR assays for gene expression normalisation in feline tissues. BMC Mol Biol 2009, 10: 106. doi:10.1186/1471-2199-10-106 10.1186/1471-2199-10-106

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors wish to thanks laboratory personnel at Virology Lab, Faculty of Veterinary Medicine, UPM and Laboratory of Vaccines and Immunotherapeutics, Institute of Bioscience, UPM. This project is funded by Fundamental Research Project No:01-11-08-6390FR, Ministry of Higher Education Malaysia. The funding source has no role in this study.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Abdul Rahman Omar.

Additional information

Competing interests

The authors declare that we have no competing interests.

Authors’ contributions

MSRH & COK - Cell culture works; virus inoculation; RNA extraction, purification and quantification; RT-qPCR assays; bioinformatic analysis; data analysis and interpretation; wrote the manuscript. GTS – Diagnose cats; collect blood samples from cats with FIP symptoms; perform FCoV, FIV and FeLV kit tests. TSW – helps with PCR; primers design; RT-qPCR assays and data analysis. SSA, MHB & ARO – secure & manage fundings; coordinated the project; designed experiments and collaborated in analyzing data and writing the manuscript. All the authors have read and approved the final manuscript.

Electronic supplementary material

12985_2013_2316_MOESM1_ESM.docx

Additional file 1: Table S1: List of 76 transcripts from 44 up regulated genes with proportions fold change of 2 or more (Kal’s Z test, FDR < 0.05) with their BLAST results, NCBI accession number and gene product function. Table S2. List of 20 transcripts from 17 down regulated genes with proportions fold change of -2 or more (Kal’s Z test, FDR < 0.05) with their BLAST results, NCBI accession number and gene product function. (DOCX 47 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Harun, M.S.R., Kuan, C.O., Selvarajah, G.T. et al. Transcriptional profiling of feline infectious peritonitis virus infection in CRFK cells and in PBMCs from FIP diagnosed cats. Virol J 10, 329 (2013). https://doi.org/10.1186/1743-422X-10-329

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1743-422X-10-329

Keywords