Identifying putative candidate genes and pathways involved in immune responses to porcine reproductive and respiratory syndrome virus (PRRSV) infection

Research progress on the N protein of porcine reproductive and respiratory syndrome virus

Porcine reproductive and respiratory syndrome (PRRS) is a highly contagious disease caused by the porcine reproductive and respiratory syndrome virus (PRRSV). PRRSV exhibits genetic diversity and complexity in terms of immune responses, posing challenges for eradication. The nucleocapsid (N) protein of PRRSV, an alkaline phosphoprotein, is important for various biological functions. This review summarizes the structural characteristics, genetic evolution, impact on PRRSV replication and virulence, interactions between viral and host proteins, modulation of host immunity, detection techniques targeting the N protein, and progress in vaccine development. The discussion provides a theoretical foundation for understanding the pathogenic mechanisms underlying PRRSV virulence, developing diagnostic techniques, and designing effective vaccines.

A New Long Noncoding RNA, MAHAT, Inhibits Replication of Porcine Reproductive and Respiratory Syndrome Virus by Recruiting DDX6 To Bind to ZNF34 and Promote an Innate Immune Response

Long noncoding RNAs (lncRNAs) have increasingly been recognized as being integral to cellular processes, including the antiviral immune response. Porcine reproductive and respiratory syndrome virus (PRRSV) is costly to the global swine industry. To identify PRRSV-related lncRNAs, we performed RNA deep sequencing and compared the profiles of lncRNAs in PRRSV-infected and uninfected Marc-145 cells. We identified a novel lncRNA called MAHAT (maintaining cell morphology-associated and highly conserved antiviral transcript; LTCON_00080558) that inhibits PRRSV replication. MAHAT binds and negatively regulates ZNF34 expression by recruiting and binding DDX6, an RNA helicase forming a complex with ZNF34. Inhibition of ZNF34 expression results in increased type I interferon expression and decreased PRRSV replication. This finding reveals a novel mechanism by which PRRSV evades the host antiviral innate immune response by downregulating the MAHAT-DDX6-ZNF34 pathway. MAHAT could be a host factor target for antiviral therapies against PRRSV infection. IMPORTANCE Long noncoding RNAs (lncRNAs) play important roles in viral infection by regulating the transcription and expression of host genes, and interferon signaling pathways. Porcine reproductive and respiratory syndrome virus (PRRSV) causes huge economic losses in the swine industry worldwide, but the mechanisms of its pathogenesis and immunology are not fully understood. Here, a new lncRNA, designated MAHAT, was identified as a regulator of host innate immune responses. MAHAT negatively regulates the expression of its target gene, ZNF34, by recruiting and binding DDX6, an RNA helicase, forming a complex with ZNF34. Inhibition of ZNF34 expression increases type I interferon expression and decreases PRRSV replication. This finding suggests that MAHAT has potential as a new target for developing antiviral drugs against PRRSV infection.

Evasion of Antiviral Innate Immunity by Porcine Reproductive and Respiratory Syndrome Virus

Innate immunity is the front line for antiviral immune responses and bridges adaptive immunity against viral infections. However, various viruses have evolved many strategies to evade host innate immunity. A typical virus is the porcine reproductive and respiratory syndrome virus (PRRSV), one of the most globally devastating viruses threatening the swine industry worldwide. PRRSV engages several strategies to evade the porcine innate immune responses. This review focus on the underlying mechanisms employed by PRRSV to evade pattern recognition receptors signaling pathways, type I interferon (IFN-α/β) receptor (IFNAR)-JAK-STAT signaling pathway, and interferon-stimulated genes. Deciphering the antiviral immune evasion mechanisms by PRRSV will enhance our understanding of PRRSV’s pathogenesis and help us to develop more effective methods to control and eliminate PRRSV.

Investigating the genetic architecture of disease resilience in pigs by genome-wide association studies of complete blood count traits collected from a natural disease challenge model

Background

Genetic improvement for disease resilience is anticipated to be a practical method to improve efficiency and profitability of the pig industry, as resilient pigs maintain a relatively undepressed level of performance in the face of infection. However, multiple biological functions are known to be involved in disease resilience and this complexity means that the genetic architecture of disease resilience remains largely unknown. Here, we conducted genome-wide association studies (GWAS) of 465,910 autosomal SNPs for complete blood count (CBC) traits that are important in an animal’s disease response. The aim was to identify the genetic control of disease resilience.

Results

Univariate and multivariate single-step GWAS were performed on 15 CBC traits measured from the blood samples of 2743 crossbred (Landrace × Yorkshire) barrows drawn at 2-weeks before, and at 2 and 6-weeks after exposure to a polymicrobial infectious challenge. Overall, at a genome-wise false discovery rate of 0.05, five genomic regions located on Sus scrofa chromosome (SSC) 2, SSC4, SSC9, SSC10, and SSC12, were significantly associated with white blood cell traits in response to the polymicrobial challenge, and nine genomic regions on multiple chromosomes (SSC1, SSC4, SSC5, SSC6, SSC8, SSC9, SSC11, SSC12, SSC17) were significantly associated with red blood cell and platelet traits collected before and after exposure to the challenge. By functional enrichment analyses using Ingenuity Pathway Analysis (IPA) and literature review of previous CBC studies, candidate genes located nearby significant single-nucleotide polymorphisms were found to be involved in immune response, hematopoiesis, red blood cell morphology, and platelet aggregation.

Conclusions

This study helps to improve our understanding of the genetic basis of CBC traits collected before and after exposure to a polymicrobial infectious challenge and provides a step forward to improve disease resilience.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07835-4.

Resilience Effects of SGK1 and TAP1 DNA Markers during PRRSV Outbreaks in Reproductive Sows

