Integrated time-serial transcriptome networks reveal common innate and tissue-specific adaptive immune responses to PRRSV infection

Cell deconvolution-based integrated time-series network of whole blood transcriptome reveals systemic antiviral activities and cell-specific immunological changes against PRRSV infection

Porcine reproductive and respiratory syndrome (PRRS) causes significant economic losses in the swine industry. However, the molecular mechanisms behind the common and cell type-specific systemic responses during PRRS virus (PRRSV) infection are not well understood. In this study, we collected viremia data, antibody levels, and whole-blood RNA-seq data obtained from eight PRRSV-infected piglets. We utilised a cell deconvolution approach to calculate cell type enrichment, constructed a time-serial gene co-expression network with differentially expressed genes, and conducted functional annotations. Three significant modules were identified within the network. The changes associated with viremia revealed an upregulated expression of genes related to antiviral activity. In the T-cell- and NK-cell-specific modules, infection led to an increased T-cell population and upregulation of genes related to T-cell defence responses. Conversely, in the monocyte- and neutrophil-specific module, genes involved in inflammatory responses were downregulated due to a decrease in monocyte proportion. This study highlights the time-series antiviral activities associated with viremia and the transcriptomic changes associated with immune responses in specific cell types. The findings provide comprehensive insights into host responses to PRRSV infection, including diagnostic biomarkers.

Comparison of pathogenicity and host responses of emerging porcine reproductive and respiratory syndrome virus variants in piglets

Porcine reproductive and respiratory syndrome virus (PRRSV) is a highly variable virus with genetic diversity. This study comparatively examines the pathogenicity and immunological impact of two emergent PRRSV strains, SD53 and HuN4, in piglets. Our results indicate that SD53 strain induces milder clinical syndromes and less severe tissue damage than HuN4, despite similar replication rates. Hematological tests showed less perturbations in peripheral blood cell profiles after SD53 infection, suggesting a less systemic impact. The neutrophil-to-lymphocyte ratio was notably lower in SD53-infected piglets, suggesting a less intense inflammatory reaction. Moreover, SD53 infection led to lower levels of pro-inflammatory cytokines, further supporting a less pronounced inflammatory profile. Both strains induced the production of PRRSV-specific antibodies. However, transcriptomic analysis of lung and lymph node tissues from infected piglets disclosed a more moderate up-regulation of core genes, including ISGs, in the SD53 group. Further analysis indicated that SD53 primarily enhanced immune-related signaling, particularly in T cell response modules, while HuN4 caused a more robust pro-inflammatory reaction and a dampening of T cell functionality. Flow cytometry analyses confirmed these findings, showing higher CD4/CD8 ratios and increased CD4+ T cell percentages in SD53-infected piglets, implying a more robust T cell response. Collectively, these findings broaden our comprehension of PRRSV pathogenesis and may inform the development of future therapeutic or prophylactic strategies for controlling PRRSV infections more effectively.

IMPORTANCE

The high mutation rate of porcine reproductive and respiratory syndrome virus (PRRSV) poses significant challenges to its accurate diagnosis and the implementation of effective control measures. This research explores the pathogenic profiles of two emerging PRRSV stains: the NADC30-like strain SD53 and the highly pathogenic strain HuN4. Our investigation reveals that SD53 initiates distinct immunopathological responses in vivo compared with those provoked by HuN4. By conducting a transcriptome analysis of differential gene expression in the lungs and lymph nodes of infected piglets, we unveil the intricate molecular mechanisms underlying the contrasting pathogenicity of these two strains. The comprehensive insights yielded by this study are instrumental in advancing our understanding of the dominant NADC30-like PRRSV strain, which has become increasingly prevalent in China's swine industry.

Veterinary systems biology for bridging the phenotype-genotype gap via computational modeling for disease epidemiology and animal welfare

