Skip to main content

Transcriptional profiling of the host cell response to feline immunodeficiency virus infection

Abstract

Background

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.

Results

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.

Discussion and conclusion

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.

Background

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[2]. FIV progresses through three clinical stages, finally leading to an acquired immunodeficiency syndrome that increases the incidence of opportunistic infections and secondary diseases[4]. 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[5]. 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[6]. 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[7]. Furthermore, microarrays have been used to analyze the consequences of the viral Orf-A expression on the cellular mRNA profile[8]. 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[9]. 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[10]. The results of this study will contribute gaining deeper insights into the complex network of virus-host interactions in FIV pathogenesis.

Results

Infection of T-cells and transcriptome sequencing

FeT-J cells, a feline T-lymphocyte cell line, were infected with the FIV Petaluma strain. A high multiplicity of infection (MOI) of 10 virions per cell was used in order to ensure infection of the majority of target cells. FIV- and mock-infected cells were harvested at 8 and 24 hpi. Four replicates were analyzed for each time-point. The amounts of FIV provirus DNA were determined by quantitative PCR (qPCR) to estimate the infection efficiencies. At both time points similar provirus loads (7 FIV copies per cell) were detected, indicating successful infection of all replicate samples (Figure 1). The small decrease observed from 8 to 24 hpi can be explained by continuous cell division, while the production of new virus particles is expected to take a minimum of 24 h. High quality total RNAs (RNA integrity numbers of 10) of FIV (24 hpi) and mock infected cells were used for poly-A mRNA purification and the subsequent preparation of cDNA libraries for transcriptome sequencing. Next generation sequencing analysis on an Illumina platform generated a total of 42–57 million 37-bp reads per replicate sample (Table 1). Out of these, 40–56 million reads passed the quality filtering and were mapped against the cat reference genome. 60% of the filtered reads could be mapped to the cat genome. Additionally, in FIV infected cells 0.2% of reads were assigned to the FIV Petaluma genome suggesting the beginning of viral RNA transcription.

Figure 1
figure 1

Quantification of FIV provirus DNA in infected T-cells. Numbers of FIV DNA copies per cell at 8 and 24 hpi were determined by qPCR.

Table 1 Numbers of Illumina sequencing reads per replicate sample

Identification of differentially expressed genes

For quantification of transcript numbers, the filtered sequencing reads were assigned to the annotated genes present on the current version of the cat genome. Numbers of reads mapping to particular genes were counted and normalized by converting them into RPKM-values (reads per kilobase per million mapped reads)[11]. Out of 19,947 annotated genes present on the reference genome, sample reads could be assigned to 17,998 genes, leaving 1,949 genes without any mapping reads. For differential expression analysis, fold-changes and p-values were calculated in between the mean read counts of FIV-and mock-infected cells. The p-values were adjusted for multiple testing, according to the procedure proposed by Benjamini and Hochberg[12]. This statistical procedure was implemented to account for the multiple comparison problem that occurs during the analysis of high-throughput data. High numbers of hypotheses (e.g. differential expression of multiple genes) that are tested in one data set increase the probability that one of them is incorrectly accepted by chance. Correction methods, like Benjamin-Hochberg, are frequently used for large genomic data and call for adjustments, so that the probability of finding at least one significant result due to chance remains below the implemented significance level[12, 13]. In order to ensure sufficient coverage of all analyzed genes, the threshold for transcription was set to a minimum RPKM-value of 0.5 in all replicate samples. Thus, genes with RPKM < 0.5 were excluded from the analysis. Genes exhibiting fold changes > 2 and adjusted p-values < 0.05 were considered as differentially expressed. This analysis strategy resulted in a total number of 69 differentially expressed genes (DEGs) at 24 h after FIV infection. The complete list of DEGs can be found in Additional file1: Table S1. Out of these genes, 44 were up- and 25 down-regulated. The observed fold changes range from 2- to 5-fold, with the majority of genes being moderately (2- to 3-fold) up-regulated. Out of the 20 most significantly regulated genes, only 3 genes were down- and 17 were up-regulated (Table 2).

Table 2 List of top 20 most significantly regulated genes

Validation of differential expressions by RT-qPCR

