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

Fuzhengjiedu San inhibits porcine reproductive and respiratory syndrome virus by activating the PI3K/AKT pathway

BACKGROUND

Traditional Chinese medicine (TCM) has been garnering ever-increasing worldwide attention as the herbal extracts and formulas prove to have potency against disease. Fuzhengjiedu San (FZJDS), has been extensively used to treat viral diseases in pigs, but its bioactive components and therapeutic mechanisms remain unclear.

METHODS

In this study, we conducted an integrative approach of network pharmacology and experimental study to elucidate the mechanisms underlying FZJDS's action in treating porcine reproductive and respiratory syndrome virus (PRRSV). We constructed PPI network and screened the core targets according to their degree of value. GO and KEGG enrichment analyses were also carried out to identify relevant pathways. Lastly, qRT-PCR, flow cytometry and western blotting were used to determine the effects of FZJDS on core gene expression in PRRSV-infected monkey kidney (MARC-145) cells to further expand the results of network pharmacological analysis.

RESULTS

Network pharmacology data revealed that quercetin, kaempferol, and luteolin were the main active compounds of FZJDS. The phosphatidylinositol-3-kinase (PI3K)/Akt pathway was deemed the cellular target as it has been shown to participate most in PRRSV replication and other PRRSV-related functions. Analysis by qRT-PCR and western blotting demonstrated that FZJDS significantly reduced the expression of P65, JNK, TLR4, N protein, Bax and IĸBa in MARC-145 cells, and increased the expression of Bcl-2, consistent with network pharmacology results. This study provides that FZJDS has significant antiviral activity through its effects on the PI3K/AKT signaling pathway.

CONCLUSION

We conclude that FZJDS is a promising candidate herbal formulation for treating PRRSV and deserves further investigation.

Differences in Immune Characteristics and Related Gene Expression in Spleen among Ningxiang, Berkshire Breeds and Their Hybrid Pigs

To investigate the differential immunology in Ningxiang and Berkshire pigs and their F1 offspring (F1 offspring), physiological and biochemical indicators in the plasma and spleen were analyzed. Then, transcriptomic analysis of the spleen identified 1348, 408, and 207 differentially expressed genes (DEGs) in comparisons of Ningxiang vs. Berkshire, Berkshire vs. F1 offspring, and Ningxiang vs. F1 offspring, respectively. In Ningxiang vs. Berkshire pigs, the gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that the DEGs included CD163, MARCO, CXCL14, CCL19, and PPBP, which are associated with immunity. GO and KEGG analyses were also conducted comparing F1 offspring and their parents. The DEGs, including BPIFB1, HAVCR2, CD163, DDX3X, CCR5, and ITGB3, were enriched in immune-related pathways. Protein-protein interaction (PPI) analysis indicated that the EGFR and ITGA2 genes were key hub genes. In conclusion, this study identifies significant immune DEGs in different pig breeds, providing data to support the exploration of breeding strategies for disease resistance in local and crossbred pig populations.

PRRSV alters m6A methylation and alternative splicing to regulate immune, extracellular matrix-associated function

The alternative splicing and N6-methyladenosine (m6A) modifications occurring during porcine reproductive and respiratory syndrome virus (PRRSV) infections remain poorly understood. Transcriptome and MeRIP-seq analyses were performed to identify the gene expression changes, splicing and m6A modifications in the lungs of PRRSV-infected pigs. In total, 1624 differentially expressed genes (DEGs) were observed between PRRSV-infected and uninfected pigs. We observed significant alterations in alternative splicing (54,367 events) and m6A modifications (2265 DASEs) in numerous genes, including LMO7, SLC25A27, ZNF185, and ECM1, during PRRSV infection. LMO7 and ZNF185 exhibited alternative splicing variants and reduced mRNA expression levels following PRRSV infection. Notably, LMO7 inhibited c-JUN, SMAD3, and FAK expression, whereas ZNF185 affected the expression of FAK, CDH1, and GSK3β downstream. Additionally, ECM1 influenced FAK expression by targeting ITGB3 and AKT2, suggesting its involvement in extracellular matrix accumulation through the ITGB3-AKT2/FAK pathway. These changes may facilitate viral invasion and replication by modulating the expression of genes and proteins participating in crucial cellular processes associated with immunity and the extracellular matrix. We highlight the importance of these genes and their associated pathways in PRRSV infections and suggest that targeting these may be a promising therapeutic approach for treating viral infections.

Integrative time-serial networks for genome-wide lncRNA-mRNA interactions reveal interferon-inducible antiviral and T-cell receptor regulations against PRRSV infection

Porcine reproductive and respiratory syndrome virus (PRRSV) infection severely affects the swine industry each year. Although the host mechanisms against PRRSV infection have been identified in key target tissues through whole transcriptome sequencing, specific molecular regulators have not been elucidated. Long non-coding RNA (lncRNA) expression is highly specific and could thus be used to effectively identify PRRSV-specific candidates. Here, we identified novel lncRNAs in lungs, bronchial lymph nodes, and tonsils after PRRSV infection and constructed phenotype-based integrative co-expression networks using time-series differentially expressed (DE) lncRNAs and mRNAs. After the analyses, a total of 309 lncRNA-mRNA interactions were identified. During early host innate signalling, interferon-inducible and interferon genes were positively regulated by specific lncRNA. Moreover, T-cell receptor genes in lung adaptive immune signalling were negatively regulated by specific lncRNA. Collectively, our findings provide insights into the genome-wide lncRNA-mRNA interactions and dynamic regulation of lncRNA-mediated mechanisms against PRRSV infection.

MARCO Inhibits Porcine Reproductive and Respiratory Syndrome Virus Infection through Intensifying Viral GP5-Induced Apoptosis

