key: cord-012909-o6t2srim
authors: Chaudhari, Jayeshbhai; Liew, Chia-Sin; Workman, Aspen M.; Riethoven, Jean-Jack M.; Steffen, David; Sillman, Sarah; Vu, Hiep L. X.
title: Host Transcriptional Response to Persistent Infection with a Live-Attenuated Porcine Reproductive and Respiratory Syndrome Virus Strain
date: 2020-07-28
journal: Viruses
DOI: 10.3390/v12080817
sha: 
doc_id: 12909
cord_uid: o6t2srim

Both virulent and live-attenuated porcine reproductive and respiratory syndrome virus (PRRSV) strains can establish persistent infection in lymphoid tissues of pigs. To investigate the mechanisms of PRRSV persistence, we performed a transcriptional analysis of inguinal lymphoid tissue collected from pigs experimentally infected with an attenuated PRRSV strain at 46 days post infection. A total of 6404 differentially expressed genes (DEGs) were detected of which 3960 DEGs were upregulated and 2444 DEGs were downregulated. Specifically, genes involved in innate immune responses and chemokines and receptors associated with T-cell homing to lymphoid tissues were down regulated. As a result, homing of virus-specific T-cells to lymphoid tissues seems to be ineffective, evidenced by the lower frequencies of virus-specific T-cell in lymphoid tissue than in peripheral blood. Genes associated with T-cell exhaustion were upregulated. Likewise, genes involved in the anti-apoptotic pathway were upregulated. Collectively, the data suggested that the live-attenuated PRRSV strain establishes a pro-survival microenvironment in lymphoid tissue by suppressing innate immune responses, T-cell homing, and preventing cell apoptosis.

Porcine reproductive and respiratory syndrome virus (PRRSV) is a positive sense, single stranded RNA virus that belongs to the family Arteriviridae, under the order Nidovirales [1] . Based on phylogenetic analysis, PRRSV was originally classified into two types: type 1 or PRRSV-1, which originated in Europe, and type 2 or PRRSV-2, which originated in North America. The International Committee on Taxonomy of Viruses (ICTV) recently updated the arterivirus taxonomic structure in which PRRSV-1 and PRRSV-2 are now respectively classified as two species: Betaarterivirus Suid 1 and Betaarterivirus Suid 2 [2] . PRRSV infects pigs of all ages; however, clinical manifestations are more severe when the virus infects pregnant sows and young pigs, causing reproductive failure and respiratory distress, respectively (reviewed in [3] ). PRRSV is endemic in most swine producing countries worldwide, causing significant economic losses to swine producers [4] .

Genes involved in innate immune response and genes encoding chemokines and receptors important for T-cell homing and trafficking were downregulated. On the other hand, genes involved in the anti-apoptotic pathway and T-cell exhaustion were upregulated. Functional studies revealed that the frequencies of virus-specific IFN-γ-secreting cells are lower in lymphoid tissue than in peripheral blood mononucleated cells (PBMCs). Collectively, the results shed important insight into the mechanisms of PRRSV persistence in the host.