RT-qPCR was performed in order to verify the results from transcriptome sequencing for 10 DEGs. The 10 genes were selected out of the 20 most significantly DEGs detected by RNA-seq (Table 2). Expression changes compared to mock-infected cells were additionally analyzed for an earlier time point, at 8 hpi. The FIV induced expression changes determined by RNA-seq could be confirmed by RT-qPCR for all investigated genes at 24 hpi (Figure 2). The highest virus induced changes were detected for TGM2 and OASL, analogue to the sequencing results. Compared to 24 hpi, the observed expression changes were generally lower at 8 hpi.

Figure 2
figure 2

Validation of 10 selected DEGs by RT-qPCR and comparison with RNA-sequencing. RT-qPCR gene expression levels at 8 and 24 hpi were compared to mock infected cells. Fold changes were obtained by merging the mean values of 4 replicates for each time point, previously normalized to the expression of the 3 reference genes ABL, B2M and RSP7. Stars indicate significant (p < 0.05) expression changes.

Discussion

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 file1: Table S1). As the validity of transcriptome data is dependent on sufficient numbers of sequencing reads mapping to all individual genes[14], 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[8]. 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[8]. Moreover, these results may also indicate differences in the cellular responses to individual viral proteins compared to complete viruses.

FIV induced gene expression changes: Similarities to HIV

Out of the 69 FIV induced genes (Additional file1: 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[17]. 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[18]. We detected TGM2 to be highly up-regulated (up 5.4-fold) after FIV infection, similar to previous findings for HIV[19]. 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[20]. HMGN2 (down 2.1-fold) is presumed to be involved in HIV replication due to interactions with the viral integrase[21]. 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[22]. 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[23]. 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[24]. 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)[25]. 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[26]. 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[29]. 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[30].

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[31], BCL6[32], IFI44[27], C1QTNF5, EEFCC1[33] and COX4I1[34]. For CDKN1A, additional interactions have been described with HIV Vpr[35]. 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.

Conclusion

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.

Materials and methods

Cells and viruses

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[36]. Parts of the virus genome, gag and orf-A genes (Genbank: JF411740, JF411742) have been sequenced for previous studies[37]. Herewith, we detected 99% sequence homology with the wild-type Petaluma isolate (Genbank: M25381).

Virus preparation and quantification

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[38]. Primer sequences can be found in Additional file2: 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[38]. 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.

Infection

A spinoculation protocol was performed in order to increase the numbers of virus particles binding to the target cells[39]. 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.