Studying viral glycoprotein-host membrane protein interactions contributes to the discovery of novel cell receptors or entry facilitators for viruses. Glycoprotein 5 (GP5), which is a major envelope protein of porcine reproductive and respiratory syndrome virus (PRRSV) virions, is a key target for the control of the virus. Here, the macrophage receptor with collagenous structure (MARCO), which is a member of the scavenger receptor family, was identified as one of the host interactors of GP5 through a DUALmembrane yeast two-hybrid screening. MARCO was specifically expressed on porcine alveolar macrophages (PAMs), and PRRSV infection downregulated MARCO expression both in vitro and in vivo. MARCO was not involved in viral adsorption and internalization processes, indicating that MARCO may not be a PRRSV-entry facilitator. Contrarily, MARCO served as a host restriction factor for PRRSV. The knockdown of MARCO in PAMs enhanced PRRSV proliferation, whereas overexpression suppressed viral proliferation. The N-terminal cytoplasmic region of MARCO was responsible for its inhibitory effect on PRRSV. Further, we found that MARCO was a proapoptotic factor in PRRSV-infected PAMs. MARCO knockdown weakened virus-induced apoptosis, whereas overexpression aggravated apoptosis. MARCO aggravated GP5-induced apoptosis, which may result in its proapoptotic function in PAMs. The interaction between MARCO and GP5 may contribute to the intensified apoptosis induced by GP5. Additionally, the inhibition of apoptosis during PRRSV infection weakened the antiviral function of MARCO, suggesting that MARCO inhibits PRRSV through the regulation of apoptosis. Taken together, the results of this study reveal a novel antiviral mechanism of MARCO and suggest a molecular basis for the potential development of therapeutics against PRRSV. IMPORTANCE Porcine reproductive and respiratory syndrome virus (PRRSV) has been one of the most serious threats to the global swine industry. Glycoprotein 5 (GP5) exposed on the surface of PRRSV virions is a major glycoprotein, and it is involved in viral entry into host cells. A macrophage receptor with collagenous structure (MARCO), which is a member of the scavenger receptor family, was identified to interact with PRRSV GP5 in a DUALmembrane yeast two-hybrid screening. Further investigation demonstrated that MARCO may not serve as a potential receptor to mediate PRRSV entry. Instead, MARCO was a host restriction factor for the virus, and the N-terminal cytoplasmic region of MARCO was responsible for its anti-PRRSV effect. Mechanistically, MARCO inhibited PRRSV infection through intensifying virus-induced apoptosis in PAMs. The interaction between MARCO and GP5 may contribute to GP5-induced apoptosis. Our work reveals a novel antiviral mechanism of MARCO and advances the development of control strategies for the virus.

Genome wide association study of thyroid hormone levels following challenge with porcine reproductive and respiratory syndrome virus

Introduction: Porcine reproductive and respiratory syndrome virus (PRRSV) causes respiratory disease in piglets and reproductive disease in sows. Piglet and fetal serum thyroid hormone (i.e., T3 and T4) levels decrease rapidly in response to Porcine reproductive and respiratory syndrome virus infection. However, the genetic control of T3 and T4 levels during infection is not completely understood. Our objective was to estimate genetic parameters and identify quantitative trait loci (QTL) for absolute T3 and/or T4 levels of piglets and fetuses challenged with Porcine reproductive and respiratory syndrome virus. Methods: Sera from 5-week-old pigs (N = 1792) at 11 days post inoculation (DPI) with Porcine reproductive and respiratory syndrome virus were assayed for T3 levels (piglet_T3). Sera from fetuses (N = 1,267) at 12 or 21 days post maternal inoculation (DPMI) with Porcine reproductive and respiratory syndrome virus of sows (N = 145) in late gestation were assayed for T3 (fetal_T3) and T4 (fetal_T4) levels. Animals were genotyped using 60 K Illumina or 650 K Affymetrix single nucleotide polymorphism (SNP) panels. Heritabilities, phenotypic correlations, and genetic correlations were estimated using ASREML; genome wide association studies were performed for each trait separately using Julia for Whole-genome Analysis Software (JWAS). Results: All three traits were low to moderately heritable (10%-16%). Phenotypic and genetic correlations of piglet_T3 levels with weight gain (0-42 DPI) were 0.26 ± 0.03 and 0.67 ± 0.14, respectively. Nine significant quantitative trait loci were identified for piglet_T3, on Sus scrofa chromosomes (SSC) 3, 4, 5, 6, 7, 14, 15, and 17, and collectively explaining 30% of the genetic variation (GV), with the largest quantitative trait loci identified on SSC5, explaining 15% of the genetic variation. Three significant quantitative trait loci were identified for fetal_T3 on SSC1 and SSC4, which collectively explained 10% of the genetic variation. Five significant quantitative trait loci were identified for fetal_T4 on SSC1, 6, 10, 13, and 15, which collectively explained 14% of the genetic variation. Several putative immune-related candidate genes were identified, including CD247, IRF8, and MAPK8. Discussion: Thyroid hormone levels following Porcine reproductive and respiratory syndrome virus infection were heritable and had positive genetic correlations with growth rate. Multiple quantitative trait loci with moderate effects were identified for T3 and T4 levels during challenge with Porcine reproductive and respiratory syndrome virus and candidate genes were identified, including several immune-related genes. These results advance our understanding of growth effects of both piglet and fetal response to Porcine reproductive and respiratory syndrome virus infection, revealing factors associated with genomic control of host resilience.

Signaling Lymphocytic Activation Molecule Family Member 1 Inhibits Porcine Reproductive and Respiratory Syndrome Virus Replication