The porcine reproductive and respiratory syndrome virus (PRRSV) is a major infectious stressor that causes serious health problems and productivity drops. Based on previous genome-wide analyses, we selected SGK1 and TAP1 as candidate genes for resilience, and genotyped three mutations, including a 3'UTR variant SGK1_rs338508371 and two synonymous variants TAP1_rs1109026889 and TAP1_rs80928141 in 305 Landrace × Large White sows. All polymorphisms affected the reproductive performance in the outbreak, but not during the endemic phase, thereby indicating a potential use of these markers for resilience. Moreover, some genotypes were associated with a stable performance across PRRSV phases. Thus, in the outbreak, the SGK1_rs338508371 AA sows had less piglets born alive (p < 0.0001) and more stillborns (p < 0.05) while other sows were able to keep their productivity. During the outbreak, TAP1_rs80928141 GG sows had less piglets born alive (p < 0.05) and both TAP1 polymorphisms influenced the number of mummies in an additive manner (p < 0.05). Remarkably, TAP1_rs80928141 AA sows had around one mummy more than GG sows (p < 0.01). Resilience to PRRSV could be improved by including the SGK1 and TAP1 markers in crossbreeding and/or selection schemes, as they contribute to maintaining a stable number of piglets born alive and lost, particularly mummies, despite the outbreak.

iTRAQ-based high-throughput proteomics analysis reveals alterations of plasma proteins in patients infected with human bocavirus

Human bocavirus (HBoV) is a member of the genus Bocavirus, family Parvoviridae, and subfamily Parvovirus and was first identified in nasopharyngeal aspirates of Swedish children with acute respiratory tract infection (ARTI) in 2005. It is the causative agent of nasopharyngeal aspirate disease and death in children. The HboV genomic structure is a linear single-stranded DNA (ssDNA). Its clinical pathogenic characteristics have been extensively studied, however, at present the molecular mechanism underlying the pathogenesis of HBoV infection is not completely clear. In this study, a total of 293 differentially expressed proteins (DEPs) between ARTI cases and healthy plasma samples were characterized using isobaric tags for relative and absolute quantitation (iTRAQ)-coupled bioinformatics analysis, among which 148 were up-regulated and 135 were down-regulated. Gene Ontology (GO) and Cluster of Orthologous Groups of proteins (COG) annotated an enrichment of DEPs in complement activation and biological processes like immunity, inflammation, signal transduction, substance synthesis, and metabolism. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis enriched DEPs mainly in the Wnt signaling pathway (ko04310), PPAR signaling pathway (ko03320), intestinal immune network for IgA production (ko04672), complement and coagulation cascades (ko04610), Toll-like receptor signaling pathway (ko04620) and B cell receptor signaling pathway (ko04662). Further, expression levels of three candidate proteins (upregulated PPP2R1A and CUL1, and downregulated CETP) were validated using western blotting. Our investigation is the first analysis of the proteomic profile of HBoV-infected ARTI cases using the iTRAQ approach, providing a foundation for a better molecular understanding of the pathogenesis of ARTI in children.

Identification of three site mutations in nonstructural protein 1β, glycoprotein 3 and glycoprotein 5 that correlate with increased interferon α resistance of porcine reproductive and respiratory syndrome virus

Porcine reproductive and respiratory syndrome virus (PRRSV) is an economically significant pathogen that has been recognized for its genetic variation, rapid evolution, and immune suppression. Type I interferons (IFNs) play an important role in host defense against viral infection by inducing many antiviral effectors, which might be a selective pressure driving viral evolution towards IFN resistance. To investigate the IFN resistance-related variation of PRRSV genome under IFN selective pressure and explore the molecular mechanism of IFN sensitivity changes, PRRSV strain JXwn06 was serially propagated in porcine pulmonary alveolar macrophages (PAMs) with IFNα treatment for 45 passages and 3 rounds of purification. Four mutant strains named JX-αP51n (n = 1, 2, 3 and 4) with reduced IFNα sensitivity were selected; the strains showed a 100-fold higher titer than the passaging-control strain JX-P51 in IFNα-treated PAMs. IFNα-resistant strains were found to antagonize the IFNα-activated JAK-STAT signaling pathway to a greater extent than the nonresistant strain by down-regulating the expression level of IFNα-activated pJAK1 through interfering with phosphatase. Furthermore, the PRRSV genetic variations interacting with IFNα were identified by full genomic sequencing and alignment. Among these mutations, amino acid substitutions in nsp1β (E87 G), GP3 (F143 L) and GP5 (Y136 H) were found to correlate with increased IFNα resistance by enhancing the suppression effect on pJAK1, which could be further increased if these three substitution sites were combined. These findings provide some novel evidence for understanding PRRSV genetic variation under host selective pressure and viral evolution strategies to evade the host innate immune response.

Transcriptional profiles of PBMCs from pigs infected with three genetically diverse porcine reproductive and respiratory syndrome virus strains

Porcine reproductive and respiratory syndrome virus is the cause of reproductive failure in sows and respiratory disease in young pigs, which has been considered as one of the most costly diseases to the worldwide pig industry for almost 30 years. This study used microarray-based transcriptomic analysis of PBMCs from experimentally infected pigs to explore the patterns of immune dysregulation after infection with two East European PRRSV strains from subtype 2 (BOR and ILI) in comparison to a Danish subtype 1 strain (DAN). Transcriptional profiles were determined at day 7 post infection in three tested groups of pigs and analysed in comparison with the expression profile of control group. Microarray analysis revealed differential regulation (> 1.5-fold change) of 4253 and 7335 genes in groups infected with BOR and ILI strains, respectively, and of 12518 genes in pigs infected with Danish strain. Subtype 2 PRRSV strains showed greater induction of many genes, especially those involved in innate immunity, such as interferon stimulated antiviral genes and inflammatory markers. Functional analysis of the microarray data revealed a significant up-regulation of genes involved in processes such as acute phase response, granulocyte and agranulocyte adhesion and diapedesis, as well as down-regulation of genes enrolled in pathways engaged in protein synthesis, cell division, as well as B and T cell signaling. This study provided an insight into the host response to three different PRRSV strains at a molecular level and demonstrated variability between strains of different pathogenicity level.

