- Research article
- Open Access
Transcription analysis of the porcine alveolar macrophage response to porcine circovirus type 2
BMC Genomics volume 14, Article number: 353 (2013)
Porcine circovirus type 2 (PCV2) is the causal agent of postweaning multisystemic wasting syndrome (PMWS), which has severely impacted the swine industry worldwide. PCV2 triggers a weak and atypical innate immune response, but the key genes and mechanisms by which the virus interferes with host innate immunity have not yet been elucidated. In this study, genes that control the response of primary porcine alveolar macrophages (PAMs), the main target of PCV2, were profiled in vitro.
PAMs were successfully infected by PCV2-WH strain, as evidenced quantitative real-time polymerase chain reaction (qPCR) and immunofluorescence assay (IFA) results. Infection-related differential gene expression was investigated using pig microarrays from the US Pig Genome Coordination Program and validated by real-time PCR and enzyme-linked immunosorbent assay (ELISA). Microarray analysis at 24 and 48 hours post-infection (HPI) revealed 266 and 175 unique genes, respectively, that were differentially expressed (false discovery rate <0.05; fold-change >2). Only six genes were differentially expressed between 24 and 48 HPI. The up-regulated genes were principally related to immune response, cytokine activity, locomotion, regulation of cell proliferation, apoptosis, cell growth arrest, and antigen procession and presentation. The down-regulated genes were mainly involved in terpenoid biosynthesis, carbohydrate metabolism, translation, proteasome degradation, signal transducer activity, and ribosomal proteins, which were representative of the reduced vital activity of PCV2-infected cells.
PCV2 infection of PAMs causes up-regulation of genes related to inflammation, indicating that PCV2 may induce systematic inflammation. PCV2 persistently induced cytokines, mainly through the Toll-like receptor (TLR) 1 and TLR9 pathways, which may promote high levels of cytokine secretion. PCV2 may prevent apoptosis in PAMs by up-regulating SERPINB9 expression, possibly to lengthen the duration of PCV2 replication-permissive conditions. The observed gene expression profile may provide insights into the underlying immunological response and pathological changes that occur in pigs following PCV2 infection.
Porcine circovirus 2 (PCV2) is a small, non-enveloped, single-stranded and closed-circular DNA virus that has been identified as the primary cause of postweaning multisystemic wasting syndrome (PMWS) in pigs. PMWS primarily affects pigs between the ages of five and 18 weeks. The clinical signs include progressive weight loss, jaundice, wasting, and respiratory disease, accompanied by increased mortality. The mortality rates of PMWS vary from 1-2%, but may reach up to 30% in complicated cases. PMWS has become endemic in many swine-producing countries, including China, Canada, the United States, Venezuela, Australia, and the United Kingdom [1–3].
DNA microarray technology, in combination with bioinformatics, has emerged as an efficient high-throughput tool that offers great advantages in the study of the genomic expression profiles of cells and host-microbe interactions . This technology provides the ability to determine gene expression levels of thousands of different genes simultaneously. Indeed, high density gene arrays have already revealed the status of host gene expression following infection by several viruses and bacteria, including influenza A virus , monkeypox virus , Epstein-Barr virus , Streptococcus suis, and Haemophilus parasuis.
The major target cells of PCV2 are the porcine monocyte/macrophage lineage cells (MLCs) [10–13]. One type of MLC, porcine alveolar macrophages (PAMs), are the primary responders of the pulmonary innate immune system and defend against various pathogens . While PAMs are known to be involved in PCV2 defense, the underlying molecular mechanisms of PCV2 pathogenicity remain to be completely elucidated. Furthermore, an immunosuppressive role has been characterized for PCV2 in porcine circovirus diseases (PCVD) [15, 16], but little is known about the genes that are involved in this process. The presence of granulomatous inflammation and multinucleated giant cells in PMWS pigs  suggests a significant role for secreted cytokines and other immune-related signaling factors. Previous studies have investigated the cytokine and chemokine gene expression profile of PMWS by transcriptomic and genomic analyses. For example, the cytokine mRNA expression profile in lymph nodes of PMWS pigs was reported , as was the global transcriptome profile of colostrums-deprived piglets experimentally infected with PCV2 [18, 19]. The genomic expression profile in lymph nodes of lymphoid-depleted PCV2-infected pigs was also reported . However, these results can not fully illustrate the pathogenesis of PCV2.
Therefore, the objective of the present study was to analyze the PCV2-induced gene expression profiles in isolated PAMs in vitro in order to reveal specific mechanisms of the host innate immune response to PCV2. The temporal gene expression profiling using a two-channel microarray carrying thousands of genes was carried out with PAMs at 24 and 48 hours post-infection (HPI) to discriminate PCV2-induced early and late changes in gene expression.
Confirmation of PCV2 infection in PAMs
Immunofluorescence assay (IFA) revealed green fluorescent cells in the infected PAMs, indicating that the PAMs were successfully infected by PCV2 at 24 and 48 hours after inoculation (Figure 1). There were no significant differences in the levels of viral DNA or virus particles between the two time points (mean viral loads: 5.67±0.35 log10 and 5.72±0.49 log10). Neither PCV2 particles nor the PCV2 nucleotide signal were detected in the mock-inoculated PAMs at any of the experimental time points.
Transcriptome analysis of the PAM response to PCV2 infection
The concentrations of RNA ranged from 0.67 to 1.97 μg/μL. All RNA samples were considered to be of high quality, based on RNA integrity number (RIN) values (ranging from 8.0 to 9.4). The isolated mRNA samples from 24 HPI and 48 HPI PAMs hybridized to 6345 and 5856 probe sets on the transcriptome microarray, respectively corresponding to 34.25% and 31.65% of all probe sets (Additional file 1: Table S1). After quantile normalization and statistical analysis, 1277 transcripts were identified (false discovery rate (FDR)-adjusted p-value <0.05). Furthermore, the criteria of a two-fold or greater change in differential expression and a FDR of 5% were chosen to determine up-regulated and down-regulated genes. According to these criteria, 266 and 175 genes were differentially expressed (DE) from baseline at 24 HPI and 48 HPI (Additional file 2: Table S2 and Additional file 3: Table S3). According to the corresponding biological process Gene Ontology (GO) terms, or potential biological processes reported in recent publications, the majority of these genes were related to Immune Response, Inflammatory Response, Antigen Processing and Presentation of Peptide or Polysaccharide Antigen via MHC Class II, Anti-viral Immune response, Cell Motility, and Negative Regulation of Apoptosis (Table 1). However, the differential expression between 24 HPI and 48 HPI reached the threshold of statistical significance for only six of those genes (Additional file 4: Table S4 and Additional file 5: Table S5). A total of 17 differentially expressed genes were selected to validate the microarray results by quantitative reverse transcription PCR (qRT-PCR) using the primers in Table 2. The qRT-PCR results confirmed, in general, the patterns of up-regulated and down-regulated expression of the 17 genes that was indicated by the microarray analysis (Table 3).
Functional GO analysis of DE genes in PCV2-infected PAMs
The Ingenuity pathway analysis (IPA) program was used to determine the relevance of PAM gene functionalities and gene networks associated with PCV2 infection. While 266 (24 HPI) and 175 (48 HPI) DE genes were imputed into IPA, only 173 and 124 were characterized with specific cellular functions.
The 173 DE genes from 24 HPI PAMs were placed into 19 functional groups (Additional file 6: Table S6 and Additional file 7: Figure S8). The most frequently represented category of gene functionality was Inflammatory Response (Cellular Movement, and Immune Cell Trafficking) with 22 genes, followed by Cell-To-Cell Signaling and Interaction, Cellular Growth and Proliferation, and Cellular Compromise.
The 124 DE genes from 48 HPI PAMs were placed into 16 functional groups (Additional file 8: Table S7 and Additional file 9: Figure S9). The most frequently represented category of gene functionality was also Inflammatory Response with 25 genes, followed by Antigen Presentation, Cellular Movement, Lipid Metabolism, Cell Cycle, and Cellular Growth and Proliferation.
qRT-PCR validated DE genes in PCV2-infected PAMs
In order to validate the differential expression of some of the DE genes identified by microarray, eleven up-regulated genes and six down-regulated genes were selected for qRT-PCR analysis. All the selected genes were successfully amplified and the expression patterns corresponded to those observed by microarray. Although the extent of fold-change (FC) varied between the qRT-PCR and microarray results (Table 3 and Figure 2), the differential expression patterns were coincident, indicating the reliability of the microarray analysis.
Toll-like receptors pathway analysis
Activation of the innate immune response is controlled in large part by the TLR family of pattern-recognition receptors (PRRs). A previous study showed that the DNA of PCV2 could impair the plasmacytoid- and monocyte-derived dendritic cells ability to induce interferon (IFN)-α and TNF-α by interacting with TLR7 and TLR9 . In the current study, we verified the expression of TLR1, TLR2, TLR3, TLR4, CD14, TLR5, TLR6, TLR8, and TLR9 (the primers are listed in Table 2. Surprisingly, the two time points examined showed significant differences in up-regulation of TLR1 (2.29-fold vs. 4.89 fold), TLR4 cofactor CD14 (2.35-fold vs. 2.71-fold) and TLR9 (4.66-fold vs. 2.24-fold), and down-regulation of TLR8 (3.01-fold vs. 1.72-fold).
Profile of PAM secreted cytokines in response to PCV2-infection
When compared with the supernatant of mock-inoculated PAMs, PCV2-infected PAMs secreted significantly higher levels of interleukin (IL)-8 (340.09 ± 53.64 vs 77.43 ± 13.15, by 4~6-fold; Figure 3D). Tumor necrosis factor-alpha (TNFα) (228.13 ± 49.36 vs 77.12 ± 9.10) and IL-6 (1827.02 ± 157.32 vs 542.70 ± 137.05) levels were significantly higher in supernatants of PCV2-inoculated PAMs at 24 HPI (2~3-fold vs. mock-inoculated PAMs; Figure 3B-C). Interferon-gamma (IFNγ), however, appeared to be significantly down-regulated (552.87 ± 83.25 vs 889.99 ± 141.03, by 1.70-fold) at 24 HPI but significantly up-regulated (1848.10 ± 190.19 vs 892.37 ± 89.66, by 2-fold) at 48 HPI (vs. mock-inoculated; Figure 3A). The relative mRNA levels of porcine IL-8, TNFα, IL-6, and IFNγ showed the same DE trends as the secreted protein levels (Table 3).
PCV2 is an important pathogen of swine that has effects on both the pork industry and animal welfare. Despite the high mortality rates and widespread incidence of this pathogen, a definitive cure has not yet been found. An effective vaccine was created by Merial in 2005 (Circovac®; Merial, France) , and entered into widespread use in 2009. Despite the fact that the current vaccine strategy can efficiently control PCV2 infection, immunization failures exist in the field, and the molecular mechanisms underlying PMWS and porcine dermatitis and nephropathy syndrome (PDNS) caused by PCV2 remain largely unknown. The present study aimed to identify genes involved in the immune response against PCV2 in primary alveolar macrophages. Furthermore, by generating a comprehensive transcriptomic profile of the temporal PCV2 pathogenic process in the target host cells, we hoped to gain insights into the underlying molecular interactions and signaling pathways that may represent novel targets of improved preventative and therapeutic strategies.
Toll-like receptors are essential sensors of microbial infection that are involved in the recognition of a variety of microbial products [25–27]. To date, 13 TLRs have been identified and characterized for their particular cognate ligands and signaling cascades to trigger pathogen-targeted immune responses . For example, TLR7 recognizes guanosine- and uridine-rich sequences derived from ssRNA viruses , while TLR9 recognizes single-stranded DNA (ssDNA) in endosomal compartments and signals to induce the production of type 1 interferons (IFNs) and pro-inflammatory cytokines and chemokines; PCV2 DNA has been shown to modulate cytokine secretion thorough interaction with TLR9 . In our study of the circular, single-stranded DNA virus, PCV2, differential gene expression of TLRs was found. Specifically, PCV2 infection of PAMs led to increased expression of TLR1, TLR4 cofactor CD14, and TLR9. Monocytic CD14 signaling is usually activated upon exposure to bacterial LPS, and this molecule’s up-regulation may induce ineffective antibacterial immunity. In addition, a similar mechanism has been demonstrated in human monocytic cells upon exposure to the severe acute respiratory syndrome (SARS) coronavirus . These receptors’ up-regulation may contribute to the increased sensitivity of macrophages to other pathogens, thereby augmenting the inflammation and host immune response. Such an enhanced immune response may be beneficial, but an over-stimulation of the innate response could result in tissue damage. For example, the same phenomenon was observed upon macrophage exposure to tobacco smoke, wherein the expression of TLR3 was increased and severe inflammation was observed after treated with TLR3 ligand .
MHC class II molecules
MHC class II molecules play a central role in the initiation of the immune response by binding and presenting immunogenic peptides to CD4+ Th lymphocytes. However, only a few cell types express MHC class II molecules, such as the antigen presenting cell types that include macrophages, B lymphocytes, Langerhans cells, and dendritic cells. The very complex and tight regulation of MHC class II gene expression has direct implications for T lymphocyte activation. These molecules include DOA, DOB, DQA1, DQB1, DRA, and DRB1, and their molecule chaperone CD74; all of which were up-regulated in PCV2-infected PAMs. This expression pattern contradicts the host response observed for other viruses, such as African swine fever virus (CSFV) , classical swine fever virus  and influenza A virus , and that was observed in follicles and PAMs of PMWS piglets [36, 37]. Increased expression of MHC class II molecular have been shown to be correlated to increased interferon production in pigs infection infected with hog cholera virus . IFN-γ up-regulation have been shown after 48 HPI, its up-regulation could benefit PCV2 replication . The aberrant expression of class II antigens could also compromise production of immunoregulatory cytokines, compromising immune function and perhaps favor viral persistence through escape from immune recognition.
Apoptosis and growth arrest
Apoptosis is considered to be an important host mechanism that interrupts viral replication and eliminates virus-infected cells. Viruses often kill infected cells by inducing apoptosis rather than necrosis, but repress apoptosis to prolong the life of the cell and increase the yield of progeny virions. SERPINB9, also called PI-9, belongs to the large superfamily of serine proteinase inhibitors, which bind to and inactivate serine proteinases, and acts to protect against granzyme B-mediated apoptosis . PCV2-infected PAMs were found to have enhanced SERPINB9 expression, which may serve to prevent apoptosis and allow for longer-term replication and more virion production. BNIP3, an apoptosis-inducing dimeric mitochondrial protein , was also found to be increased in PCV2-infected PAMs. It is possible that this factor may represent the host’s counter response to PCV2, by which host-induced apoptosis would act to limit the amount of infected cells. It is also possible that these BNIP3 overexpressing cells represent a subset of PAMs that have reached full capacity of virus load, in which the virus itself may have manipulated the host gene expression to facilitate cell death and release of the virus progeny.
Actin alpha 2 (ACTA2) is a smooth muscle α-actin, whose expression is transformation-sensitive to growth signals in normal cells. Actively proliferating fibroblasts and smooth muscle cells have low levels of ACTA2, but the inhibition of cell proliferation by density arrest or treatment with antimitotic agents induces the smooth muscle α-actin promoter . Thus, the ACTA2 up-regulation that was detected in PCV2-infected PAMs may be the result of growth arrest induced by the virus The same phenomenon was found in Rhabdomyosarcoma cells after EV71 infection .
ENO1 is involved in the induction of cell death in fibroblasts and has been implicated in tumorigenicity of human breast carcinoma cells . Its regulation of cell growth is further suggested by the presence of functional transcriptional repression domains . PCV2-infected PAMs showed enhanced ENO1 expression. It is possible that up-regulation of this gene may enhance growth arrest in PCV2-infected cells, thereby culminating in cell death.
Serum amyloid A 1 (SAA1) is a member of the SAA family of apolipoproteins synthesized in response to cytokines released by activated macrophages, which often serve as clinical markers of acute and chronic inflammatory diseases . The pattern of increased SAA1 expression observed in PCV2-infected PAMs suggested that PCV2 induced an acute inflammatory response at the early time point (24 HPI) of infection. TNF-α and IL-1β are acute-phase proinflammatory mediators that promote inflammation and induce fever, tissue destruction, and, in some cases, shock and death . Our study showed that PCV2 could induce TNF-α and IL-1β expression (Table 1) in PAMs at 24 HPI and 48 HPI. These results are in good agreement with those from a previous study of lymphoid tissues and peripheral blood mononuclear cells (PBMCs) of PMWS pigs [17, 48]. SAA1, TNF-α, and IL-1β up-regulation after PCV2 infection suggests the occurrence of systemic inflammation following PCV2 infection.
Lymphohistiocytic to granulomatous pneumonia is a common histopathological finding in naturally-acquired PMWS and experimentally-infected PCV2 pigs [49, 50]. Chemokines, such as MCP-1/CCL-2, were found to be up-regulated in the PCV2-infected PAMs. In a previous study, MCP-1 was suggested to play an important role in the pathogenesis of granulomatous inflammation in PMWS-affected pigs [51–53]. The progressive granulomatous inflammation has been speculated to be compromise organ function and result in PMWS-associated fatalities . Furthermore, the MCP-1 cytokine is known to be an important chemokine for the recruitment of monocytes from the blood . Thus, macrophage infiltration has been suggested as an important feature of the pathogenesis and progression of PMWS .
Recruitment of inflammatory cells in injured tissue is regulated by various cytokines and chemokines. PMWS-affected pigs present with increased expression of MIP-1, CCL-2, IL-1β, IL-8, and TNF-α mRNA in the lymphoid tissues and PBMCs, suggesting that such cytokines and chemokines may contribute to the development of granulomatous lymphadenitis and recruitment of macrophages or other inflammatory cells to the PCV2-infected tissues [17, 52, 57]. The PCV2-infected PAMs in our study showed up-regulation of MIP-1, CCL-2, IL-1β, IL-8 and TNF-α mRNA expression. TNF-α is known to activate microvascular endothelium and cause a pyrexic response . Increased production of TNF-α may induce fever and respiratory distress in PCV2-infected pigs [58, 59]. IL-8 is a chemoattractant for neutrophils and other polymorphonuclear leukocytes (PMNs) that are produced following acute infection . Interaction of IL-8 and MCP-1/CCL-2 is also known to trigger tight adhesion of monocytes to vascular endothelium under flow conditions . The ability of PCV2 to induce IL-8 production in PAMs or PBMCs has been previously reported [57, 62, 63]. Similar to IL-8, the PAM-derived chemokine AMCF-II/CXCL5  is able to recruit PMNs to a lesion area . Unlike the PMN-derived chemoattractants discussed above, MCP-1 is a powerful chemoattractant for monocytes and its activities can be further enhanced by IL-8 . It is reasonable to speculate that the development of interstitial pneumonia in pigs suffering from naturally-acquired PMWS or experimental infection of PCV2 may result from the recruitment of acute phase inflammatory cells and subsequent mononuclear phagocytic cells from blood vessels to interalveolar septa and alveolar spaces by up-regulated cytokines and chemokines released from PCV2-containing PAMs.
Antiviral immune response
IFNγ is a Th1-specific cytokine produced by macrophages, NK cells, and other cell types during the onset of viral infection. It has a potent antiviral property that contributes to the control of acute viral infections, and is an important mediator of cellular responses. IFNγ mRNA was found to be slightly down-regulated at the early stage of in vitro PCV2 infection (24 HPI) of PAMs. This repressed expression may be indicative of suppressed Th1 responses, which may facilitate viral persistence and delayed viral clearance. At the late stage of in vitro infection (48 HPI), however, the IFNγ expression level increased. Such an up-regulation of expression may indicate a virus-modulated mechanism to enhance infection and replication at the late stage of PCV2 infection in vivo.
This is the first study to evaluate the gene expression profile of PCV2-infected PAMs in vitro. Microarray analysis showed that expression of 266 and 175 PAM genes was altered after 24 and 48 hours of PCV2 infection, respectively. Among these, several genes related to the inflammatory response were up-regulated, suggesting that PCV2 may induce systemic inflammation. Collectively, this work established a comprehensive differential transcription profile of early and late PCV2 infection in PAMs, which revealed that PCV2 can induce persistent cytokine production mainly through the TLR1 and TLR9 pathways.
Virus and cells
The PCV2-WH strain (GenBank accession number: FJ870967) was isolated from pooled samples of spleen and lymph nodes of a PMWS-affected pig. The virus was amplified by passaging in PCV-free porcine kidney (PK-15) cells and confirmed as PCV2 by nucleotide sequence analysis and reactivity with PCV2 Cap protein-specific monoclonal antibody . The titer of the PCV2 inoculum was 1.0×107 TCID50/mL, as determined by titration in PCV-free PK-15 cells.
All animals’ tissue collection procedures were performed according to protocols approved by the Hubei Province PR China for Biological Studies Animal Care and Use Committee. Three 3-week-old Large White piglets were obtained from a herd in Hubei province, which tested negative for both PCV2 and porcine reproductive and respiratory syndrome virus (PRRSV), and were raised in controlled lab conditions. For the first three days of containment, the pigs were administered daily doses of enrofloxacin (1 mL of 5% solution) and lincospectin/spectinomycin (1 mL of 5 or 10% solution) to eliminate any residual bacterial pathogens . Seven days later, the piglets were sacrificed. PAMs were collected by bronchoalveolar lavage and frozen in liquid nitrogen, as previously described . Prior to use, the cells were confirmed as negative for PCV1, PCV2, parvovirus, pseudorabies virus (PRV), classical swine fever virus (CSFV), and PRRSV variously by PCR and RT-PCR [69, 70]. The PAM phenotype was confirmed by flow cytometric detection of the macrophage markers SWC3, CD169, and SLAII .
RNA isolation and IFA
PAMs were thawed and cultured for 48 h before treatment, as previously described [71, 72]. One primary culture from each animal was split into two for experimental analysis. The first was infected with PCV2-WH from the 11th passage at a multiplicity of infection (MOI) of 1. The second was mock-inoculated with DMEM to serve as a control. Virus uptake by the cells was determined by IFA using a PCV2 Cap protein-specific monoclonal antibody  and by quantitative real-time PCR (qPCR) . Briefly, DNA was extracted from the cells using the E.Z.N.A.™ Viral DNA Kit (OMEGA, USA) according to the manufacturer’s instructions. The amplification was performed in a 25 μL reaction mixture containing 12.5 μL of 2× THUNDERBIRD Probe qPCR Mix (TOYOBO, Japan), 9.45 μL sterile, nuclease-free water, 10 pmol of forward primer (5′-CCAGGAGGGCGTTCTGACT-3'), 10 pmol of reverse primer (5′-CGTTACCGCTGGAGAAGGAA-3′), 4 pmol of TaqMan probe (5′-AATGGCATCTTCAACACCCGCCTCT-3′), 0.05 μL of 50× ROX reference dye, and 2 μL of the extracted DNA. The reaction was run in a real-time thermocycler (7500 PCR System; Applied BioSystems Inc.) with the following conditions: 1 cycle at 50°C for 2 min, 1 cycle at 95°C for 10 min, and 45 cycles at 95°C for 15 s and 62°C for 60 s (real time).
At 24 HPI and 48 HPI, the cells were collected for RNA extraction with the TRIzol reagent (Invitrogen Life Technologies, USA). The extracted RNA was further purified with an RNA clean-up kit (Macherey-Nagel, Germany) and applied to a Bioanalyzer 2100 spectrophotometer (Agilent Technologies, USA) to assess the integrity, quality, and quantity of the sample. The corresponding supernatants were collected for subsequent protein analysis and stored at −80°C until use.
RNA labeling, microarray hybridization, and data analysis
RNA was converted to aminoallyl-RNA with the Amino Allyl II MessageAmp™ aRNA Amplification Kit (Ambion, USA) and labeled with either Cy3 or Cy5 (Reactuve Dye Pack; Amersham, Sweden). Each Cy3- (or Cy5-) coupled PAM sample was then hybridized to its corresponding Cy5- (or Cy3-) coupled reference sample anchored on a pig microarray slide containing 20,400 oligonucleotides (US Pig Genome Coordination Program; http://www.animalgenome.org/pig/resources/array_request.html). Hybridizations and analyses were carried out according to a previously published protocol from CapitalBio Corporation (Beijing, China) . Raw data was extracted from the TIFF images using LuxScan 3.0 software (CapitalBio). A spot-exclusion method was adopted to filter out faint spots, with signal intensities in the lowest 50%, and exclude the corresponding genes from further analysis. An intensity-dependent LOWESS program in the R language package was used to normalize the two-channel ratio values. As a measure of technical replication, each experiment was conducted with the corresponding dye swap.
Gene Cluster 3.0 and Eisen’s Treeview software (Stanford University, USA) were used to compare similarities among individual samples. Cy3/Cy5 ratios were log-transformed (base 2), median centered by arrays and genes, and hierarchically clustered (average linkage correlation metric). To determine gene products with significantly up-regulated expression, the PAM 24 HPI data sets were assessed for increases in average intensity of at least 2-fold. Statistical comparisons were performed by the one-class method in the significance analysis of microarray (SAM software version 3.0, Stanford University) . All the microarray results from this study were deposited in the NCBI Gene Expression Omnibus database under the following accession numbers: platform, GPL13968; Series, GSE30918; Samples, GSM766444, GSM766445, GSM766446, GSM766447, GSM766448, GSM766449.
Functional interpretation of microarray data as well as pathway and network analysis
Differentially expressed genes were selected for pathway exploration using the IPA web-based query system (Ingenuity Systems, http://www.ingenuity.com), this system can generate the functional analysis and networks that are most significant to the data set, IPA compares proteins in the input group and displays a rank-ordered list of pathways and networks, whose activities are most likely affected. Ficher’s exact test was used to calculate a p-value determining the probability that each biofunction assigned to that data set is due to chance alone. Genes from the dataset that met the fold change cut-off of 2.0 were considered for the analysis. The score is the probability that a collection of proteins equal to or greater than the number in a network could be achieved by chance alone. A score of 3 indicates that there is a 1/1000 chance that the focus proteins are in a network due to random chance, therefore, score of 3 or higher have a 99.9% confidence of not being generated by random chance alone.
Equal amounts of RNA (1 μg) from PAMs was used to generate cDNA with oligo (dT) primers using the SuperScript II RNase H Reverse Transcriptase Kit (Invitrogen, USA). Quantitative RT-PCR reactions were set up in triplicate using reagents from the SYBR Green I PCR Kit (Toyobo) and 10 μM primers (designed by the Beacon Designer 7.0 program, http://www.premierbiosoft.com; Table 2). The amplification reactions were performed on a LightCycler 480 (Roche) with the following conditions: 30 s at 95°C, 40 cycles of 95°C for 10 s, 60°C for 5 s and 72°C for 10 s. Subsequently, melting curve analysis and quantitative analysis of the data were performed using the Roche LightCycler 480 software version 1.5.0. Each sample was run in triplicate. Samples were normalized at each respective time point, using mock-infected PAMs as calibrators and the hypoxantine phosphoribosyltransferase (HPRT-1) and glyceraldehyde-3-phosphatedehydrogenase (GAPDH) as the reference genes.. The average threshold cycles of the replicates (each gene was detected three times) were used to calculate the fold-changes between PCV2-infected or mock-inoculated PAMs by using the formula of 2-delta-delta Ct described by Livak and Schmittgen .
Cytokine detection in culture supernatant
The secreted protein levels of IFNγ, IL-6, IL-8 and TNFα were determined by testing culture supernatants with the respective commercial ELISA kits (R&D Systems, USA). All samples were tested in triplicate and read at 450 nm using an ELISA plate reader (Bio-Tek, USA).
Statistical analysis of quantitative PCR data
Statistical analyses were carried out using Microsoft Excel 2007 (Microsoft Co., USA). Differences between groups were assessed by one-way repeated measures analysis of variance followed by Tukey’s multiple comparison tests. P-values less than 0.05 were considered to indicate statistical significance.
Classical swine fever virus
False discovery rate
Hours post infection
Ingenuity pathway analysis
Open reading frame
Porcine alveolar macrophages
Peripheral blood mononuclear cell
Polymerase chain reaction
Porcine reproductive and respiratory syndrome virus
Quantitative reverse transcription PCR
Tumor necrosis factor-alpha.
Cano J, Rodr¨ªguez-Cari o C, Sogbe E, Utrera V, D¨ªaz C, Segal¨¦s J, Olvera A, Calsamiglia M: Postweaning multisystemic wasting syndrome in pigs in Venezuela. Vet Rec. 2005, 156 (19): 620-
Muhling J, Raye WS, Buddle JR, Wilcox GE: Genetic characterisation of Australian strains of porcine circovirus types 1 and 2. Aust Vet J. 2006, 84 (12): 421-425. 10.1111/j.1751-0813.2006.00081.x.
Done S, Gresham A, Potter R, Chennells D: PMWS and PDNS-two recently recognised diseases of pigs in the UK. In Pract. 2001, 23 (1): 14-10.1136/inpract.23.1.14.
Solis M, Wilkinson P, Romieu R, Hernandez E, Wainberg MA, Hiscott J: Gene expression profiling of the host response to HIV-1 B, C, or A/E infection in monocyte-derived dendritic cells. Virology. 2006, 352 (1): 86-99. 10.1016/j.virol.2006.04.010.
Chakrabarti AK, Vipat VC, Mukherjee S, Singh R, Pawar SD, Mishra AC: Host gene expression profiling in influenza A virus-infected lung epithelial (A549) cells: a comparative analysis between highly pathogenic and modified H5N1 viruses. Virol J. 2010, 7: 219-10.1186/1743-422X-7-219.
Alkhalil A, Hammamieh R, Hardick J, Ichou MA, Jett M, Ibrahim S: Gene expression profiling of monkeypox virus-infected cells reveals novel interfaces for host-virus interactions. Virol J. 2010, 7: 173-10.1186/1743-422X-7-173.
Bernasconi M, Berger C, Sigrist JA, Bonanomi A, Sobek J, Niggli FK, Nadal D: Quantitative profiling of housekeeping and Epstein-Barr virus gene transcription in Burkitt lymphoma cell lines using an oligonucleotide microarray. Virol J. 2006, 3: 43-10.1186/1743-422X-3-43.
Li R, Zhang A, Chen B, Teng L, Wang Y, Chen H, Jin M: Response of swine spleen to Streptococcus suis infection revealed by transcription analysis. BMC Genomics. 2010, 11: 556-10.1186/1471-2164-11-556.
Chen H, Li C, Fang M, Zhu M, Li X, Zhou R, Li K, Zhao S: Understanding Haemophilus parasuis infection in porcine spleen through a transcriptomics approach. BMC Genomics. 2009, 10: 64-10.1186/1471-2164-10-64.
Gilpin DF, McCullough K, Meehan BM, McNeilly F, McNair I, Stevenson LS, Foster JC, Ellis JA, Krakowka S, Adair BM, et al: In vitro studies on the infection and replication of porcine circovirus type 2 in cells of the porcine immune system. Vet Immunol Immunop. 2003, 94 (3–4): 149-161.
Chang HW, Pang VF, Chen LJ, Chia MY, Tsai YC, Jeng CR: Bacterial lipopolysaccharide induces porcine circovirus type 2 replication in swine alveolar macrophages. Vet Microbiol. 2006, 115 (4): 311-319. 10.1016/j.vetmic.2006.03.010.
Allan GM, Ellis JA: Porcine circoviruses: a review. J Vet Diagn Invest. 2000, 12 (1): 3-14. 10.1177/104063870001200102.
Choi C, Chae C: Distribution of porcine parvovirus in porcine circovirus 2-infected pigs with postweaning multisystemic wasting syndrome as shown by in-situ hybridization. J Comp Pathol. 2000, 123 (4): 302-305. 10.1053/jcpa.2000.0421.
McGuire K, Glass EJ: The expanding role of microarrays in the investigation of macrophage responses to pathogens. Vet Immunol Immunop. 2005, 105 (3): 259-275.
Chae C: Postweaning multisystemic wasting syndrome: a review of aetiology, diagnosis and pathology. Vet J. 2004, 168 (1): 41-49.
Chae C: A review of porcine circovirus 2-associated syndromes and diseases. Vet J. 2005, 169 (3): 326-336. 10.1016/j.tvjl.2004.01.012.
Darwich L, Pie S, Rovira A, Segales J, Domingo M, Oswald IP, Mateu E: Cytokine mRNA expression profiles in lymphoid tissues of pigs naturally affected by postweaning multisystemic wasting syndrome. J Gen Virol. 2003, 84: 2117-2125. 10.1099/vir.0.19124-0.
Tomas A, Fernandes LT, Sanchez A, Segales J: Time course differential gene expression in response to porcine circovirus type 2 subclinical infection. Vet Res. 2010, 41 (1): 12-10.1051/vetres/2009060.
Fernandes LT, Tomas A, Bensaid A, Perez-Enciso M, Sibila M, Sanchez A, Segales J: Exploratory Study on the Transcriptional Profile of Pigs Subclinically Infected with Porcine Circovirus Type 2. Anim Biotechnol. 2009, 20 (3): 96-109. 10.1080/10495390902885785.
Lee G, Han D, Song JY, Lee YS, Kang KS, Yoon S: Genomic expression profiling in lymph nodes with lymphoid depletion from porcine circovirus 2-infected pigs. J Gen Virol. 2010, 91 (Pt 10): 2585-2591.
Tomás A, Fernandes LT, Sánchez A, Segalés J: Time course differential gene expression in response to porcine circovirus type 2 subclinical infection. Vet Res. 2009, 41 (1): 12-12.
Cheng S, Zhang M, Li W, Wang Y, Liu Y, He Q: Proteomic analysis of porcine alveolar macrophages infected with porcine circovirus type 2. J Proteomics. 2012, 75 (11): 3258-3269. 10.1016/j.jprot.2012.03.039.
Vincent IE, Balmelli C, Meehan B, Allan G, Summerfield A, McCullough KC: Silencing of natural interferon producing cell activation by porcine circovirus type 2 DNA. Immunology. 2007, 120 (1): 47-56.
Charreyre C, Bésème S, Brun A, Bublot M, Joisel F, Lapostolle B, Sierra P, Vaganay A: Protection of piglets against a PCV2 experimental challenge by vaccinating gilts with inactivated adjuvanted PCV2 vaccine. Proceedings of the 20th IPVS. 2005, 26-30.
Janeway CA, Medzhitov R: Innate immune recognition. Annu Rev Immunol. 2002, 20: 197-216. 10.1146/annurev.immunol.20.083001.084359.
Kawai T, Akira S: TLR signaling. Cell Death Differ. 2006, 13 (5): 816-825. 10.1038/sj.cdd.4401850.
Takeda K, Kaisho T, Akira S: Toll-like receptors. Annu Rev Immunol. 2003, 21: 335-376. 10.1146/annurev.immunol.21.120601.141126.
Kawai T, Akira S: Toll-like receptors and their crosstalk with other innate receptors in infection and immunity. Immunity. 2011, 34 (5): 637-650. 10.1016/j.immuni.2011.05.006.
Lund J, Alexopoulou L, Sato A, Karow M, Adams N, Gale N, Iwasaki A, Flavell R: Recognition of single-stranded RNA viruses by Toll-like receptor 7. P Natl Acad Sci USA. 2004, 101 (15): 5598-10.1073/pnas.0400937101.
Balmelli C, Steiner E, Moulin H, Peduto N, Herrmann B, Summerfield A, McCullough K: Porcine circovirus type 2 DNA influences cytoskeleton rearrangements in plasmacytoid and monocyte‒derived dendritic cells. Immunology. 2011, 132 (1): 57-65. 10.1111/j.1365-2567.2010.03339.x.
Hu W, Yen Y-T, Singh S, Kao C-L, Wu-Hsieh BA: SARS-CoV Regulates Immune Function-Related Gene Expression in Human Monocytic Cells. Viral Immunol. 2012, 25 (4): 277-288. 10.1089/vim.2011.0099.
Koarai A, Yanagisawa S, Sugiura H, Ichikawa T, Akamatsu K, Hirano T, Nakanishi M, Matsunaga K, Minakata Y, Ichinose M: Cigarette smoke augments the expression and responses of toll‒like receptor 3 in human macrophages. Respirology. 2012, 17 (6): 1018-1025. 10.1111/j.1440-1843.2012.02198.x.
Juarrero MG, Mebus CA, Pan R, Revilla Y, Alonso JM, Lunney JK: Swine leukocyte antigen and macrophage marker expression on both African swine fever virus-infected and non-infected primary porcine macrophage cultures. Vet Immunol Immunop. 1992, 32 (3–4): 243-259.
Li J, Yu YJ, Feng L, Cai XB, Tang HB, Sun SK, Zhang HY, Liang JJ, Luo TR: Global transcriptional profiles in peripheral blood mononuclear cell during classical swine fever virus infection. Virus Res. 2010, 148 (1–2): 60-70.
Liang Q, Luo J, Zhou K, Dong J, He H: Immune-related gene expression in response to H5N1 avian influenza virus infection in chicken and duck embryonic fibroblasts. Mol Immunol. 2011, 48 (6-7): 924-930. 10.1016/j.molimm.2010.12.011.
Mandrioli L, Sarli G, Zengarini M, Panarese S, Marcato P: Immunohistochemical MHC-II and interleukin 2-R (CD25) expression in lymph nodes of pigs with spontaneous postweaning multisystemic wasting syndrome (PMWS). Vet Pathol. 2006, 43 (6): 993-997. 10.1354/vp.43-6-993.
McNeilly F, Allan G, Foster J, Adair B, McNulty M: Effect of porcine circovirus infection on porcine alveolar macrophage function. Vet Immunol Immunop. 1996, 49 (4): 295-306. 10.1016/0165-2427(95)05476-6.
Torlone V, Titoli F, Gialletti L: Circulating interferon production in pigs infected with hog cholera virus. Life Sci. 1965, 4 (17): 1707-10.1016/0024-3205(65)90218-3.
Meerts P, Misinzo G, Nauwynck HJ: Enhancement of porcine circovirus 2 replication in porcine cell lines by IFN-gamma before and after treatment and by IFN-alpha after treatment. J Interf Cytok Res. 2005, 25 (11): 684-693. 10.1089/jir.2005.25.684.
Ekert PG, Silke J, Vaux DL: Caspase inhibitors. Cell Death Differ. 1999, 6 (11): 1081-1086. 10.1038/sj.cdd.4400594.
Chen G, Ray R, Dubik D, Shi L, Cizeau J, Bleackley RC, Saxena S, Gietz RD, Greenberg AH: The E1B 19K/Bcl-2–binding Protein Nip3 is a Dimeric Mitochondrial Protein that Activates Apoptosis. J Exp Med. 1997, 186 (12): 1975-1983. 10.1084/jem.186.12.1975.
Kumar CC, Kim JH, Bushel P, Armstrong L, Catino JJ: Activation of smooth muscle alpha-actin promoter in ras-transformed cells by treatments with antimitotic agents: correlation with stimulation of SRF:SRE mediated gene transcription. J Biochem. 1995, 118 (6): 1285-1292.
Leong WF, Chow VTK: Transcriptomic and proteomic analyses of rhabdomyosarcoma cells reveal differential cellular gene expression in response to enterovirus 71 infection. Cell Microbiol. 2006, 8 (4): 565-580. 10.1111/j.1462-5822.2005.00644.x.
Ray RB, Steele R, Seftor E, Hendrix M: Human breast carcinoma cells transfected with the gene encoding a c-myc promoter-binding protein (MBP-1) inhibits tumors in nude mice. Cancer Res. 1995, 55 (17): 3747-3751.
Ghosh AK, Steele R, Ray RB: Functional domains of c-myc promoter binding protein 1 involved in transcriptional repression and cell growth regulation. Mol Cell Biol. 1999, 19 (4): 2880-2886.
Malle E, De Beer FC: Human serum amyloid A (SAA) protein: a prominent acute-phase reactant for clinical practice. Eur J Clin Invest. 1996, 26 (6): 427-435. 10.1046/j.1365-2362.1996.159291.x.
Dinarello CA: Proinflammatory cytokines. Chest. 2000, 118 (2): 503-508. 10.1378/chest.118.2.503.
Sipos W, Duvigneau JC, Willheim M, Schilcher F, Hartl RT, Hofbauer G, Exel B, Pietschmann P, Schmoll F: Systemic cytokine profile in feeder pigs suffering from natural postweaning multisystemic wasting syndrome (PMWS) as determined by semiquantitative RT-PCR and flow cytometric intracellular cytokine detection. Vet Immunol Immunop. 2004, 99 (1–2): 63-71.
Opriessnig T, Halbur PG: Concurrent infections are important for expression of porcine circovirus associated disease. Virus Res. 2012, 164 (1-2): 20-32. 10.1016/j.virusres.2011.09.014.
Magar R, Larochelle R, Thibault S, Lamontagne L: Experimental transmission of porcine circovirus type 2 (PCV2) in weaned pigs: a sequential study. J Comp Pathol. 2000, 123 (4): 258-269. 10.1053/jcpa.2000.0413.
Kim J, Chae C: Expression of monocyte chemoattractant protein-1 but not interleukin-8 in granulomatous lesions in lymph nodes from pigs with naturally occurring postweaning multisystemic wasting syndrome. Vet Pathol. 2003, 40 (2): 181-186. 10.1354/vp.40-2-181.
Kim J, Chae C: Expression of monocyte chemoattractant protein-1 and macrophage inflammatory protein-1 in porcine circovirus 2-induced granulomatous inflammation. J Comp Pathol. 2004, 131 (2–3): 121-126.
Tsai Y, Jeng C, Hsiao S, Chang H, Liu J, Chang C, Lin C, Chia M, Pang V: Porcine circovirus type 2 (PCV2) induces cell proliferation, fusion, and chemokine expression in swine monocytic cells in vitro. Vet Res. 2010, 41 (5): 60-10.1051/vetres/2010032.
Krakowka S, Ellis J, McNeilly F, Gilpin D, Meehan B, McCullough K, Allan G: Immunologic features of porcine circovirus type 2 infection. Viral Immunol. 2002, 15 (4): 567-582. 10.1089/088282402320914511.
Mantovani A, Bonecchi R, Locati M: Tuning inflammation and immunity by chemokine sequestration: decoys and more. Nat Rev Immunol. 2006, 6 (12): 907-918. 10.1038/nri1964.
Segalés J, Domingo M, Chianini F, Majó N, Domı́nguez J, Darwich L, Mateu E: Immunosuppression in postweaning multisystemic wasting syndrome affected pigs. Vet Microbiol. 2004, 98 (2): 151-158. 10.1016/j.vetmic.2003.10.007.
Darwich L, Balasch M, Plana-Duran J, Segales J, Domingo M, Mateu E: Cytokine profiles of peripheral blood mononuclear cells from pigs with postweaning multisystemic wasting syndrome in response to mitogen, superantigen or recall viral antigens. J Gen Virol. 2003, 84: 3453-3457. 10.1099/vir.0.19364-0.
Albina E, Truong C, Hutet E, Blanchard P, Cariolet R, L’Hospitalier R, Mahe D, Allee C, Morvan H, Amenna N, et al: An experimental model for post-weaning multisystemic wasting syndrome (PMWS) in growing piglets. J Comp Pathol. 2001, 125 (4): 292-303. 10.1053/jcpa.2001.0508.
Rovira A, Balasch M, Segales J, Garcia L, Plana-Duran J, Rosell C, Ellerbrok H, Mankertz A, Domingo M: Experimental inoculation of conventional pigs with porcine reproductive and respiratory syndrome virus and porcine circovirus 2. J Virol. 2002, 76 (7): 3232-3239. 10.1128/JVI.76.7.3232-3239.2002.
Baggiolini M, Dewald B, Moser B: Interleukin-8 and related chemotactic cytokines-CXC and CC chemokines. Adv Immunol. 1994, 55: 97-179.
Gerszten RE, Garcia-Zepeda EA, Lim YC, Yoshida M, Ding HA, Gimbrone MA, Luster AD, Luscinskas FW, Rosenzweig A: MCP-1 and IL-8 trigger firm adhesion of monocytes to vascular endothelium under flow conditions. Nature. 1999, 398 (6729): 718-723. 10.1038/19546.
Chang HW, Jeng CR, Lin TL, Liu JJ, Chiou MT, Tsai YC, Chia MY, Jan TR, Pang VF: Immunopathological effects of porcine circovirus type 2 (PCV2) on swine alveolar macrophages by in vitro inoculation. Vet Immunol Immunop. 2006, 110 (3–4): 207-219.
Forta M, Fernandes LT, Nofrarias M, Diaz I, Sibila M, Pujols J, Mateu E, Segales J: Development of cell-mediated immunity to porcine circovirus type 2 (PCV2) in caesarean-derived, colostrum-deprived piglets. Vet Immunol Immunop. 2009, 129 (1–2): 101-107.
Smith JB, Rovai LE, Herschman HR: Sequence similarities of a subgroup of CXC chemokines related to murine LIX: implications for the interpretation of evolutionary relationships among chemokines. J Leukoc Biol. 1997, 62 (5): 598-603.
Goodman RB, Forstrom JW, Osborn SG, Chi EY, Martin TR: Identification of two neutrophil chemotactic peptides produced by porcine alveolar macrophages. J Biol Chem. 1991, 266 (13): 8455-8463.
Cao S, Chen H, Zhao J, Lu J, Xiao S, Jin M, Guo A, Wu B, He Q: Detection of porcine circovirus type 2, porcine parvovirus and porcine pseudorabies virus from pigs with postweaning multisystemic wasting syndrome by multiplex PCR. Vet Res Commun. 2005, 29 (3): 263-269.
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 (Pt 10): 2550-2564.
Wensvoort G, Terpstra C, Pol JM, Ter Laak EA, Bloemraad M, De Kluyver EP, Kragten C, Van Buiten L, Den Besten A, Wagenaar F, et al: Mystery swine disease in The Netherlands: the isolation of Lelystad virus. Vet Q. 1991, 13 (3): 121-130. 10.1080/01652176.1991.9694296.
Ogawa H, Taira O, Hirai T, Takeuchi H, Nagao A, Ishikawa Y, Tuchiya K, Nunoya T, Ueda S: Multiplex PCR and multiplex RT-PCR for inclusive detection of major swine DNA and RNA viruses in pigs with multiple infections. J Virol Methods. 2009, 160 (1–2): 210-214.
Choi KS, Chae JS: Genetic characterization of porcine circovirus type 2 in Republic of Korea. Res Vet Sci. 2008, 84 (3): 497-501. 10.1016/j.rvsc.2007.05.017.
Delputte PL, Nauwynck HJ: Porcine arterivirus infection of alveolar macrophages is mediated by sialic acid on the virus. J Virol. 2004, 78 (15): 8094-8101. 10.1128/JVI.78.15.8094-8101.2004.
Chueh LL, Lee KH, Wang FI, Pang VF, Weng CN: Sequence analysis of the nucleocapsid protein gene of the porcine reproductive and respiratory syndrome virus Taiwan MD-001 strain. Adv Exp Med Biol. 1998, 440: 795-799.
Yang K, Li W, Niu H, Yan W, Liu X, Wang Y, Cheng S, Ku X, He Q: Efficacy of single dose of an inactivated porcine circovirus type 2 (PCV2) whole-virus vaccine with oil adjuvant in piglets. Acta Vet Scand. 2012, 54 (1): 67-10.1186/1751-0147-54-67.
Patterson TA, Lobenhofer EK, Fulmer-Smentek SB, Collins PJ, Chu TM, Bao WJ, Fang H, Kawasaki ES, Hager J, Tikhonova IR, et al: Performance comparison of one-color and two-color platforms within the MicroArray Quality Control (MAQC) project. Nat Biotechnol. 2006, 24 (9): 1140-1150. 10.1038/nbt1242.
Livak K, Schmittgen T: Analysis of relative gene expression data using real-time quantitative PCR and the 2-[Delta][Delta] CT method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
This work was supported by grants from China Agricultural Research System (No. CARS-36) and Merial Limited.
The authors declare that they have no competing interests.
LWT carried out most of the experiments and wrote the manuscript. HQG, CHC, CC, and JCA critically revised the manuscript and the experimental design. LSQ, WY, DF, YWD, and YK helped with the experiments. All of the authors read and approved the final version of the manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Li, W., Liu, S., Wang, Y. et al. Transcription analysis of the porcine alveolar macrophage response to porcine circovirus type 2. BMC Genomics 14, 353 (2013). https://doi.org/10.1186/1471-2164-14-353
- Differentially Express
- Ingenuity Pathway Analysis
- Differentially Express Gene
- Classical Swine Fever Virus
- Severe Acute Respiratory Syndrome