The porcine reproductive and respiratory syndrome virus (PRRSV) causes a highly contagious disease in domestic swine. Signaling lymphocytic activation molecule family member 1 (SLAMF1) is a costimulatory factor that is involved in innate immunity, inflammation, and infection. Here, we demonstrate that overexpression of the SLAMF1 gene inhibited PRRSV replication significantly and reduced the levels of key signaling pathways, including MyD88, RIG-I, TLR2, TRIF, and inflammatory factors IL-6, IL-1β, IL-8, TNF-β, TNF-α, and IFN-α in vitro. However, the knockdown of the SLAMF1 gene could enhance replication of the PRRSV and the levels of key signaling pathways and inflammatory factors. Overall, our results identify a new, to our knowledge, antagonist of the PRRSV, as well as a novel antagonistic mechanism evolved by inhibiting innate immunity and inflammation, providing a new reference and direction for PRRSV disease resistance breeding.

Integrated Proteotranscriptomics Reveals Differences in Molecular Immunity between Min and Large White Pig Breeds

Long-term selection or evolution is an important factor governing the development of disease resistance in pigs. To better clarify the molecular mechanisms underlying different levels of disease resistance, we used transcriptomics and proteomics analysis to characterize differences in the immunities between six resistant (Min pig) and six susceptible (Large White, LW) pigs which were raised in the same environment. A total of 135 proteins and 791 genes were identified as being differentially expressed between the Large White and Min pig groups. Protein expression clustering and functional analysis revealed that proteins related to immune system process, humoral immune response, the B cell receptor signaling pathway, lymphocyte-mediated immunity, and innate immune responses were more highly expressed in Min pigs. Transcriptome gene set enrichment analysis was used to reveal that pathways of cell adhesion molecules and antigen processing and presentation are significantly enriched in Min pigs. Integrated proteomics and transcriptomics data analysis identified 16 genes that are differentially expressed at both the mRNA and protein levels. In addition, 13 out of these 16 genes were related to the quantitative trait loci of immune diseases, including neural EGFL-like 2 (NELL2) and lactate dehydrogenase B (LDHB), which are involved in innate immunity. Correlation analysis between the genes/proteins and cytokines shows upregulated proteins in LW pigs in association with immunosuppressive/pro-inflammatory cytokines, such as interleukin (IL) 10, IL6, and tumor necrosis factor alpha. This was further validated using parallel reaction monitoring analysis. In summary, we discovered several potential candidate pathways and key genes/proteins involved in determining differences in disease resistance between the two studied pig breeds, which could provide new insights into the breeding of pigs for disease resistance.

The transcriptional characteristics of NADC34-like PRRSV in porcine alveolar macrophages

The widespread and endemic circulation of porcine reproductive and respiratory syndrome virus (PRRSV) cause persistent financial losses to the swine industry worldwide. In 2017, NADC34-like PRRSV-2 emerged in northeastern China and spread rapidly. The dynamics analysis of immune perturbations associated with novel PRRSV lineage is still incomplete. This study performed a time-course transcriptome sequencing of NADC34-like PRRSV strain YC-2020-infected porcine alveolar macrophages (PAMs) and compared them with JXA1-infected PAMs. The results illustrated dramatic changes in the host’s differentially expressed genes (DEGs) presented at different timepoints after PRRSV infection, and the expression profile of YC-2020 group is distinct from that of JXA1 group. Functional enrichment analysis showed that the expression of many inflammatory cytokines was up-regulated following YC-2020 infection but at a significantly lower magnitude than JXA1 group, in line with the trends for most interferon-stimulated genes (ISGs) and their regulators. Meanwhile, numerous components of histocompatibility complex (MHC) class II and phagosome presented a stronger transcription suppression after the YC-2020 infection. All results imply that YC-2020 may induce milder inflammatory responses, weaker antiviral processes, and more severe disturbance of antigen processing and presentation compared with HP-PRRSV. Additionally, LAPTM4A, GLMP, and LITAF, which were selected from weighted gene co-expression network analysis (WGCNA), could significantly inhibit PRRSV proliferation. This study provides fundamental data for understanding the biological characteristics of NADC34-like PRRSV and new insights into PRRSV evolution and prevention.

Genetic diversity and population structure of Tongcheng pigs in China using whole-genome SNP chip

Tongcheng (TC) pigs, distinguished by their superior meat quality, are a Chinese indigenous pig breed. Recently, the genetic resources of TC pigs are under tremendous threat due to the introduction of cosmopolitan pig breeds and African swine fever disease. To promote their management and conservation, the present study assessed genetic diversity and population structure of TC pigs using single nucleotide polymorphism (SNP) markers. A total of 26, 999 SNPs were screened from 51, 315 SNPs in 68 TC pigs. The multi-dimensional scaling (MDS) analysis and neighbor-joining tree revealed that all 68 pigs were from a purebred population. The effective population size decreased over time, and it was 96 prior to generation 20. Both linkage disequilibrium (LD) and neutrality test indicated a low selection of TC pigs with average LD value of 0.15 ± 0.23. Genetic diversity results exhibited a minor allele frequency (MAF) of 0.23, observed heterozygosity (HO) of 0.32, expected heterozygosity (He) of 0.31, and nucleotide diversity (Pi) of 0.31. All these parameters indicated a remarkably high genetic diversity of TC pigs. Additionally, 184 runs of homozygosity (ROH) segments were detected from the whole genome of TC pigs with an average ROH length of 23.71Mb, ranging from 11.26Mb to 69.02 Mb. The highest ROH coverage was found on chromosome 1 (10.12%), while the lowest was on chromosome 18 (1.49%). The average inbreeding coefficients based on ROH (FROH) was 0.04%. Fourteen ROH islands containing 240 genes were detected on 9 different autosomes. Some of these 240 genes were overlapped with the genes related to biological processes such as immune function, reproduction, muscular development, and fat deposition, including FFAR2, FFAR4, MAPK8, NPY5R, KISS1, and these genes might be associated with such traits as meat quality and disease resistance in TC pigs. Taken together, population structure and genetic diversity results suggested that the TC pig represented a valuable genetic resource. However, TC pig breed conservation program remains to be further optimized to ensure adequate genetic diversity and avoid inbreeding depression. Our findings provide theoretical basis for formulating management and conservation strategies for TC pigs.