The attenuated PRRSV strain designated CON90 used in this study was recovered from an infectious cDNA clone constructed using viral RNA extracted from the attenuated CON-P90 [27] . MARC-145 cells, a monkey kidney cell line [28] , were cultured in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and 1 × penicillin-streptomycin (100 units/mL of penicillin, and 100 µg/mL of streptomycin (Life Technologies, Grand Island, NY, USA). PBMCs and lymph node-derived cells were cultured in complete RPMI-1640 medium (cRPMI) supplemented with 10% FBS, 1 × of GlutaMax-1 (Life Technologies, Grand Island, NY, USA), 100 units/mL of penicillin, and 100 µg/mL of streptomycin (Sigma-Aldrich, St. Louis, MO, USA). Mouse monoclonal antibody specific to PRRSV N protein SDOW17 was purchased from the National Veterinary Services Laboratories (Ames, IA, USA). Anti-porcine CD3ε (clone BB23-8E6-8C8; FITC-conjugated), anti-porcine CD4 (clone 74-12-4; PECy7 -conjugated), anti-porcine CD8 (clone 76-2-11; Alexa Flour 647 -conjugated), and anti-porcine IFN-γ (clone P2G10; PE-conjugated), were purchased from BD Biosciences (San Diego, CA, USA). Alexa Flour 488-conjugated goat anti-mouse IgG was purchased from Invitrogen (Eugen, OR, USA).

The pig experiment conducted in this study was approved by the University of Nebraska-Lincoln (UNL) Institutional Animal Care and Use Committee under the protocol number 1360. Ten 3-week-old PRRSV, porcine circovirus type 2 (PCV2), and SIV negative pigs were purchased from the Midwest Research Swine (Glencoe, MN, USA). The pigs were randomly assigned to two groups of five pigs, which were accommodated in two separate rooms in the animal biosecurity level 2 (ABSL-2) research facilities at UNL. After 1 week of acclimation, pigs in group 1 were injected with DMEM to serve as negative controls, whereas pigs in groups 2 were inoculated intramuscularly with 10 5.0 TCID50 of live-attenuated PRRSV strain CON90. Whole blood samples with anticoagulant EDTA were collected to obtain plasma and PBMCs. Plasma samples were stored at −70 • C for measurement of viremia and antibody response. A portion of freshly isolated PBMCs were used for measurement of T cell responses while the rest of PBMCs were cryopreserved for future analysis. At 46 days post-infection (dpi), sample of inguinal lymph node (ILN) was aseptically collected from the pigs under anesthesia. One half of the INL was stored in cRPMI for lymphocyte isolation while the other half of the INL was minced and stored in TRIzol Reagent (Life Technologies, Carlsbad, CA, USA) for RNA extraction.

Viral loads in plasma and tissues were measured by a commercial RT-qPCR kit (Tetracore Inc., Rockville, MD, USA). Viral loads in plasma was reported as log 10 copies per ml, whereas viral loads in tissues were reported as log 10 copy per µg of total RNA used in the RT-qPCR reaction. For statistical purposes, samples that had no detectable levels of viral RNA were assigned a value of 0 log 10 copies.

PBMCs were isolated from EDTA-whole blood as previously described [29, 30] . Single cell suspension was isolated from ILN immediately after collection as follows. The tissue was processed to Viruses 2020, 12, 817 4 of 22 remove connecting tissue and cut into small pieces which were placed in a 70-µm nylon cell strainer (Corning, Durham, NC, USA) in the presence of cRPMI. The tissue pieces were pressed against the nylon mesh by using the plunger of a 3 mL syringe. The resulting cell suspension was collected into a 50 mL conical tube which was passed through a 70-µm cell strainer one more time to remove large tissue debris. Cells were pelleted by centrifugation at 700× g for 10 min at room temperature and treated with a 5 mL red blood cell (RBC) lysis buffer (Life Technologies, Carlsbad, CA, USA). RBC lysis reaction was ceased by adding ice cold PBS containing 4% FBS, followed by centrifugation at 700× g for 10 min at room temperature. Cell pellet was resuspended in cRPMI. To determine cell concentration and viability, samples of PBMC and ILN-derived cells were stained with acridine orange and propidium iodide (ViaStainTM AOPI Staining Solution, Nexcelom, Lawrence, MA, USA) and counted using an automatic cell counter (Cellometer Auto 2000, Nexcelom, Lawrence, MA, USA). Freshly isolated cells were used for Elispot and flow cytometric analysis. The remaining cells were cryopreserved in 10% DMSO, 40% FBS, and 50% RPMI-1640 and stored in liquid nitrogen.

PRRSV antibody levels in plasma were measured at the Veterinary Diagnostic Center of the University of Nebraska by using the commercial ELISA IDEXX PRRS X3 Ab Test (IDEXX Laboratories, Westbrook, ME, USA) following the manufacturer's instructions. The serum-virus neutralization (SVN) assay was performed as previously described [31] using plasma rather than serum. Results were expressed as the log 2 of the reciprocal of the highest dilution that showed a ≥90% reduction in the number of fluorescent foci presenting in the control wells.

The frequencies of IFN-γ-secreting cells (IFN-γ SCs) in PBMCs and ILN-derived cells were measured by using an IFN-γ Elispot assay as previously described [32, 33] . Briefly two replicates of 250,000 PBMCs or ILN cells freshly collected from each pig were plated into two wells of a 96-well plate with PVDF membrane that were coated with anti-porcine IFN-γ antibody. The cells were stimulated with 100 µL of cRPMI containing 2.5 × 10 4 TCID50 of CON90. For positive control, cells were plated at 5000 cells per well, followed by stimulation with a 100 µL cRPMI containing 10 ng/mL of phorbol 12 myristate 13-acetate (PMA) and 1 µg/mL of ionomycin. For negative control, cells were simply cultured in cRPMI. At 18 h post stimulation, the plate was washed with PBS-containing 0.05% Tween 20 (PBS-T) followed by incubation with biotin-labeled antibody against porcine IFN-γ (clone P2C11, BD Biosciences Pharmingen, San Diego, CA, USA). Spots were developed by using alkaline phosphatase-conjugated streptavidin (Southern Biotech, Birmingham, AL, USA) in conjunction with alkaline phosphatase substrate (Vector laboratories, Burlingame, CA, USA). Spots were counted and analyzed using an AID Elispot Reader Version 7.0 (AID GmbH, Strassberg, Germany).

Freshly isolated PBMCs and ILN cells were ex vivo stimulated with 1 × 10 5 TCID50 of CON90 as described above. A Cocktail solution consisting of 10 ng/mL of PMA and 1 µg/mL of ionomycin was included as a positive control while cRPMI served as a negative control. At 12 hrs post-stimulation, 100 µL of cRPMI containing 1 µg/mL of Golgi-plug brefeldin A (BD Bioscience, San Jose, CA, USA) was added and further incubated for 6 hrs. Samples were centrifuged at 500× g for 10 min at room temperature. Cells were resuspended in flow cytometry staining buffer (FACS buffer) (PBS with 4% FBS and 0.1% sodium azide), and stained with anti-porcine CD3ε, CD4, and CD8 monoclonal antibodies and incubated on ice for 30 min in the dark. Cells were washed thrice using FACS buffer. The cells were fixed and permeabilized with 4% paraformaldehyde and 0.1% triton X-100, respectively, followed by intracellular staining with an IFN-γ antibody. Cells were analyzed by using a Cytek DxP10 cytometer (Cytek Biosciences, Fremont, CA, USA) and acquired data were analyzed using the FlowJo software (BD Biosciences, San Jose, CA, USA) with gating based upon fluorescence minus one (FMO) control. For each sample, 100,000 events were acquired. Relevant cell population was gated on CD3 + prior to analyzing the CD4 + and CD8 + cells ( Figure S1a ). IFN-γ positive cells were counted from individual CD4 + , CD8 + , and CD4 + CD8 + double positive (DP) cells ( Figure S1b ).

Total RNA was isolated from ILN collected from the infected and non-infected using the TRIzol reagent according to the manufacturer's protocol. RNA was re-purified using RNeasy Mini Kit (Qiagen, Hilden, NRW, Germany) according to the manufacturer instructions, followed by a DNase treatment using Turbo DNA-free TM kit (Life Technologies, V.A. Graiciuno 8, Vilnius, Lithuania) to remove DNA contamination. Purified RNA was submitted to Novogene Bioinformatics Technology Co.Ltd for RNA sequencing.

Sequencing libraries were generated using NEBNext ® UltraTM RNA Library Prep Kit for Illumina ® (NEB, Ipswich, MA, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from 1 mg of total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase (RNase H). Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3 ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, MA, USA). Then 3 µL USER enzyme (NEB, Ipswich, MA, USA) was used with size-selected, adaptor-ligated cDNA at 37 • C for 15 min followed by 5 min at 95 • C before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. Finally, PCR products were purified using AMPure XP system and library quality was assessed on the Agilent Bioanalyzer 2100 system. The libraries were sequenced on an Illumina platform Hiseq 2000 and approximately 40 million raw 125 bp/150 bp paired-end) reads were generated.

RNA-Seq reads were first trimmed with Trim Galore version 0.6.4 [34, 35] , which include filtering reads shorter than 30 bp and low-quality ends from reads (Phred score: <20) and the option of removing reads with Ns. The quality of reads before and after trimming was checked with FastQC version 0.11.7 [36] . The filtered reads were aligned to the Sus scrofa genome (Sscrofa11.1; GCF_000003025.6) using TopHat version 2.0.14 [37] [38] [39] with a read mismatch of 0 and a maximum multi hits of 5 and further default parameters. Alignment files generated by TopHat were then used to generate gene counts using HTSeq version 0.9.1 (htseq-count), with a minimum alignment quality of 10 [40] . A data matrix of gene counts for all samples was created using a custom Python script and the data matrix was used to run the differential gene expression analysis in DESeq2 version 1.22.1 [41] . Genes with an adjusted p value smaller than 0.05 and a log 2 fold change larger than 1 were considered as differentially expressed. The differential expression result table generated by DESeq2 was annotated with gene information obtained from Ensembl BioMart for Sus scrofa, supplemented by annotation from the GTF file. To visualize the level of PRRSV RNA genome, reads mapped to the PRRSV CON90 genome were counted base by base, normalized to counts per million and subsequently used to generate a coverage track using deepTools version 3.4.3 [42] .