MicroRNA expression profiling in alveolar macrophages of indigenous Chinese Tongcheng pigs infected with PRRSV in vivo

Porcine respiratory and reproductive syndrome (PRRS), caused by PRRS virus (PRRSV), is one of the most serious infectious diseases in the swine industry worldwide. Indigenous Chinese Tongcheng (TC) pigs reportedly show strong resistance to PRRSV infection. The miRNA expression profiles of porcine alveolar macrophages (PAMs) of control TC pigs and those infected with PRRSV in vivo were analyzed by high-throughput sequencing to explore changes induced by infection. A total of 182 known miRNAs including 101 miRNA-5p and 81 miRNA-3p were identified with 23 up-regulated differentially expressed miRNAs (DEmiRNAs) and 25 down-regulated DEmiRNAs. Gene Ontology analysis showed that predicted target genes for the DEmiRNAs were enriched in immune response, transcription regulation, and cell death. The integrative analysis of mRNA and miRNA expression revealed that down-regulated methylation-related genes (DNMT1 and DNMT3b) were targeted by five up-regulated DEmiRNAs. Furthermore, 35 pairs of miRNAs (70 miRNAs) were co-expressed after PRRSV infection and six pairs were co-expressed differently. Our results describe miRNA expression profiles of TC pigs in response to PRRSV infection and lay a strong foundation for developing novel therapies to control PRRS in pigs.

Antagonizing cytokine-mediated JAK-STAT signaling by porcine reproductive and respiratory syndrome virus

Janus kinase (JAK)-signal transducer and activator of transcription (STAT) signaling pathway is activated by myriad cytokines, which are involved in regulation of cell growth, proliferation, differentiation, apoptosis, angiogenesis, immunity and inflammatory response. Because of its significance in immune response, JAK-STAT pathway is often targeted by pathogens, including porcine reproductive and respiratory syndrome virus (PRRSV). PRRSV causes reproductive failure in sows and respiratory disease in pigs of all ages. A typical feature of the immune response to PRRSV infection in pigs is delayed production and low titer of virus neutralizing antibodies, and weak cell-mediated immune response. One of the possible reasons for the weak protective immune response is that PRRSV interferes with cytokine-mediated JAK-STAT signaling. PRRSV inhibits interferon-activated JAK-STAT signaling by blocking nuclear translocation of STAT1 and STAT2. The mechanism is that PRRSV non-structural protein 1β (nsp1β) induces degradation of karyopherin α1 (KPNA1), a critical adaptor in nucleo-cytoplasmic transport. PRRSV also antagonizes IL6-activated JAK-STAT3 signaling via inducing degradation of STAT3. In this review, we briefly introduce JAK-STAT signaling, summarize the PRRSV interference with it, and provide perspective on the perturbation in the context of PRRSV-elicited immune response.

Effect of single nucleotide polymorphism on the total number of piglets born per parity of three different pig breeds

Objective

To determine the effects of genomic breeding values (GBV) and single nucleotide polymorphisms (SNP) on the total number of piglets born (TNB) in 3 pig breeds (Berkshire, Landrace, and Yorkshire).

Methods

After collecting genomic information (Porcine SNP BeadChip) and phenotypic TNB records for each breed, the effects of GBV and SNP were estimated by using single step best linear unbiased prediction (ssBLUP) method.

Results

The heritability estimates for TNB in Berkshire, Landrace, and Yorkshire breeds were 0.078, 0.107, and 0.121, respectively. The breeding value estimates for TNB in Berkshire, Landrace, and Yorkshire breeds were in the range of −1.34 to 1.47 heads, −1.79 to 1.87 heads, and −2.60 to 2.94 heads, respectively. Of sows having records for TNB, the reliability of breeding value for individuals with SNP information was higher than that for individuals without SNP information. Distributions of the SNP effects on TNB did not follow gamma distribution. Most SNP effects were near zero. Only a few SNPs had large effects. The numbers of SNPs with absolute value of more than 4 standard deviations in Berkshire, Landrace, and Yorkshire breeds were 11, 8, and 19, respectively. There was no SNP with absolute value of more than 5 standard deviations in Berkshire or Landrace. However, in Yorkshire, four SNPs (ASGA 0089457, ASGA0103374, ALGA0111816, and ALGA0098882) had absolute values of more than 5 standard deviations.

Conclusion

There was no common SNP with large effect among breeds. This might be due to the large genetic composition differences and the small size of reference population. For the precise evaluation of genetic performance of individuals using a genomic selection method, it may be necessary to establish the appropriate size of reference population.

Transcriptome Differences in Porcine Alveolar Macrophages from Tongcheng and Large White Pigs in Response to Highly Pathogenic Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Infection

Porcine reproductive and respiratory syndrome virus (PRRSV) is a single-stranded positive-sense RNA virus that can cause devastating reproductive failure and respiratory tract lesions, which has led to serious damage to the swine industry worldwide. Our previous studies have indicated that Tongcheng (TC) pigs, a Chinese local breed, have stronger resistance or tolerance to PRRSV infection than Large White (LW) pigs. This study aims to investigate their host transcriptome differences in porcine alveolar macrophages (PAMs) at 7 days post challenge. Transcriptome profiling of PAMs from PRRSV infected and control pigs of these two breeds were performed using RNA-sequencing. For both breeds, there were 1257 common differentially expressed genes (DEGs) in response to PRRSV infection, involving hepatic fibrosis/hepatic stellate cell activation, phospholipase C, and granulocyte adhesion and diapedesis pathways. For TC pig, 549 specific DEGs were identified, including VAV2, BCL2 and BAX, which were enriched in activation of leukocyte extravasation and suppression of apoptosis. While, 898 specific DEGs were identified in LW pigs, including GNAQ, GNB5, GNG2, CALM4 and RHOQ, which were involved in suppression of Gαq and PI3K-AKT signaling. This study provides an insight into the transcriptomic comparison of resistant and susceptible pigs to PRRSV infection. TC pigs may promote the extravasation and migration of leukocytes to defend against PRRSV infections and suppress apoptosis of the infected macrophages to increase antigen presentation, thereby reducing the lung lesions.