MiR-142-5p/FAM134B Axis Manipulates ER-Phagy to Control PRRSV Replication

Porcine reproductive and respiratory syndrome virus (PRRSV) can replicate its RNA genome in endoplasmic reticulum (ER) and utilize ER to facilitate its assembly and maturation. To maintain ER homeostasis, host cells initiate reticulophagy (known as ER-phagy) to effectively digest the stressed ER. In this study, we found that PRRSV infection subverted ER-phagy by downregulating ER-phagy receptor FAM134B. PRRSV-induced miR-142-5p directly targeted FAM134B and significantly promoted PRRSV replication. Meanwhile, siRNA-mediated depletion of FAM134B protein and overexpression of FAM134B mutant protein significantly disrupted ER-phagy and facilitated PRRSV replication. Furthermore, our results showed that FAM134B-mediated ER-phagy activated type I interferon signaling to inhibit PRRSV replication. Overall, this study reveals the important role of ER-phagy in PRRSV replication in a FAM134B-dependent manner. Our findings provide an insight into the pathogenesis of PRRSV and offer a theoretical basis for further development of antiviral therapeutic targets.

Integration of Transcriptome and Proteome in Lymph Nodes Reveal the Different Immune Responses to PRRSV Between PRRSV-Resistant Tongcheng Pigs and PRRSV-Susceptible Large White Pigs

Porcine reproductive and respiratory syndrome (PRRS) is an infectious disease that seriously affects the swine industry worldwide. Understanding the interaction between the host immune response and PRRS virus (PRRSV) can provide insight into the PRRSV pathogenesis, as well as potential clues to control PRRSV infection. Here, we examined the transcriptome and proteome differences of lymph nodes between PRRSV-resistant Tongcheng (TC) pigs and PRRSV-susceptible Large White (LW) pigs in response to PRRSV infection. 2245 and 1839 differentially expressed genes (DEGs) were detected in TC and LW pigs upon PRRSV infection, respectively. Transcriptome analysis revealed genetic differences in antigen presentation and metabolism between TC pigs and LW pigs, which may lead to different immune responses to PRRSV infection. Furthermore, 678 and 1000 differentially expressed proteins (DEPs) were identified in TC and LW pigs, and DEPs were mainly enriched in the metabolism pathways. Integrated analysis of transcriptome and proteome datasets revealed antigen recognition capacity, immune activation, cell cycles, and cell metabolism are important for PRRSV clearance. In conclusion, this study provides important resources on transcriptomic and proteomic levels in lymph nodes for further revealing the interaction between the host immune response and PRRSV, which would give us new insight into molecular mechanisms related to genetic complexity against PRRSV.

Transcriptome profiling identifies immune response genes against porcine reproductive and respiratory syndrome virus and Haemophilus parasuis co-infection in the lungs of piglets

Background

Co-infections of the porcine reproductive and respiratory syndrome virus (PRRSV) and the Haemophilus parasuis (HPS) are severe in Chinese pigs, but the immune response genes against co-infected with 2 pathogens in the lungs have not been reported.

Objectives

To understand the effect of PRRSV and/or HPS infection on the genes expression associated with lung immune function.

Methods

The expression of the immune-related genes was analyzed using RNA-sequencing and bioinformatics. Differentially expressed genes (DEGs) were detected and identified by quantitative real-time polymerase chain reaction (qRT-PCR), immunohistochemistry (IHC) and western blotting assays.

Results

All experimental pigs showed clinical symptoms and lung lesions. RNA-seq analysis showed that 922 DEGs in co-challenged pigs were more than in the HPS group (709 DEGs) and the PRRSV group (676 DEGs). Eleven DEGs validated by qRT-PCR were consistent with the RNA sequencing results. Eleven common Kyoto Encyclopedia of Genes and Genomes pathways related to infection and immune were found in single-infected and co-challenged pigs, including autophagy, cytokine-cytokine receptor interaction, and antigen processing and presentation, involving different DEGs. A model of immune response to infection with PRRSV and HPS was predicted among the DEGs in the co-challenged pigs. Dual oxidase 1 (DUOX1) and interleukin-21 (IL21) were detected by IHC and western blot and showed significant differences between the co-challenged pigs and the controls.

Conclusions

These findings elucidated the transcriptome changes in the lungs after PRRSV and/or HPS infections, providing ideas for further study to inhibit ROS production and promote pulmonary fibrosis caused by co-challenging with PRRSV and HPS.

Transcriptome sequencing analysis of porcine alveolar macrophages infected with PRRSV strains to elucidate virus pathogenicity and immune evasion strategies

Highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV) causes a serious disease to the swine industry worldwide. To understand the mechanisms of HP-PRRSV infection, RNA-seq-based transcriptome analyses were performed on porcine alveolar macrophages (PAMs) infected with a HP-PRRSV strain (TJ), a less virulent strain of a classical lineage (CH-1a), and a vaccine strain TJM-F92. Gene ontology, Kyoto Encyclopedia of Genes and Genomes analyses indicate that TJM-F92 led to significant up-regulation of gene expression for proteins associated with membrane-bound organelles. The differentially expressed genes of HP-PRRSV TJ-infected PAM cells were up-regulated in the special G-protein coupled receptor. The six cytokines were tested by real time Reverse Transcription-Polymerase Chain Reaction (RT-PCR). The relative expression levels showed the same trend of expression difference. Significant up-regulation of TMEM173 plays an important role in the cytosolic DNA-sensing pathway and the RIG-I-like receptor signaling pathway in TJM-F92 infected PAM cells. These data provide new insight into PRRSV pathogenicity and immune evasion strategies.