Gene Ontology (GO) in the database (http://www.geneontology.org/) is an international standardized classification system for gene function, and it supplies a set of controlled vocabulary to comprehensively describe the properties of genes and gene products. GO enrichment analysis of DEGs was implemented by the GOseq R package [43] , in which gene length bias was corrected. GO terms with Viruses 2020, 12, 817 6 of 22 corrected p-value less than 0.05 were considered significantly enriched by differentially expressed genes. Kyoto Encyclopedia of Genes and Genomes (KEGG) is a database resource for understanding high-level functions and utilities of the biological system, such as the cell, the organism, and the ecosystem, from molecular-level information, especially large-scale molecular datasets generated by genome sequencing and other high-throughput experimental technologies (http://www.genome.jp/kegg/). KOBAS version 3.0 [44] was used to test the statistical enrichment of differential expression genes in KEGG pathway.

Virus-specific T-cell data were analyzed by unpaired t-test analysis using GraphPad prism software version 8.3.1 (GraphPad Software, LLC, San Diego, CA, USA). A p < 0.05 was considered significant.

Sequencing data from RNA-Seq were deposited in the NCBI GEO and are available under accession number GEO: GSE153174

All pigs inoculated with CON90 became viremic starting from 2 dpi. The viremia level increased and peaked at 7 dpi ( Figure 1a ). Pigs inoculated with CON90 did not displayed any clinical signs throughout the course of this study. At 46 dpi, ILN was collected and RNA was extracted. RT-PCR analysis revealed that viral RNA genome was detected in all five CON90-infected pigs but none of the sham-inoculated pigs ( Figure 1b) . Subsequently, the RNA samples were subjected to RNA-seq. An average of 32 million paired-end reads were obtained for each RNA sample (Table S1 ). To confirm the presence of viral RNA during persistent infection in pigs, the reads were mapped to the CON90 genome. Viral RNA reads were only detected from ILN of CON90-infected pigs which mapped throughout the viral genome ( Figure 1c ). The results demonstrate that all five pigs in this study were persistently infected with the attenuated PRRSV strain CON90.

To examine the host responses to CON90 persistent infection, RNA reads were mapped to the reference pig genome (Sscrofa11.1; GCF_000003025.6). Principal component analysis (PCA) indicated that control and CON90-infected pigs formed separated clusters indicating that they had distinct transcriptional profiles ( Figure 2a ). The infected pigs showed more variability than the control pigs. In addition, one pig in the CON90-infected group (CON_429) showed a unique transcriptional profile, as it associated closer to the control pigs. This association is likely due to the different genetic makeup of this pig than the other four animals within this group.

Direct comparison of RNA transcripts between control and CON90-infected animals revealed a profound change in the transcriptional profiles in the ILN tissue of CON90-infected pigs. Out of 17553 genes in the annotated porcine genome, there were 6404 DEGs (FDR < 5%, |log 2 | fold change ≥ 1), of which 3960 DEGs were upregulated and 2444 DEGs were downregulated ( Figure 2b and Table S2 ). 

To examine the host responses to CON90 persistent infection, RNA reads were mapped to the reference pig genome (Sscrofa11.1; GCF_000003025.6). Principal component analysis (PCA) indicated that control and CON90-infected pigs formed separated clusters indicating that they had distinct transcriptional profiles (Figure 2a ). The infected pigs showed more variability than the control pigs. In addition, one pig in the CON90-infected group (CON_429) showed a unique transcriptional profile, as it associated closer to the control pigs. This association is likely due to the different genetic makeup of this pig than the other four animals within this group. Direct comparison of RNA transcripts between control and CON90-infected animals revealed a profound change in the transcriptional profiles in the ILN tissue of CON90-infected pigs. Out of 17553 genes in the annotated porcine genome, there were 6404 DEGs (FDR < 5%, |log2| fold change ≥ 1), of which 3960 DEGs were upregulated and 2444 DEGs were downregulated ( Figure 2b and Table S2 ). The control animals are indicated by red dots labeled "DMEM_ animal number," while the CON90-infected animals are indicated by cyan dots labeled "CON_animal number"; (b) the volcano plot of differentially expressed genes between CON90-infected and control pigs. Red dots denote significantly down-regulated genes (p < 0.05, log2 fold change ≤−1), green dots indicated significantly up-regulated genes (p < 0.05, log2 fold change ≥1), and grey dots represent genes that were not differentially expressed.

To understand their biological functions, DEGs were subjected to GO enrichment analysis. GO terms having a corrected value of p < 0.05 were considered statistically significant (Table S3) . Significantly enriched GO terms were grouped under three categories: biological processes, molecular function, and cellular components. Thirty highly enriched GO terms in CON90-infected pigs are shown (Figure 3a) . A large number of DEGs are enriched into the category of biological processes (cytoskeleton organization, multicellular organism development, regulation of intracellular signal transduction, cell communication, cell surface receptor signaling, regulation of programmed cell death etc.), followed by molecular function (cellular macromolecular catabolic processes, kinase activity, phosphotransferase activity, DNA binding, metal ion binding) and cellular components (ribonucleoprotein complexes, bounding membrane of organelles). KEGG pathway analysis revealed multiple enriched pathways including apoptosis, chemokine signaling, cellular, senescence, mitophagy, lysosome, endocytosis, and MAPK ( Figure 3b and Table S4 ). Detailed analysis of selected enriched pathways is presented in the respective sections below. To understand their biological functions, DEGs were subjected to GO enrichment analysis. GO terms having a corrected value of p < 0.05 were considered statistically significant (Table S3) . Significantly enriched GO terms were grouped under three categories: biological processes, molecular function, and cellular components. Thirty highly enriched GO terms in CON90-infected pigs are shown (Figure 3a) . A large number of DEGs are enriched into the category of biological processes (cytoskeleton organization, multicellular organism development, regulation of intracellular signal transduction, cell communication, cell surface receptor signaling, regulation of programmed cell death etc.), followed by molecular function (cellular macromolecular catabolic processes, kinase activity, phosphotransferase activity, DNA binding, metal ion binding) and cellular components (ribonucleoprotein complexes, bounding membrane of organelles). KEGG pathway analysis revealed multiple enriched pathways including apoptosis, chemokine signaling, cellular, senescence, mitophagy, lysosome, endocytosis, and MAPK ( Figure 3b and Table S4 ). Detailed analysis of selected enriched pathways is presented in the respective sections below. 

The innate immune system is the first line of defense against invading pathogens, with the major influence on the development of strong adaptive immune responses. CON90 is capable of inducing type I IFNs both in vitro and in vivo [27] . In the current study, we did not detect differential expression of type I IFN or interferon stimulated genes (ISG) RNA transcripts, suggesting that the type-I IFN was not induced in ILN of CON-infected pigs at 46 dpi. We observed an upregulation of RNA-sensing molecules including TLR3 (dsRNA), TLR7, and TLR8 (ssRNA) (Figure 4a ). On the other hand, no differential expression of signaling molecules downstream of TLR pathways including interferon regulatory factors (IRF) IRF3 and IRF7, was observed (Table S2 ). In addition, we observed an upregulation of TLR10 mRNA in CON90-infected pigs. While the ligands and signaling pathway involving TLR10 remain poorly understood, it has been demonstrated that TLR10 can act as an anti-inflammatory receptor which can also suppress TLR3-induced IFN-β production [45] . Interestingly, we observed an increased expression of CD200 (OX2) and its receptor CD200R (Figure 4a ). CD200R is an inhibitory immune receptor that is expressed on myeloid cells and B-and T-lymphocytes [46] , while CD200 is widely expressed on multiple cell types including endothelial cells, neurons, and lymphocytes [47] . CD200/CD200R interaction suppresses cytokine production, and inflammatory responses [48] . The innate immune system is the first line of defense against invading pathogens, with the major influence on the development of strong adaptive immune responses. CON90 is capable of inducing type I IFNs both in vitro and in vivo [27] . In the current study, we did not detect differential expression of type I IFN or interferon stimulated genes (ISG) RNA transcripts, suggesting that the type-I IFN was not induced in ILN of CON-infected pigs at 46 dpi. We observed an upregulation of RNA-sensing molecules including TLR3 (dsRNA), TLR7, and TLR8 (ssRNA) (Figure 4a ). On the other hand, no differential expression of signaling molecules downstream of TLR pathways including interferon regulatory factors (IRF) IRF3 and IRF7, was observed (Table S2 ). In addition, we observed an upregulation of TLR10 mRNA in CON90-infected pigs. While the ligands and signaling pathway involving TLR10 remain poorly understood, it has been demonstrated that TLR10 can act as an antiinflammatory receptor which can also suppress TLR3-induced IFN-β production [45] . Interestingly, we observed an increased expression of CD200 (OX2) and its receptor CD200R (Figure 4a ). CD200R is an inhibitory immune receptor that is expressed on myeloid cells and B-and T-lymphocytes [46] , while CD200 is widely expressed on multiple cell types including endothelial cells, neurons, and lymphocytes [47] . CD200/CD200R interaction suppresses cytokine production, and inflammatory responses [48] .

The complement system is another crucial component of the innate immunity. It acts as an important link between innate and adaptive immune system. Thus, viruses have developed a mechanism to modulate complement responses by down-regulating activation proteins and upregulation of regulatory proteins [49] . In the current study, increased expression of regulatory proteins (CFH, CFI, and CD55) was observed from CON90-infected animals (Figure 4b) . Conversely, C1q, a classical complement component involved in increasing virus neutralizing ability of antibodies [50] was down-regulated ( Figure 4b) . Collectively, the data suggest that the CON90 virus suppresses both complement activation and inflammatory responses at the site of persistent infection. 

Apoptosis is a powerful innate immunity mechanism to curtail viral spread through eliminating virally infected cells. Apoptosis can be triggered by both extrinsic and intrinsic stimuli [51] . It is well documented that PRRSV induces apoptosis in tissues of infected pigs, especially during an acute phase of the infection [15, 16, 52] . In this study, expression of TNFRSF1A (TNFα receptor 1) gene that The complement system is another crucial component of the innate immunity. It acts as an important link between innate and adaptive immune system. Thus, viruses have developed a mechanism to modulate complement responses by down-regulating activation proteins and up-regulation of regulatory proteins [49] . In the current study, increased expression of regulatory proteins (CFH, CFI, and CD55) was observed from CON90-infected animals (Figure 4b) . Conversely, C1q, a classical complement component involved in increasing virus neutralizing ability of antibodies [50] was down-regulated ( Figure 4b) . Collectively, the data suggest that the CON90 virus suppresses both complement activation and inflammatory responses at the site of persistent infection.

Apoptosis is a powerful innate immunity mechanism to curtail viral spread through eliminating virally infected cells. Apoptosis can be triggered by both extrinsic and intrinsic stimuli [51] . It is well documented that PRRSV induces apoptosis in tissues of infected pigs, especially during an acute phase of the infection [15, 16, 52] . In this study, expression of TNFRSF1A (TNFα receptor 1) gene that contains a death domain [53] , was downregulated in the ILN of infected pigs ( Figure 5) . Similarly, expression of pro-apoptotic genes including AIFM2, CHAC1, and OSR1, was downregulated ( Figure 5 ). On the other hand, expression of BIRC3 and Bfl-1/A1 (BCL2A1), two apoptosis inhibitors that interfere with caspase activation [54] , was upregulated. Furthermore, expression of several negative regulators of apoptotic genes including BCL-2-associated killer 1 (BAK1), damage induced apoptosis suppressor (DDIAS), X-linked inhibitor of apoptosis protein (XIAP), MCL1, API, BNIP2, and FAIM was upregulated. Collectively, these results suggest that the pro-apoptotic signaling pathway was suppressed in ILN tissue of pigs persistently infected with CON90. contains a death domain [53] , was downregulated in the ILN of infected pigs ( Figure 5) . Similarly, expression of pro-apoptotic genes including AIFM2, CHAC1, and OSR1, was downregulated ( Figure  5 ). On the other hand, expression of BIRC3 and Bfl-1/A1 (BCL2A1), two apoptosis inhibitors that interfere with caspase activation [54] , was upregulated. Furthermore, expression of several negative regulators of apoptotic genes including BCL-2-associated killer 1 (BAK1), damage induced apoptosis suppressor (DDIAS), X-linked inhibitor of apoptosis protein (XIAP), MCL1, API, BNIP2, and FAIM was upregulated. Collectively, these results suggest that the pro-apoptotic signaling pathway was suppressed in ILN tissue of pigs persistently infected with CON90. 

Lymphocyte activation occurs in the secondary lymphoid organs. Activated T-cells then migrate to the sites of infection where they exert their effector functions to eliminate infected cells [55] . Celladhesion molecules, chemokines, and receptors play a central role in regulating T-cell migration [56] . PRRSV mainly persists in lymphoid tissue of infected pigs [17, 18] . In this study, expression of several important chemokine ligands (CCL19, CCL21, CCL24, CCL22, CX3CL1, and CCL14) and chemokine receptors (CCR6 and CCR10) which play an essential role in migration and localization of lymphocytes and antigen-presenting cells (APCs) to the lymphoid tissues was down regulated in ILN of CON90-infected pigs (Figure 6a ). On the other hand, expression of CD274 (PD-1), a marker of Tcell exhaustion, was upregulated ( Figure 6b) . Likewise, expression of inhibitory receptors HAVCR2 (also known as TIM3) and TGIT, which transmit the inhibitory signals for T-cell differentiation and effector activities [57] , was upregulated. We also found increased expression of other co-inhibitory molecules (BTLA, FASLG, FAS, and IDO1) that are associated with the regulation of T-cell exhaustion during chronic viral infection (Figure 6b) . Together, the results suggest that T-cell migration to ILN, one of the sites of PRRSV persistence, might be affected because of the down regulation of important chemokines and receptors, and that the T-cells in ILN might be exhausted. Interestingly, markers of 

Lymphocyte activation occurs in the secondary lymphoid organs. Activated T-cells then migrate to the sites of infection where they exert their effector functions to eliminate infected cells [55] . Cell-adhesion molecules, chemokines, and receptors play a central role in regulating T-cell migration [56] . PRRSV mainly persists in lymphoid tissue of infected pigs [17, 18] . In this study, expression of several important chemokine ligands (CCL19, CCL21, CCL24, CCL22, CX3CL1, and CCL14) and chemokine receptors (CCR6 and CCR10) which play an essential role in migration and localization of lymphocytes and antigen-presenting cells (APCs) to the lymphoid tissues was down regulated in ILN of CON90-infected pigs (Figure 6a ). On the other hand, expression of CD274 (PD-1), a marker of T-cell exhaustion, was upregulated ( Figure 6b) . Likewise, expression of inhibitory receptors HAVCR2 (also known as TIM3) and TGIT, which transmit the inhibitory signals for T-cell differentiation and effector activities [57] , was upregulated. We also found increased expression of other co-inhibitory molecules (BTLA, FASLG, FAS, and IDO1) that are associated with the regulation of T-cell exhaustion during chronic viral infection (Figure 6b) . Together, the results suggest that T-cell migration to ILN, one of the sites of PRRSV persistence, might be affected because of the down regulation of important chemokines and receptors, and that the T-cells in ILN might be exhausted. Interestingly, markers of regulatory T-cells (T reg ) were downregulated in ILN of infected pigs (Figure 6c) , suggesting that T reg might not be present in ILN at 46 dpi. regulatory T-cells (Treg) were downregulated in ILN of infected pigs (Figure 6c ), suggesting that Treg might not be present in ILN at 46 dpi. 

Since the expression of chemokines and receptors important for T-cell migration was downregulated in ILN, we sought to compare the frequencies of virus-specific T-cells in PBMCs and in ILN using the IFN-γ SC Elispot assay. The number of spots was similar when PBMCs and ILN cells were stimulated with non-specific T-cell activator PHA/Ionomycin. However, the number of spots was significantly lower in ILN than in PBMCs when the cells were stimulated with whole PRRSV antigen (Figure 7a ). The results clearly indicate that the frequencies of PRRSV-specific IFN-γ SCs were significantly lower in ILN than in PBMCs. To further elucidate the phenotypes of PRRSVspecific T-cells, we used a multi-color flow cytometric assay to identify subsets of T-cells that secrete IFN-γ in response to PRRSV activation ex vivo. CD4 + CD8 + DP cells were the major T-cell population secreting IFN-γ, both in PBMCs and ILN. Swine have a significant population of extrathymic CD4 + CD8 + DP T cells that represent memory T-cells [58] . Similar to the Elispot assay, the flow cytometric assay also indicated that the percentage of T-cell secreting IFN-γ was comparatively lower in ILN than in PBMCs (Figure 7b ). Cells collected from DMEM-inoculated pigs did not show any significant T-cell reactivities after ex vivo stimulation with CON90, both in Elispot and in flow cytometry assays (data not shown). 

Since the expression of chemokines and receptors important for T-cell migration was downregulated in ILN, we sought to compare the frequencies of virus-specific T-cells in PBMCs and in ILN using the IFN-γ SC Elispot assay. The number of spots was similar when PBMCs and ILN cells were stimulated with non-specific T-cell activator PHA/Ionomycin. However, the number of spots was significantly lower in ILN than in PBMCs when the cells were stimulated with whole PRRSV antigen (Figure 7a ). The results clearly indicate that the frequencies of PRRSV-specific IFN-γ SCs were significantly lower in ILN than in PBMCs. To further elucidate the phenotypes of PRRSV-specific T-cells, we used a multi-color flow cytometric assay to identify subsets of T-cells that secrete IFN-γ in response to PRRSV activation ex vivo. CD4 + CD8 + DP cells were the major T-cell population secreting IFN-γ, both in PBMCs and ILN. Swine have a significant population of extrathymic CD4 + CD8 + DP T cells that represent memory T-cells [58] . Similar to the Elispot assay, the flow cytometric assay also indicated that the percentage of T-cell secreting IFN-γ was comparatively lower in ILN than in PBMCs (Figure 7b ). Cells collected from DMEM-inoculated pigs did not show any significant T-cell reactivities after ex vivo stimulation with CON90, both in Elispot and in flow cytometry assays (data not shown). Viruses 2020, 12, x FOR PEER REVIEW 13 of 22 

A number of genes associated with Th2, or humoral response, were upregulated in the ILN of CON90-infected pigs (Figure 8a ). IL-21, RGS13, and NUGGC are involved in the development of germinal center (GC) and activation of B-cell follicles [59] . TNFSF13B and TNFSF8 are potent activators of B-cell lineage and Ig class switching, respectively [60] . B-cell surface antigens MS4A1 (CD20), activation-induced cytidine deaminase (AICDA) [61] , and rafting family member 2 (RFTN2) are associated with B-cell receptor signaling. The upregulated expression of these genes in ILN of infected animals suggested that the humoral immune response to CON90-infection was not affected. To corroborate transcriptome data, we measured both virus-specific antibody levels at various time points post-infection. High levels of non-neutralizing antibodies specific to viral N protein (measured by a commercial ELISA) were detected in all pigs (Figure 8b ). However, only minimal levels of virusneutralizing antibodies (titer 1:2) were detected in the serum of infected pigs at 46 dpi (Figure 8c ). Together, both transcriptomic and serological data indicate that B-cell development and antibody production are not affected by CON90-infection. However, low levels of virus-neutralizing antibodies are likely due to the virus ability to escape antibody neutralization (see below). 

A number of genes associated with Th2, or humoral response, were upregulated in the ILN of CON90-infected pigs (Figure 8a ). IL-21, RGS13, and NUGGC are involved in the development of germinal center (GC) and activation of B-cell follicles [59] . TNFSF13B and TNFSF8 are potent activators of B-cell lineage and Ig class switching, respectively [60] . B-cell surface antigens MS4A1 (CD20), activation-induced cytidine deaminase (AICDA) [61] , and rafting family member 2 (RFTN2) are associated with B-cell receptor signaling. The upregulated expression of these genes in ILN of infected animals suggested that the humoral immune response to CON90-infection was not affected. To corroborate transcriptome data, we measured both virus-specific antibody levels at various time points post-infection. High levels of non-neutralizing antibodies specific to viral N protein (measured by a commercial ELISA) were detected in all pigs (Figure 8b ). However, only minimal levels of virus-neutralizing antibodies (titer 1:2) were detected in the serum of infected pigs at 46 dpi (Figure 8c ). Together, both transcriptomic and serological data indicate that B-cell development and antibody production are not affected by CON90-infection. However, low levels of virus-neutralizing antibodies are likely due to the virus ability to escape antibody neutralization (see below). Viruses 2020, 12, x FOR PEER REVIEW 14 of 22 

PRRSV persists in lymphoid tissue of infected pigs for several months [17, 18] . The mechanism of PRRSV persistence is not fully understood. We performed genome-wide transcriptome analysis of lymphoid tissue collected from pigs persistently infected with an attenuated PRRSV strain using RNA-seq technology that detects both host and viral RNA. Viral RNA reads were detected in ILN of all five infected pigs. It was reported previously that PRRSV genome mainly exists in dsRNA forms in lymphoid tissues during persistent infection [26] . Since only mRNA (purified by using poly-T oligo-attached magnetic beads) was used for library construction and RNA sequencing, the viral RNA reads detected in this study must be derived from either viral genomic RNA or sub-genomic mRNA, not from dsRNA. The viral RNA reads map throughout the viral genome. However, we are not able to discern whether these reads are derived from genomic or sub-genomic mRNA because we used a short-read RNA sequencing platform. It would be interesting to use long-read RNA sequencing to study the viral transcriptome at different states of infection in pigs [62] .

It was reported previously that no significant DEGs were observed in lymphoid tissue of persistently infected animals [26] . In the current study, we identified a large number of DEGs in persistently infected animals (Figure 2b ). This might be due to the difference in the experimental setup. In the previous study, pigs were infected with a wild-type PRRSV-1 whereas in this current 

PRRSV persists in lymphoid tissue of infected pigs for several months [17, 18] . The mechanism of PRRSV persistence is not fully understood. We performed genome-wide transcriptome analysis of lymphoid tissue collected from pigs persistently infected with an attenuated PRRSV strain using RNA-seq technology that detects both host and viral RNA. Viral RNA reads were detected in ILN of all five infected pigs. It was reported previously that PRRSV genome mainly exists in dsRNA forms in lymphoid tissues during persistent infection [26] . Since only mRNA (purified by using poly-T oligo-attached magnetic beads) was used for library construction and RNA sequencing, the viral RNA reads detected in this study must be derived from either viral genomic RNA or sub-genomic mRNA, not from dsRNA. The viral RNA reads map throughout the viral genome. However, we are not able to discern whether these reads are derived from genomic or sub-genomic mRNA because we used a short-read RNA sequencing platform. It would be interesting to use long-read RNA sequencing to study the viral transcriptome at different states of infection in pigs [62] .

It was reported previously that no significant DEGs were observed in lymphoid tissue of persistently infected animals [26] . In the current study, we identified a large number of DEGs in persistently infected animals (Figure 2b ). This might be due to the difference in the experimental setup. In the previous study, pigs were infected with a wild-type PRRSV-1 whereas in this current study, pigs were infected with an attenuated synthetic PRRSV strain [27, 63] . Besides, the previous study looked at a small set of selected genes while in this study we look at the genome-wide RNA transcriptome.

GO terms and KEGG analysis revealed that genes involved in the innate immune (complement and TLR) pathways, apoptosis, cytokine-chemokine signaling, T-cell exhaustion, and humoral responses are highly differentially expressed in PRRSV-infected animals compared to control. The complement system is a constituent of innate immunity that serves by neutralizing cell-free viruses, lysing virus-infected cells, and boosting virus specific responses [64] . It also links the innate and adaptive immune responses, enhances humoral immunity, regulates antibody effector mechanisms, and modulates T-cell function [65] . Many viruses have developed a strategy to evade the complement pathway by recruiting or enhancing the production of host transcription regulatory components. During acute infection (7 dpi), PRRSV significantly represses the expression of complement regulatory components (CD55 and C4BPB) in lung tissue [16] . Early complement activation during acute infection facilitates the release of newly formed virions from infected cells, yet, chronic viral infection is reported to suppress the activity of complement activation proteins and increases the activity of regulatory proteins (Reviewed [49] ).

In the current study, we found overexpression of CFH and CFI along with another regulatory factor CD55 (decay-accelerating factor). Overexpression of CFH, a major soluble regulator of the alternative pathway, results in inhibition of C3 and C5 convertase enzymes. CFI suppresses the complement active proteins C3b (opsonin) and C4b via mediating cleavage to their inactive form [66, 67] . Regulatory factor CD55 is an inhibitor of C3 convertase which prevents C3b deposition on the cell surface [68] . Apart from that, we also see the downregulation of C1q, which increases neutralizing and hemagglutination inhibition activity of anti-influenza antibodies [50] . Available data from previous studies and the current study suggest that activation of complement pathway during acute infection helps to disseminate virus, while suppression of complement components like C1q, C1r, and C5, and upregulation in regulatory components during persistent infection facilitates the virus to escape from the complement system to maintain persistence in lymphoid tissues.

Most naturally occurring PRRSV strains suppress type I IFN production by inhibiting the activation and nuclear translocation of IRF3/IRF7 and NF-kB [69, 70] . Deficiency of IRF3/7 results in severe mortality to infection with West Nile virus (WNV), Chikungunya virus (CHIKV), and Ross River virus infection, and promotes viral persistence in the infected hosts [71] [72] [73] . It has been reported that PRRSV nsp1β inhibits IRF3 phosphorylation and nuclear translocation; thus, inhibiting IFN production [69, 74, 75] . Interestingly, the synthetic PRRSV-CON and its attenuated form CON90 were able to induce type I IFNs [27, 76] . Contrary to its nature to induce type I IFN, in the current study no upregulation of canonical type I IFN signaling pathway-associated genes was observed. However, expression of RNA sensing molecules including TLR3, TLR7, and TLR8 was upregulated in CON90-infected animals. On the basis of available data, we hypothesize that reduced levels of viral replication and sequestration of viral RNA in infected cells during persistent infection might limit further activation of type I IFN signaling pathway by preventing interaction with cytoplasmic pattern recognition receptors (PRRs).

Apoptosis or programmed cell death is a potential host immune mechanism against virally infected cells to curtail the spread of newly formed viral progenies. Apoptosis is induced by two distinct, yet inter-connected signaling pathways, the extrinsic and intrinsic pathways [51] . In order to successfully establish persistent infection, viruses have developed mechanisms to inhibit apoptosis. For instance, adenovirus (E1B-19K), human cytomegalovirus (UL37), poxviruses (F1L), and myxoma virus (M11L) have an inhibitory effect on proapoptotic proteins Bak/Bax [77] . It appears that PRRSV can also modulate apoptosis. Studies of pulmonary alveolar macrophages (PAM) infected with PRRSV in vitro reveal that the virus stimulates anti-apoptotic pathways early in infection while it induces apoptosis late in infection [78] . On the other hand, the virus induces apoptosis in the tissues of infected animals during acute infection, but the frequencies of apoptotic cells reduced to normal levels observed in control, non-infected pigs from day 20 post-infection [52] . In this study, expression of multiple anti-apoptotic genes including XIAP, Bfl-1/A1 (BCL2A1), and BIRC3 was upregulated while expression of pro-apoptotic genes (AIFM2, CHAC1, and OSR1) was downregulated. XIAP is an endogenous caspase 9 and 3 inhibitor [79] . Bfl-1/A1 (BCL2A1) is a transcriptional target of nuclear factor-kB (NF-kB) that suppresses caspase activation and apoptosis in response to death-inducing stimuli like TNFα [54] . BIRC3 is an interaction partner to TNFRSF1B and inhibits apoptosis by interfering with the caspase activation [53] . The pro-apoptotic BCL-2 family member like BCL-2-associated killer 1 (BAK1), which allows the release of cytochrome C via formation of homo-oligomers and stable insertion into outer mitochondrial membrane was significantly suppressed. Thus, the data suggest that apoptosis was suppressed in the ILN of CON90-infected animals during persistent infection.

Chemokine molecules CCL19 and CCL21 play a crucial role in migration, activation, expansion, and survival of antiviral T-cells. Suppression of the CCR7 and CCL19/CCL21 axis results in dysfunction of T-cells during viral infection [80, 81] . CCL19-CCR7 axis plays role in protection against multiple viruses, such as HIV [82] , herpes simplex virus (HSV-1) [83] , HSV-2, hepatitis C virus [84] , and pseudorabies [85] . Acute PRRSV infection increases the expression of chemokines, leading to the infiltration of immune cells toward the sites of infection [16, 86] . In this study, expression of both CCL19 and CCL21 was significantly downregulated in CON90-infected lymph node (Figure 6a ). Additionally, expression of different chemoattractant molecules (fractalkine/CX3CL1, CCL14, CCL16, CCL17, CCL22, and CCL8) and receptors (CCR6 and CCR10) was also downregulated. It is possible that downregulation of these chemokines and receptors might impair T-cells trafficking to inguinal lymph node, the site of PRRSV persistence. Our transcriptional data are supported by the functional data which demonstrate that the frequencies of IFN-γ SC in ILN were significantly lower than in PBMCs (Figure 7) .

Chronic or persistent viral stimulation results in hierarchical loss of effector functions of T-lymphocytes, including proliferation, cytokine production (e.g., IFN-γ and IL-2), and cytolytic responses [57, 87, 88] . Although the molecular signatures involved in T-cell exhaustion are not completely understood, the overexpression of cell surface inhibitory receptors (e.g., PD-1, CTLA-4, and others) primarily mediates CD8+ T-cell dysfunction [87] [88] [89] . T-cell exhaustion seems to be a common phenomenon caused by arteriviruses, as T-cell exhaustion was also reported in horses persistently infected with EAV [24] . Recently, it was shown that PRRSV infection alone or co-infection with porcine circovirus type 2 (PCV2) can significantly upregulate surface expression of PD-L1 (CD274) on porcine monocytes-derived dendritic cells (MoDCs) [90] . Increased expression of PD-L1 on the surface of antigen presenting cells (APCs) possibly contributes to the ineffective T-cell responses during the PRRSV infection. In the present study, expression of several markers for T-cells exhaustion such as PD-1 (CD279) and its ligand, CD274 (PD-L1) was upregulated in CON90-infected animals, suggesting that the T-cells in lymph nodes of persistently infected animals might be exhausted (Figure 6b) . However, additional studies (e.g., T-cell cytotoxicity and cytokine production) are required to fully assess the functionality of T-cells during PRRSV persistence.

Germinal centers (GCs) are the specialized locations in the secondary lymphoid tissues where B-cell maturation, differentiation, somatic hypermutation, and class switching of isotypes take place [91] . In the current study, the GC development-associated genes, which includes IL-21 produced by follicular helper T-cells [61] , the regulator of G protein signaling (RGS13) [92] , and GC-associated nuclear GTPase (NUGGC) [59] , were upregulated in the lymph node of CON90-infected animals. Thus, PRRSV infection does not seem to impair GC formation and B-cell activation. This is supported by the fact that all five CON90-infected pigs developed a robust antibody response, as measured by the IDEXX ELISA ( Figure 8b) . However, only minimal titers of virus-neutralizing antibodies were detected at 46 dpi. Perhaps, PRRSV does not suppress the humoral immune response. Instead, the virus evades from antibody-mediated neutralization through different mechanisms such as glycan shielding and decoy epitopes [93] [94] [95] .

In summary, we identified a robust host transcriptional response in the inguinal lymphoid tissue of pigs persistently infected with the attenuated PRRSV strain CON90. Genes involved in the innate immune responses are downregulated. Similarly, chemokines and receptors associated with T-cell homing to the lymphoid tissue are downregulated. This might lead to the lower frequencies of virus-specific T-cells in lymphoid tissue than in peripheral blood. Additionally, the genes involved in the anti-apoptotic pathways are upregulated. Collectively, the data suggest that PRRSV can create a pro-survival microenvironment at the lymphoid tissue which allows the virus to persist for an extended period.

Supplementary Materials: The following are available online at http://www.mdpi.com/1999-4915/12/8/817/s1. Figure S1 . Representative flow cytometry gating strategy used to identify IFN-γ secreting CD4α + , CD8α + , and CD4α + CD8α + (DP) T cells from PBMCs and ILN. Table S1 . Summary of sequencing quality control and mapping data of samples. 

Arterivirus molecular biology and pathogenesis

ICTV-International Committee on Taxonomy of Viruses. Virus Taxonomy: 2019 Release. EC 51

Porcine reproductive and respiratory syndrome

Assessment of the economic impact of porcine reproductive and respiratory syndrome virus on United States pork producers

Effects of origin and state of differentiation and activation of monocytes/macrophages on their susceptibility to porcine reproductive and respiratory syndrome virus (PRRSV)

CD163 expression confers susceptibility to porcine reproductive and respiratory syndrome viruses

In vivo detection of porcine reproductive and respiratory syndrome virus RNA by in situ hybridization at different times postinfection

Apoptosis induced in vivo during acute infection by porcine reproductive and respiratory syndrome virus

Type 2 porcine reproductive and respiratory syndrome virus infection increases apoptosis at the maternal-fetal interface in late gestation pregnant gilts

Open reading frame 5 of porcine reproductive and respiratory syndrome virus as a cause of virus-induced apoptosis

Characterization of the porcine reproductive and respiratory syndrome virus glycoprotein 5 (GP5) in stably expressing cells

Induction of Apoptosis by the Nonstructural Protein 4 and 10 of Porcine Reproductive and Respiratory Syndrome Virus

Pig immune response to general stimulus and to porcine reproductive and respiratory syndrome virus infection: A meta-analysis approach

Expression of Toll-like receptor mRNA and cytokines in pigs infected with porcine reproductive and respiratory syndrome virus

Analysis of the swine tracheobronchial lymph node transcriptomic response to infection with a Chinese highly pathogenic strain of porcine reproductive and respiratory syndrome virus

Understanding PRRSV infection in porcine lung based on genome-wide transcriptome response identified by deep sequencing

Porcine reproductive and respiratory syndrome virus: Description of persistence in individual pigs upon experimental infection

Duration of infection and proportion of pigs persistently infected with porcine reproductive and respiratory syndrome virus

Humoral immune responses and viral shedding following vaccination with modified live porcine reproductive and respiratory syndrome virus vaccines

Simian hemorrhagic fever virus: Recent advances

Lactate dehydrogenase-elevating virus

Equine viral arteritis

Genome-wide association study among four horse breeds identifies a common haplotype associated with in vitro CD3+ T cell susceptibility/resistance to equine arteritis virus infection

Equine arteritis virus long-term persistence is orchestrated by CD8+ T lymphocyte transcription factors, inhibitory receptors, and the CXCL16/CXCR6 axis

Identification of factors associated with virus level in tonsils of pigs experimentally infected with porcine reproductive and respiratory syndrome virus

Double-stranded viral RNA persists in vitro and in vivo during prolonged infection of porcine reproductive and respiratory syndrome virus

Development of a broadly protective modified-live virus vaccine candidate against porcine reproductive and respiratory syndrome virus

Enhanced replication of porcine reproductive and respiratory syndrome (PRRS) virus in a homogeneous subpopulation of MA-104 cell line

Gradual development of the interferon-gamma response of swine to porcine reproductive and respiratory syndrome virus infection or vaccination

Location of T-cell epitopes in nonstructural proteins 9 and 10 of type-II porcine reproductive and respiratory syndrome virus

A 10-kDa structural protein of porcine reproductive and respiratory syndrome virus encoded by ORF2b

Use of interleukin 12 to enhance the cellular immune response of swine to an inactivated herpesvirus vaccine

Cross reactivity of immune responses to porcine reproductive and respiratory syndrome virus infection

A Wrapper Script to Automate Quality and Adapter Trimming

Cutadapt removes adapter sequences from high-throughput sequencing reads

A Quality Control Tool for High throughput Sequence Data

Ultrafast and memory-efficient alignment of short DNA sequences to the human genome

Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks

Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions

HTSeq-a Python framework to work with high-throughput sequencing data

Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

Manke, T. deepTools2: A next generation web server for deep-sequencing data analysis

Gene ontology analysis for RNA-seq: Accounting for selection bias

KOBAS 2.0: A web server for annotation and identification of enriched pathways and diseases

TLR10 Is a Negative Regulator of Both MyD88-Dependent and -Independent TLR Signaling

The inhibitory CD200R is differentially expressed on human and mouse T and B lymphocytes

CD200 and membrane protein interactions in the control of myeloid cells

CD200/CD200R paired potent inhibitory molecules regulating immune and inflammatory responses; Part II: CD200/CD200R potential clinical applications

Complement Evasion Strategies of Viruses: An Overview. Front

Influenza A virus M1 blocks the classical complement pathway through interacting with C1qA

Apoptosis: A review of programmed cell death

Apoptosis in the lungs of pigs infected with porcine reproductive and respiratory syndrome virus and associations with the production of apoptogenic cytokines

NF-kappaB antiapoptosis: Induction of TRAF1 and TRAF2 and c-IAP1 and c-IAP2 to suppress caspase-8 activation

Bfl-1/A1 functions, similar to Mcl-1, as a selective tBid and Bak antagonist

Activation, differentiation, and migration of naive virus-specific CD8+ T cells during pulmonary influenza virus infection

Regulation of T cell migration during viral infection: Role of adhesion molecules and chemokines

Molecular and cellular insights into T cell exhaustion

Functional and phenotypic analysis of porcine peripheral blood CD4/CD8 double-positive T cells

Speckled-like pattern in the germinal center (SLIP-GC), a nuclear GTPase expressed in activation-induced deaminase-expressing lymphomas and germinal center B cells

Member of the tumor necrosis factor family and B lymphocyte stimulator

IL-21 and IL-4 Collaborate to Shape T-Dependent Antibody Responses

The Architecture of SARS-CoV-2 Transcriptome

A Synthetic Porcine Reproductive and Respiratory Syndrome Virus Strain Confers Unprecedented Levels of Heterologous Protection

Viral mimicry of the complement system

The complement system in regulation of adaptive immunity

Control of the amplification convertase of complement by the plasma protein beta1H

Human complement C3b inactivator: Isolation, characterization, and demonstration of an absolute requirement for the serum protein beta1H for cleavage of C3b and C4b in solution

Decay accelerating factor (CD55): A target for cancer vaccines?

Porcine reproductive and respiratory syndrome virus nonstructural protein 1beta modulates host innate immune response by antagonizing IRF3 activation

The cysteine protease domain of porcine reproductive and respiratory syndrome virus nonstructural protein 2 possesses deubiquitinating and interferon antagonism functions

Regulation of interferon regulatory factor-3 by the hepatitis C virus serine protease

Infectious Chikungunya Virus in the Saliva of Mice, Monkeys and Humans

Induction of IFN-beta and the innate antiviral response in myeloid cells occurs through an IPS-1-dependent signal that does not require IRF-3 and IRF-7

Interplay between interferon-mediated innate immunity and porcine reproductive and respiratory syndrome virus

The cysteine protease domain of porcine reproductive and respiratory syndrome virus non-structural protein 2 antagonizes interferon regulatory factor 3 activation

Identification of viral genes associated with the interferon-inducing phenotype of a synthetic porcine reproductive and respiratory syndrome virus strain

BAX insertion, oligomerization, and outer membrane permeabilization in brain mitochondria: Role of permeability transition and SH-redox regulation

Porcine reproductive and respiratory syndrome virus modulates apoptosis during replication in alveolar macrophages

Structural biology. Controlling the caspases

CCR7 and its ligands: Balancing immunity and tolerance

CCR7 coordinates the primary immune response by establishing functional microenvironments in secondary lymphoid organs

CCL28 induces mucosal homing of HIV-1-specific IgA-secreting plasma cells in mice immunized with HIV-1 virus-like particles

Influence of CCR7 ligand DNA preexposure on the magnitude and duration of immunity

Enhancement of Immune Responses by Co-delivery of CCL19/MIP-3beta Chemokine Plasmid with HCV Core DNA/Protein Immunization

Genetic co-transfer of CCR7 ligands enhances immunity and prolongs survival against virulent challenge of pseudorabies virus

Molecular characterization of transcriptome-wide interactions between highly pathogenic porcine reproductive and respiratory syndrome virus and porcine alveolar macrophages in vivo

Molecular signature of CD8+ T cell exhaustion during chronic viral infection

T cell exhaustion

T cells and viral persistence: Lessons from diverse infections

PD-L1 expression is increased in monocyte derived dendritic cells in response to porcine circovirus type 2 and porcine reproductive and respiratory syndrome virus infections

Role of NF-kappa B in cell survival and transcription of latent membrane protein 1-expressing or Epstein-Barr virus latency III-infected cells

Influence of N-linked glycosylation of porcine reproductive and respiratory syndrome virus GP5 on virus infectivity, antigenicity, and ability to induce neutralizing antibodies

Immune evasion of porcine reproductive and respiratory syndrome virus through glycan shielding involves both glycoprotein 5 as well as glycoprotein 3

Identification of neutralizing and nonneutralizing epitopes in the porcine reproductive and respiratory syndrome virus GP5 ectodomain

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license

We would like to thank Dirk Anderson at Flow Cytometry Core Facility, Nebraska Center for Biotechnology for assistance on flow cytometric analysis and the staff members of UNL Life Sciences Annex and Veterinary Diagnostic Center for the care of animals.

The use of product and company names is necessary to accurately report the methods and results; however, the United States Department of Agriculture (USDA) neither guarantees nor warrants the standard of the products. The use of names by the USDA implies no approval of the product to the exclusion of others that may also be suitable. The USDA is an equal opportunity provider and employer.