PBMC transcriptome profiles identifies potential candidate genes and functional networks controlling the innate and the adaptive immune response to PRRSV vaccine in Pietrain pig

The porcine reproductive and respiratory syndrome (PRRS) is a devastating viral disease affecting swine production, health and welfare throughout the world. A synergistic action of the innate and the adaptive immune system of the host is essential for mounting a durable protective immunity through vaccination. Therefore, the current study aimed to investigate the transcriptome profiles of peripheral blood mononuclear cells (PBMCs) to characterize the innate and the adaptive immune response to PRRS Virus (PRRSV) vaccination in Pietrain pigs. The Affymetrix gene chip porcine gene 1.0 ST array was used for the transcriptome profiling of PBMCs collected at immediately before (D0), at one (D1) and 28 days (D28) post PRRSV vaccination with three biological replications. With FDR <0.05 and log2 fold change ±1.5 as cutoff criteria, 295 and 115 transcripts were found to be differentially expressed in PBMCs during the stage of innate and adaptive response, respectively. The microarray expression results were technically validated by qRT-PCR. The gene ontology terms such as viral life cycle, regulation of lymphocyte activation, cytokine activity and inflammatory response were enriched during the innate immunity; cytolysis, T cell mediated cytotoxicity, immunoglobulin production were enriched during adaptive immunity to PRRSV vaccination. Significant enrichment of cytokine-cytokine receptor interaction, signaling by interleukins, signaling by the B cell receptor (BCR), viral mRNA translation, IFN-gamma pathway and AP-1 transcription factor network pathways were indicating the involvement of altered genes in the antiviral defense. Network analysis revealed that four network modules were functionally involved with the transcriptional network of innate immunity, and five modules were linked to adaptive immunity in PBMCs. The innate immune transcriptional network was found to be regulated by LCK, STAT3, ATP5B, UBB and RSP17. While TGFß1, IL7R, RAD21, SP1 and GZMB are likely to be predictive for the adaptive immune transcriptional response to PRRSV vaccine in PBMCs. Results of the current immunogenomics study advances our understanding of PRRS in term of host-vaccine interaction, and thereby contribute to design a rationale for disease control strategy.

Genetic resistance - an alternative for controlling PRRS?

PRRS is one of the most challenging diseases for world-wide pig production. Attempts for a sustainable control of this scourge by vaccination have not yet fully satisfied. With an increasing knowledge and methodology in disease resistance, a new world-wide endeavour has been started to support the combat of animal diseases, based on the existence of valuable gene variants with regard to any host-pathogen interaction. Several groups have produced a wealth of evidence for natural variability in resistance/susceptibility to PRRS in our commercial breeding lines. However, up to now, exploiting existing variation has failed because of the difficulty to detect the carriers of favourable and unfavourable alleles, especially with regard to such complex polygenic traits like resistance to PRRS. New hope comes from new genomic tools like next generation sequencing which have become extremely fast and low priced. Thus, research is booming world-wide and the jigsaw puzzle is filling up – slowly but steadily. On the other hand, knowledge from virological and biomedical basic research has opened the way for an “intervening way”, i.e. the modification of identified key genes that occupy key positions in PRRS pathogenesis, like CD163. CD163 was identified as the striking receptor in PRRSV entry and its knockout from the genome by gene editing has led to the production of pigs that were completely resistant to PRRSV – a milestone in modern pig breeding. However, at this early step, concerns remain about the acceptance of societies for gene edited products and regulation still awaits upgrading to the new technology. Further questions arise with regard to upcoming patents from an ethical and legal point of view. Eventually, the importance of CD163 for homeostasis, defence and immunity demands for more insight before its complete or partial silencing can be answered. Whatever path will be followed, even a partial abolishment of PRRSV replication will lead to a significant improvement of the disastrous herd situation, with a significant impact on welfare, performance, antimicrobial consumption and consumer protection. Genetics will be part of a future solution.

Genome-wide analysis of the transcriptional response to porcine reproductive and respiratory syndrome virus infection at the maternal/fetal interface and in the fetus

Background

Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) infection of pregnant pigs can result in congenital infection and ultimately fetal death. Little is known about immune responses to infection at the maternal-fetal interface and in the fetus itself, or the molecular events behind virus transmission and disease progression in the fetus. To investigate these processes, RNA-sequencing of two sites, uterine endothelium with adherent placental tissue and fetal thymus, was performed 21 days post-challenge on four groups of fetuses selected from a large PRRSV challenge experiment of pregnant gilts: control (CON), uninfected (UNINF), infected (INF), and meconium-stained (MEC) (n = 12/group). Transcriptional analyses consisted of multiple contrasts between groups using two approaches: differential gene expression analysis and weighted gene co-expression network analysis (WGCNA). Biological functions, pathways, and regulators enriched for differentially expressed genes or module members were identified through functional annotation analyses. Expression data were validated by reverse transcription quantitative polymerase chain reaction (RTqPCR) carried out for 16 genes of interest.

Results

The immune response to infection in endometrium was mainly adaptive in nature, with the most upregulated genes functioning in either humoral or cell-mediated immunity. In contrast, the expression profile of infected fetal thymus revealed a predominantly innate immune response to infection, featuring the upregulation of genes regulated by type I interferon and pro-inflammatory cytokines. Fetal infection was associated with an increase in viral load coupled with a reduction in T cell signaling in the endometrium that could be due to PRRSV-controlled apoptosis of uninfected bystander cells. There was also evidence for a reduction in TWIST1 activity, a transcription factor involved in placental implantation and maturation, which could facilitate virus transmission or fetal pathology through dysregulation of placental function. Finally, results suggested that events within the fetus could also drive fetal pathology. Thymus samples of meconium-stained fetuses exhibited an increase in the expression of pro-inflammatory cytokine and granulocyte genes previously implicated in swine infectious disease pathology.