Alternative polyadenylation events differ dramatically between Tongcheng and Large White pigs in response to PRRSV infection

Alternative polyadenylation (APA) is a widespread post-transcriptional regulation mechanism that increases the biological complexity of transcriptome and proteome. However, it is unclear whether APA regulation plays a role in genetic resistance to porcine reproductive and respiratory syndrome virus (PRRSV). Here, we reported genome-wide APA regulation of porcine alveolar macrophages in PRRSV-resistant Tongcheng (TC) pigs and PRRSV-susceptible Large White (LW) pigs upon PRRSV infection. Using 3' mRNA sequencing strategy, we detected 75 981 high-quality APA sites in porcine alveolar macrophages of TC and LW pigs. Furthermore, 1202 and 1089 differentially expressed APA sites, as well as 79 and 117 untranslated region-APA switching genes were identified in TC pigs and LW pigs upon PRRSV infection respectively. The APA events in TC pigs and LW pigs were involved in different biological pathways, while APA events in TC pigs are directly associated with the immune response to PRRSV infection. In addition, we identified genetic variations affecting polyadenylation signal between TC pigs and LW pigs. These findings would provide helpful information on APA regulation for further understanding of genetic resistance to PRRSV.

Comprehensive Transcriptomic Comparison between Porcine CD8- and CD8+ Gamma Delta T Cells Revealed Distinct Immune Phenotype

Simple Summary

This study was conducted to comprehensively understand the functional mechanisms of CD8^+/− porcine gamma delta (γδ) T cells related to the immune system using RNA-sequencing technology. In total, 646 upregulated and 561 downregulated differentially expressed genes (DEGs) for CD8⁺ were identified and functional annotation was performed. A cytokine–cytokine receptor interaction and T cell receptor (TCR) signaling pathway were enriched in the upregulated DEG group, whereas the B cell receptor signaling pathway was enriched in the downregulated DEG group. Chemokine-related genes (CXCR3, CCR5, CCL4, CCL5), interferon gamma (IFNG), and CD40 ligand (CD40LG) identified in the cytokine–cytokine receptor interaction and TCR signaling pathway may affect the inter-regulation of immune signaling. Our results are expected to contribute to the understanding of mechanisms of porcine γδ T cells.

Abstract

We aimed to comprehensively understand the functional mechanisms of immunity, especially of the CD8^+/− subsets of gamma delta (γδ) T cells, using an RNA-sequencing analysis. Herein, γδ T cells were obtained from bronchial lymph node tissues of 38-day-old (after weaning 10-day: D10) and 56-day-old (after weaning 28-day: D28) weaned pigs and sorted into CD8⁺ and CD8⁻ groups. Differentially expressed genes (DEGs) were identified based on the CD8 groups at D10 and D28 time points. We confirmed 1699 DEGs between D10 CD8⁺ versus D10 CD8⁻ groups and 1784 DEGs between D28 CD8⁺ versus D28 CD8⁻ groups; 646 upregulated and 561 downregulated DEGs were common. The common upregulated DEGs were enriched in the cytokine–cytokine receptor interaction and T cell receptor (TCR) signaling pathway, and the common downregulated DEGs were enriched in the B cell receptor signaling pathway. Further, chemokine-related genes, interferon gamma, and CD40 ligand were involved in the cytokine–cytokine receptor interaction and TCR signaling pathway, which are associated with inter-regulation in immunity. We expect our results to form the basic data required for understanding the mechanisms of γδ T cells in pigs; however, further studies are required in order to reveal the dynamic changes in γδ T cells under pathogenic infections, such as those by viruses.

Ovary-derived circular RNAs profile analysis during the onset of puberty in gilts

Background

In mammals, the ovary is the essential system of female reproduction for the onset of puberty, and the abnormal puberty has negative outcomes on health. CircRNA is a non-coding RNA produced by non-canonical alternative splicing (AS). Several studies have reported that circRNA is involved in the gene regulation and plays an important role in some human diseases. However, the contribution of circRNA has received little known within the onset of puberty in ovary.

Results

Here, the profiles of ovarian circRNAs across pre-, in- and post-pubertal stages were established by RNA-sEq. In total, 972 circRNAs were identified, including 631 stage-specific circRNAs and 8 tissue-specific circRNAs. The biological functions of parental genes of circRNAs were enriched in steroid biosynthesis, autophagy-animal, MAPK signaling pathway, progesterone-mediated oocyte maturation and ras signaling pathway. Moreover, 5 circRNAs derived from 4 puberty-related genes (ESR1, JAK2, NF1 and ARNT) were found in this study. The A3SS events were the most alternative splicing, but IR events were likely to be arose in post-pubertal ovaries. Besides, the circRNA-miRNA-gene networks were explored for 10 differentially expressed circRNAs. Furthermore, the head-to-tail exon as well as the expressions of 10 circRNAs were validated by the divergent RT-qPCR and sanger sequencing.

Conclusions

In summary, the profiles of ovarian circRNAs were provided during pubertal transition in gilts, and these results provided useful information for the investigation on the onset of puberty at the ovarian-circRNAs-level in mammals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07786-w.

Structure and function of the porcine TAP protein and its inhibition by the viral immune evasion protein ICP47