Veterinary systems biology is an innovative approach that integrates biological data at the molecular and cellular levels, allowing for a more extensive understanding of the interactions and functions of complex biological systems in livestock and veterinary science. It has tremendous potential to integrate multi-omics data with the support of vetinformatics resources for bridging the phenotype-genotype gap via computational modeling. To understand the dynamic behaviors of complex systems, computational models are frequently used. It facilitates a comprehensive understanding of how a host system defends itself against a pathogen attack or operates when the pathogen compromises the host's immune system. In this context, various approaches, such as systems immunology, network pharmacology, vaccinology and immunoinformatics, can be employed to effectively investigate vaccines and drugs. By utilizing this approach, we can ensure the health of livestock. This is beneficial not only for animal welfare but also for human health and environmental well-being. Therefore, the current review offers a detailed summary of systems biology advancements utilized in veterinary sciences, demonstrating the potential of the holistic approach in disease epidemiology, animal welfare and productivity.

Effect of dietary betaine supplementation on the liver transcriptome profile in broiler chickens under heat stress conditions

OBJECTIVE

The objective of the present study was to investigate the effect of dietary betaine (BT) supplementation on the hepatic transcriptome profiles in broiler chickens raised under heat stress (HS) conditions.

METHODS

A total of 180 (21-d-old) Ross 308 male broiler chicks were allotted to 1 of 3 treatment groups with 6 replicated cages in a completely randomized design. One group was kept under thermoneutral conditions at all times and was fed a basal diet (PC). Other 2 groups were exposed to a cyclic heat stress condition. One of the 2 groups under heat stress conditions was fed the basal diet as a negative control (NC), whereas the other group was fed the basal diet supplemented with 0.2% BT. All chickens were provided with diets and water ad libitum for 21 d. Following the experiment, the liver samples were collected for RNA sequencing analysis.

RESULTS

Broiler chickens in NC and BT group had decreased (p<0.05) growth performance. In the transcriptome analysis, the number of differentially expressed genes were identified in the liver by HS conditions and dietary BT supplementation. In the comparison between NC and PC treatments, genes related to energy and nucleic acid metabolism, amino acid metabolism, and immune system were altered by HS, which support the reason why heat-stressed poultry had decreased growth performance. In the comparison between NC and BT treatments, genes related to lipid metabolism, carbohydrate metabolism, and immune system were differently expressed under HS conditions.

CONCLUSION

HS negatively impacts various physiological processes, including DNA replication, metabolism of amino acids, lipids, and carbohydrates, and cell cycle progression in broiler chickens. Dietary BT supplementation, however, offers potential counteractive effects by modulating liver function, facilitating gluconeogenesis, and enhancing immune systems. These findings provide a basis for understanding molecular responses by HS and the possible benefits of dietary BT supplementation in broiler chickens exposed to HS.

Comparative transcriptome analysis between long- and short-term survival after pig-to-monkey cardiac xenotransplantation reveals differential heart failure development

Cardiac xenotransplantation is the potential treatment for end-stage heart failure, but the allogenic organ supply needs to catch up to clinical demand. Therefore, genetically-modified porcine heart xenotransplantation could be a potential alternative. So far, pig-to-monkey heart xenografts have been studied using multi-transgenic pigs, indicating various survival periods. However, functional mechanisms based on survival period-related gene expression are unclear. This study aimed to identify the differential mechanisms between pig-to-monkey post-xenotransplantation long- and short-term survivals. Heterotopic abdominal transplantation was performed using a donor CD46-expressing GTKO pig and a recipient cynomolgus monkey. RNA-seq was performed using samples from POD60 XH from monkey and NH from age-matched pigs, D35 and D95. Gene-annotated DEGs for POD60 XH were compared with those for POD9 XH (Park et al. 2021). DEGs were identified by comparing gene expression levels in POD60 XH versus either D35 or D95 NH. 1,804 and 1,655 DEGs were identified in POD60 XH versus D35 NH and POD60 XH versus D95 NH, respectively. Overlapped 1,148 DEGs were annotated and compared with 1,348 DEGs for POD9 XH. Transcriptomic features for heart failure and inhibition of T cell activation were observed in both long (POD60)- and short (POD9)-term survived monkeys. Only short-term survived monkey showed heart remodeling and regeneration features, while long-term survived monkey indicated multi-organ failure by neural and hormonal signaling as well as suppression of B cell activation. Our results reveal differential heart failure development and survival at the transcriptome level and suggest candidate genes for specific signals to control adverse cardiac xenotransplantation effects.

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.