Conclusions

This study identified major differences in the response to PRRSV infection in the uterine endometrium and fetus at the gene expression level, and provides insight into the molecular basis of virus transmission and disease progression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2720-4) contains supplementary material, which is available to authorized users.

Differences in Whole Blood Gene Expression Associated with Infection Time-Course and Extent of Fetal Mortality in a Reproductive Model of Type 2 Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Infection

Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) infection of pregnant females causes fetal death and increased piglet mortality, but there is substantial variation in the extent of reproductive pathology between individual dams. This study used RNA-sequencing to characterize the whole blood transcriptional response to type 2 PRRSV in pregnant gilts during the first week of infection (at 0, 2, and 6 days post-inoculation), and attempted to identify gene expression signatures associated with a low or high level of fetal mortality rates (LFM and HFM; n = 8/group) at necropsy, 21 days post-inoculation. The initial response to infection measured at 2 days post-inoculation saw an upregulation of genes involved in innate immunity, such as interferon-stimulated antiviral genes and inflammatory markers, and apoptosis. A concomitant decrease in expression of protein synthesis and T lymphocyte markers was observed. By day 6 the pattern had reversed, with a drop in innate immune signaling and an increase in the expression of genes involved in cell division and T cell signaling. Differentially expressed genes (DEGs) associated with extremes of litter mortality rate were identified at all three time-points. Among the 15 DEGs upregulated in LFM gilts on all three days were several genes involved in platelet function, including integrins ITGA2B and ITGB3, and the chemokine PF4 (CXCL4). LFM gilts exhibited a higher baseline expression of interferon-stimulated and pro-inflammatory genes prior to infection, and of T cell markers two days post-infection, indicative of a more rapid progression of the immune response to PRRSV. This study has increased our knowledge of the early response to PRRSV in the blood of pregnant gilts, and could ultimately lead to the development of a biomarker panel that can be used to predict PRRSV-associated reproductive pathology.

Porcine Reproductive and Respiratory Syndrome Virus (PRRSV): Pathogenesis and Interaction with the Immune System

This review addresses important issues of porcine reproductive and respiratory syndrome virus (PRRSV) infection, immunity, pathogenesis, and control. Worldwide, PRRS is the most economically important infectious disease of pigs. We highlight the latest information on viral genome structure, pathogenic mechanisms, and host immunity, with a special focus on immune factors that modulate PRRSV infections during the acute and chronic/persistent disease phases. We address genetic control of host resistance and probe effects of PRRSV infection on reproductive traits. A major goal is to identify cellular/viral targets and pathways for designing more effective vaccines and therapeutics. Based on progress in viral reverse genetics, host transcriptomics and genomics, and vaccinology and adjuvant technologies, we have identified new areas for PRRS control and prevention. Finally, we highlight the gaps in our knowledge base and the need for advanced molecular and immune tools to stimulate PRRS research and field applications.

Whole blood microarray analysis of pigs showing extreme phenotypes after a porcine reproductive and respiratory syndrome virus infection

Background

The presence of variability in the response of pigs to Porcine Reproductive and Respiratory Syndrome virus (PRRSv) infection, and recent demonstration of significant genetic control of such responses, leads us to believe that selection towards more disease resistant pigs could be a valid strategy to reduce its economic impact on the swine industry. To find underlying molecular differences in PRRS susceptible versus more resistant pigs, 100 animals with extremely different growth rates and viremia levels after PRRSv infection were selected from a total of 600 infected pigs. A microarray experiment was conducted on whole blood RNA samples taken at 0, 4 and 7 days post infection (dpi) from these pigs. From these data, we examined associations of gene expression with weight gain and viral load phenotypes. The single nucleotide polymorphism (SNP) marker WUR10000125 (WUR) on the porcine 60 K SNP chip was shown to be associated with viral load and weight gain after PRRSv infection, and so the effect of the WUR10000125 (WUR) genotype on expression in whole blood was also examined.

Results

Limited information was obtained through linear modeling of blood gene differential expression (DE) that contrasted pigs with extreme phenotypes, for growth or viral load or between animals with different WUR genotype. However, using network-based approaches, molecular pathway differences between extreme phenotypic classes could be identified. Several gene clusters of interest were found when Weighted Gene Co-expression Network Analysis (WGCNA) was applied to 4dpi contrasted with 0dpi data. The expression pattern of one such cluster of genes correlated with weight gain and WUR genotype, contained numerous immune response genes such as cytokines, chemokines, interferon type I stimulated genes, apoptotic genes and genes regulating complement activation. In addition, Partial Correlation and Information Theory (PCIT) identified differentially hubbed (DH) genes between the phenotypically divergent groups. GO enrichment revealed that the target genes of these DH genes are enriched in adaptive immune pathways.

Conclusion

There are molecular differences in blood RNA patterns between pigs with extreme phenotypes or with a different WUR genotype in early responses to PRRSv infection, though they can be quite subtle and more difficult to discover with conventional DE expression analyses. Co-expression analyses such as WGCNA and PCIT can be used to reveal network differences between such extreme response groups.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1741-8) contains supplementary material, which is available to authorized users.

Comparative analysis of cytokine transcript profiles within mediastinal lymph node compartments of pigs after infection with porcine reproductive and respiratory syndrome genotype 1 strains differing in pathogenicity