The binding mode to TAP (i.e., the peptide transporter associated with antigen processing) from a viral peptide thus far has been unknown in the field of antiviral immunity, but an interfering mode from a virus-encoded TAP inhibitor has been well documented with respect to blocking the TAP function. In the current study, we predicted the structure of the pig TAP transporter and its inhibition complex by the small viral protein ICP47 of the herpes simplex virus (HSV) encoded by the TAP inhibitor to exploit inhibition of the TAP transporter as the host's immune evasion strategy. We found that the hot spots (residues Leu5, Tyr22, and Leu51) on the ICP47 inhibitor interface tended to prevail over the favored Leu and Tyr, which contributed to significant functional binding at the C-termini recognition principle of the TAP. We further characterized the specificity determinants of the peptide transporter from the pig TAP by the ICP47 inhibitor effects and multidrug TmrAB transporter from the Thermus thermophillus and its immunity regarding its structural homolog of the pig TAP. The specialized structure-function relationship from the pig TAP exporter could provide insight into substrate specificity of the unique immunological properties from the host organism. The TAP disarming capacity from all five viral inhibitors (i.e., the five virus-encoded TAP inhibitors of ICP47, UL49.5, U6, BNLF2a, and CPXV012 proteins) was linked to the infiltration of the TAP functional structure in an unstable conformation and the mounting susceptibility caused by the host's TAP polymorphism. It is anticipated that the functional characterization of the pig TAP transporter based on the pig genomic variants will lead to additional insights into the genotype and single nucleotide polymorphism (SNP) in relation to antiviral resistance and disease susceptibility.

Role of IFNLR1 gene in PRRSV infection of PAM cells

Background

Interferon lambda receptor 1 (IFNLR1) is a type II cytokine receptor that clings to interleukins IL-28A, IL29B, and IL-29 referred to as type III IFNs (IFN-λs). IFN-λs act through the JAK-STAT signaling pathway to exert antiviral effects related to preventing and curing an infection. Although the immune function of IFN-λs in virus invasion has been described, the molecular mechanism of IFNLR1 in that process is unclear.

Objectives

The purpose of this study was to elucidate the role of IFNLR1 in the pathogenesis and treatment of porcine reproductive and respiratory syndrome virus (PRRSV).

Methods

The effects of IFNLR1 on the proliferation of porcine alveolar macrophages (PAMs) during PRRSV infection were investigated using interference and overexpression methods.

Results

In this study, the expressions of the IFNLR1 gene in the liver, large intestine, small intestine, kidney, and lung tissues of Dapulian pigs were significantly higher than those in Landrace pigs. It was determined that porcine IFNLR1 overexpression suppresses PRRSV replication. The qRT-PCR results revealed that overexpression of IFNLR1 upregulated antiviral and IFN-stimulated genes. IFNLR1 overexpression inhibits the proliferation of PAMs and upregulation of p-STAT1. By contrast, knockdown of IFNLR1 expression promotes PAMs proliferation. The G0/G1 phase proportion in IFNLR1-overexpressing cells increased, and the opposite change was observed in IFNLR1-underexpressing cells. After inhibition of the JAK/STAT signaling pathway, the G2/M phase proportion in the IFNLR1-overexpressing cells showed a significant increasing trend. In conclusion, overexpression of IFNLR1 induces activation of the JAK/STAT pathway, thereby inhibiting the proliferation of PAMs infected with PRRSV.

Conclusion

Expression of the IFNLR1 gene has an important regulatory role in PRRSV-infected PAMs, indicating it has potential as a molecular target in developing a new strategy for the treatment of PRRSV.

Whole-genome sequencing reveals breed-differential CNVs between Tongcheng and Large White pigs

Large phenotypic differences have been observed between Tongcheng and Large White pigs. However, little is known about their genetic basis. This study performed a genome-wide comparison of CNVs between Tongcheng and Large White pigs using genome sequencing data. By combining the advantages of three different strategies (read depth, paired-end mapping and split read), we detected in total 18 687 CNVs that covered approximately 3.5% of the pig genome length for Tongcheng and Large White pigs. We identified 1864 breed-stratified CNVs (top 10%) by performing V_ST statistics. Functional enrichment analyses for genes located in breed-stratified CNVs were found to be involved in pigmentation, behavior, immune system and reproductive processes, which coincide with phenotypic differences between the two breeds. Using a systematic analysis of the genome and transcriptome data, we further identified four novel breed-differential CNVs on the functional genes (disease-resistant, DCUN1D2 and SPARCL1; lipid metabolism, PLEKHA2 and SLCO1A2). Subsequent PCR validation confirmed their accurate breakpoint positions in 33 Tongcheng pigs and 33 Large White pigs. This study provides essential information on differential CNVs for further research on the genetic basis of phenotypic differences between Tongcheng and Large White pigs.

Improvement of the Enterotoxigenic Escherichia coli Infection Model for Post-Weaning Diarrhea by Controlling for Bacterial Adhesion, Pig Breed and MUC4 Genotype

Enterotoxigenic Escherichia coli (ETEC) is a major cause of post-weaning diarrhea (PWD) in pigs and causes significant damage to the swine industry worldwide. In recent years, there has been increased regulation against the use of antibacterial agents in swine due to their health risks. Utilizing experimental models that consistently recapitulate PWD is important for the development of non-antibacterial agents against PWD in pigs. In this study, we established a highly reproducible PWD infection model by examining differences in adhesion of ETEC to the intestinal tissue as well as the association between MUC4 polymorphisms and sensitivity to PWD. Post-weaning diarrhea differences between pig breeds were also examined. The adhesion to enterocytes varied from 10^4.0 to 10^6.4 CFU/mL even among the F4 ETEC strains. Experimental infection revealed that PWD can be induced in all MUC4 genotypes after infection with 10¹⁰ CFU/pig of highly adherent ETEC, although there were variable sensitivities between the genotypes. Lowly adherent ETEC did not cause PWD as efficiently as did highly adherent ETEC. The incidence of PWD was confirmed for all pigs with the ETEC-susceptible MUC4 genotypes in all of the breeds. These results indicate that high-precision and reproducible experimental infection is possible regardless of pig breeds by controlling factors on the pig-end (MUC4 genotype) and the bacterial-end (adhesion ability).