The Genetic Variation of Porcine Reproductive and Respiratory Syndrome Virus Replicase Protein nsp2 Modulates Viral Virulence and Persistence

Fast evolution in the field of the replicase nsp2 represents a most prominent feature of porcine reproductive and respiratory syndrome virus (PRRSV). Here, we determined its biological significance in viral pathogenesis by constructing interlineage chimeric mutants between the Chinese highly pathogenic PRRSV (HP-PRRSV) strain JXwn06 (lineage 8) and the low-virulent NADC30-like strain CHsx1401 (lineage 1). Replacement with nsp2 from JXwn06 was surprisingly lethal to the backbone virus CHsx1401, but combined substitution with the structural protein-coding region (SP) gave rise to viable virus CHsx1401-SPnsp2JX. Meanwhile, a derivative carrying only the SP region (CHsx1401-SPJX) served as a control. Subsequent animal experiments revealed that acquisition of SP alone (CHsx1401-SPJX) did not allow CHsx1401 to gain much virulence, but additional swapping of HP-PRRSV nsp2 (CHsx1401-SPnsp2JX) enabled CHsx1401 to acquire some properties of HP-PRRSV, exemplified by prolonged high fever, microscopic lung hemorrhage, and a significant increase in proinflammatory cytokines in the acute stage. Consistent with this was the transcriptomic analysis of persistently infected secondary lymphoid tissues that revealed a much stronger induction of host cellular immune responses in this group and identified several core immune genes (e.g., TLR4, IL-1β, MPO, etc.) regulated by HP-PRRSV nsp2. Interestingly, immune activation status in the individual groups correlated well with the rate of viremia clearance and viral tissue load reduction. Overall, the above results suggest that the Chinese HP-PRRSV nsp2 is a critical virulence regulator and highlight the importance of nsp2 genetic variation in modulating PRRSV virulence and persistence via immune modulation. IMPORTANCE Porcine reproductive and respiratory syndrome virus (PRRSV) has been a major threat to the world swine industry. In the field, rapid genetic variations (e.g., deletion, mutation, recombination, etc.) within the nsp2 region present an intriguing conundrum to PRRSV biology and pathogenesis. By making chimeric mutants, here, we show that the Chinese highly pathogenic PRRSV (HP-PRRSV) nsp2 is a virulence factor and a much stronger inducer of host immune responses (e.g., inflammation) than its counterpart, currently epidemic, NADC30-like strains. Differences in the ability to modulate host immunity provide insight into the mechanisms of why NADC30-like strains and their derivatives are rising to be the dominant viruses, whereas the Chinese HP-PRRSV strains gradually give away center stage in the field. Our results have important implications in understanding PRRSV evolution, interlineage recombination, and persistence.

Integrated transcriptome analysis for the hepatic and jejunal mucosa tissues of broiler chickens raised under heat stress conditions

Background

Heat stress (HS) is one of the most important threats for the current poultry industry. Therefore, many efforts have been made to ameliorate the adverse effect of HS on poultry production; however, physiological and molecular mechanisms pertaining to HS are still limited in poultry. Therefore, the objective of the current study was to investigate functional alterations based on individual and integrated transcriptomes in the liver and jejunal mucosa tissues of broiler chickens exposed to HS conditions.

Results

Broiler chickens exposed to HS showed decreased growth performance and increased corticosterone concentrations in the feather. In the transcriptome analysis, the number of differentially expressed genes (DEGs) were identified in the liver and jejunal mucosa by HS conditions. In the liver, genes related to amino acid oxidation, tryptophan metabolism, lipid metabolism, oxidative phosphorylation, and immune responses were altered by HS, which support the reason why heat-stressed poultry had decreased productive performance. In the jejunal mucosa, genes related to defense systems, glutathione metabolism, detoxification of xenobiotics, and immune responses were differently expressed by HS conditions. The integrated transcriptome analysis with DEGs found in the liver and jejunal mucosa showed a considerable connectivity between the core nodes in the constructed networks, which includes glutathione metabolism, xenobiotic metabolism, carbon metabolism, and several amino acid metabolisms.

