Figures
Abstract
Porcine reproductive and respiratory syndrome (PRRS) caused by PRRS virus (PRRSV) is an infectious disease characterized by severe reproductive deficiency in pregnant sows, typical respiratory symptoms in piglets, and high mortality rate of piglets. In this study, we employed an Affymetrix microarray chip to compare the gene expression profiles of lung tissue samples from Dapulian (DPL) pigs (a Chinese indigenous pig breed) and Duroc×Landrace×Yorkshire (DLY) pigs after infection with PRRSV. During infection with PRRSV, the DLY pigs exhibited a range of clinical features that typify the disease, whereas the DPL pigs showed only mild signs of the disease. Overall, the DPL group had a lower percentage of CD4+ cells and lower CD4+/CD8+ratios than the DLY group (p<0.05). For both IL-10 and TNF-α, the DLY pigs had significantly higher levels than the DPL pigs (p<0.01). The DLY pigs have lower serum IFN-γ levels than the DPL pigs (p<0.01). The serum IgG levels increased slightly from 0 dpi to 7 dpi, and peaked at 14 dpi (p<0.0001). Microarray data analysis revealed 16 differentially expressed (DE) genes in the lung tissue samples from the DLY and DPL pigs (q≤5%), of which LOC100516029 and LOC100523005 were up-regulated in the PRRSV-infected DPL pigs, while the other 14 genes were down-regulated in the PRRSV-infected DPL pigs compared with the PRRSV-infected DLY pigs. The mRNA expression levels of 10 out of the 16 DE genes were validated by real-time quantitative RT-PCR and their fold change was consistent with the result of microarray data analysis. We further analyzed the mRNA expression level of 8 differentially expressed genes between the DPL and DLY pigs for both uninfected and infected groups, and found that TF and USP18 genes were important in underlying porcine resistance or susceptibility to PRRSV.
Citation: Xing J, Xing F, Zhang C, Zhang Y, Wang N, Li Y, et al. (2014) Genome-Wide Gene Expression Profiles in Lung Tissues of Pig Breeds Differing in Resistance to Porcine Reproductive and Respiratory Syndrome Virus. PLoS ONE 9(1): e86101. https://doi.org/10.1371/journal.pone.0086101
Editor: Shuhong Zhao, Huazhong Agricultural University, China
Received: May 17, 2013; Accepted: December 5, 2013; Published: January 23, 2014
Copyright: © 2014 Xing et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the Natural Science Foundation of Shandong Province, China (www.sdnsf.gov.cn) (ZR2012CZ002) and the NSF, USA (www.nsf.gov) (NSF #EPS–0903787). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Porcine reproductive and respiratory syndrome (PRRS), also known as blue-ear pig disease, is a widespread infectious disease caused by porcine reproductive and respiratory syndrome virus (PRRSV), which is characterized by severe reproductive deficiency in pregnant sows, typical respiratory symptoms in pigs of all ages, increased susceptibility to secondary bacterial infection, and high mortality rate of piglets [1], [2], [3]. It was first reported in 1987 in North America and Central Europe, and is now one of the most severe diseases in the swine industry, costing approximately $560 million each year in the United States alone [4], [5]. The PRRSV is a member of the enveloped, positive-sense, single-stranded RNA virus family Arteriviridae [6] and exhibits a very narrow host cell tropism. PRRSV-infected pigs are susceptible to pneumonia and growth reduction, while infected sows have increased rates of abortion, stillbirth, and mummification, and also give birth to weak piglets with chronic respiratory problems [7].
Breed is one of the factors that determines resistance or susceptibility to PRRSV in pigs [8]. Experimental PRRSV challenge infections in Duroc, Hampshire and Meishan pigs revealed that the latter breed had significantly less PRRSV antigen in their lungs than the other breeds, suggesting that differences in the severity and distribution of PRRSV-induced lesions and normal serum antibody responses to PRRSV may be influenced, at least in part, by genetic factors [9]. Another two independent studies indicated that different pig breeds or lines responded differently to PRRSV infection, while analysis of the serum cytokine levels confirmed gene expression differences in the spleens and bronchial lymph nodes of these pigs [10], [11]. Our previous study revealed differences in the PRRSV copy number and NMMHC-IIA and CD163 mRNA expression levels between virus infected Dapulian (DPL) and Duroc×Landrace×Yorkshire (DLY) pigs [12].
Since the early 1990s, cDNA microarray genomic technology has been used to simultaneously assess mRNA transcription patterns for thousands of genes simultaneously, and has been commonly employed for determining patterns of differential gene expression in a wide range of organisms, thereby attracting considerable attention from the biomedical research community [13], [14]. Microarrays are particularly useful for studying whether cellular mRNAs, differentially regulated by each viral protein, play a crucial role for virus multiplication in the cell and interactions between the host and virus [15], [16]. Several microarray studies have reported that many genes are involved in porcine responses to infection with PRRSV [16], [17], [18].
DPL pig, a breed of indigenous Chinese pigs which is native in Jining City, Shandong Province, China, has some desirable traits, such as strong disease-resistance, high prolificacy and good meat quality. During the summer of 2006, a large-scale devastating outbreak of PRRS overwhelmed about 10 provinces in China, resulting in more than 2 million infected pigs and a death toll of at least 0.4 million animals [19]. However, the DPL pigs showed strong resistance to PRRSV infection and only a few individuals died. The objective of this study was to examine the gene expression differences between lung tissues from DPL and DLY pigs infected with PRRSV (strain JXA1), a highly pathogenic virus isolated from DLY commercial pig stocks, to identify genes associated with PRRSV-resistance. To achieve this objective, we used Affymetrix microarray technology to examine gene expression profiles in the PRRSV infected DPL- and DLY- pigs.
Materials and Methods
Animals and Tissue Collection
All animal procedures were performed according to protocols approved by the Animal Care and Use Committee of Shandong Agricultural University. Sixteen healthy 30-day-old weaned DPL pigs were selected from Jiaxiang Dapulian Farm, Jining City, China, and 15 healthy 30-day-old weaned DLY pigs were obtained from a commercial farm with high standards of animal health. These pigs were free from PRRSV, porcine circovirus type 2 (PCV2), pseudorabies virus (PRV), and classical swine fever virus (CSFV) as determined by ELISA tests for serum antibodies; the absence of PRRSV was also confirmed by real-time quantitative reverse transcription PCR (qRT-PCR). Pigs were randomly assigned into two groups and reared in separate places: the PRRSV-infected group consisted of eleven DPL and ten DLY pigs, and the control group consisted of five DPL and five DLY pigs. Infection in the pigs was administrated via inoculation with 2 ml of a viral suspension of PRRSV (at a tissue culture infectious dose of 105) by dripping the solution into the nasal cavity of each pig. The control group was treated with an identical volume of phosphate buffer saline (PBS) by the same method. Rectal temperatures and clinical examinations on the pigs were recorded daily during the experiment. Anticoagulant-treated blood and untreated blood samples were collected separately at 0, 7, 14, and 21 days post-infection (dpi) from the infected and control groups for assaying CD4+, CD8+, cytokines (interleukin (IL) 1 beta (IL-1β), IL-2, IL-10), interferon (IFN)-gamma (IFN-γ), tumor necrosis factor-alpha (TNF-α), and immunoglobulin G (IgG) protein levels. Lung samples for microarray analysis and real-time qRT-PCR analysis were collected from six infected DLY and DPL pigs (three pigs for each breed) immediately post-slaughter at 28 dpi.
Flow Cytometry and ELISA Assays
Peripheral blood mononuclear cells (PBMC) were separated from whole blood at 0, 7, 14 and 21 dpi by density-gradient centrifugation with Histopaque 1077 (Sigma, Steinheim, Germany). Phenotypic analysis of the PBMC subsets was done by flow cytometry using anti-CD4-FITC and anti-CD8-RPE monoclonal antibodies (AbD Serotec, Kidlington, UK). Analyses were conducted using an Easy Cyte Mini Cytometer System (Guava, Billerica, MA, USA) at an excitation wavelength of 488 nm. Sera collected from the infected pigs at 0, 7 and 14 dpi were used for cytokine detection (IL-1β, IL-2, IL-10, IFN-γ, TNF-α) and IgG levels by ELISA (R & D Systems, Minneapolis, MN, USA) according to manufacturer’s instructions.
Microarray Hybridization
Utilizing the porcine Affymetrix GeneChip to compare the mRNA expression profiles of lung tissue samples from the PRRSV-infected DPL pigs with those from DLY pigs was carried out at Capital Biochip Limited Company, Beijing, China. In brief, total RNA was isolated from lung tissue samples and purified using an RNeasy Mini kit (Qiagen, Shanghai, China) according to the manufacturer’s protocol. RNA was prepared using the GeneChip (AFF-900623) one cycle target for the labeling and control reagents, and the labeled RNA was hybridized in an Affymetrix Hybridization Oven 640. After hybridization, the microarray slides were washed and dyed with the Affymetrix Fluidics Station 450. The probe arrays were scanned with an Affymetrix GeneChip Scanner 3000. The signal intensities of the spots on each image were quantified by the Affymetrix GeneChip Operating Software Version 1.4 (GCOS 1.4), and the data were adjusted and normalized by the dChip software.
Microarray Data Analysis
The collected microarray data were MIAME compliant. The raw data (tag sequences and counts) have been deposited in the NCBI Gene Expression Omnibus (GEO) with accession number GSE49306.
Microarray data was analyzed with BRB-ArrayTools software developed by Dr. Richard Simon and BRB-Array Tools Development Team and available at http://linus.nci.nih.gov/BRB-ArrayTools.html. The class comparisons for identifying differentially expressed (DE) genes between the PRRSV-infected DPL and DLY pigs were performed with two-sample t-test with q-value ≤5% to control false discovery rate and a fold change of 1.5 as a threshold [20]. We also utilized BRB-ArrayTools to obtain hierarchical clustering results and to graphically represent DE genes between the DPL and the DLY PRRSV-infected pigs. Gene ontology (GO) analyses were conducted using the NetAffx (http://www.affymetrix.com) and AgBase (http://www.agbase.msstate.edu/cgi-bin/tools/goslimviewer_select.pl) tools for comparative gene ontology categories, including molecular function (MF), biological process (BP) and cellular component (CC). The DE genes were also imported into the online software KEGG (http://www.genome.jp/kegg/pathway.html) for biological pathway mapping.
Real-time Quantitative Reverse Transcription PCR
Real-time qRT-PCR was employed to validate the microarray data analysis results and to compare mRNA expression levels of DE genes. Genes were selected based on the magnitude of changes in gene expression and its relevance to immune function, the following 10 genes were selected: adenosine monophosphate deaminase 3 (AMPD3), torsinA interacting protein 2 (TOR1AIP2), thymosin beta-15A-like (LOC100525848, THYMOSIN), PRVE (LOC100151993), C6ORF52 (LOC100270682, chromosome 6 open reading frame 52-like), HSPCB (LOC100396742), cytochrome P450 2J2-like (CYP2J2), coiled-coil domain-containing protein 84-like (LOC100516029, CCDC84), ACOX3 (LOC100513192), and the inactive progesterone receptor (PTGES3, LOC100155956). Furthermore, the expression levels of eight more genes AMPD3, ACOX3, CYP3A88, LPL (Lipoprotein lipase), MRC1 (mannose receptor, C type 1), USP18 (ubiquitin sqecific peptidase 18), TF (transferrin) and ENPEP (glutamyl aminopeptidase (aminopeptidase A)), were compared in both uninfected and infected DPL pigs with those in the uninfected and infected DLY pigs. The RNA samples prepared for microarray analysis were also used for real-time qRT-PCR validation. cDNA synthesis was conducted according to the manufacturer’s instruction and the primer sequences for these reactions are shown in Table 1. qRT-PCR was carried out with the Brilliant SYBR Green qRT-PCR Master Mix (Stratagene, La Jolla, CA, USA). All PCR reactions were performed in triplicate with negative controls. Using the GAPDH gene for normalization, the relative expression levels of the DE genes were analyzed by the 2−ΔΔCT method [21].
Statistical Analysis
CD4+ and CD8+ T cell numbers, and their ratios as well as the cytokine and IgG levels in both groups were repeated measures on the same pigs. These data sets are typical longitudinal data sets. We analyzed these data sets with linear mixed-effects model fit by restricted maximum likelihood procedure. This was done with lme function in R package nlme. The significance level of 0.05 was chosen for testing the difference between different pig breeds. Two-sample t-test and Wilcoxon rank test were employed to compare the mRNA expression levels of DE genes.
Results
Clinical Features of PRRSV-infected Pigs
The PRRSV-infected DLY pigs showed the following clinical symptoms within 3 days of infection: anorexia, depression, diarrhea, reluctance to drink, and wasting. At 14 dpi, most of the PRRSV-infected DLY pigs had signs of coughing, sneezing, conjunctivitis, eye discharge, increased respiratory rates, and constipation. Some DLY pigs exhibited shivering, spasms, cyanosis, and reddening of the ear skin after 20 dpi. On the contrary, the PRRSV-infected DPL pigs exhibited only mild depression and did not show any obvious clinical characteristics of the disease. The rectal temperatures of the DLY pigs revealed a persistent high fever (above 39.6°C) after 18 dpi, which was higher than in the DPL pigs.
CD4+ and CD8+ T lymphocytes
The percentage of CD4+ cells in the peripheral blood samples from both DPL and DLY pigs significantly increased from 0 dpi to 7 dpi (p<0.05), then decreased to lower levels at 14 and 21 dpi, which were not significantly different with that at 0 dpi (p>0.05) (Table S1). Overall, the DPL group had a lower percentage of CD4+ cells than the DLY group (p<0.05). In contrast, the percentage of CD8+ T cells in peripheral blood samples from both DPL and DLY pigs did not show significant difference (p>0.05); and it was significantly lower at 7 and 21 dpi than at 0 and 14 dpi (p<0.05) (Table S1). Furthermore, the DPL pigs had relatively lower CD4+/CD8+ratios than the DLY pigs (p<0.05), the CD4+/CD8+ratio at 0, 7 and 21 dpi did not show significant difference (p>0.05) and was significantly lower than that at 14 dpi (p<0.01) (Table S1).
Cytokines and IgG
During 0–14 days after PRRSV infection, the change of IL-1β and IL-2 was not significantly different between the DPL and the DLY pigs (p>0.05) although both IL-1β and IL-2 levels were much higher at 7 dpi than those at 0 (p<0.01) and 14 dpi (Table S2). For both IL-10 and TNF-α, the DLY pigs had significantly higher levels than the DPL pigs (p<0.01); for both the DPL and the DLY pigs, their IL-10 and TNF-α levels increased from 0 dpi to 7 dpi (p<0.001), decreased from 7 dpi to 14 dpi, at which the levels were significantly lower than those at 0 dpi (p<0.01) (Table S2). In contrast, serum IFN-γ levels fell from 0 dpi to 7 dpi (p<0.01), then increased at 14 dpi; significant differences between DPL and DLY pigs were detected across the whole experiment (p<0.01), more specifically, the DLY pigs have lower serum IFN-γ levels than the DPL pigs (p<0.01) (Table S2). Overall, the serum IgG levels increased slightly from 0 dpi to 7 dpi, which is not statistically significant (p>0.05), and peaked at 14 dpi (p<0.0001) (Table S2). In addition, the PRRSV-infected DLY pigs had slightly higher IgG levels than the PRRSV-infected DPL pigs during 0–14 dpi (Table S2); however, no statistically significant differences were detected between the DPL and the DLY pigs (p>0.05).
Microarray Gene Expression Analysis
Lung samples collected from six of the PRRSV-infected pigs (three DPL and three DLY) at 28 dpi were analyzed by using a porcine Affymetrix GeneChip. First, the gene probe sets were filtered to eliminate those with very low expression summary values and low variability across all of the samples. The remaining 1419 transcripts were used for subsequent analysis (Table S3). After quantile normalization, 16 genes showed differential expression at the q≤5% level (Table 2, Figure 1), among which two genes (LOC100516029 and LOC100523005) were up-regulated by a fold change of 2.63 and 2.47 in the infected DPL pigs compared to those in the infected DLY pigs, respectively. The other fourteen genes were down-regulated in the infected DPL pigs compared with the infected DLY pigs.
Each point represents a gene, and colors correspond to the ranges of the negative log10 P and log2 fold change values. Blue circles: differentially expressed genes. Fold changes >0 indicate up-regulated genes, whereas fold changes <0 indicate down-regulated genes. Black circles: no statistically significant expressed genes.
Cluster Analyses of DE Genes
Clustering was used to define sets of genes with similar responses to viral infection at the q≤0.005 level. Clustering analysis of gene expression profiles for DPL versus DLY pigs post-infection showed that most of the 16 genes from the DPL pigs had lower expression levels when compared with those from the DLY pigs (Figure 2).
Each column represents one pig, and each horizontal line refers to one gene. The three columns on the left and right represent the DPL and the DLY pigs, respectively. The color bar at the bottom of the figure indicates the expression level of the genes, those in the lightest blue have lower expression relative to the geometrical means, while dark blue indicates genes with higher expression relative to the geometrical means.
GO Annotation of DE Genes
GO analysis revealed 39 classes of genes, many of which share common transcripts and are grouped into 20 BP, 13 MF and 6 CC (Table 3, Table S4). Most genes in the BP class are related to metabolic processes, but catalytic activity and binding are only represented in the MF group.
Pathway Analysis
Pathway analysis of the DE genes was performed using the KEGG database (http://www.genome.jp/kegg/pathway.html). The results revealed eight main pathway categories (Figure S1): ERBB signaling, endocytosis, JAK-STAT signaling, T cell receptor signaling, regulation of actin cytoskeleton, insulin signaling, bacterial invasion of epithelial cells and chronic myeloid leukemia.
Validation of the Microarray Data by qRT-PCR
We selected 10 of the DE genes identified by microarray analysis for real-time qRT-PCR analysis, including one up-regulated gene (CCDC84, which encodes coiled-coil domain-containing 84) and nine down-regulated genes (C6ORF52, THYMOSIN, PRVE, HSPCB, CYP2J2, AMPD3, TOR1AIP2, PTGES3, ACOX3). The fold change in the mRNA expression levels of the 10 genes between the PRRSV-infected DPL and DLY pigs as determined by qRT-PCR was consistent with the microarray results (Figure 3), suggesting that microarray hybridization is a reliable method for genome-wide analysis of expression profiles in the PRRSV-infected DPL and DLY pigs.
The p-values of the qRT-PCR data are as follows: 0.0067 (C6ORF52), 0.3920 (CCDC84), 0.3661 (THYMOSIN), 0.8195 (PRVE), 0.3154 (HSPCB), 0.6534 (CYP2J2), 0.002 (AMPD3), 0.1622 (TOR1AIP2), 0.8329 (PTGES3), 0.8808 (ACOX3).
Comparison of DE Gene mRNA Expressions in Uninfected and Infected DPL and DLY Pigs
In addition, we also analyzed the mRNA expression level of 8 out of 67 genes which are differentially expressed at significance level 0.01 between the DPL and DLY pigs for both uninfected and infected groups. The results showed that AMPD3 gene expression level for both uninfected and infected DPL pigs was significantly lower compared with the uninfected and infected DLY pigs (p<0.05, p<0.01, Figure 4A). For ACOX3 gene, it did not show significant difference between uninfected DPL and DLY pigs or between infected DPL and DLY pigs (p>0.05, Figure 4B). Moreover, three down-regulated genes (CYP3A88, LPL and MRC1) behaved similarly in uninfected pigs, that is, their gene expression levels of uninfected DLY pigs were significantly higher than those of uninfected DPL pigs (p<0.01, Figure 4C, D and E), and no significant difference of their mRNA expression levels was detected between infected DLY and DPL pigs (p>0.05, Figure 4C, D and E). In uninfected pigs, the mRNA expression level of TF gene was significantly higher in DLY pigs than that in DPL pigs (p<0.05, Figure 4F); whereas in infected pigs, it was significantly lower in DLY than that in DPL pigs (p<0.05, Figure 4F). USP18 mRNA levels were not significantly different in uninfected groups (p>0.05) and were significantly lower in infected DLY pigs than those in infected DPL pigs (p<0.05, Figure 4G); ENPEP mRNA expression levels did not show significant difference in uninfected/infected groups of DPL and DLY pigs (p>0.05, Figure 4H).
* or ** indicates a statistically significant difference at p<0.05 or 0.01, respectively.
Discussion
PRRS is pandemic in pigs and causes significant losses to the pig industry via reproductive disorders and growth retardation. In particular, a highly pathogenic disease caused by PRRSV emerged in swine farms all over China from 2006 to 2010, which caused continuous high fever and a high proportion of deaths in vaccinated pigs of all ages [19], [22], [23]. Studies have shown that responses of different pig breeds and lines to PRRSV infection differ greatly [8], [9], [10], [11], [12], [24], [25], [26]. Nevertheless, the molecular mechanisms underlying PRRSV remain largely unknown. In the present study, using Affymetrix GeneChip technology, we have reported the results of a genome-wide comparison of gene expression profiles in lung tissues from two pig breeds differing in resistance to PRRSV. This is the first report of differential expression at the q≤5% level for 16 genes that differ between two breeds (Table 2, Figure 1). Observations of interest are the differential responses to PRRSV infection in gene expression, clinical symptoms and biochemical measurements of CD4+, CD8+, cytokine and IgG profiles between the DPL and the DLY pigs.
The PRRSV-infected DPL and DLY lines showed marked differences in body weight, rectal temperature and serum cytokine levels [10], [11], [27]. Vincent et al. (2006) reported that the serum viral titers at 21 dpi differed among different pig lines after infection with PRRSV [26]. Variation in resistance to PRRSV in Pietrain and Miniature pigs was reported by Reiner et al. [28]. In the present study, we also found differences between DPL and DLY pigs in their clinical features of PRRSV infection: DLY pigs had more severe clinical symptoms and lesions than the DPL pigs, which implied that DPL pigs had a stronger resistance phenotype than DLY pigs when challenged with PRRSV. These results indicate that underlying genetic variation exists in disease susceptibility to PRRSV.
CD4+ is a co-receptor that assists the T cell receptor (TCR) to communicate with an antigen-presenting cell [29], while CD8+ is a transmembrane glycoprotein that serves as a co-receptor for the TCR [30]. Consequently, CD4+ and CD8+ T cells play important roles in protection against viral infections. Shimizu et al. observed a remarkable decrease in CD4+ T cells after 3 days infection with PRRSV [31]; and this study also reported slight decreases in CD8+ T cells at 3 dpi, followed by substantially increased levels [31], while at the same time, the ratios of CD4+/CD8+ T cells were significantly lower between day 3 and 28 post-inoculation compared with day 0 [31]. A different study revealed that the percentages of CD4+ and CD8+ T cells were significantly decreased for a few days shortly after PRRSV infection, but returned to pre-infection levels on days 8–10 post-inoculation [32]. Comparing with the results of Shimizu et al. [31] and Nielsen and Bøtner [32], we observed that CD4+ T cell numbers increased at 7 dpi, then decreased to its normal range; and the CD8+ T cell numbers was lower at 7, 21 dpi than that at 0, 14 dpi. Furthermore, the CD4+/CD8+ratio was relatively stable from 0 dpi to 21 dpi although it had low value at 14 dpi (Table S1).
Cytokines, including IL-1β, IL-2, IL-8, IL-10, IFN-γ, and TNF-α, have been shown to play important roles in defense against PRRSV [33]. Increased levels of IL-1β, IL-8, IL-10 and TNF-α were present in the porcine PBMC culture supernatant of the PRRSV-infected group as measured by ELISA [34], [35]. Similar findings were reported for IL10, IFN-γ and TNF-α, the levels of which were elevated in the blood of PRRSV-infected pigs [36], [37]. We also observed remarkable increases in IL-1β, IL-2, IL-10 and TNF-α at 7 dpi after which IL-1β and IL-2 dropped to pre-infection levels, IL-10 and TNF-α decreased to lower levels than pre-infection levels. Overall, the IL-10, TNF-α and IFN-γ cytokines in DLY pigs were different with those in DPL pigs (Table S2). They are probably caused by fast PRRSV replication in the DLY pigs and more cytokines production to resist PRRSV replication for DPL pigs, a finding that requires further investigation.
Genome-wide comparisons of the microarray gene expression profiles revealed many complex biologic processes involved in PRRSV infection. By using Affymetrix microarrays, 1409 DE transcripts were identified by analysis of variance, of which 2, 5, 25, 16 and 100 differed between the two breeds by a minimum of 1.5-fold at 1, 3, 6, 9 and 12 h pi, respectively [17]. A study investigating interactions between PRRSV and porcine alveolar macrophages (PAMs) showed that highly pathogenic PRRSV infection affected PAM expression of important genes [18]. The present study identified 16 DE genes using a 1.5-fold change filter, and the number of down-regulated genes (14) was greater than that of the up-regulated genes (2) (Table 2). Cluster analysis of DE genes revealed that most of the 16 genes from the PRRSV-infected DPL pigs had low expression values in comparison with the PRRSV-infected DLY pigs, indicating that they had distinctly different responses to PRRSV infection, thus implying that DPL pigs are more resistant to PRRSV than DLY pigs.
In mice, remote reperfusion lung injury is associated with AMP deaminase 3 activation and attenuation by inosine monophosphate [38], showing that AMPD3 plays a critical role in remote reperfusion lung injury in mice [39]. In this study, the lower AMPD3 expression in DPL pigs compared with their DLY counterparts indicates that the DPL pig lungs were less affected by PRRSV. TorsinA interacting protein 2 (TOR1AIP2), also called LULL1, is a single-pass membrane protein located on the endoplasmic reticulum membrane and nuclear membrane and functions to regulate the distribution of TOR1A between the endoplasmic reticulum and the nuclear envelope (http://www.genecards.org/). Functional defects in the mutant protein are responsible for DYT1 dystonia [40]; however, the role of TOR1AIP2 in lung development is still unknown. Coiledcoil domains typically function in homodimerization and are present in a number of proteins involved in intracellular transport [41]. In this study, we found that expression of the coiled-coil domain-containing protein 84-like (CCDC84) mRNA was 2.63-fold higher in PRRSV-infected DPL pigs than DLY pigs. The role of CCDC84 in resistance to PRRSV requires further study. The other 13 DE genes identified are neither annotated nor encode hypothetical proteins; hence, their functions remain unknown.
DNA microarrays are widely used for measuring the expression levels of large numbers of genes simultaneously; however, microarray results for any given gene are often noisy or ambiguous [42]. Many external factors can significantly affect the accuracy of array data, including RNA preparation, image acquisition, methods for normalization and background subtraction, data processing and standardization, use of visualization tools, and probes used for labeling [43]. At present, real-time qRT-PCR is thought to be the best technique for validating microarray data. However, not all microarray data are verified by this technique. It has been reported that real-time qRT-PCR was able to confirm the change in expression for 17 out of 24 genes identified by high-density filter arrays, and genes showing less than a four-fold difference on the array were not verified by real-time qRT-PCR [44]. A later study revealed that 10 of the genes (10 versus 33) tested that had differential expression in the microarray experiment were not confirmed by the real-time qRT-PCR [45]. Similarly, in our work, 10 out of 16 genes differentially expressed between the PRRSV-infected DPL and DLY pigs as determined by qRT-PCR were consistent with the microarray results. These contradictions between the microarray data and qRT-PCR results probably arise from the greater accuracy of the latter technique compared with microarrays or the lack of specificity in the primer design.
To investigate if the differential expression between the two breeds of infected pigs was confounded with genetic background, 8 genes (6 down-regulated genes and 2 up-regulated genes) were further selected for real-time qRT-PCR analysis in both uninfected and infected groups. The results indicated that the mRNA expression of AMPD3, CYP3A88, LPL, MRC1 and TF genes was significantly different between uninfected DPL and DLY pigs (Figure 4A, C, D, E and F), suggesting that these differential gene expression levels may be also caused by the genetic differences between the two breeds. However, there were genes showed significant changes at expression level between uninfected and infected pigs within breed. In both uninfected and infected groups, the expression differences of AMPD3 between DLY and DPL pigs are similar, suggesting that the difference is not caused by PRRSV infection. The expression of TF gene is significantly lower in DPL pigs than in DLY pigs in uninfected group; however, in infected group, it is much higher in DPL pigs than in DLY pigs (Figure 4F), suggesting that TF gene is likely involved in PRRSV resistance. Similarly, one study indicated that the expression of TF gene was down regulated in high (HR) pigs, with high viremia, low/no weight gain, and many lung lesions [46]. The increase in the mRNA expression in infected DPL pigs compared to uninfected ones suggests a role of TF gene in PRRSV resistance. The expression of CYP3A88, LPL and MRC1 is not significantly different between infected DPL and DLY pigs, indicating that their response to PRRSV infection is not different and they may not control the resistance of DPL pigs to PRRSV infection. The expression of ACOX3 and ENPEP genes were not significantly different between either uninfected or infected DPL and DLY pigs, suggesting that they are not related to PRRSV infection. Studies demonstrated that overexpression of the porcine USP18 leads to reduced replication and/or growth of PRRSV [47] and USP18 restricts PRRSV growth through alteration of nuclear translocation of NF-κB p65 and p50 in MARC-145 cells [48]. In this study, we also noted a slight increase of USP18 mRNA expression level in uninfected DPL pigs compared to uninfected DLY pigs and significantly higher expression level of USP18 gene in infected DPL than that in infected DLY pigs, implying that it also plays an important role in the resistance of DPL pigs to PRRSV infection (Figure 4G).
Conclusion
In this study, we found that PRRSV-infected DPL pigs showed mild clinical symptoms and had a lower percentage of CD4+ T cells. Moreover, the levels of IL-10, TNF-α and IFN-γ also differed between the PRRSV-infected DPL and DLY pigs. Using Affymetrix microarray chip technology, we compared the gene expression profiles of lung tissues in DPL and DLY pigs after infection with PRRSV and identified sixteen DE genes. We further analyzed the mRNA expression level of 8 out of 67 genes which are differentially expressed at significance level 0.01 between the DPL and DLY pigs for both uninfected and infected groups, and found that TF and USP18 genes were important in underlying porcine resistance or susceptibility to PRRSV.
Supporting Information
Figure S1.
Signaling pathways of DE genes. Pathway analysis was mainly based on the KEGG database. A p-value of <0.05 and an FDR of <0.05 in the two-side Fisher’s exact test were selected as the statistical significance criteria.
https://doi.org/10.1371/journal.pone.0086101.s001
(PDF)
Table S1.
Statistical analysis for the percentage of CD4+, CD8+ cells and CD4+/CD8+ ratio.
https://doi.org/10.1371/journal.pone.0086101.s002
(DOC)
Table S2.
Statistical analysis for Cytokines and IgG.
https://doi.org/10.1371/journal.pone.0086101.s003
(DOC)
Table S3.
1419 probe set after filtered by 5% discovery rate and a fold change of 1.5.
https://doi.org/10.1371/journal.pone.0086101.s004
(XLS)
Table S4.
GO annotation of differentially expressed genes.
https://doi.org/10.1371/journal.pone.0086101.s005
(XLS)
Author Contributions
Conceived and designed the experiments: JX FX YJ. Performed the experiments: FX YZ YL CJ JX. Analyzed the data: JX FX LY NW CH. Zhang CY. Zhang CW. Contributed reagents/materials/analysis tools: FX YL CJ LY CH. Zhang NW. Wrote the paper: JX YJ CH. Zhang.
References
- 1. Collins JE, Benfield DA, Christianson WT, Harris L, Hennings JC, et al. (1992) Isolation of swine infertility and respiratory syndrome virus (isolate ATCC VR-2332) in North America and experimental reproduction of the disease in gnotobiotic pigs. J Vet Diagn Invest 4: 117–126.
- 2. Botner A (1997) Diagnosis of PRRS. Vet Microbiol 55: 295–301.
- 3. Done SH, Paton DJ (1995) Porcine reproductive and respiratory syndrome: clinical disease, pathology and immunosuppression. Vet Rec 136: 32–35.
- 4. Dee SA, Joo HS, Polson DD, Marsh WE (1997) Evaluation of the effects of nursery depopulation on the profitability of 34 pig farms. Vet Rec 140: 498–500.
- 5. Neumann EJ, Kliebenstein JB, Johnson CD, Mabry JW, Bush EJ, et al. (2005) Assessment of the economic impact of porcine reproductive and respiratory syndrome on swine production in the United States. J Am Vet Med Assoc 227: 385–392.
- 6. Cavanagh D (1997) Nidovirales: a new order comprising Coronaviridae and Arteriviridae. Arch Virol 142: 629–633.
- 7. Lunney JK, Steibel JP, Reecy JM, Fritz E, Rothschild MF, et al. (2011) Probing genetic control of swine responses to PRRSV infection: current progress of the PRRS host genetics consortium. BMC Proc 5 Suppl 4S30.
- 8. Lunney JK, Chen H (2010) Genetic control of host resistance to porcine reproductive and respiratory syndrome virus (PRRSV) infection. Virus Res 154: 161–169.
- 9. Halbur PG, Rothschild MF, Thacker BJ, Meng XJ, Paul PS, et al. (1998) Differences in susceptibility of Duroc, Hampshire, and Meishan pigs to infection with a high virulence strain (VR2385) of porcine reproductive and respiratory syndrome virus (PRRSV). J Anim Breed Genet 115: 181–189.
- 10. Petry DB, Holl JW, Weber JS, Doster AR, Osorio FA, et al. (2005) Biological responses to porcine respiratory and reproductive syndrome virus in pigs of two genetic populations. J Anim Sci 83: 1494–1502.
- 11. Petry DB, Lunney J, Boyd P, Kuhar D, Blankenship E, et al. (2007) Differential immunity in pigs with high and low responses to porcine reproductive and respiratory syndrome virus infection. J Anim Sci 85: 2075–2092.
- 12. Jiang C, Xing F, Xing J, Jiang Y, Zhou E (2013) Different expression patterns of PRRSV mediator genes in the lung tissues of PRRSV resistant and susceptible pigs. Dev Comp Immunol 39: 127–131.
- 13. Schena M, Shalon D, Davis RW, Brown PO (1995) Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science 270: 467–470.
- 14. Lockhart DJ, Winzeler EA (2000) Genomics, gene expression and DNA arrays. Nature 405: 827–836.
- 15. Lee C, Bachand A, Murtaugh MP, Yoo D (2004) Differential host cell gene expression regulated by the porcine reproductive and respiratory syndrome virus GP4 and GP5 glycoproteins. Vet Immunol Immunopathol 102: 189–198.
- 16. Bates JS, Petry DB, Eudy J, Bough L, Johnson RK (2008) Differential expression in lung and bronchial lymph node of pigs with high and low responses to infection with porcine reproductive and respiratory syndrome virus. J Anim Sci 86: 3279–3289.
- 17. Genini S, Delputte PL, Malinverni R, Cecere M, Stella A, et al. (2008) Genome-wide transcriptional response of primary alveolar macrophages following infection with porcine reproductive and respiratory syndrome virus. J Gen Virol 89: 2550–2564.
- 18. Zhou P, Zhai S, Zhou X, Lin P, Jiang T, et al. (2011) Molecular characterization of transcriptome-wide interactions between highly pathogenic porcine reproductive and respiratory syndrome virus and porcine alveolar macrophages in vivo. Int J Biol Sci 7: 947–959.
- 19. Tian K, Yu X, Zhao T, Feng Y, Cao Z, et al. (2007) Emergence of fatal PRRSV variants: unparalleled outbreaks of atypical PRRS in China and molecular dissection of the unique hallmark. PLoS One 2: e526.
- 20. Patterson TA, Lobenhofer EK, Fulmer-Smentek SB, Collins PJ, Chu TM, et al. (2006) Performance comparison of one-color and two-color platforms within the MicroArray Quality Control (MAQC) project. Nat Biotechnol 24: 1140–1150.
- 21. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods 25: 402–408.
- 22. Tong GZ, Zhou YJ, Hao XF, Tian ZJ, An TQ, et al. (2007) Highly pathogenic porcine reproductive and respiratory syndrome, China. Emerg Infect Dis 13: 1434–1436.
- 23. Li B, Fang L, Guo X, Gao J, Song T, et al. (2011) Epidemiology and evolutionary characteristics of the porcine reproductive and respiratory syndrome virus in China between 2006 and 2010. J Clin Microbiol 49: 3175–3183.
- 24. Lowe JE, Husmann R, Firkins LD, Zuckermann FA, Goldberg TL (2005) Correlation of cell-mediated immunity against porcine reproductive and respiratory syndrome virus with protection against reproductive failure in sows during outbreaks of porcine reproductive and respiratory syndrome in commercial herds. J Am Vet Med Assoc 226: 1707–1711.
- 25. Vincent AL, Thacker BJ, Halbur PG, Rothschild MF, Thacker EL (2005) In vitro susceptibility of macrophages to porcine reproductive and respiratory syndrome virus varies between genetically diverse lines of pigs. Viral Immunol 18: 506–512.
- 26. Vincent AL, Thacker BJ, Halbur PG, Rothschild MF, Thacker EL (2006) An investigation of susceptibility to porcine reproductive and respiratory syndrome virus between two genetically diverse commercial lines of pigs. J Anim Sci 84: 49–57.
- 27. Doeschl-Wilson AB, Kyriazakis I, Vincent A, Rothschild MF, Thacker E, et al. (2009) Clinical and pathological responses of pigs from two genetically diverse commercial lines to porcine reproductive and respiratory syndrome virus infection. J Anim Sci 87: 1638–1647.
- 28. Reiner G, Willems H, Pesch S, Ohlinger VF (2010) Variation in resistance to the porcine reproductive and respiratory syndrome virus (PRRSV) in Pietrain and Miniature pigs. J Anim Breed Genet 127: 100–106.
- 29. Brady RL, Dodson EJ, Dodson GG, Lange G, Davis SJ, et al. (1993) Crystal structure of domains 3 and 4 of rat CD4: relation to the NH2-terminal domains. Science 260: 979–983.
- 30. Gao GF, Jakobsen BK (2000) Molecular interactions of coreceptor CD8 and MHC class I: the molecular basis for functional coordination with the T-cell receptor. Immunol Today 21: 630–636.
- 31. Shimizu M, Yamada S, Kawashima K, Ohashi S, Shimizu S, et al. (1996) Changes of lymphocyte subpopulations in pigs infected with porcine reproductive and respiratory syndrome (PRRS) virus. Vet Immunol Immunopathol 50: 19–27.
- 32. Nielsen J, Botner A (1997) Hematological and immunological parameters of 4 1/2-month old pigs infected with PRRS virus. Vet Microbiol 55: 289–294.
- 33. Murtaugh MP, Xiao Z, Zuckermann F (2002) Immunological responses of swine to porcine reproductive and respiratory syndrome virus infection. Viral Immunol 15: 533–547.
- 34. Thanawongnuwech R, Young TF, Thacker BJ, Thacker EL (2001) Differential production of proinflammatory cytokines: in vitro PRRSV and Mycoplasma hyopneumoniae co-infection model. Vet Immunol Immunopathol 79: 115–127.
- 35. Zhuge ZY, Zhu YH, Liu PQ, Yan XD, Yue Y, et al. (2012) Effects of Astragalus polysaccharide on immune responses of porcine PBMC stimulated with PRRSV or CSFV. PLoS One 7: e29320.
- 36. Xiao S, Mo D, Wang Q, Jia J, Qin L, et al. (2010) Aberrant host immune response induced by highly virulent PRRSV identified by digital gene expression tag profiling. BMC Genomics 11: 544.
- 37. Dwivedi V, Manickam C, Binjawadagi B, Linhares D, Murtaugh MP, et al. (2012) Evaluation of immune responses to porcine reproductive and respiratory syndrome virus in pigs during early stage of infection under farm conditions. Virol J 9: 45.
- 38. Li P, Ogino K, Hoshikawa Y, Morisaki H, Cheng J, et al. (2007) Remote reperfusion lung injury is associated with AMP deaminase 3 activation and attenuated by inosine monophosphate. Circ J 71: 591–596.
- 39. Li P, Ogino K, Hoshikawa Y, Morisaki H, Toyama K, et al. (2013) AMP deaminase 3 plays a critical role in remote reperfusion lung injury. Biochem Biophys Res Commun 434: 131–136.
- 40. Vander Heyden AB, Naismith TV, Snapp EL, Hodzic D, Hanson PI (2009) LULL1 retargets TorsinA to the nuclear envelope revealing an activity that is impaired by the DYT1 dystonia mutation. Mol Biol Cell 20: 2661–2672.
- 41. Burkhard P, Stetefeld J, Strelkov SV (2001) Coiled coils: a highly versatile protein folding motif. Trends Cell Biol 11: 82–88.
- 42.
Huang TT, Deoghare HV, Smith BK, Beaver TM, Baker HV, et al.. (2011) Gene expression changes in the human diaphragm after cardiothoracic surgery. J Thorac Cardiovasc Surg 142: 1214–1222, 1222 e1211–1220.
- 43. Chuaqui RF, Bonner RF, Best CJ, Gillespie JW, Flaig MJ, et al. (2002) Post-analysis follow-up and validation of microarray experiments. Nat Genet 32 Suppl: 509–514
- 44. Rajeevan MS, Vernon SD, Taysavang N, Unger ER (2001) Validation of array-based gene expression profiles by real-time (kinetic) RT-PCR. J Mol Diagn 3: 26–31.
- 45. Zhao SH, Kuhar D, Lunney JK, Dawson H, Guidry C, et al. (2006) Gene expression profiling in Salmonella Choleraesuis-infected porcine lung using a long oligonucleotide microarray. Mamm Genome 17: 777–789.
- 46. Wysocki M, Chen H, Steibel JP, Kuhar D, Petry D, Bates J, et al. (2012) Identifying putative candidate genes and pathways involved in immune responses to porcine reproductive andrespiratory syndrome virus (PRRSV) infection. Anim Genet 43: 328–332.
- 47. Ait-Ali T, Wilson AW, Finlayson H, Carré W, Ramaiahgari SC, et al. (2009) Functional analysis of the porcine USP18 and its role during porcine arterivirus replication. Gene 439: 35–42.
- 48. Xu D, Lillico SG, Barnett MW, Whitelaw CB, Archibald AL, et al. (2012) USP18 restricts PRRSV growth through alteration of nuclear translocation of NF-κB p65 and p50 in MARC-145 cells. Virus Res 169: 264–267.