Two coupled mutations abolished the binding of CEBPB to the promoter of CXCL14 that displayed an antiviral effect on PRRSV by activating IFN signaling

Porcine reproductive and respiratory syndrome (PRRS) is the most economically important infectious disease of pigs worldwide. Our previous study revealed that Tongcheng (TC) pigs display higher resistance to PRRS than Largewhite (LW) pigs, but the genetic mechanism remains unknown. Here, we first confirmed that CXCL14 was downregulated in lungs and porcine alveolar macrophages (PAMs) responding to PRRS virus (PRRSV) infection, but the decline in LW pigs was more obvious than that in TC pigs. Then, we found that the overexpression of CXCL14 activated type-I interferon (IFN-I) signaling by upregulating interferon beta (IFNB), which plays a major role in the antiviral effect. To further decipher the mechanism underlying its differential expression, we characterized the core promoter of CXCL14 as being located from -145 to 276 bp of the transcription start site (TSS) and identified two main haplotypes that displayed significant differential transcriptional activities. We further identified two coupled point mutations that altered the binding status of CEBPB and were responsible for the differential expression in TC and LW pigs. The regulatory effect of CEBPB on CXCL14 was further confirmed by RNA interference (RNAi) and chromatin immunoprecipitation (ChIP), providing crucial clues for deciphering the mechanism of CXCL14 downregulation in unusual conditions. The present study revealed the potential antiviral effect of CXCL14, occurring via activation of interferon signaling, and suggested that CXCL14 contributes to the PRRS resistance of TC pigs.

Adaptive changes in the genomes of wild rabbits after 16 years of viral epidemics

Since its introduction to control overabundant invasive European rabbits (Oryctolagus cuniculus), the highly virulent rabbit haemorrhagic disease virus (RHDV) has caused regular annual disease outbreaks in Australian rabbit populations. Although initially reducing rabbit abundance by 60%, continent-wide, experimental evidence has since indicated increased genetic resistance in wild rabbits that have experienced RHDV-driven selection. To identify genetic adaptations, which explain the increased resistance to this biocontrol virus, we investigated genome-wide SNP (single nucleotide polymorphism) allele frequency changes in a South Australian rabbit population that was sampled in 1996 (pre-RHD genomes) and after 16 years of RHDV outbreaks. We identified several SNPs with changed allele frequencies within or close to genes potentially important for increased RHD resistance. The identified genes are known to be involved in virus infections and immune reactions or had previously been identified as being differentially expressed in healthy versus acutely RHDV-infected rabbits. Furthermore, we show in a simulation study that the allele/genotype frequency changes cannot be explained by drift alone and that several candidate genes had also been identified as being associated with surviving RHD in a different Australian rabbit population. Our unique data set allowed us to identify candidate genes for RHDV resistance that have evolved under natural conditions, and over a time span that would not have been feasible in an experimental setting. Moreover, it provides a rare example of host genetic adaptations to virus-driven selection in response to a suddenly emerging infectious disease.

Global Analysis of Alternative Splicing Difference in Peripheral Immune Organs between Tongcheng Pigs and Large White Pigs Artificially Infected with PRRSV In Vivo

Alternative splicing (AS) plays a significant role in regulating gene expression at the transcriptional level in eukaryotes. Flexibility and diversity of transcriptome and proteome can be significantly increased through alternative splicing of genes. In the present study, transcriptome data of peripheral immune organs including spleen and inguinal lymph nodes (ILN) were used to identify AS difference between PRRSV-resistant Tongcheng (TC) pigs and PRRSV-susceptible Large White (LW) pigs artificially infected with porcine reproductive and respiratory syndrome virus (PRRSV) in vivo. The results showed that PRRSV infection induced global alternative splicing events (ASEs) with different modes. Among them, 373 genes and 595 genes in the spleen and ILN of TC pigs, while 458 genes and 560 genes in the spleen and ILN of LW pigs had significantly differential ASEs. Alternative splicing was subject to tissue-specific and lineage-specific regulation in response to PRRSV infection. Enriched GO terms and pathways showed that genes with differential ASEs played important roles in transcriptional regulation, immune response, metabolism, and apoptosis. Furthermore, a splicing factor associated with apoptosis, SRSF4, was significantly upregulated in LW pigs. Functional analysis on apoptosis associated genes was validated by RT-PCR and DNA sequencing. These findings revealed different response to PRRSV between PRRSV-resistant TC pigs and PRRSV-susceptible LW pigs at the level of alternative splicing, suggesting the potential relationship between AS and disease resistance to PRRSV.

Cellular Innate Immunity against PRRSV and Swine Influenza Viruses

Porcine respiratory disease complex (PRDC) is a polymicrobial syndrome that results from a combination of infectious agents, such as environmental stressors, population size, management strategies, age, and genetics. PRDC results in reduced performance as well as increased mortality rates and production costs in the pig industry worldwide. This review focuses on the interactions of two enveloped RNA viruses—porcine reproductive and respiratory syndrome virus (PRRSV) and swine influenza virus (SwIV)—as major etiological agents that contribute to PRDC within the porcine cellular innate immunity during infection. The innate immune system of the porcine lung includes alveolar and parenchymal/interstitial macrophages, neutrophils (PMN), conventional dendritic cells (DC) and plasmacytoid DC, natural killer cells, and γδ T cells, thus the in vitro and in vivo interactions between those cells and PRRSV and SwIV are reviewed. Likewise, the few studies regarding PRRSV-SwIV co-infection are illustrated together with the different modulation mechanisms that are induced by the two viruses. Alterations in responses by natural killer (NK), PMN, or γδ T cells have not received much attention within the scientific community as their counterpart antigen-presenting cells and there are numerous gaps in the knowledge regarding the role of those cells in both infections. This review will help in paving the way for future directions in PRRSV and SwIV research and enhancing the understanding of the innate mechanisms that are involved during infection with these viruses.