Conclusions

The core network analysis may indicate that increased requirement of energy and amino acids in the jejunal mucosa of broiler chickens exposed to HS conditions is likely compromised by increased oxidation and synthesis of amino acids in the liver. Therefore, our results may provide comprehensive insights for molecular and metabolic alterations of broiler chickens raised under HS conditions, which can aid in the development of the novel strategies to ameliorate the negative effect of HS on poultry productivity and health.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00734-y.

Integrated time-series transcriptomic and metabolomic analyses reveal different inflammatory and adaptive immune responses contributing to host resistance to PRRSV

Porcine reproductive and respiratory syndrome virus (PRRSV) is a highly contagious disease that affects the global pig industry. To understand mechanisms of susceptibility/resistance to PRRSV, this study profiled the time-serial white blood cells transcriptomic and serum metabolomic responses to PRRSV in piglets from a crossbred population of PRRSV-resistant Tongcheng pigs and PRRSV-susceptible Large White pigs. Gene set enrichment analysis (GSEA) illustrated that PRRSV infection up-regulated the expression levels of marker genes of dendritic cells, monocytes and neutrophils and inflammatory response, but down-regulated T cells, B cells and NK cells markers. CIBERSORT analysis confirmed the higher T cells proportion in resistant pigs during PRRSV infection. Resistant pigs showed a significantly higher level of T cell activation and lower expression levels of monocyte surface signatures post infection than susceptible pigs, corresponding to more severe suppression of T cell immunity and inflammatory response in susceptible pigs. Differentially expressed genes between resistant/susceptible pigs during the course of infection were significantly enriched in oxidative stress, innate immunity and humoral immunity, cell cycle, biotic stimulated cellular response, wounding response and behavior related pathways. Fourteen of these genes were distributed in 5 different QTL regions associated with PRRSV-related traits. Chemokine CXCL10 levels post PRRSV infection were differentially expressed between resistant pigs and susceptible pigs and can be a promising marker for susceptibility/resistance to PRRSV. Furthermore, the metabolomics dataset indicated differences in amino acid pathways and lipid metabolism between pre-infection/post-infection and resistant/susceptible pigs. The majority of metabolites levels were also down-regulated after PRRSV infection and were significantly positively correlated to the expression levels of marker genes in adaptive immune response. The integration of transcriptome and metabolome revealed concerted molecular events triggered by the infection, notably involving inflammatory response, adaptive immunity and G protein-coupled receptor downstream signaling. This study has increased our knowledge of the immune response differences induced by PRRSV infection and susceptibility differences at the transcriptomic and metabolomic levels, providing the basis for the PRRSV resistance mechanism and effective PRRS control.

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.

Hepatic transcriptome profiling according to growth rate reveals acclimation in metabolic regulatory mechanisms to cyclic heat stress in broiler chickens

Climate change has numerous effects on poultry that result in welfare concerns and economic losses in agricultural industries. However, the mechanisms underlying the acclimation to heat stress in poultry have not been comprehensively defined. Therefore, identifying associated patterns of gene regulation and understanding the molecular mechanisms of acclimation to a warmer environment will provide insights into the acclimation system of broiler chickens. We profiled differentially expressed genes (DEGs) associated with differences in growth performance under heat stress conditions in the liver tissues of broilers based on RNA sequencing data. The DEGs were identified by comparison to the gene expression levels of broilers exhibiting average growth at 28 d of age (D28A) and D36A relative to those at D21A. In D36A, 507 and 312 DEGs were up- and downregulated, respectively, whereas 400 and 156 DEGs were up- and downregulated in D28A, respectively. Pathway enrichment analysis further revealed that “fatty acid degradation” and “heat shock protein expression” were upregulated in broilers exhibiting a higher growth and weight, whereas “cell cycle arrest” and “amino acid metabolism” were downregulated. Transcriptome profiling revealed that the acclimatized group supplied fat and energy from the liver to tissues through the breakdown of fatty acids. Furthermore, homeostasis was maintained via heat shock proteins and antioxidant enzymes. The characterized candidate genes and mechanisms associated with the response to heat stress might serve as a foundation for improving the ability of broilers to acclimatize under heat stress conditions.