Porcine reproductive and respiratory syndrome virus (PRRSV) induces a weak immune response enabling it to persist in different organs of infected pigs. This has been attributed to the ability of PRRSV to influence the induction of cytokine responses. In this study, we investigated the cytokine transcriptional profiles in different compartments of the mediastinal lymph node of pigs infected with three genotype 1 PRRSV strains of differing pathogenicity: the low virulence prototype Lelystad virus (LV), and UK field strain 215–06 and the highly virulent subtype 3 SU1-Bel isolate from Belarus. We have used a combination of laser capture micro-dissection (LCM) followed by real time quantitative PCR (RT-qPCR) and immunohistochemical (IHC) detection of immune cell markers (CD3, CD79a and MAC387) and RT-qPCR quantification of PRRSV and cytokine transcripts. Compared to mock infected pigs, we found a significant downregulation of TNF-α and IFN-α in follicular and interfollicular areas of the mediastinal lymph node from 3 days post-infection (dpi) in animals infected with all three strains. This was accompanied by a transient B cell depletion and T cell and macrophage infiltration in the follicles together with T cell depletion in the interfollicular areas. A delayed upregulation of IFN-γ and IL-23p19 was observed mainly in the follicles. The PRRSV load was higher in all areas and time-points studied in the animals infected with the SU1-Bel strain. This paper describes the first application of LCM to study the cytokine transcript profiles and virus distribution in different compartments of the lymph node of pigs.

Current transcriptomics in pig immunity research

Swine performance in the face of disease challenge is becoming progressively more important. To improve the pig’s robustness and resilience against pathogens through selection, a better understanding of the genetic and epigenetic factors in the immune response is required. This review highlights results from the most recent transcriptome research, and the meta-analyses performed, in the context of pig immunity. A technological overview is given including wholegenome microarrays, immune-specific arrays, small-scale high-throughput expression methods, high-density tiling arrays, and next generation sequencing (NGS). Although whole genome microarray techniques will remain complementary to NGS for some time in domestic species, research will transition to sequencing-based methods due to cost-effectiveness and the extra information that such methods provide. Furthermore, upcoming high-throughput epigenomic studies, which will add greatly to our knowledge concerning the impact of epigenetic modifications on pig immune response, are listed in this review. With emphasis on the insights obtained from transcriptomic analyses for porcine immunity, we also discuss the experimental design in pig immunity research and the value of the newly published porcine genome assembly in using the pig as a model for human immune response. We conclude by discussing the importance of establishing community standards to maximize the possibility of integrative computational analyses, such as was clearly beneficial for the human ENCODE project.

Applications of omics approaches to the development of microbiological risk assessment using RNA virus dose-response models as a case study

T e in the amount of ‘omics’ data available and in our ability to interpret those data. The aim of this paper was to consider how omics techniques can be used to improve and refine microbiological risk assessment, using dose–response models for RNA viruses, with particular reference to norovirus through the oral route as the case study. The dose–response model for initial infection in the gastrointestinal tract is broken down into the component steps at the molecular level and the feasibility of assigning probabilities to each step assessed. The molecular mechanisms are not sufficiently well understood at present to enable quantitative estimation of probabilities on the basis of omics data. At present, the great strength of gene sequence data appears to be in giving information on the distribution and proportion of susceptible genotypes (for example due to the presence of the appropriate pathogen‐binding receptor) in the host population rather than in predicting specificities from the amino acid sequences concurrently obtained. The nature of the mutant spectrum in RNA viruses greatly complicates the application of omics approaches to the development of mechanistic dose–response models and prevents prediction of risks of disease progression (given infection has occurred) at the level of the individual host. However, molecular markers in the host and virus may enable more broad predictions to be made about the consequences of exposure in a population. In an alternative approach, comparing the results of deep sequencing of RNA viruses in the faeces/vomitus from donor humans with those from their infected recipients may enable direct estimates of the average probability of infection per virion to be made.

Genetic variations of TAP1 gene exon 3 affects gene expression and Escherichia coli F18 resistance in piglets

Firstly, our research group identified Sutai pigs’ phenotypes that exhibited extreme resistance and susceptibility to the Escherichia coli F18 respectively, and then eight ETEC (Enterotoxigenic Escherichia coli) F18-resistant piglets and eight ETEC F18-sensitive piglets were selected. Then, the TAP1 (Transporter associated with antigen processing) mRNA relative expression levels were analyzed in 11 tissues of the resistant and susceptible phenotypes. Simultaneously, we detected the genetic variations in exon 3 of the TAP1 gene and evaluated the TAP1 mRNA expression levels among the different genotype pigs to study the effects of the genetic variation on gene expression, and the E. coli F18 resistance. The results revealed higher expression levels in the resistant genotypes than that in the susceptible genotypes in 11 tissues, with significant differences in the spleen, lymph node, lung, thymus, duodenum and jejunum. Furthermore, a G729A mutation was identified in the TAP1 gene exon 3, and this mutation deviates from Hardy-Weinberg equilibrium (p < 0.01). The TAP1 mRNA levels in GG genotype were significantly higher than that in the other two genotypes, with significant differences in the liver, lung, kidney, thymus, lymph node, duodenum and jejunum tissues. We speculated that high expression of the TAP1 gene might confer resistance against the E. coli F18, the G729A mutation had a significant effect on the mRNA expression, and individuals with the GG genotype possessed a stronger ability to resist the E. coli F18 infection.

Footprints of directional selection in wild Atlantic salmon populations: evidence for parasite-driven evolution?

Mechanisms of host-parasite co-adaptation have long been of interest in evolutionary biology; however, determining the genetic basis of parasite resistance has been challenging. Current advances in genome technologies provide new opportunities for obtaining a genome-scale view of the action of parasite-driven natural selection in wild populations and thus facilitate the search for specific genomic regions underlying inter-population differences in pathogen response. European populations of Atlantic salmon (Salmo salar L.) exhibit natural variance in susceptibility levels to the ectoparasite Gyrodactylus salaris Malmberg 1957, ranging from resistance to extreme susceptibility, and are therefore a good model for studying the evolution of virulence and resistance. However, distinguishing the molecular signatures of genetic drift and environment-associated selection in small populations such as land-locked Atlantic salmon populations presents a challenge, specifically in the search for pathogen-driven selection. We used a novel genome-scan analysis approach that enabled us to i) identify signals of selection in salmon populations affected by varying levels of genetic drift and ii) separate potentially selected loci into the categories of pathogen (G. salaris)-driven selection and selection acting upon other environmental characteristics. A total of 4631 single nucleotide polymorphisms (SNPs) were screened in Atlantic salmon from 12 different northern European populations. We identified three genomic regions potentially affected by parasite-driven selection, as well as three regions presumably affected by salinity-driven directional selection. Functional annotation of candidate SNPs is consistent with the role of the detected genomic regions in immune defence and, implicitly, in osmoregulation. These results provide new insights into the genetic basis of pathogen susceptibility in Atlantic salmon and will enable future searches for the specific genes involved.