Quantification of FIV DNA

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[38] and the cat RPP30 gene, 30 kDa subunit of ribonuclease P (Additional file2: 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[40], 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[37].

Transcriptome sequencing

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.

Sequencing data analysis

The Trimmomatic filtering tool was used for quality control and trimming of the sequencing reads[41]. 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[12]. 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[11]. 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.

RT-qPCR analysis of gene expression

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[42]. All primers are described in Additional file2: 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[43]. Relative expression changes of FIV- vs. mock-infected cells were estimated using the comparative 2^-∆∆CT method[44].

Data deposition

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.

Abbreviations

DEGs:

Differentially expressed genes

FIV:

Feline immunodeficiency virus

HIV:

Human immunodeficiency virus

hpi:

Hours post infection

IFN:

Interferon

MOI:

Multiplicity of infection

qPCR:

Quantitative PCR

RNA-seq:

RNA-sequencing

RPKM:

Reads per kilobase per million mapped reads

RT-qPCR:

Reverse transcription quantitative PCR.

References

  1. Pedersen NC, Ho EW, Brown ML, Yamamoto JK: Isolation of a T-lymphotropic virus from domestic cats with an immunodeficiency-like syndrome. Science 1987, 235: 790-793. 10.1126/science.3643650

    Article  PubMed  Google Scholar 

  2. Elder JH, Lin YC, Fink E, Grant CK: Feline immunodeficiency virus (FIV) as a model for study of lentivirus infections: parallels with HIV. Curr HIV Res 2010, 8: 73-80. 10.2174/157016210790416389

    Article  PubMed  PubMed Central  Google Scholar 

  3. Medeiros SDO, Martins AN, Dias CGA, Tanuri A, Brindeiro RDM: Natural transmission of feline immunodeficiency virus from infected queen to kitten. Virol J 2012, 9: 99. 10.1186/1743-422X-9-99

    Article  PubMed Central  Google Scholar 

  4. Hartmann K: Clinical aspects of feline immunodeficiency and feline leukemia virus infection. Vet Immunol Immunopathol 2011, 143: 190-201. 10.1016/j.vetimm.2011.06.003

    Article  PubMed  Google Scholar 

  5. Kenyon JC, Lever AM: The molecular biology of feline immunodeficiency virus (FIV). Viruses 2011, 3: 2192-2213. 10.3390/v3112192

    Article  PubMed  PubMed Central  Google Scholar 

  6. Giri MS, Nebozhyn M, Showe L, Montaner LJ: Microarray data on gene modulation by HIV-1 in immune cells: 2000–2006. J Leukocyte Biol 2006, 80: 1031-1043. 10.1189/jlb.0306157

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  8. Sundstrom M, Chatterji U, Schaffer L, de Rozières S, Elder JH: Feline immunodeficiency virus OrfA alters gene expression of splicing factors and proteasome-ubiquitination proteins. Virology 2008, 371: 394-404. 10.1016/j.virol.2007.09.039

    Article  PubMed  PubMed Central  Google Scholar 

  9. Radford AD, Chapman D, Dixon L, Chantrey J, Darby AC, Hall N: Application of next-generation sequencing technologies in virology. J Gen Virol 2012, 2012: 1853-1868.

    Article  Google Scholar 

  10. Klein D: Quantification using real-time PCR technology: Applications and limitations. Trends Mol Med 2002, 8: 257-260. 10.1016/S1471-4914(02)02355-9

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  12. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol 1995, 57: 289-300.

    Google Scholar 

  13. Asomaning N, Archer KJ: High-throughput DNA methylation datasets for evaluating false discovery rate methodologies. Comput Stat Data Anal 2012, 56: 1748-1756. 10.1016/j.csda.2011.10.020

    Article  PubMed  PubMed Central  Google Scholar 

  14. Tarazona S, García-Alcalde F, Dopazo J, Ferrer A, Conesa A: Differential expression in RNA-seq: A matter of depth. Genome Res 2011, 21: 2213-2223. 10.1101/gr.124321.111

    Article  PubMed  PubMed Central  Google Scholar 

  15. Ishibashi M, Wakita T, Esumi M: 2′,5′-Oligoadenylate synthetase-like gene highly induced by hepatitis C virus infection in human liver is inhibitory to viral replication in vitro. Biochem Biophys Res Commun 2010, 392: 397-402. 10.1016/j.bbrc.2010.01.034

    Article  PubMed  Google Scholar 

  16. Melchjorsen J, Kristiansen H, Christiansen R, Rintahaka J, Matikainen S, Paludan SR, Hartmann R: Differential regulation of the OASL and OAS1 genes in response to viral infections. J Interferon Cytokine Res 2009, 29: 199-207. 10.1089/jir.2008.0050

    Article  PubMed  Google Scholar 

  17. Steimer L, Klostermeier D: RNA helicases in infection and disease. RNA Biol 2012, 9: 751-771. 10.4161/rna.20090

    Article  PubMed  Google Scholar 

  18. Ali SA, Huang MB, Campbell PE, Roth WW, Campbell T, Khan M, Newman G, Villinger F, Powell MD, Bond VC: Genetic characterization of HIV type 1 nef-induced vesicle secretion. AIDS Res Hum Retroviruses 2010, 26: 173-192. 10.1089/aid.2009.0068

    Article  PubMed  PubMed Central  Google Scholar 

  19. Amendola A, Gougeon ML, Poccia F, Bondurand A, Fesus L, Piacentini M: Induction of "tissue" transglutaminase in HIV pathogenesis: evidence for high rate of apoptosis of CD4+ T lymphocytes and accessory cells in lymphoid tissues. Proc Natl Acad Sci U S A 1996, 93: 11057-11062. 10.1073/pnas.93.20.11057

    Article  PubMed  PubMed Central  Google Scholar 

  20. Amendola A, Rodolfo C, Caro A, Ciccosanti F, Falasca L, Piacentini M: “Tissue” Transglutaminase Expression in HIV‐Infected Cells. Ann N Y Acad Sci 2001, 946: 108-120.

    Article  PubMed  Google Scholar 

  21. Allouch A, Cereseto A: Identification of cellular factors binding to acetylated HIV-1 integrase. Amino Acids 2011, 41: 1137-1145. 10.1007/s00726-009-0444-3

    Article  PubMed  Google Scholar 

  22. Amini S, Khalili K, Sawaya BE: Effect of HIV-1 Vpr on cell cycle regulators. DNA Cell Biol 2004, 23: 249-260. 10.1089/104454904773819833

    Article  PubMed  Google Scholar 

  23. Goh WC, Rogel ME, Kinsey CM, Michael SF, Fultz PN, Nowak MA, Hahn BH, Emerman M: HIV-1 Vpr increases viral expression by manipulation of the cell cycle: a mechanism for selection of Vpr in vivo. Nature Med 1998, 4: 65-71. 10.1038/nm0198-065

    Article  PubMed  Google Scholar 

  24. Maeda M, Sawa H, Tobiume M, Tokunaga K, Hasegawa H, Ichinohe T, Sata T, Moriyama M, Hall WW, Kurata T, Takahashi H: Tristetraprolin inhibits HIV-1 production by binding to genomic RNA. Microb Infect 2006, 8: 2647-2656. 10.1016/j.micinf.2006.07.010

    Article  Google Scholar 

  25. Baron BW, Desai M, Julia Baber L, Paras L, Zhang Q, Sadhu A, Duguay S, Nucifora G, McKeithan TW, Zeleznik-Le N: BCL6 can repress transcription from the human immunodeficiency virus type I promoter/enhancer region. Genes Chromosomes Cancer 1997, 19: 14-21. 10.1002/(SICI)1098-2264(199705)19:1<14::AID-GCC3>3.0.CO;2-3

    Article  PubMed  Google Scholar 

  26. Swingler S, Mann AM, Zhou J, Swingler C, Stevenson M: Apoptotic killing of HIV-1-infected macrophages is subverted by the viral envelope glycoprotein. PLoS Pathog 2007, 3: 1281-1290.

    Article  PubMed  Google Scholar 

  27. Izmailova E, Bertley FM, Huang Q, Makori N, Miller CJ, Young RA, Aldovini A: HIV-1 Tat reprograms immature dendritic cells to express chemoattractants for activated T cells and macrophages. Nat Med 2003, 9: 191-197. 10.1038/nm822

    Article  PubMed  Google Scholar 

  28. Hallen LC, Burki Y, Ebeling M, Broger C, Siegrist F, Oroszlan-Szovik K, Bohrmann B, Certa U, Foser S: Antiproliferative activity of the human IFN-α-inducible protein IFI44. J Interferon Cytokine Res 2007, 27: 675-680. 10.1089/jir.2007.0021

    Article  PubMed  Google Scholar 

  29. Rawat P, Mitra D: Cellular heat shock factor 1 positively regulates human immunodeficiency virus-1 gene expression and replication by two distinct pathways. Nucleic Acids Res 2011, 39: 5879-5892. 10.1093/nar/gkr198

    Article  PubMed  PubMed Central  Google Scholar 

  30. Foley JF, Yu CR, Solow R, Yacobucci M, Peden KW, Farber JM: Roles for CXC chemokine ligands 10 and 11 in recruiting CD4+ T cells to HIV-1-infected monocyte-derived macrophages, dendritic cells, and lymph nodes. J Immunol 2005, 174: 4892-4900.

    Article  PubMed  Google Scholar 

  31. Thakur BK, Chandra A, Dittrich T, Welte K, Chandra P: Inhibition of SIRT1 by HIV-1 viral protein Tat results in activation of p53 pathway. Biochem Biophys Res Commun 2012, 424: 245-25. 10.1016/j.bbrc.2012.06.084

    Article  PubMed  Google Scholar 

  32. Deregibus MC, Cantaluppi V, Doublier S, Brizzi MF, Deambrosis I, Albini A, Camussi G: HIV-1-Tat protein activates phosphatidylinositol 3-kinase/AKT-dependent survival pathways in Kaposi's sarcoma cells. J Biol Chem 2002, 277: 25195-25202. 10.1074/jbc.M200921200

    Article  PubMed  Google Scholar 

  33. Lin CW, Kuo JHS, Jan MS: The global gene-expression profiles of U-937 human macrophages treated with Tat peptide and Tat-FITC conjugate. J Drug Target 2012, 20: 515-523. 10.3109/1061186X.2012.693496

    Article  PubMed  Google Scholar 

  34. Lecoeur H, Borgne-Sanchez A, Chaloin O, El-Khoury R, Brabant M, Langonné A, Porceddu M, Brier̀e JJ, Buron N, Rebouillat D, Péchoux C, Deniaud A, Brenner C, Briand JP, Muller S, Rustin P, Jacotot E: HIV-1 Tat protein directly induces mitochondrial membrane permeabilization and inactivates cytochrome c oxidase. Cell Death Dis 2012, 3: e282. 10.1038/cddis.2012.21

    Article  PubMed  PubMed Central  Google Scholar 

  35. Chowdhury IH, Wang XF, Landau NR, Robb ML, Polonis VR, Birx DL, Kim JH: HIV-1 Vpr activates cell cycle inhibitor p21/Waf1/Cip1: A potential mechanism of G2/M cell cycle arrest. Virology 2003, 305: 371-377. 10.1006/viro.2002.1777

    Article  PubMed  Google Scholar 

  36. Olmsted RA, Barnes AK, Yamamoto JK, Hirsch VM, Purcell RH, Johnson PR: Molecular cloning of feline immunodeficiency virus. Proc Natl Acad Sci U S A 1989, 86: 2448-2452. 10.1073/pnas.86.7.2448

    Article  PubMed  PubMed Central  Google Scholar 

  37. 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: 314-324. 10.1016/j.vetimm.2011.06.010

    Article  PubMed  Google Scholar 

  38. Klein D, Leutenegger CM, Bahula C, Gold P, Hofmann-Lehmann R, Salmons B, Lutz H, Gunzburg WH: Influence of preassay and sequence variations on viral load determination by a multiplex read-time reverse transcriptase-polymerase chain reaction for feline immunodeficiency virus. J Acquir Immune Defic Syndr 2001, 26: 8-20. 10.1097/00126334-200101010-00002

    Article  PubMed  Google Scholar 

  39. O'Doherty U, Swiggard WJ, Malim MH: Human immunodeficiency virus type 1 spinoculation enhances infection through virus binding. J Virol 2000, 74: 10074-10080. 10.1128/JVI.74.21.10074-10080.2000

    Article  PubMed  PubMed Central  Google Scholar 

  40. Baer M, Nilsen TW, Costigan C, Altman S: Structure and transcription of a human gene for H1 RNA, the RNA component of human RNase P. Nucleic Acids Res 1990, 18: 97-103. 10.1093/nar/18.1.97

    Article  PubMed  PubMed Central  Google Scholar 

  41. Lohse M, Bolger AM, Nagel A, Fernie AR, Lunn JE, Stitt M, Usadel B: RobiNA: A user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic Acids Res 2012, 40: W622-W627. 10.1093/nar/gks540

    Article  PubMed  PubMed Central  Google Scholar 

  42. Klein D, Janda P, Steinborn R, Müller M, Salmons B, Günzburg WH: Proviral load determination of different feline immunodeficiency virus isolates using real-time polymerase chain reaction: Influence of mismatches on quantification. Electrophoresis 1999, 20: 291-299. 10.1002/(SICI)1522-2683(19990201)20:2<291::AID-ELPS291>3.0.CO;2-R

    Article  PubMed  Google Scholar 

  43. Kessler Y, Helfer-Hungerbuehler AK, Cattori V, Meli M, 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. 10.1186/1471-2199-10-106

    Article  PubMed  PubMed Central  Google Scholar 

  44. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2^-ΔΔCT method. Methods 2001, 25: 402-408. 10.1006/meth.2001.1262

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We wish to thank Werner Rust for providing reagents for Illumina sequencing.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Reinhard Ertl.

Additional information

Competing interest

The authors declare that they have no competing interests.

Authors’ contributions

RE and DK conceived and designed the experiments. RE performed the experiments, analyzed the data and wrote the manuscript. All authors edited and approved the final manuscript.

Electronic supplementary material

12985_2014_2397_MOESM1_ESM.xlsx

Additional file 1: Table S1: List of differentially expressed genes in FIV infected T-cells at 24 hpi as detected by RNA-seq analysis. Fold changes refer to mock-infected cells. Read counts (log2) represent the mean values of four replicate samples. P-values were adjusted according to the Benjamini-Hochberg method. (XLSX 16 KB)

Additional file 2: Table S2: Primers used for quantitative PCR. (PDF 127 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

Authors’ original file for figure 2

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Ertl, R., Klein, D. Transcriptional profiling of the host cell response to feline immunodeficiency virus infection. Virol J 11, 52 (2014). https://doi.org/10.1186/1743-422X-11-52

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1743-422X-11-52

Keywords