Molecular cloning, characterization, and functional analysis of the uncharacterized C11orf96 gene

Background

The mammalian genome encodes millions of proteins. Although many proteins have been discovered and identified, a large part of proteins encoded by genes are yet to be discovered or fully characterized. In the present study, we successfully identified a host protein C11orf96 that was significantly upregulated after viral infection.

Results

First, we successfully cloned the coding sequence (CDS) region of the cat, human, and mouse C11orf96 gene. The CDS region of the C11orf96 gene is 372 bp long, encodes 124 amino acids, and is relatively conserved in different mammals. From bioinformatics analysis, we found that C11orf96 is rich in Ser and has multiple predicted phosphorylation sites. Moreover, protein interaction prediction analysis revealed that the protein is associated with several transmembrane family proteins and zinc finger proteins. Subsequently, we found that C11orf96 is strictly distributed in the cytoplasm. According to the tissue distribution characteristics, C11orf96 is distributed in all tissues and organs, with the highest expression levels in the kidney. These results indicate that C11orf96 may play a specific biological role in the kidney.

Conclusions

Summarizing, these data lay the foundation for studying the biological functions of C11orf96 and for exploring its role in viral replication.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-022-03224-5.

The Antiviral Effect of Isatis Root Polysaccharide against NADC30-like PRRSV by Transcriptome and Proteome Analysis

(1) Background: In recent years, the porcine reproductive and respiratory syndrome virus (PRRSV) has become a virulent pathogen that has caused devastating diseases and economic losses worldwide in the swine industry. IRPS has attracted extensive attention in the field of virology. However, it is not clear that IRPS has an antiviral effect on PRRSV at gene and protein levels. (2) Methods: We used transcriptomic and proteomic analysis to investigate the antiviral effect of IRPS against PRRSV. Additionally, a microbiome was used to explore the effects of IRPS on gut microbes. (3) Results: IRPS significantly extenuated the pulmonary pathological lesions and inflammatory response. We used transcriptomic and proteomic analysis to investigate the antiviral effect of IRPS against PRRSV. In the porcine model, 1669 differentially expressed genes (DEGs) and 370 differentially expressed proteins (DEPs) were identified. Analysis of the DEG/DEP-related pathways indicated immune-system and infectious-disease (viral) pathways, such as the NOD-like receptor (NLR) signaling pathway, toll-like receptor (TLR) signaling pathway, and Influenza A-associated signaling pathways. It is noteworthy that IRPS can inhibit NLR-dependent gene expression, then reduce the inflammatory damage. IRPS could exert beneficial effects on the host by regulating the structure of intestinal flora. (4) Conclusions: The antiviral effect of IRPS on PRRSV can be directly achieved by omics techniques. Specifically, the antiviral mechanism of IPRS can be better elucidated by screening target genes and proteins using transcriptome and proteome sequencing, and then performing enrichment and classification according to DEGs and DEPs.

Ribosome profiling of porcine reproductive and respiratory syndrome virus reveals novel features of viral gene expression

The arterivirus porcine reproductive and respiratory syndrome virus (PRRSV) causes significant economic losses to the swine industry worldwide. Here we apply ribosome profiling (RiboSeq) and parallel RNA sequencing (RNASeq) to characterise the transcriptome and translatome of both species of PRRSV and to analyse the host response to infection. We calculated programmed ribosomal frameshift (PRF) efficiency at both sites on the viral genome. This revealed the nsp2 PRF site as the second known example where temporally regulated frameshifting occurs, with increasing -2 PRF efficiency likely facilitated by accumulation of the PRF-stimulatory viral protein, nsp1β. Surprisingly, we find that PRF efficiency at the canonical ORF1ab frameshift site also increases over time, in contradiction of the common assumption that RNA structure-directed frameshift sites operate at a fixed efficiency. This has potential implications for the numerous other viruses with canonical PRF sites. Furthermore, we discovered several highly translated additional viral ORFs, the translation of which may be facilitated by multiple novel viral transcripts. For example, we found a highly expressed 125-codon ORF overlapping nsp12, which is likely translated from novel subgenomic RNA transcripts that overlap the 3' end of ORF1b. Similar transcripts were discovered for both PRRSV-1 and PRRSV-2, suggesting a potential conserved mechanism for temporally regulating expression of the 3'-proximal region of ORF1b. We also identified a highly translated, short upstream ORF in the 5' UTR, the presence of which is highly conserved amongst PRRSV-2 isolates. These findings reveal hidden complexity in the gene expression programmes of these important nidoviruses.

