Bioinformatic analyses in early host response to Porcine Reproductive and Respiratory Syndrome virus (PRRSV) reveals pathway differences between pigs with alternate genotypes for a major host response QTL

Background

A region on Sus scrofa chromosome 4 (SSC4) surrounding single nucleotide polymorphism (SNP) marker WUR10000125 (WUR) has been reported to be strongly associated with both weight gain and serum viremia in pigs after infection with PRRS virus (PRRSV). A proposed causal mutation in the guanylate binding protein 5 gene (GBP5) is predicted to truncate the encoded protein. To investigate transcriptional differences between WUR genotypes in early host response to PRRSV infection, an RNA-seq experiment was performed on globin depleted whole blood RNA collected on 0, 4, 7, 10 and 14 days post-infection (dpi) from eight littermate pairs with one AB (favorable) and one AA (unfavorable) WUR genotype animal per litter.

Results

Gene Ontology (GO) enrichment analysis of transcripts that were differentially expressed (DE) between dpi across both genotypes revealed an inflammatory response for all dpi when compared to day 0. However, at the early time points of 4 and 7dpi, several GO terms had higher enrichment scores compared to later dpi, including inflammatory response (p < 10^-7), specifically regulation of NFkappaB (p < 0.01), cytokine, and chemokine activity (p < 0.01). At 10 and 14dpi, GO term enrichment indicated a switch to DNA damage response, cell cycle checkpoints, and DNA replication. Few transcripts were DE between WUR genotypes on individual dpi or averaged over all dpi, and little enrichment of any GO term was found. However, there were differences in expression patterns over time between AA and AB animals, which was confirmed by genotype-specific expression patterns of several modules that were identified in weighted gene co-expression network analyses (WGCNA). Minor differences between AA and AB animals were observed in immune response and DNA damage response (p = 0.64 and p = 0.11, respectively), but a significant effect between genotypes pointed to a difference in ion transport/homeostasis and the participation of G-coupled protein receptors (p = 8e-4), which was reinforced by results from regulatory and phenotypic impact factor analyses between genotypes.

Conclusion

We propose these pathway differences between WUR genotypes are the result of the inability of the truncated GBP5 of the AA genotyped pigs to inhibit viral entry and replication as quickly as the intact GBP5 protein of the AB genotyped pigs.

Electronic supplementary material

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

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.