Expression of Immune Regulatory Genes in the Porcine Internal Genital Tract Is Differentially Triggered by Spermatozoa and Seminal Plasma

Mating or cervical deposition of spermatozoa or seminal plasma (SP) modifies the expression of genes affecting local immune defense processes at the oviductal sperm reservoir in animals with internal fertilization, frequently by down-regulation. Such responses may occur alongside sperm transport to or even beyond the reservoir. Here, immune-related gene expression was explored with cDNA microarrays on porcine cervix-to-infundibulum tissues, pre-/peri-ovulation. Samples were collected 24 h post-mating or cervical deposition of sperm-peak spermatozoa or SP (from the sperm-peak fraction or the whole ejaculate). All treatments of this interventional study affected gene expression. The concerted action of spermatozoa and SP down-regulated chemokine and cytokine (P00031), interferon-gamma signaling (P00035), and JAK/STAT (P00038) pathways in segments up to the sperm reservoir (utero-tubal junction (UTJ)/isthmus). Spermatozoa in the vanguard sperm-peak fraction (P1-AI), uniquely displayed an up-regulatory effect on these pathways in the ampulla and infundibulum. Sperm-free SP, on the other hand, did not lead to major effects on gene expression, despite the clinical notion that SP mitigates reactivity by the female immune system after mating or artificial insemination.

Identification of Differentially Expressed Non-coding RNA in Porcine Alveolar Macrophages from Tongcheng and Large White Pigs Responded to PRRSV

Porcine reproductive and respiratory syndrome (PRRS) is one of the most ruinous diseases in pig production. Our previous work showed that Tongcheng pigs (TC) were less susceptible to PRRS virus (PRRSV) than Large White (LW) pigs. To elucidate the difference in PRRSV resistance between the two breeds, small RNA-seq and ribo-zero RNA-seq were used to identify differentially expressed non-coding RNAs (including miRNAs and lincRNAs) responded to PRRSV in porcine alveolar macrophages (PAMs) from TC and LW pigs. Totally, 250 known mature miRNAs were detected. For LW pigs, there were 44 down-regulated and 67 up-regulated miRNAs in infection group; while for TC pigs, 12 down-regulated and 23 up-regulated miRNAs in TC infection group were identified. The target genes of the common differentially expressed miRNAs (DEmiRNAs) in these two breeds were enriched in immune-related processes, including apoptosis process, inflammatory response, T cell receptor signaling pathway and so on. In addition, 5 shared DEmiRNAs (miR-181, miR-1343, miR-296-3p, miR-199a-3p and miR-34c) were predicted to target PRRSV receptors, of which miR-199a-3p was validated to inhibit the expression of CD151. Interestingly, miR-378 and miR-10a-5p, which could inhibit PRRSV replication, displayed higher expression level in TC control group than that in LW control group. Contrarily, miR-145-5p and miR-328, which were specifically down-regulated in LW pigs, could target inhibitory immunoreceptors and may involve in immunosuppression caused by PRRSV. This indicates that DEmiRNAs are involved in the regulation of the immunosuppression and immune escape of the two breeds. Furthermore, we identified 616 lincRNA transcripts, of which 48 and 30 lincRNAs were differentially expressed in LW and TC pigs, respectively. LincRNA TCONS_00125566 may play an important role in the entire regulatory network, and was predicted to regulate the expression of immune-related genes through binding with miR-1343 competitively. In conclusion, this study provides an important resource for further revealing the interaction between host and virus, which will specify a new direction for anti-PRRSV research.

Transcriptome Analysis Reveals Dynamic Gene Expression Profiles in Porcine Alveolar Macrophages in Response to the Chinese Highly Pathogenic Porcine Reproductive and Respiratory Syndrome Virus

Porcine reproductive and respiratory syndrome virus (PRRSV) is one of the most economically important swine pathogens and causes reproductive failure in sows and respiratory disease in growing pigs. PRRSV mainly infects porcine alveolar macrophages (PAMs), leading to the subversion of innate and adaptive immunity of pigs. The transcriptome analysis of gene expression profiles in PRRSV-infected PAMs is essential for understanding the pathogenesis of PRRSV. Here we performed next-generation RNA sequencing and a comprehensive bioinformatics analysis to characterize the dynamic transcriptome landscapes in PAMs following PRRSV infection. Totally 38222 annotated mRNAs, 12987 annotated long noncoding RNAs (lncRNAs), and 17624 novel lncRNAs in PRRSV-infected PAMs were identified through a transcripts computational identification pipeline. The differentially expressed mRNAs and lncRNAs during PRRSV infection were characterized. Several differentially expressed transcripts were validated using qRT-PCR. Analyses on dynamic overrepresented GO terms and KEGG pathways in PRRSV-infected PAMs at different time points were performed. Meanwhile the genes involved in IFN-related signaling pathways, proinflammatory cytokines and chemokines, phagocytosis, and antigen presentation and processing were significantly downregulated, indicating the aberrant function of PAMs during PRRSV infection. Moreover, the differentially and highly expressed lncRNA XR_297549.1 was predicted to both cis-regulate and trans-regulate its neighboring gene, prostaglandin-endoperoxide synthase 2 (PTGS2), indicating its role in inflammatory response. Our findings reveal the transcriptome profiles and differentially expressed mRNAs and lncRNAs in PRRSV-infected PAMs in vitro, providing valuable information for further exploration of PRRSV pathogenesis.