Structural insights into inhibition of PRRSV Nsp4 revealed by structure-based virtual screening, molecular dynamics, and MM-PBSA studies

BACKGROUND

Porcine reproductive and respiratory syndrome respiratory sickness in weaned and growing pigs, as well as sow reproductive failure, and its infection is regarded as one of the most serious swine illnesses worldwide. Given the current lack of an effective treatment, in this study, we identified natural compounds capable of inhibiting non-structural protein 4 (Nsp4) of the virus, which is involved in their replication and pathogenesis.

RESULTS

We screened natural compounds (n = 97,999) obtained from the ZINC database against Nsp4 and selected the top 10 compounds for analysing protein-ligand interactions and physicochemical properties. The five compounds demonstrating strong binding affinity were then subjected to molecular dynamics simulations (100 ns) and binding free energy calculations. Based on analysis, we identified four possible lead compounds that represent potentially effective drug-like inhibitors.

CONCLUSIONS

These methods identified that these natural compounds are capable of inhibiting Nsp4 and possibly effective as antiviral therapeutics against PRRSV.

Time Series Ovarian Transcriptome Analyses of the Porcine Estrous Cycle Reveals Gene Expression Changes during Steroid Metabolism and Corpus Luteum Development

Simple Summary

The estrous cycle, which is divided into follicular and luteal phases based on ovulation, is influenced by reproductive hormones which affect reproduction and cause changes in the reproductive system of the pig. As the main reproductive organ, the ovary is involved in ovulation and changes in the corpus luteum. We aimed to identify dynamic changes in gene expression through differentially expressed gene profiling and to provide a comprehensive understanding of the molecular mechanisms that occur in the pig ovary during the estrous cycle. The transcriptome analysis revealed that the dynamic change in gene expression was more activated in the luteal phase than in the follicular phase. Functional analysis revealed that the metestrus and diestrus periods are important in preparation for pregnancy or the next estrous cycle after ovulation.

Abstract

The porcine estrous cycle is influenced by reproductive hormones, which affect porcine reproduction and result in physiological changes in the reproductive organs. The ovary is involved in ovulation, luteinization, corpus luteum development, and luteolysis. Here, we aimed to provide a comprehensive understanding of the gene expression patterns in porcine ovarian transcriptomes during the estrous cycle through differentially expressed genes profiling and description of molecular mechanisms. The transcriptomes of porcine ovary were obtained during the estrous cycle at three-day intervals from day 0 to day 18 using RNA-seq. At seven time points of the estrous cycle, 4414 DEG were identified; these were classified into three clusters according to their expression patterns. During the late metestrus and diestrus periods, the expression in cluster 1 increased rapidly, and steroid biosynthesis was significant in the pathway. Cluster 2 gene expression patterns represented the cytokine–cytokine receptor interaction in significant pathways. In cluster 3, the hedgehog signaling pathway was selected as the significant pathway. Our study exhibited dynamic gene expression changes with these three different patterns of cluster 1, 2, and 3. The results helped identify the functions and related significant genes especially during the late metestrus and diestrus periods in the estrous cycle.

Time-series transcriptomic analysis of bronchoalveolar lavage cells from virulent and low virulent PRRSV-1-infected piglets

