Aberrant host immune response induced by highly virulent PRRSV identified by digital gene expression tag profiling
© Xiao et al; licensee BioMed Central Ltd. 2010
Received: 6 April 2010
Accepted: 7 October 2010
Published: 7 October 2010
There was a large scale outbreak of the highly pathogenic porcine reproductive and respiratory syndrome (PRRS) in China and Vietnam during 2006 and 2007 that resulted in unusually high morbidity and mortality among pigs of all ages. The mechanisms underlying the molecular pathogenesis of the highly virulent PRRS virus (H-PRRSV) remains unknown. Therefore, the relationship between pulmonary gene expression profiles after H-PRRSV infection and infection pathology were analyzed in this study using high-throughput deep sequencing and histopathology.
H-PRRSV infection resulted in severe lung pathology. The results indicate that aberrant host innate immune responses to H-PRRSV and induction of an anti-apoptotic state could be responsible for the aggressive replication and dissemination of H-PRRSV. Prolific rapid replication of H-PRRSV could have triggered aberrant sustained expression of pro-inflammatory cytokines and chemokines leading to a markedly robust inflammatory response compounded by significant cell death and increased oxidative damage. The end result was severe tissue damage and high pathogenicity.
The systems analysis utilized in this study provides a comprehensive basis for better understanding the pathogenesis of H-PRRSV. Furthermore, it allows the genetic components involved in H-PRRSV resistance/susceptibility in swine populations to be identified.
In 2006 there was an unparalleled large-scale outbreak of so-called 'high fever' disease in many areas of China that affected more than 2,000,000 pigs. There were approximately 400,000 fatal cases in 2006 and 243,000 fatalities in 2007, leading to concerns within the global swine industry and in relation to public health . In March 2007 the disease was identified in the Hai Duong Province of Vietnam and it spread country-wide affecting more than 65,000 pigs . The outbreaks caused extensive concern worldwide . Affected pigs of all ages presented with clinical signs including continuous high fever (40.5°C-42°C), depression, anorexia, dyspnoea, reddening of the skin, edema of the eyelids, conjunctivitis, mild diarrhea, shivering, lamping, and unusually high morbidity (50%-100%) and mortality (20%-100%). Studies demonstrated that highly virulent porcine reproductive and respiratory syndrome virus (H-PRRSV) was the major causative pathogen of the so-called "high fever" disease. Genetic analysis indicated that the H-PRRSVs isolated from China and Vietnam shared a discontinuous deletion of 30 aa in non-structural protein 2 (NSP2), as compared with the North American type PRRSV strains (NA PRRSV). However, the mechanisms underlying the molecular pathogenesis of the H-PRRSV that emerged in China and Vietnam have not been elucidated.
Preliminary results indicated that PRRSV modulates the host immune responses and alters host gene expression. PRRSV infection up-regulated expression of mRNA for interleukin-10 (IL10), interferon gamma (IFN-γ), tumor necrosis factor-alpha (TNF-α), myxovirus resistance 1 (MX1), ubiquitin specific proteases (USP) and toll-like receptors (TLR), and inhibited expression of type I interferons [4–7]. A study concerning the genome-wide transcriptional response of primary alveolar macrophages (PAMs) following infection with the Lelystad PRRSV strain (European type, EU PRRSV) reported that the expression of beta interferon 1 (IFN-β1) was strongly up-regulated while expression of IL-10 and TNF-α was up-regulated slightly . A further study concerning the effect of the VR-2332 PRRSV strain (NA PRRSV) on PAM function utilized serial analysis of gene expression and demonstrated that expression of MX1 and USP were significantly up-regulated 24 hours post-infection (h pi) . These studies have provided a genome-wide gene expression profile of PAMs in vitro following infection with EU PRRSV or NA PRRSV. However, in vitro studies have significant limitations owing to disparities between the in vitro and in vivo environments. Therefore, characterization of host immune responses to PRRSV in vivo is required. PRRSV infection causes widespread apoptosis in pulmonary and lymphoid tissues of infected pigs , but the cause of the increased severity of the symptoms and the unusually high mortality of pigs infected with H-PRRSV (a mutant North American highly virulent porcine reproductive and respiratory syndrome virus) remains unknown.
High-throughput sequencing technology has been adapted for transcriptome analysis . The technology developed by Illumina (formerly Solexa sequencing) , also referred to as Digital Gene Expression tag profiling (DGE) , allows millions of short RNAs and differentially expressed genes to be identified in a sample without the need for prior annotations. DGE has many advantages including greater sequencing depth, detection of unknown transcripts, practical implementation of digital tags, generation of absolute rather than relative gene expression measurements, detection of high levels of differential polyadenylation, detection of low-abundance transcripts and small changes in gene expression, that make it particularly attractive for measuring mRNA expression and identifying differentially expressed genes [14, 15]. The hippocampal expression profiles of wild-type mice and δC-doublecortin-like kinase transgenic mice have been compared using Solexa sequencing technology , as have differences in gene expression between the liver and kidney . Furthermore, the Illumina Genome Analyzer II platform has been used to perform DGE analysis of the zebra fish transcriptome response to mycobacterium infection . However, DGE analysis has not been carried out on H-PRRSV infected pigs.
Herein histopathology, high-throughput deep sequencing and bioinformatics were utilized to analyze the relationship between pulmonary gene expression profiles after H-PRRSV infection and infection pathology. Comprehensive analysis of the global host response induced by H-PRRSV demonstrated that aggressive replication and dissemination of H-PRRSV resulted in an excessively vigorous immune and inflammatory response, contributing to severe tissue damage and high pathogenicity. This systems analysis could lead to a better understanding of the pathogenesis of H-PRRSV and to the identification of genetic components associated with H-PRRSV resistance/susceptibility in swine populations.
Clinical and pathological features of H-PRRSV-infected pigs
H-PRRSV-infected pigs exhibited signs of 'high fever' disease within 3 days post-infection (d pi). They developed a persistent high fever of 41.0°C-41.7°C between 3d pi and 7d pi, presenting with reddening of the skin, dyspnoea, depression, anorexia, edema of the eyelids, conjunctivitis, mild diarrhea, rough hair coats, shivering and lamping.
Quantitative PCR (QPCR) demonstrated H-PRRSV virus 4 and 7d pi in all tissues tested, namely serum, heart, liver, spleen, lung, kidney, lymph and brain (Table S1 in Additional file 1). Moreover, the H-PRRSV virus was successfully recovered from each of the eight tissues investigated in the affected pigs. Higher levels of H-PRRSV virus were detected in serum, lung, spleen and lymph than in other tissues. Uninfected negative control (C) pigs had no clinical signs of disease, and H-PRRSV pathogens and viral re-isolates were absent.
Analysis of DGE libraries
Major characteristics of DGE libraries
Low Quality Tag
Tag CopyNum = 1
CopyNum > 1
CopyNum > 5
CopyNum > 10
CopyNum > 20
CopyNum > 50
CopyNum > 100
Analysis of tag mapping
Identification of differentially expressed (DE) genes and signaling pathway analysis
To identify global transcriptional changes in H-PRRSV infected porcine lungs, a previously described method  modified properly was utilized to identify DE genes from normalized DGE data via pairwise comparisons between differential time points (H96/C, H168/C, H168/H96) during infection: 4520 genes had p values < 0.005, the false discovery rate (FDR) was < 0.01 and the estimated absolute log2-fold change was > 0.5 in at least one of the pairwise comparisons, which were declared to be differentially expressed during the course of infection (Additional file 2).
Validation of DGE data using qPCR and serum cytokine analysis
TNFα expression was elevated 2.27- to 6.29-fold in the sera of H-PRRSV-infected pigs on days 4 (p = 0.000) and 7 (p = 0.000) post infection, respectively, compared with C levels (Figure 4B). In H-PRRSV-infected animals, IFN-γ expression increased 1.2-fold by 3 d pi (p = 0.216), and 4.3-fold by 7 d pi (p = 0.000) (Figure 4C).
STC (Series Test of Cluster) and STC-GO (Series Test Cluster of Gene Ontology) analysis
In order to profile the gene expression time series and search for the most probable set of clusters generating the observed time series, the STC algorithm of gene expression dynamics was used, which took into account the dynamic nature of temporal gene expression profiles during clustering and identified the number of distinct clusters. DE genes exhibited eight types of temporal expression pattern (Figure S4 in Additional file 1; Additional file 2) with four significant cluster profiles (profile 6,0,1,7), which have significantly more genes assigned under the true ordering of time points than the average number assigned to the model profile in the permutation runs (Figure S4 in Additional file 1). One striking observation was the relative constancy in gene expression profiles of four (6, 0, 1, 7) significant cluster profiles, particularly profiles 6 and 1. The sustained host response in H-PRRSV-infected pigs indicated that critical decisions influencing the outcome of infection occur very early after infection. Gene Ontology (GO) based on biological process (BP) enrichment analyses for sets of DE genes having significant cluster profiles was performed using the two-sided Fisher's exact test (Additional file 4; Additional file 5; Figure S5 and S6 in Additional file 1). Significant GO categories that had a P-value of < 0.05 were used.
The most prominently over-represented GO terms of significant cluster profile 6 (0,1,1, those are up-regulated genes) included pathogenesis, release of cytochrome c from mitochondria, complement activation, regulation of ubiquitin-protein ligase activity, mitochondrial transport, inflammatory response, innate immune response, defense response, response to stress, oxidation reduction, proteolysis, lipid metabolic process, cell communication, signal transduction, endocytosis, apoptosis, regulation of caspase activity, and regulation of interferon biosynthesis and production (Additional file 4). These results are consistent with these genes and their associated processes having important roles in H-PRRSV replication and pathogenesis.
The most prominently over-represented GO terms of significant cluster profile 1 (0,-1,-1, those are down-regulated genes) included epithelial cell differentiation, sterol, steroid, cholesterol, lipid biosynthetic and metabolic process, actin cytoskeleton reorganization, regulation of transport, cell proliferation and adhesion, and cellular biosynthetic process (Additional file 4). These results suggest that H-PRRSV infection could inhibit epithelial cell differentiation. Impaired regulation associated with the biosynthesis and metabolism of steroids, cholesterol and lipids indicated that they could be involved in H-PRRSV pathogenesis.
Cytotoxic T-lymphocytes (CTLs) are the major cell-mediated immune response to viral infections and are MHC restricted. Clones of CTL cells recognize a specific antigen when it is presented to the T-cell receptor (TCR)/CD3 complex on the surface of the CTL by MHC I on the surface of the target cell. CT activity requires 'help' (i.e. cytokine production) from T-helper (Th) lymphocytes. TCRs of Th lymphocytes recognize specific antigens presented by MHC II molecules on antigen presenting cells (APCs). T cell activation requires TCR signals and co-stimulators. Co-receptor molecules and CAMs ensure that APCs are in contact with T cells for a substantial time, enhancing the interactions of APCs and T cells. Gene expression of TCR signals (CD3, IFNG, TNF, NF-kappa-B inhibitor alpha (NFKBIA) and IL10) and co-stimulators (inducible T-cell co-stimulator (ICOS) and CD86) were significantly up-regulated in H-PRRSV infected lungs (Figure 6). Furthermore, gene expression of co-receptor molecules (CD8A) and CAMs (CAM-1, VCAM1, CD40, CD2, ITGAL, ITGAV, ITGB2 and SELL) increased significantly. Collaborative action of TCR signals, co-stimulators, co-receptor molecules and CAMs leads to activation of Th cells. Activated Th cells produced cytokines (IFN-γ, TNF and IL10) and expressed CD40L, which bound to CD40 on APCs to activate them; activated APCs are more efficient in stimulating the differentiation of CD8+ T cells. Through recognition of peptide-class I MHC complexes by the TCR and involvement of the CD8 co-receptor, co-stimulator molecules and Th cells, naïve CD8+ T cells differentiated into functional CTLs capable of recognizing and 'killing' target cells bearing the same epitope on their MHC class Imolecules. Activated CTLs release perforin (PFR1) and granzymes to kill target cells. Gene expression for PRF1 and granzymes B, A and H were significantly up-regulated in H-PRRSV infected lungs, relative to C (Figure 6).
H-PRRSV induced up-regulation of anti-apoptotic genes in infected lungs including the BCL2-related myeloid cell leukemia sequence 1 (MCL1), Bcl-2-related protein A1 (protein BFL-1), putative inhibitor of apoptosis, adrenomedullin (ADM) and IL10, and the down-regulation of pro-apoptotic genes including p53 protein, apoptosis-inducing TAF9-like domain 1, apoptosis-related protein 1 (APR-1), secreted apoptosis-related protein 3 (SARP-3) and nucleoside diphosphate kinase homolog 5 (NDK-H 5). These actions of H-PRRSV serve to inhibit apoptosis, possibly prolonging the life-span of the cell and thereby increasing the yield of progeny virions.
The results from the present study are in agreement with previous research that demonstrated that H-PRRSV infected pigs exhibit severe clinical symptoms including persistent high fever, reddening of the skin, conjunctivitis, dyspnoea and severe diffuse pulmonary consolidation lesions [2, 22]. Histopathological examination demonstrated robust interstitial pneumonia in the lungs with thickening of alveolar septa accompanied by extensive infiltration of immune cells (Figure 1). The H-PRRSV virus replicates prolifically in the lungs, spleen and lymphoid organs. During infection an invading virus is recognized by PRRs that engage PAMPs and trigger signaling pathways within infected cells that are involved in innate immune (Figure 5) and adaptive immune (Figure 6) responses. Host immune responses are normally protective but if numerous cells are infected before immune induction, immune-mediated destruction can result in severe or fatal pathological consequences. Global profiling of transcriptional changes occurring in host lungs during H-PRRSV viral infection, analyzed using high-throughput Solexa sequencing, has provided important information regarding how H-PRRSV viruses trigger and regulate host immune responses and cause disease.
QPCR assays demonstrated that the H-PRRSV virus replicated rapidly and persisted in infected cells (Table S1 in Additional file 1). Substantial viral antigen was detected in alveolar cells and bronchiolar epithelial cells (Figure 1). The ability of a virus to induce and sustain an infection depends partly on its ability to block host innate immune responses or to modulate the activity of antiviral effector proteins. Production of type I IFN (IFN-α/β, SPI IFN) is an innate antiviral immune reaction in virus-infected cells that prevents viral replication and restricts the spread of the virus to neighboring cells. However, the present study demonstrated that H-PRRSV infection suppressed production of SPI IFN and down-regulated expression of IFN-α (Figure 5E). Previous in vitro and in vivo studies [8, 23, 24] have demonstrated that PRRSV infection results in minimal IFN-α production or suppresses its production, and IFN-α has been shown to inhibit PRRSV replication. During H-RRRSV viral infection, blocking SPI IFN production and particularly production of IFN-α could result in rapid spread of the virus and a high rate of viral replication. Other viral infections including the 1918 influenza virus , hepatitis C virus (HCV)  and Ebola virus  suppress type I IFN gene expression, leading to extensive viral replication and increased pathogenesis. IRF3 plays an important role in typeI IFN gene expression and the present study demonstrated that IRF3 gene expression was suppressed during H-PRRSV infection (Figure 5D). This result is in agreement with a previous study reporting that PRRSV NSP1β inhibited IRF3 and NF-κB transactivation, and down-regulated IFN-β gene expression. This suggested that NSP1β mediates subversion of the host innate immune response and plays an important role in PRRSV pathogenesis. Furthermore, influenza A NSP1 can suppress innate immunity by inhibiting activation of IRF3, and subsequently disrupting the induction of α/β-interferon .
Many viruses induce apoptosis in infected cells but some can block the apoptosis pathway, leading to prolonged life of the cell and an increase in the yield of progeny virions. H-PRRSV up-regulated expression of anti-apoptotic genes and down-regulated expression of some pro-apoptotic genes in H-PRRSV infected lungs (Figure 7C). MCL1, BFL-1, putative inhibitor of apoptosis, ADM and IL10 were up-regulated. MCL1 and BFL-1 belong to the BCL-2 subfamily, which negatively regulates apoptosis and blocks the apoptosis pathway; ADM is an anti-apoptotic peptide ; and IL10 protects cells against apoptosis . The pro-apoptotic genes APR-1, p53 protein, SARP-3, and NDK-H 5 were down-regulated to prevent the occurrence of apoptosis. These findings indicate that H-PRRSV could induce an anti-apoptotic state to prolong the life-span of infected cells and increase the yield of progeny virions.
IL10 could have an important role in the regulation of the immune response to PRRSV. Up-regulation of IL10 gene expression has been demonstrated in PRRSV-infected porcine leukocytes, alveolar macrophages, dendritic cells, and in vivo in PRRSV infected pigs [8, 31]. Incubation of freshly isolated CD14 positive cells with IL10 during differentiation increased susceptibility to PRRSV infection and was correlated with up-regulation of CD163 on the cell surface . This suggests that IL10 plays an important role in CD163 up-regulation and susceptibility to PRRSV during differentiation of macrophages in vivo. CD163 alone can confer PRRSV replication on a non-permissive pig cell line and its expression on macrophages in vivo could determine the efficiency of replication and subsequent pathogenicity of PRRSV . It is possible that internalization of H-PRRSV via CD163 on the target cells could induce expression of IL10 and subsequently induce the expression of CD163 on neighboring undifferentiated monocytes, increasing overall susceptibility to PRRSV.
Taken together, the above findings suggested that the H-PRRSV virus aggressively replicated and disseminated by subverting the host innate immune response, inducing an anti-apoptotic state and up-regulating expression of CD163.
Prolific replication and rapid spread of H-PRRSV virus caused severe lung damage, hemorrhage and extensive infiltration of immune cells throughout the course of infection. Accordingly, significant increases in the expression of a number of genes involved in phagocytic cell activation were observed including CAMs, and several pro-inflammatory cytokines and chemokines such as IFN-γ, TNF, SELL, ICAM, integrin, C-type lectin, IL2RG, IL8, CSF2, IRG6, macrophage inflammatory protein 3 (MIP-3), CXCL2, CXCL9, CXCL10, CCL2 and CCR5 (Figure 5H, J, K). Up-regulated expression of these genes resulted in recruitment of neutrophils, macrophages and other immune cells to sites of infection, and excessive infiltration resulted in destruction of tissues . Moreover, H-PRRSV infection resulted in the activation of CD4 and CD8 T lymphocytes specific for H-PRRSV antigens, and these secreted vasoactive cytokines including TNFα and IFN-γ. This cytokine 'storm' increased capillary fragility (with associated hemorrhages) and permeability. H-PRRSV infection activated complement proteins, which enhanced vascular permeability and were associated with sequestration of thrombocytes. The sustained induction of pro-inflammatory cytokines and chemokines contributed to a robust inflammatory response in the lung.
Fever is frequently the initial response to infection and it is triggered by PRR-PAMP interactions that activate a signaling cascade that causes the production of inflammatory cytokines responsible for fever including CASP1, the IL1-converting enzyme responsible for cleaving the IL-1β precursor and resulting in production of the mature form . TLR2, 4, 6, 7, 9 (Figure 5A) and CASP1 (Figure 7E) were significantly up-regulated in H-PRRSV infected lungs. Heat shock proteins, referred to as stress proteins, are induced in cells exposed to a wide range of environmental stressors including infection and extreme temperature. Gene expression levels of heat shock genes including HSPA5, HSP27, HSP90, HSP90B1, HSPCB and HSPD1 were significantly elevated in H-PRRSV infected lungs relative to C (Figure S7 in Additional file 1).
During H-RRRSV virus infection, activated CTLs and NK cells release perforin and granzymes to kill target cells. Gene expression of PRF1 and granzyme B, A and H were significantly up-regulated in H-PRRSV infected lungs (Figure 6H). Perforin is exocytosed and polymerizes in the target cell plasma membrane to form pores. Granzymes enter target cells through the perforin pores and induce target cell apoptosis. The perforin pores also allow the release of intracellular calcium from the target cell, which acts to trigger apoptotic pathways. The induction of a CTL response results in the release of various cytokines from Th cells, some of which result in clonal proliferation of antigen-specific CTLs, and others that have direct antiviral effects. Diffusion of perforin and local cytokine production frequently results in inflammation and bystander cell damage.
There was up-regulation of genes involved in TNF signaling including TNFα and TNFRSF1A (Figure 7A) and the up-regulation of TNF during PRRSV infection has been reported to have an important role in pathogenicity. It has been reported that PRRSV infection is a potent inducer of TNFα in PAMs . In the present study, continuously up-regulated expression of TNFα (at mRNA and protein levels) from 96 h pi to 168 h pi was observed (Figure 4B). Interestingly, infection with H-PRRSV led to up-regulation of NFKBIA, an inhibitor of the TNF receptor activated transcription factor NF-κB. Loss of NF-κB activity has been reported to increase the cytotoxic effects of TNF and result in increased cell death . TNF and NFKBIA could act synergistically to cause significant alveolar and bronchial epithelial cell necrosis during H-PRRSV infections.
This study has indicated that H-PRRSV could induce apoptosis through a mitochondria-mediated pathway, and previous research provided evidence that PRRSV induces apoptosis in MARC-145 cells through an intrinsic mitochondria-mediated pathway . Pro-apoptotic genes (BAX, BAK, BID, PIK3C3), cytochrome C, and caspases (CASP-10, CASP1, CASP4, CASP15, CASP3) were up-regulated (Figure 7). These results indicate that up-regulation of pro-apoptotic genes resulted in disruption of the mitochondrial transmembrane potential, thereby inducing release of cytochrome c, AIF-like mitochondrion-associated inducer of death and CASP3 from mitochondrial membranes, leading to induction of apoptosis and secondary necrosis. The release of cytochrome c can also induce necrosis through a slower non-apoptotic mechanism due to the electrochemical gradient across the inner membrane, production of reactive oxygen species (ROS) and declining ATP production . The production of ROS, particularly superoxide radicals, and the subsequent oxidative damage to cells and tissues, are recognized as key contributors to viral pathogenesis [36, 39]. ROS-mediated oxidative stress could also contribute to PRRSV-induced apoptosis . In the current study, continuous up-regulated expression of cytochrome b245 heavy chain (GP91-PHOX), a critical component of the membrane-bound oxidase of phagocytes (macrophages and neutrophils) that generates superoxide radicals, was observed from 96 h pi to 168 h pi (Figure 7F). Increased expression of cytochrome b245 in H-PRRSV infected lungs implies the increased production of oxygen radicals and the activation of phagocytic cells. Taken together, these data suggested that the severe pulmonary pathology caused by H-PRRSV infection was induced by significant production of TNF, PRF1, granzymes, cytochrome c and oxygen radicals.
Experimental animals and tissue collection
All animal procedures were performed according to guidelines developed by the China Council on Animal Care and protocols were approved by the Animal Care and Use Committee of Guangdong Province, China. Nine conventionally-reared healthy 6-week-old crossbred weaned pigs (Landrace × Yorkshire) were selected from a high-health commercial farm that has historically been free from all major pig diseases including PRRSV, porcine circovirus type 2, classical swine fever virus, porcine parvovirus, pseudorabies virus, swine influenza virus and Mycoplasma hyopneumoniae infections. All pigs were PRRSV-seronegative as determined by ELISA (HerdChek PRRS 2XR; IDEXX Laboratories) and absence of PRRSV was confirmed by RT-PCR. Pigs were randomly assigned to three groups and raised in isolation rooms. Six pigs were inoculated with 6 ml viral suspension (4 ml intranasally and 2 ml intramuscularly) of H-PRRSV (H-PRRSV gift from Dr. Zhang Guihong, South China Agricultural University) at a dose of 106.0 TCID50 ml-1 on day 0. Three C pigs were treated with an identical volume of DMEM culture medium from uninfected MARC-145 cells 1 day prior to experimental infection, and were immediately necropsied. H-PRRSV-inoculated pigs were clinically examined daily and rectal body temperatures were recorded from day 2 to 7 pi. Viral re-isolates and H-PRRSV specific RT-PCR were performed after the pigs were sacrificed. Three infected pigs randomly chosen were necropsied at 96 h pi and 168 h pi. Lung samples were collected from C, H96 and H168, and frozen in liquid nitrogen for RNA isolation or fixed in 10% neutralized buffered formalin for histological processing.
Virus re-isolation and QPCR detection
Heart, liver, spleen, lung, kidney, lymph and brain tissue were collected at autopsy. Samples were placed in sterilized PBS buffer, homogenized and centrifuged to harvest supernatants for virus re-isolation and detection of H-PRRSV using QPCR. Tissue homogenates were inoculated on monolayers of MARC-145 cells and maintained in Dulbecco's modified Eagle's medium (DMEM) supplemented with 75 μg of penicillin and 50 μg of streptomycin per ml. Cytopathic effects (CPE) were monitored daily. Cultures not displaying CPE after three passages were considered negative. Tissue homogenates were examined by H-PRRSV-specific QPCR. The oligonucleotide primers used were NSP2F2 (5'-GTGGGTCGGCACCAGTT-3') and NSP2R2 (5'- GACGCAGACAAATCCAGAGG-3'), designed from the gene segment encoding for NSP2 with a deletion of 87 bases in the fixed site as compared to the NA PRRSV. The TaqMan probe, 5'FAM-CGCGTAGAACTGTGACAACAACGCTGA-TAMRA3', was synthesized according to previous protocols .
Lungs of C and experimentally infected pigs were processed for haematoxylin and eosin (H&E) and immunohistochemistry staining, as described previously .
RNA extraction, library construction and sequencing
Total RNA was extracted from frozen lungs using standard protocols (Trizol) and treated with DNase to remove potential genomic DNA contamination, according to the manufacturer's protocol. RNA integrity and concentration were evaluated using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA).
For RNA library construction and deep sequencing, RNA samples were prepared as follows: for each time point equal quantities of RNA isolated from three individual lungs were pooled from the H-PRRSV-inoculated group and the C group. A 6 μg sample of RNA from each group was submitted to Solexa (now Illumina Inc.) for sequencing.
Sequence tag preparation was carried out using Illumina's Digital Gene Expression Tag Profiling Kit according to the manufacturer's protocol. In brief, mRNA was isolated from 6 μg of total RNA by binding the mRNA to a magnetic oligo bead. First- and second-strand cDNA were synthesized while the mRNA was attached to the beads. Double stranded cDNA was digested with NlaIII to remove all fragments other than the 3' most CATG fragment attached to the oligobead. GEX NlaIII Adapter 1 was ligated at the site of NlaIII cleavage. GEX NlaIII Adapter 1 contains the sequence for the restriction enzyme MmeI, and the restriction enzyme MmeI was used to create the 17 bp tag. The GEX Adapter 2 was ligated at the site of MmeI cleavage. A 12 cycle PCR was performed with two primers that anneal to the ends of the adapters to enrich the adapter-ligated cDNA construct. The amplified cDNA construct was purified from a 6% Novex TBE PAGE gel. The purified cDNA tags were sequenced on the Illumina Cluster Station and Genome Analyzer. Image recognition and base calling were performed using the Illumina Pipeline.
Analysis of sequencing data
For the raw data adaptor tags, low quality tags and tags of copy number = 1 were filtered to produce clean tags. The raw data (tag sequences and counts) have been submitted to Gene Expression Omnibus (GEO) under series GSE19456. The clean tags were classified according their copy number in the library and their percentages in the total clean tags were provided. Saturation of the library was also analyzed.
The pre-processed database of all possible CATG+17-nt tag sequences was created using sus scrofa UniGene (http://www.ncbi.nlm.nih.gov/UniGene/UGOrg.cgi?TAXID=9823, UniGene Build #35 Sus scrofa, Nov, 7th, 2008) from NCBI. Clean tags were aligned to the reference sequences, and unambiguous tags were annotated. The clean tag number corresponding to each gene was counted.
Differential expression detection
To compare the differential expression of genes across samples (H96/C, H168/C, H168/H96), the number of raw clean tags in each sample was normalized to Tags Per Million (TPM) to obtain normalized gene expression levels. Differential expression detection of genes or tags across samples was performed . Genes were classed as significantly differentially expressed if they had a P-value < 0.005, a false discovery rate (FDR) < 0.01 and an estimated absolute log2-fold change > 0.5 in sequence counts across libraries.
qPCR and serum cytokine analysis
The RNA samples used for the qPCR assays were the same as those used for the DGE experiments and independent RNA extractions from biological replicates. qPCR was carried out on the Lightcycler480 (Roche) with SYBR-Green detection (SYBR PrimeScript RT-PCR Kit, TaKaRa Biotechnology Co., Ltd.), according to the manufacturer's instructions. Each cDNA was analyzed in triplicate, and the average threshold cycle was calculated. Relative expression levels were calculated using the 2-ΔΔCt method. The results were normalized to the expression level of HPRT1 and relative to the C sample. Levels of cytokines (TNFα and IFN-γ) from serum were assayed using swine commercial ELISA kits from R&D Systems according to the manufacturer's instructions.
STC and STC-GO analysis
STC is implemented entirely in java. The clustering algorithm selects a set of distinct and representative temporal expression profiles. These model profiles are selected independently of the data. The clustering algorithm assigns each gene passing the filtering criteria to the model profile that most closely matches the gene's expression profile as determined by the correlation coefficient. Since the model profiles are selected independently of the data, the algorithm can determine which profiles have a statistically significant higher number of genes assigned using a permutation test. This test determines an assignment of genes to model profiles using a large number of permutations of the time points. It uses standard hypothesis testing to determine which model profiles have significantly more genes assigned under the true ordering of time points compared to the average number assigned to the model profile in the permutation runs. Significant model profiles can be either analyzed independently or grouped together on the basis of similarity to form clusters of significant profiles [42, 43].
STC-GO supports Gene Ontology enrichment analyses for sets of genes having the same significant temporal expression pattern. Random samples of S a (S a is the number of genes assigned to the same model temporal expression profile r) were selected and genes at each iteration and Fisher's exact test p-values for the selected genes in all GO biological categories were calculated . The two-sided Fisher's exact test p-value for a category reflects a test of the null hypothesis that the category is enriched in genes assigned to profile r with respect to what would have been expected by chance alone. To decide whether to investigate a category that appears enriched in these genes further, the statistical reliability of the apparent enrichment would be calculated. To assess the significance of a particular category, the distribution of p-values that would occur by random chance must be known. The percentage of false positives to be tolerated will generally depend on the relative costs of false positives and false negatives in whatever follow-up study is to be carried out. This way of framing the question leads us to specify the false discovery rate (FDR) for a set of categories, rather than the significance level (p-value) for each category. With the significance at the 0.05 level for a given category, the enrichment Re is given by Re = (i/m)/(Sa/N) where i is the number of genes assigned to profile r within the GO category of interest, m is the total number of genes within the GO category of interest, and N is total number of unique genes in the gene reference database list.
Pathway analysis was predominantly based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The two-side Fisher's exact test with multiple testing and the χ2 test were used to classify the pathway category. The false discovery rate (FDR) was used to correct the P-value. Only pathway categories that had a P < 0.05 were chosen. Within the significant category, the enrichment Re was given by: (Re = ENRICHMENT), where n f is the number of flagged proteins within the particular category, n is the total number of proteins within the same category, N f is the number of flagged proteins in the protein reference database list, and N is the total number of proteins in the gene reference database list.
highly virulent porcine reproductive and respiratory syndrome virus
Digital Gene Expression tag profiling
We thank the Beijing Genomics Institute (BGI) Shenzhen and Genminix Informatics Co. Ltd. for providing technical assistance in DGE and bioinformatics analysis. This research was supported by the National Natural Science Foundation of China (Grant No. U0731003), National Key Basic Research Plan (973 Project) (Grant No. 2006CB102101) and the Ministry of National Agricultural Special Public Sector Research (Grant No. nyhyzx07-034).
- Tian K, Yu X, Zhao T, Feng Y, Cao Z, Wang C, Hu Y, Chen X, Hu D, Tian X: Emergence of fatal PRRSV variants: unparalleled outbreaks of atypical PRRS in China and molecular dissection of the unique hallmark. PLoS ONE. 2007, 2: e526-10.1371/journal.pone.0000526.PubMed CentralPubMedView ArticleGoogle Scholar
- Feng Y, Zhao T, Nguyen T, Inui K, Ma Y, Nguyen TH, Nguyen VC, Liu D, Bui QA, To LT: Porcine respiratory and reproductive syndrome virus variants, Vietnam and China, 2007. Emerg Infect Dis. 2008, 14: 1774-1776. 10.3201/eid1411.071676.PubMed CentralPubMedView ArticleGoogle Scholar
- Normile D: Virology. China, Vietnam grapple with 'rapidly evolving' pig virus. Science. 2007, 317: 1017-10.1126/science.317.5841.1017.PubMedView ArticleGoogle Scholar
- Choi C, Cho WS, Kim B, Chae C: Expression of Interferon-gamma and tumour necrosis factor-alpha in pigs experimentally infected with Porcine Reproductive and Respiratory Syndrome Virus (PRRSV). J Comp Pathol. 2002, 127: 106-113. 10.1053/jcpa.2002.0566.PubMedView ArticleGoogle Scholar
- Zhang X, Shin J, Molitor TW, Schook LB, Rutherford MS: Molecular responses of macrophages to porcine reproductive and respiratory syndrome virus infection. Virology. 1999, 262: 152-162. 10.1006/viro.1999.9914.PubMedView ArticleGoogle Scholar
- Liu CH, Chaung HC, Chang HL, Peng YT, Chung WB: Expression of Toll-like receptor mRNA and cytokines in pigs infected with porcine reproductive and respiratory syndrome virus. Vet Microbiol. 2009, 136: 266-276. 10.1016/j.vetmic.2008.11.016.PubMedView ArticleGoogle Scholar
- Kimman TG, Cornelissen LA, Moormann RJ, Rebel JM, Stockhofe-Zurwieden N: Challenges for porcine reproductive and respiratory syndrome virus (PRRSV) vaccinology. Vaccine. 2009, 27: 3704-3718. 10.1016/j.vaccine.2009.04.022.PubMedView ArticleGoogle Scholar
- Genini S, Delputte PL, Malinverni R, Cecere M, Stella A, Nauwynck HJ, Giuffra E: Genome-wide transcriptional response of primary alveolar macrophages following infection with porcine reproductive and respiratory syndrome virus. J Gen Virol. 2008, 89: 2550-2564. 10.1099/vir.0.2008/003244-0.PubMed CentralPubMedView ArticleGoogle Scholar
- Miller LC, Harhay GP, Lager KM, Smith TP, Neill JD: Effect of porcine reproductive and respiratory syndrome virus on porcine alveolar macrophage function as determined using serial analysis of gene expression (SAGE). Dev Biol (Basel). 2008, 132: 169-174.Google Scholar
- Labarque G, Van Gucht S, Nauwynck H, Van Reeth K, Pensaert M: Apoptosis in the lungs of pigs infected with porcine reproductive and respiratory syndrome virus and associations with the production of apoptogenic cytokines. Vet Res. 2003, 34: 249-260. 10.1051/vetres:2003001.PubMedView ArticleGoogle Scholar
- Velculescu VE, Kinzler KW: Gene expression analysis goes digital. Nat Biotechnol. 2007, 25: 878-880. 10.1038/nbt0807-878.PubMedView ArticleGoogle Scholar
- Bennett ST, Barnes C, Cox A, Davies L, Brown C: Toward the 1,000 dollars human genome. Pharmacogenomics. 2005, 6: 373-382. 10.1517/146224220.127.116.113.PubMedView ArticleGoogle Scholar
- Saliki JT, Rodgers SJ, Eskew G: Serosurvey of selected viral and bacterial diseases in wild swine from Oklahoma. J Wildl Dis. 1998, 34: 834-838.PubMedView ArticleGoogle Scholar
- Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008, 18: 1509-1517. 10.1101/gr.079558.108.PubMed CentralPubMedView ArticleGoogle Scholar
- t Hoen PA, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, de Menezes RX, Boer JM, van Ommen GJ, den Dunnen JT: Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36: e141-10.1093/nar/gkn705.View ArticleGoogle Scholar
- Hegedus Z, Zakrzewska A, Agoston VC, Ordas A, Racz P, Mink M, Spaink HP, Meijer AH: Deep sequencing of the zebrafish transcriptome response to mycobacterium infection. Mol Immunol. 2009, 46: 2918-2930. 10.1016/j.molimm.2009.07.002.PubMedView ArticleGoogle Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
- Katze MG, Fornek JL, Palermo RE, Walters KA, Korth MJ: Innate immune modulation by RNA viruses: emerging insights from functional genomics. Nat Rev Immunol. 2008, 8: 644-654. 10.1038/nri2377.PubMedView ArticleGoogle Scholar
- Van Gorp H, Van Breedam W, Delputte PL, Nauwynck HJ: Sialoadhesin and CD163 join forces during entry of the porcine reproductive and respiratory syndrome virus. J Gen Virol. 2008, 89: 2943-2953. 10.1099/vir.0.2008/005009-0.PubMedView ArticleGoogle Scholar
- Lau DT, Fish PM, Sinha M, Owen DM, Lemon SM, Gale M: Interferon regulatory factor-3 activation, hepatic interferon-stimulated gene expression, and immune cell infiltration in hepatitis C virus patients. Hepatology. 2008, 47: 799-809. 10.1002/hep.22076.PubMedView ArticleGoogle Scholar
- Tsichlis PN, Donehower L, Hager G, Zeller N, Malavarca R, Astrin S, Skalka AM: Sequence comparison in the crossover region of an oncogenic avian retrovirus recombinant and its nononcogenic parent: genetic regions that control growth rate and oncogenic potential. Mol Cell Biol. 1982, 2: 1331-1338.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhou YJ, Hao XF, Tian ZJ, Tong GZ, Yoo D, An TQ, Zhou T, Li GX, Qiu HJ, Wei TC, Yuan XF: Highly virulent porcine reproductive and respiratory syndrome virus emerged in China. Transbound Emerg Dis. 2008, 55: 152-164. 10.1111/j.1865-1682.2008.01020.x.PubMedView ArticleGoogle Scholar
- Buddaert W, Van Reeth K, Pensaert M: In vivo and in vitro interferon (IFN) studies with the porcine reproductive and respiratory syndrome virus (PRRSV). Adv Exp Med Biol. 1998, 440: 461-467.PubMedView ArticleGoogle Scholar
- Van Reeth K, Labarque G, Nauwynck H, Pensaert M: Differential production of proinflammatory cytokines in the pig lung during different respiratory virus infections: correlations with pathogenicity. Res Vet Sci. 1999, 67: 47-52. 10.1053/rvsc.1998.0277.PubMedView ArticleGoogle Scholar
- Kobasa D, Jones SM, Shinya K, Kash JC, Copps J, Ebihara H, Hatta Y, Kim JH, Halfmann P, Hatta M: Aberrant innate immune response in lethal infection of macaques with the 1918 influenza virus. Nature. 2007, 445: 319-323. 10.1038/nature05495.PubMedView ArticleGoogle Scholar
- Kash JC, Muhlberger E, Carter V, Grosch M, Perwitasari O, Proll SC, Thomas MJ, Weber F, Klenk HD, Katze MG: Global suppression of the host antiviral response by Ebola- and Marburgviruses: increased antagonism of the type I interferon response is associated with enhanced virulence. J Virol. 2006, 80: 3009-3020. 10.1128/JVI.80.6.3009-3020.2006.PubMed CentralPubMedView ArticleGoogle Scholar
- Beura LK, Sarkar SN, Kwon B, Subramaniam S, Jones C, Pattnaik AK, Osorio FA: Porcine reproductive and respiratory syndrome virus nonstructural protein 1beta modulates host innate immune response by antagonizing IRF3 activation. J Virol. 2010, 84: 1574-1584. 10.1128/JVI.01326-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Mibayashi M, Martinez-Sobrido L, Loo YM, Cardenas WB, Gale M, Garcia-Sastre A: Inhibition of retinoic acid-inducible gene I-mediated induction of beta interferon by the NS1 protein of influenza A virus. J Virol. 2007, 81: 514-524. 10.1128/JVI.01265-06.PubMed CentralPubMedView ArticleGoogle Scholar
- Uzan B, Ea HK, Launay JM, Garel JM, Champy R, Cressent M, Liote F: A critical role for adrenomedullin-calcitonin receptor-like receptor in regulating rheumatoid fibroblast-like synoviocyte apoptosis. J Immunol. 2006, 176: 5548-5558.PubMedView ArticleGoogle Scholar
- Zhou JH, Broussard SR, Strle K, Freund GG, Johnson RW, Dantzer R, Kelley KW: IL-10 inhibits apoptosis of promyeloid cells by activating insulin receptor substrate-2 and phosphatidylinositol 3'-kinase. J Immunol. 2001, 167: 4436-4442.PubMedView ArticleGoogle Scholar
- Flores-Mendoza L, Silva-Campa E, Resendiz M, Osorio FA, Hernandez J: Porcine reproductive and respiratory syndrome virus infects mature porcine dendritic cells and up-regulates interleukin-10 production. Clin Vaccine Immunol. 2008, 15: 720-725. 10.1128/CVI.00224-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Patton JB, Rowland RR, Yoo D, Chang KO: Modulation of CD163 receptor expression and replication of porcine reproductive and respiratory syndrome virus in porcine macrophages. Virus Res. 2009, 140: 161-171. 10.1016/j.virusres.2008.12.002.PubMedView ArticleGoogle Scholar
- Wang T, Town T, Alexopoulou L, Anderson JF, Fikrig E, Flavell RA: Toll-like receptor 3 mediates West Nile virus entry into the brain causing lethal encephalitis. Nat Med. 2004, 10: 1366-1373. 10.1038/nm1140.PubMedView ArticleGoogle Scholar
- Yu HB, Finlay BB: The caspase-1 inflammasome: a pilot of innate immune responses. Cell Host Microbe. 2008, 4: 198-208. 10.1016/j.chom.2008.08.007.PubMedView ArticleGoogle Scholar
- Chang HW, Jeng CR, Liu JJ, Lin TL, Chang CC, Chia MY, Tsai YC, Pang VF: Reduction of porcine reproductive and respiratory syndrome virus (PRRSV) infection in swine alveolar macrophages by porcine circovirus 2 (PCV2)-induced interferon-alpha. Vet Microbiol. 2005, 108: 167-177. 10.1016/j.vetmic.2005.03.010.PubMedView ArticleGoogle Scholar
- Kash JC, Basler CF, Garcia-Sastre A, Carter V, Billharz R, Swayne DE, Przygodzki RM, Taubenberger JK, Katze MG, Tumpey TM: Global host immune response: pathogenesis and transcriptional profiling of type A influenza viruses expressing the hemagglutinin and neuraminidase genes from the 1918 pandemic virus. J Virol. 2004, 78: 9499-9511. 10.1128/JVI.78.17.9499-9511.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee SM, Kleiboeker SB: Porcine reproductive and respiratory syndrome virus induces apoptosis through a mitochondria-mediated pathway. Virology. 2007, 365: 419-434. 10.1016/j.virol.2007.04.001.PubMedView ArticleGoogle Scholar
- Green DR, Reed JC: Mitochondria and apoptosis. Science. 1998, 281: 1309-1312. 10.1126/science.281.5381.1309.PubMedView ArticleGoogle Scholar
- Akaike T, Noguchi Y, Ijiri S, Setoguchi K, Suga M, Zheng YM, Dietzschold B, Maeda H: Pathogenesis of influenza virus-induced pneumonia: involvement of both nitric oxide and oxygen radicals. Proc Natl Acad Sci USA. 1996, 93: 2448-2453. 10.1073/pnas.93.6.2448.PubMed CentralPubMedView ArticleGoogle Scholar
- Xiao XL, Wu H, Yu YG, Cheng BZ, Yang XQ, Chen G, Liu DM, Li XF: Rapid detection of a highly virulent Chinese-type isolate of Porcine Reproductive and Respiratory Syndrome Virus by real-time reverse transcriptase PCR. J Virol Methods. 2008, 149: 49-55. 10.1016/j.jviromet.2008.01.009.PubMedView ArticleGoogle Scholar
- Qiu HJ, Tian ZJ, Tong GZ, Zhou YJ, Ni JQ, Luo YZ, Cai XH: Protective immunity induced by a recombinant pseudorabies virus expressing the GP5 of porcine reproductive and respiratory syndrome virus in piglets. Vet Immunol Immunopathol. 2005, 106: 309-319. 10.1016/j.vetimm.2005.03.008.PubMedView ArticleGoogle Scholar
- Ramoni MF, Sebastiani P, Kohane IS: Cluster analysis of gene expression dynamics. Proc Natl Acad Sci USA. 2002, 99: 9121-9126. 10.1073/pnas.132656399.PubMed CentralPubMedView ArticleGoogle Scholar
- Miller LD, Long PM, Wong L, Mukherjee S, McShane LM, Liu ET: Optimal gene expression analysis by microarrays. Cancer Cell. 2002, 2: 353-361. 10.1016/S1535-6108(02)00181-2.PubMedView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralPubMedView ArticleGoogle Scholar
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 cited.