Genome-wide analysis of antiviral signature genes in porcine macrophages at different activation statuses

Macrophages (MФs) can be polarized to various activation statuses, including classical (M1), alternative (M2), and antiviral states. To study the antiviral activation status of porcine MФs during porcine reproductive and respiratory syndrome virus (PRRSV) infection, we used RNA Sequencing (RNA-Seq) for transcriptomic analysis of differentially expressed genes (DEGs). Sequencing assessment and quality evaluation showed that our RNA-Seq data met the criteria for genome-wide transcriptomic analysis. Comparisons of any two activation statuses revealed more than 20,000 DEGs that were normalized to filter out 153–5,303 significant DEGs [false discovery rate (FDR) ≤0.001, fold change ≥2] in each comparison. The highest 5,303 significant DEGs were found between lipopolysaccharide- (LPS) and interferon (IFN)γ-stimulated M1 cells, whereas only 153 significant DEGs were detected between interleukin (IL)-10-polarized M2 cells and control mock-activated cells. To identify signature genes for antiviral regulation pertaining to each activation status, we identified a set of DEGs that showed significant up-regulation in only one activation state. In addition, pathway analyses defined the top 20–50 significantly regulated pathways at each activation status, and we further analyzed DEGs pertinent to pathways mediated by AMP kinase (AMPK) and epigenetic mechanisms. For the first time in porcine macrophages, our transcriptomic analyses not only compared family-wide differential expression of most known immune genes at different activation statuses, but also revealed transcription evidence of multiple gene families. These findings show that using RNA-Seq transcriptomic analyses in virus-infected and status-synchronized macrophages effectively profiled signature genes and gene response pathways for antiviral regulation, which may provide a framework for optimizing antiviral immunity and immune homeostasis.

Genome-wide gene expression profiles in lung tissues of pig breeds differing in resistance to porcine reproductive and respiratory syndrome virus

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.

In vitro and ex vivo analyses of co-infections with swine influenza and porcine reproductive and respiratory syndrome viruses

Viral respiratory diseases remain problematic in swine. Among viruses, porcine reproductive and respiratory syndrome virus (PRRSV) and swine influenza virus (SIV), alone or in combination, are the two main known contributors to lung infectious diseases. Previous studies demonstrated that experimental dual infections of pigs with PRRSV followed by SIV can cause more severe disease than the single viral infections. However, our understanding of the impact of one virus on the other at the molecular level is still extremely limited. Thus, the aim of the current study was to determine the influence of dual infections, compared to single infections, in porcine alveolar macrophages (PAMs) and precision cut lung slices (PCLS). PAMs were isolated and PCLS were acquired from the lungs of healthy 8-week-old pigs. Then, PRRSV (ATCC VR-2385) and a local SIV strain of H1N1 subtype (A/Sw/Saskatchewan/18789/02) were applied simultaneously or with 3 h apart on PAMs and PCLS for a total of 18 h. Immuno-staining for both viruses and beta-tubulin, real-time quantitative PCR and ELISA assays targeting various genes (pathogen recognition receptors, interferons (IFN) type I, cytokines, and IFN-inducible genes) and proteins were performed to analyze the cell and the tissue responses. Interference caused by the first virus on replication of the second virus was observed, though limited. On the host side, a synergistic effect between PRRSV and SIV co-infections was observed for some transcripts such as TLR3, RIG-I, and IFNβ in PCLS. The PRRSV infection 3 h prior to SIV infection reduced the response to SIV while the SIV infection prior to PRRSV infection had limited impact on the second infection. This study is the first to show an impact of PRRSV/SIV co-infection and superinfections in the cellular and tissue immune response at the molecular level. It opens the door to further research in this exciting and intriguing field.

Characterization of the microRNAome in porcine reproductive and respiratory syndrome virus infected macrophages

Porcine Reproductive and Respiratory Syndrome Virus (PRRSV), a member of the arterivirus family, is the causative agent of Porcine Reproductive and Respiratory Syndrome (PRRS). PRRS is characterized by late term abortions and respiratory disease, particularly in young pigs. Small regulatory RNAs termed microRNA (miRNA) are associated with gene regulation at the post-transcriptional level. MiRNAs are known to play many diverse and complex roles in viral infections. To discover the impact of PRRSV infections on the cellular miRNAome, Illumina deep sequencing was used to construct small RNA expression profiles from in vitro cultured PRRSV-infected porcine alveolar macrophages (PAMs). A total of forty cellular miRNAs were significantly differentially expressed within the first 48 hours post infection (hpi). The expression of six miRNAs, miR-30a-3p, miR-132, miR-27b*, miR-29b, miR-146a and miR-9-2, were altered at more than one time point. Target gene identification suggests that these miRNAs are involved in regulating immune signaling pathways, cytokine, and transcription factor production. The most highly repressed miRNA at 24 hpi was miR-147. A miR-147 mimic was utilized to maintain miR-147 levels in PRRSV-infected PAMs. PRRSV replication was negatively impacted by high levels of miR-147. Whether down-regulation of miR-147 is directly induced by PRRSV or if it is part of the cellular response and PRRSV indirectly benefits remains to be determined. No evidence could be found of PRRSV-encoded miRNAs. Overall, the present study has revealed that a large and diverse group of miRNAs are expressed in swine alveolar macrophages and that the expression of a subset of these miRNAs is altered in PRRSV infected macrophages.

Functional genetic polymorphisms in PP2A subunit genes confer increased risks of lung cancer in southern and eastern Chinese

