Transcriptional profiling of the host cell response to feline immunodeficiency virus infection
© Ertl and Klein; licensee BioMed Central Ltd. 2014
Received: 8 January 2014
Accepted: 12 March 2014
Published: 19 March 2014
Skip to main content
© Ertl and Klein; licensee BioMed Central Ltd. 2014
Received: 8 January 2014
Accepted: 12 March 2014
Published: 19 March 2014
Feline immunodeficiency virus (FIV) is a widespread pathogen of the domestic cat and an important animal model for human immunodeficiency virus (HIV) research. In contrast to HIV, only limited information is available on the transcriptional host cell response to FIV infections. This study aims to identify FIV-induced gene expression changes in feline T-cells during the early phase of the infection. Illumina RNA-sequencing (RNA-seq) was used identify differentially expressed genes (DEGs) at 24 h after FIV infection.
After removal of low-quality reads, the remaining sequencing data were mapped against the cat genome and the numbers of mapping reads were counted for each gene. Regulated genes were identified through the comparison of FIV and mock-infected data sets. After statistical analysis and the removal of genes with insufficient coverage, we detected a total of 69 significantly DEGs (44 up- and 25 down-regulated genes) upon FIV infection. The results obtained by RNA-seq were validated by reverse transcription qPCR analysis for 10 genes.
Out of the most distinct DEGs identified in this study, several genes are already known to interact with HIV in humans, indicating comparable effects of both viruses on the host cell gene expression and furthermore, highlighting the importance of FIV as a model system for HIV. In addition, a set of new genes not previously linked to virus infections could be identified. The provided list of virus-induced genes may represent useful information for future studies focusing on the molecular mechanisms of virus-host interactions in FIV pathogenesis.
Feline immunodeficiency virus (FIV) was first isolated from domestic cats in 1986 and has since then become an important animal model for human immunodeficiency virus (HIV) research [1, 2]. Both viruses belong to the genus of lentiviruses and can be transmitted pre- and postnatally, via blood transfer and mucosal contact [2, 3]. FIV can infect a broad range of cell types, including T-cells, B-cells, macrophages, and cells of the central nervous system. The genomes of FIV and HIV share the three main structural genes: gag, pol and env, as well as the two accessory genes: rev and vif. However, HIV encodes four additional accessory genes not present in FIV: vpr, vpu, tat and nef. In contrast, the orf-A gene can be found uniquely in FIV. Despite these differences on the genome level, FIV induced pathogeneses display striking similarities to human AIDS . FIV progresses through three clinical stages, finally leading to an acquired immunodeficiency syndrome that increases the incidence of opportunistic infections and secondary diseases . After the infection of a target cell, lentiviruses parasitize the cellular machinery to complete their life cycles. Following virus entry into the cell, the viral RNA genome is reverse transcribed and subsequently integrated into the cellular genome. The host cell machinery is then used to generate viral transcripts. These transcripts will be partially spliced and used as templates for the translation of the respective viral proteins. On the contrary, unspliced transcripts are incorporated into new virus particles, assembled from the structural viral proteins . During all these steps, viral factors interact with multiple cellular proteins and hence, affect the normal course of cellular processes. These virus-induced changes of physiological processes can be detected on transcriptional and protein levels. For HIV, multiple cellular genes have been detected as differentially expressed at several stages of infection . Some of these genes directly interact with viral proteins, whereas others might be only side-effects of the virus-induced changes in the host cell environment. Currently, only limited data is available about the impact of FIV infection on the host cell transcription. However, previous studies using cDNA microarrays suggest that transcriptional changes induced by FIV differ in between different cell types . Furthermore, microarrays have been used to analyze the consequences of the viral Orf-A expression on the cellular mRNA profile . In the present study, we use for the first time next-generation RNA sequencing (RNA-seq) to investigate the transcriptional host cell response to FIV . T-cells were infected with FIV and the virus induced gene expression changes were investigated at 24 h post infection (hpi). The most significantly affected genes were additionally investigated by reverse transcriptase qPCR (RT-qPCR) at 8 and 24 hpi . The results of this study will contribute gaining deeper insights into the complex network of virus-host interactions in FIV pathogenesis.
Numbers of Illumina sequencing reads per replicate sample
Numbers of reads
Mapping cat genome
Mapping FIV genome
FIV 24 hpi
FIV 24 hpi
FIV 24 hpi
FIV 24 hpi
List of top 20 most significantly regulated genes
Log2 read counts
Validation by RT-qPCR
FIV 24 hpi
Since viruses have only a few genes, the use of the cellular machinery for the generation of viral transcripts and proteins is essential for their replication. Thus, virus replication interferes with normal cellular processes and therefore affects the expression of cellular genes. At present, only limited information is available on the effect of FIV infections on the host cell transcription. Here, we perform transcriptome analysis to investigate FIV induced gene expression changes in infected T-cells. A total of 69 DEGs could be identified at 24 hpi (Additional file 1: Table S1). As the validity of transcriptome data is dependent on sufficient numbers of sequencing reads mapping to all individual genes , we considered only genes with a minimum coverage (RPKM values > 0.5) for differential expression analysis. Due to this stringent threshold level, some lowly expressed genes were removed from the analysis. However, this approach could improve the reliability of the final list of DEGs. Furthermore, the expression patterns of 10 genes were confirmed by RT-qPCR (Figure 2), confirming that RNA-seq is a suitable method for the quantification of RNA expression levels. The comparison of RT-qPCR results for two time points (8 and 24 hpi) revealed generally lower expression changes for the 10 investigated genes at 8 hpi, thus indicating time-dependent induction of particular genes during the initial phase of infection. So far, only limited information is available about transcriptional changes modulated by FIV. However, previous studies investigated the impact of the viral orf-A gene on the host cell transcription. Here, alterations in several genes, responsible for post-transcriptional RNA modification and protein ubiquitination, were identified in T-cells expressing Orf-A, independently of the residual viral genes . In contrast to that, we could not detect changes for these genes in our data. We assume that these discrepancies can be mainly attributed to different experimental setups, especially regarding the exposure time to the stimulus, as we investigated the early antiviral response whereas the previous study used stable Orf-A expressing cell lines . Moreover, these results may also indicate differences in the cellular responses to individual viral proteins compared to complete viruses.
Out of the 69 FIV induced genes (Additional file 1: Table S1); several genes have already been described to be affected also by infection with other viruses, in particular HIV. OASL, the most significantly regulated gene upon FIV infection (up-regulated 3.4-fold; Table 2) is known as an antiviral protein, induced by a variety of viruses, such as the influenza A and hepatitis C virus [15, 16]. Other examples include DHX58 and DDX58 (up 2.6- and 2.8-folds), belonging to the group of RNA helicases. These enzymes are involved in many processes of the cellular RNA metabolism including the sensing of viral RNAs and the mediation of antiviral immune responses . ACHE (up 2.6-fold) encodes for an enzyme necessary for the termination of neuro-transmitted signals. Furthermore, the plasma secretion of this protein is accomplished via exosome-like vesicles that are also used by HIV for the release of viral Nef proteins. Thus, HIV increases ACHE expression and vesicle production in order to enhance the export of its own proteins . We detected TGM2 to be highly up-regulated (up 5.4-fold) after FIV infection, similar to previous findings for HIV . This multifunctional protein is involved in several cellular processes including apoptosis. The induction of TGM2 might be an antiviral mechanism of the host cells in order to trigger apoptosis and interfere with virus replication . HMGN2 (down 2.1-fold) is presumed to be involved in HIV replication due to interactions with the viral integrase . CDKN1A (other names: p21/waf1; up 2.1-fold) encodes a cyclin-dependent kinase inhibitor that functions as a regulator of cell cycle progression, especially in the promotion of cell cycle arrest. HIV affects the function of the encoded protein and thus, influences cell cycle progression . The induction of cell cycle arrest might benefit viruses because of the resulting suppression of immune responses and the generation of an optimized cellular environment for viral reproduction . In contrast, the increase of ZFP36 (tristetraprolin; up 2.0-fold) might be an antiviral response of the host cell. ZFP36, due to its zinc finger motifs, is able to interact with specific AU-rich regions of HIV RNAs and consequently targets them for rapid degradation . Together with other cellular defense mechanisms, this process may represent an additional antiviral barrier that aims to block virus replication at the transcription level. Similar functions, as a transcriptional repressor of HIV have also been described for BCL6 (up 2.1-fold) . Other examples of genes affected by FIV that have previously been associated with HIV pathogenesis include CSF1, IFI44, HSPE1 and CXCL11. The protein encoded by CSF1 (up 2.0-fold) is a secreted cytokine that controls the production, differentiation and function of macrophages. The induction of this pro-survival cytokine by HIV facilitates protection infected cells from apoptosis . IFI44 (MTAp44; up 2.6-fold) is inducible by interferons (IFNs) which are released by cells in response to the presence of viruses, bacteria and parasites. IFNs trigger the early antiviral defense of the immune system that aims to eradicate the pathogens. IFI44 is frequently upregulated in virus infected cells generating an antiproliferative state, impeding viral replication [27, 28]. HSPE1 (up 3.0-fold) encodes a major heat shock protein. These stress proteins are expressed in response to a broad range of stimuli such as heat and microbial infections. Additionally, interactions with HIV suggest functions in the regulation of viral gene expression and replication . Another IFN induced gene is CXCL11 (up 3.1-fold) which functions as a chemokine and is the dominant ligand for the CXC receptor-3 on immune cells. Its release initiates a chemotactic response in activated T-cells recruiting them to lymphoid organs. In the case of lentiviral infections, the migration of additional cells to infected lymph nodes increases the numbers of potential target cells and thus, contributes to the spread of infection .
Several of the DEGs identified in this study are known to interact with HIV accessory proteins that are missing in FIV. The HIV Tat protein functions as an activator of viral RNA transcription and is responsible for the induction of several cellular genes, such as CDKN1A , BCL6 , IFI44 , C1QTNF5, EEFCC1  and COX4I1 . For CDKN1A, additional interactions have been described with HIV Vpr . As we observed expression changes for all these genes also in this study, similar mechanisms can be assumed for FIV. Here, other viral genes might adopt the particular functions attributed to HIV Tat and Vpr.
Taken together, we investigated the virus induced gene expression changes in FIV infected T-cells. Using next-generation sequencing technology, a total of 69 DEGs were identified at 24 hpi. Out of those genes most significantly affected by FIV, several genes could be associated to known interactions with HIV in humans. These findings indicate a similar effect of both viruses on the host cell transcription and further underline the importance of FIV as a model system for HIV research. The virus-induced genes identified in this study represent potential cellular factors involved in FIV pathogenesis and may serve as a basis for future investigations aiming to provide new insights into the molecular mechanisms of lentiviral infections.
FeT-J cells were purchased from ATCC (CRL-11967). Cells were cultivated in RPMI 1640 medium supplemented with 10% FBS, 1 mM sodium pyruvate, 10 mM HEPES and 0.05 mM 2-mercaptoethanol. All media and supplements were purchased from PAA Laboratories (Pasching, Austria). A molecular clone of the FIV Petaluma strain was propagated in FeT-J cells . Parts of the virus genome, gag and orf-A genes (Genbank: JF411740, JF411742) have been sequenced for previous studies . Herewith, we detected 99% sequence homology with the wild-type Petaluma isolate (Genbank: M25381).
Supernatants from FIV-infected FeT-J cells were harvested by centrifugation and 0.45 μm filtered. For virus quantification, the amount of FIV RNA in the supernatant was measured by RT-qPCR. Therefore, viral RNA was isolated and DNase I digested using the Direct-zol RNA Miniprep Kit (Zymo Research, Irvine, USA) and used as a template for RT-qPCR analysis with FIV gag specific primers and a TaqMan probe . Primer sequences can be found in Additional file 2: Table S2. FIV RNA copy numbers were estimated relative to a tenfold dilution series of an in-vitro transcribed RNA standard. PCR reactions conditions and standard generation have been described earlier . The obtained FIV RNA copy numbers were then converted into virions/ml supernatant. Virus preparations were stored at -80°C until further use for infection experiments.
A spinoculation protocol was performed in order to increase the numbers of virus particles binding to the target cells . Briefly, FeT-J cells were resuspended in FIV-containing supernatant at a MOI of 10 and a cell density of 1 × 106 cells/ml. Cells were then seeded into 6-well plates and centrifuged at 1200 g for 2 h at 25°C. Afterwards, the cells were washed three times with cold medium, transferred into new culture plates and cultivated for 8 and 24 h. Control cells were mock-infected with FIV-free supernatant and treated in the same way. Four replicate samples were infected for each time point. The harvested cells were split in half and either frozen at -20°C for DNA isolation or immediately resuspended in Trizol reagent (Life Technologies, Carlsbad, USA) for subsequent RNA isolation.
FIV provirus DNA in infected cells was estimated by qPCR in order to quantify the infection rates. DNA was isolated utilizing a commercial kit (DNeasy Blood and Tissue Kit, Qiagen, Hilden, Germany) according to the manufacturers’ instructions. qPCR analysis was performed with TaqMan probe assays targeting the FIV gag and the cat RPP30 gene, 30 kDa subunit of ribonuclease P (Additional file 2: Table S2). The RPP30 sequence was used for a blast search against the cat genome, locating one hit on chromosome D2. Thus, analog to the human genome , we assume that there are two allelic copies RPP30 per cell in the cat genome and used this qPCR assay for the quantification of cell numbers (2 copies RPP30 = 1 cell). For the generation of plasmid DNA standards, parts of the FIV gag and RPP30 gene were amplified by PCR and cloned into the pCR2.1 plasmid using the TOPO TA cloning Kit (Invitrogen, Carlsbad, USA). The prepared plasmid DNA was used for the generation of standard dilutions ranging from 1 × 10° – 1 × 106 copies FIV or RPP30. Relating to both standards, FIV and RPP30 copy numbers were determined and converted into FIV copies/cell. Reaction conditions for the TaqMan probe qPCR have been described earlier .
Cellular RNA was isolated using the Direct-zol RNA Miniprep Kit (Zymo Research) including an on-column DNase I digestion for the removal of gDNA. RNA quality control was performed by capillary electrophoretic separation of the samples on the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, USA) and subsequent estimation of RNA integrity numbers. 1 μg of total RNA were used as the starting material for the preparation of cDNA libraries utilizing the TruSeq RNA Sample Prep Kit v2 (Illumina, San Diego, USA) following the recommended protocol. Sequencing was performed on a Genome Analyzer IIx instrument (Illumina) with reagents taken from the TruSeq SR Cluster Kit v5 and TruSeq SBS Kit v5 (Illumina). The lengths of the generated sequencing reads were 37 bp. The raw sequencing data are available in the ArrayExpress database [http://www.ebi.ac.uk/arrayexpress] under accession number E-MTAB-2083.
The Trimmomatic filtering tool was used for quality control and trimming of the sequencing reads . The data were scanned for contamination with adapter sequences from library preparation. Additionally, all bases below an Illumina quality score of 25 and reads shorter than 25 bp were removed. The trimmed data set was then mapped against the cat genome release 72 from Ensembl database [ftp://ftp.ensembl.org/pub/release-72/genbank/felis_catus/] and the FIV Petaluma genome (Genbank: M25381) using the QSeq v5 analysis software (DNASTAR, Madison, USA) implementing the default options for Illumina short reads (matching of a minimum of 20 bases and at least 80% of bases within each read). The numbers of raw reads mapping to annotated genes on the genome reference were counted and fold changes were calculated by comparing the mean values of four FIV infected samples to four control samples. P-values (Student’s t-test) were calculated in QSeq v5 and adjusted according to the Benjamini-Hochberg method . Genes were considered to be differentially expressed; exhibiting fold changes > 2, adjusted p-values < 0.05 and RPKM expression values greater than 0.5 in all replicate samples . For those genes, which were specified as ‘novel’ or ‘uncharacterized’ in the present version of the cat genome, a blast search in NCBI database was performed and the homologous human gene symbols were used for further analysis.
Total RNA was converted into cDNA using the High Capacity Reverse Transcription Kit (Life Technologies) following the manufacturer’s instructions. Primers were designed using the Primer Express 3.0 software (Life Technologies) and validated by the determination of PCR reaction efficiencies as previously described . All primers are described in Additional file 2: Table S2. RT-qPCR was performed in 20 μl reactions including 2 μl cDNA template, 200 nM of each primer, 0.2 mM of each dNTP, 3 mM MgCl2, 1 x Solution S, 1 x buffer B2, 0.4 x EvaGreen fluorescent dye (Biotium, Hayward, USA) and 1 unit of HOT FIREPol DNA polymerase. All reagents for RT-qPCR were purchased from Solis BioDyne (Tartu, Estonia). Reactions were run on a Viia 7 real-time PCR instrument (Life Technologies) using the following temperature protocol: 95°C for 10 min, 45 cycles of 95°C for 15 sec and 60°C for 1 min, followed by the generation of melting curves. Target gene expression levels were normalized to those of three references genes, namely ABL, B2M and RSP7 . Relative expression changes of FIV- vs. mock-infected cells were estimated using the comparative 2^-∆∆CT method .
The RNA-sequencing reads reported in this paper have been deposited in the ArrayExpress database [http://www.ebi.ac.uk/arrayexpress] under accession number E-MTAB-2083.
Differentially expressed genes
Feline immunodeficiency virus
Human immunodeficiency virus
Hours post infection
Multiplicity of infection
Reads per kilobase per million mapped reads
Reverse transcription quantitative PCR.
We wish to thank Werner Rust for providing reagents for Illumina sequencing.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.