Porcine reproductive and respiratory syndrome virus (PRRSV) has evolved to escape the immune surveillance for a survival advantage leading to a strong modulation of host’s immune responses and favoring secondary bacterial infections. However, limited data are available on how the immunological and transcriptional responses elicited by virulent and low-virulent PRRSV-1 strains are comparable and how they are conserved during the infection. To explore the kinetic transcriptional signature associated with the modulation of host immune response at lung level, a time-series transcriptomic analysis was performed in bronchoalveolar lavage cells upon experimental in vivo infection with two PRRSV-1 strains of different virulence, virulent subtype 3 Lena strain or the low-virulent subtype 1 3249 strain. The time-series analysis revealed overlapping patterns of dysregulated genes enriched in T-cell signaling pathways among both virulent and low-virulent strains, highlighting an upregulation of co-stimulatory and co-inhibitory immune checkpoints that were disclosed as Hub genes. On the other hand, virulent Lena infection induced an early and more marked “negative regulation of immune system process” with an overexpression of co-inhibitory receptors genes related to T-cell and NK cell functions, in association with more severe lung lesion, lung viral load, and BAL cell kinetics. These results underline a complex network of molecular mechanisms governing PRRSV-1 immunopathogenesis at lung level, revealing a pivotal role of co-inhibitory and co-stimulatory immune checkpoints in the pulmonary disease, which may have an impact on T-cell activation and related pathways. These immune checkpoints, together with the regulation of cytokine-signaling pathways, modulated in a virulence-dependent fashion, orchestrate an interplay among pro- and anti-inflammatory responses.

IMPORTANCE Porcine reproductive and respiratory syndrome virus (PRRSV) is one of the major threats to swine health and global production, causing substantial economic losses. We explore the mechanisms involved in the modulation of host immune response at lung level performing a time-series transcriptomic analysis upon experimental infection with two PRRSV-1 strains of different virulence. A complex network of molecular mechanisms was revealed to control the immunopathogenesis of PRRSV-1 infection, highlighting an interplay among pro- and anti-inflammatory responses as a potential mechanism to restrict inflammation-induced lung injury. Moreover, a pivotal role of co-inhibitory and co-stimulatory immune checkpoints was evidenced, which may lead to progressive dysfunction of T cells, impairing viral clearance and leading to persistent infection, favoring as well secondary bacterial infections or viral rebound. However, further studies should be conducted to evaluate the functional role of immune checkpoints in advanced stages of PRRSV infection and explore a possible T-cell exhaustion state.

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.

Porcine Reproductive and Respiratory Syndrome Virus: Immune Escape and Application of Reverse Genetics in Attenuated Live Vaccine Development

Porcine reproductive and respiratory syndrome virus (PRRSV), an RNA virus widely prevalent in pigs, results in significant economic losses worldwide. PRRSV can escape from the host immune response in several processes. Vaccines, including modified live vaccines and inactivated vaccines, are the best available countermeasures against PRRSV infection. However, challenges still exist as the vaccines are not able to induce broad protection. The reason lies in several facts, mainly the variability of PRRSV and the complexity of the interaction between PRRSV and host immune responses, and overcoming these obstacles will require more exploration. Many novel strategies have been proposed to construct more effective vaccines against this evolving and smart virus. In this review, we will describe the mechanisms of how PRRSV induces weak and delayed immune responses, the current vaccines of PRRSV, and the strategies to develop modified live vaccines using reverse genetics systems.

Multi-omics integration strategies for animal epigenetic studies - A review

Genome-wide studies provide considerable insights into the genetic background of animals; however, the inheritance of several heritable factors cannot be elucidated. Epigenetics explains these heritabilities, including those of genes influenced by environmental factors. Knowledge of the mechanisms underlying epigenetics enables understanding the processes of gene regulation through interactions with the environment. Recently developed next-generation sequencing (NGS) technologies help understand the interactional changes in epigenetic mechanisms. There are large sets of NGS data available; however, the integrative data analysis approaches still have limitations with regard to reliably interpreting the epigenetic changes. This review focuses on the epigenetic mechanisms and profiling methods and multi-omics integration methods that can provide comprehensive biological insights in animal genetic studies.