Protein phosphatase-2A (PP2A) is one of the major cellular serine-threonine phosphatases and functions as a tumor suppressor that negatively regulates the activity of some oncogenic kinases. Recent studies have reported that PP2A expression was suppressed during lung carcinogenesis, we there hypothesized that the single nucleotide polymorphisms (SNPs) in PP2A subunit genes may affect PP2A function and thus contribute to lung cancer susceptibility. In a two-stage case-control study with a total of 1559 lung cancer patients and 1679 controls, we genotyped eight putative functional SNPs and one identified functional SNP (i.e., rs11453459) in seven major PP2A subunits (i.e., PPP2R1A, PPP2R1B, PPP2CA, PPP2R2A, PPP2R2B, PPP2R5C, PPP2R5E) in southern and eastern Chinese. We found that rs11453459G (-G/GG) variant genotypes of PPP2R1A and the rs1255722AA variant genotype of PPP2R5E conferred increased risks of lung cancer (rs11453459, -G/GG vs. –: OR = 1.31, 95% CI = 1.13–1.51; rs1255722, AA vs. AG/GG: OR = 1.27, 95% CI = 1.07–1.51). After combined the two variants, the number of the adverse genotypes was positively associated with lung cancer risk in a dose-response manner (P_trend  = 5.63×10⁻⁶). Further functional assay showed that lung cancer tissues carrying rs1255722AA variant genotype had a significantly lower mRNA level of PPP2R5E compared with tissues carrying GG/GA genotypes. However, such effect was not observed for the other SNPs and other combinations. Our findings suggested that the two functional variants in PPP2R1A and PPP2R5E and their combination are associated with lung cancer risk in Chinese, which may be valuable biomarkers to predict risk of lung cancer.

Significant association between SNPs in the superoxide dismutase 3, extracellular (SOD3) gene and resistance to Aeromonas hydrophila in the freshwater mussel Hyriopsis cumingii

Extracellular superoxide dismutase (SOD3) is a major antioxidant enzyme that protects organs from damage by reactive oxygen species (ROS). In this study, the SOD3 gene was identified and characterized from the freshwater mussel Hyriopsis cumingii (Hc-SOD3). The cDNA sequence consists of 763 bp, encoding a protein of 208 amino acids. The amino acid sequence possesses two CuZnSOD signature sequences, and amino acids required for binding of Cu (His-93, -95, -110 and -169) and Zn (His-110, -118, -129 and Asp-132) were conserved in Hc-SOD3. The Hc-SOD3 genomic sequence was 9165 bp in length, containing four exons and three introns. Eighteen single nucleotide polymorphisms were detected in the Hc-SOD3 gene from resistant stock (RS) and susceptible stock (SS) of H. cumingii to Aeromonas hydrophila. The genotype and allele distribution were examined in resistant and susceptible stocks. Among them, a C/G substitution at the g.7994C>G locus and G/C substitution at the g.8087G>C locus were significantly associated with resistance/susceptibility of H. cumingii to A. hydrophila, both in genotype (P = 0.017, P = 0.004 respectively) and allele frequency (P = 0.021, P = 0.006 respectively). Linkage disequilibrium analysis revealed that g.7994C>G, g.8001A>G, g.8035G>A, g.8087G>C and g.8191T>A were in linkage disequilibrium. The results suggest that the two polymorphic loci, g.7994C>G and g.8087G>C, could be potential genetic markers for future molecular selection of strains that are resistant to diseases.

Characterizing differential individual response to porcine reproductive and respiratory syndrome virus infection through statistical and functional analysis of gene expression

We evaluated differences in gene expression in pigs from the Porcine Reproductive and Respiratory Syndrome (PRRS) Host Genetics Consortium initiative showing a range of responses to PRRS virus infection. Pigs were allocated into four phenotypic groups according to their serum viral level and weight gain. RNA obtained from blood at 0, 4, 7, 11, 14, 28, and 42 days post-infection (DPI) was hybridized to the 70-mer 20K Pigoligoarray. We used a blocked reference design for the microarray experiment. This allowed us to account for individual biological variation in gene expression, and to assess baseline effects before infection (0 DPI). Additionally, this design has the flexibility of incorporating future data for differential expression analysis. We focused on evaluating transcripts showing significant interaction of weight gain and serum viral level. We identified 491 significant comparisons [false discovery rate (FDR) = 10%] across all DPI and phenotypic groups. We corroborated the overall trend in direction and level of expression (measured as fold change) at 4 DPI using qPCR (r = 0.91, p ≤ 0.0007). At 4 and 7 DPI, network and functional analyses were performed to assess if immune related gene sets were enriched for genes differentially expressed (DE) across four phenotypic groups. We identified cell death function as being significantly associated (FDR ≤ 5%) with several networks enriched for DE transcripts. We found the genes interferon-alpha 1(IFNA1), major histocompatibility complex, class II, DQ alpha 1 (SLA-DQA1), and major histocompatibility complex, class II, DR alpha (SLA-DRA) to be DE (p ≤ 0.05) between phenotypic groups. Finally, we performed a power analysis to estimate sample size and sampling time-points for future experiments. We concluded the best scenario for investigation of early response to PRRSV infection consists of sampling at 0, 4, and 7 DPI using about 30 pigs per phenotypic group.

SOCS proteins in infectious diseases of mammals

As for most biological processes, the immune response to microbial infections has to be tightly controlled to remain beneficial for the host. Inflammation is one of the major consequences of the host's immune response. For its orchestration, this process requires a fine-tuned interplay between interleukins, endothelial cells and various types of recruited immune cells. Suppressors of cytokine signalling (SOCS) proteins are crucially involved in the complex control of the inflammatory response through their actions on various signalling pathways including the JAK/STAT and NF-κB pathways. Due to their cytokine regulatory functions, they are frequent targets for exploitation by infectious agents trying to escape the host's immune response. This review article aims to summarize our current knowledge regarding SOCS family members in the different mammalian species studied so far, and to display their complex molecular interactions with microbial pathogens.