Transcriptomic Analysis Identifies Candidate Genes and Gene Sets Controlling the Response of Porcine Peripheral Blood Mononuclear Cells to Poly I:C Stimulation

Pathogen stimulations and immune cells synergistically affect the gene expression profile characteristics of porcine peripheral blood mononuclear cells

BACKGROUND

Pigs serve as a crucial source of protein in the human diet and play a fundamental role in ensuring food security. However, infectious diseases caused by bacteria or viruses are a major threat to effective global pig farming, jeopardizing human health. Peripheral blood mononuclear cells (PBMCs) are a mixture of immune cells that play crucial roles in immunity and disease resistance in pigs. Previous studies on the gene expression regulation patterns of PBMCs have concentrated on a single immune stimulus or immune cell subpopulation, which has limited our comprehensive understanding of the mechanisms of the pig immune response.

RESULTS

Here, we integrated and re-analyzed RNA-seq data published online for porcine PBMC stimulated by lipopolysaccharide (LPS), polyinosinic acid (PolyI:C), and various unknown microorganisms (EM). The results revealed that gene expression and its functional characterization are highly specific to the pathogen, identifying 603, 254, and 882 pathogen-specific genes and 38 shared genes, respectively. Notably, LPS and PolyI:C stimulation directly triggered inflammatory and immune-response pathways, while exposure to mixed microbes (EM) enhanced metabolic processes. These pathogen-specific genes were enriched in immune trait-associated quantitative trait loci (QTL) and eGenes in porcine immune tissues and were implicated in specific cell types. Furthermore, we discussed the roles of eQTLs rs3473322705 and rs1109431654 in regulating pathogen- and cell-specific genes CD300A and CD93, using cellular experiments. Additionally, by integrating genome-wide association studies datasets from 33 complex traits and diseases in humans, we found that pathogen-specific genes were significantly enriched for immune traits and metabolic diseases.

CONCLUSIONS

We systematically analyzed the gene expression profiles of the three stimulations and demonstrated pathogen-specific and cell-specific gene regulation across different stimulations in porcine PBMCs. These findings enhance our understanding of shared and distinct regulatory mechanisms of genetic variants in pig immune traits.

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.

Development of an In Vitro Model for Inflammation Mediated Renal Toxicity Using 3D Renal Tubules and Co-Cultured Human Immune Cells

BACKGROUND

The emergence of various infectious diseases and the toxic effects of hyperinflammation by biotherapeutics have highlighted the need for in vitro preclinical models mimicking the human immune system. In vitro models studying the relationship between hyperinflammation and acute renal injury mainly rely on 2D culture systems, which have shown limitations in recapitulating kidney function. Herein, we developed an in vitro kidney toxicity model by co-culturing 3D engineered kidney proximal tubules cells (RPTEC/TERT1) with human peripheral blood mononuclear cells (PBMC).

METHODS

RPTEC/TERT1 were sandwich cultured to form 3D renal tubules for 16 days. The tubules were then co-cultured with PBMC using transwell (0.4 μm pores) for 24 h. Hyperinflammation of PBMC was induced during co-culture using polyinosinic-polycytidylic acid (polyI:C) and lipopolysaccharide (LPS) to investigate the effects of the induced hyperinflammation on the renal tubules.

RESULTS

Encapsulated RPTEC/TERT1 cells in Matrigel exhibited elevated renal function markers compared to 2D culture. The coexistence of PBMC and polyI:C induced a strong inflammatory response in the kidney cells. This hyperinflammation significantly reduced primary cilia formation and upregulated kidney injury markers along the 3D tubules. Similarly, treating co-cultured PBMC with LPS to induce hyperinflammation resulted in comparable inflammatory responses and potential kidney injury.

CONCLUSION

The model demonstrated similar changes in kidney injury markers following polyI:C and LPS treatment, indicating its suitability for detecting immune-associated kidney damage resulting from infections and biopharmaceutical applications.

Quantification of camelid cytokine mRNA expression in PBMCs by microfluidic qPCR technology

Camelids are economically and socially important in several parts of the world and might carry pathogens with epizootic or zoonotic potential. However, biological research in these species is limited due to lack of reagents. Here, we developed RT-qPCR assays to quantify a panel of camelid innate and adaptive immune response genes, which can be monitored in a single run. The assays were validated with PHA, PMA-ionomycin, and Poly I:C-stimulated PBMCs from alpaca, dromedary camel and llama, including normalization by multiple reference genes. Further, comparative gene expression analyses for the different camelid species were performed by a unique microfluidic qPCR assay. Compared to unstimulated controls, PHA and PMA-ionomycin stimulation elicited robust Th1 and Th2 responses in PBMCs from camelid species. Additional activation of type I and type III IFN signalling pathways was described exclusively in PHA-stimulated dromedary lymphocytes, in contrast to those from alpaca and llama. We also found that PolyI:C stimulation induced robust antiviral response genes in alpaca PBMCs. The proposed methodology should be useful for the measurement of immune responses to infection or vaccination in camelid species.

Transcriptional Profiling of Human Peripheral Blood Mononuclear Cells Stimulated by Mycobacterium tuberculosis PPE57 Identifies Characteristic Genes Associated With Type I Interferon Signaling

Proline-glutamic acid (PE)- and proline-proline-glutamic acid (PPE)-containing proteins are exclusive to Mycobacterium tuberculosis (MTB), the leading cause of tuberculosis (TB). In this study, we performed global transcriptome sequencing (RNA-Seq) on PPE57-stimulated peripheral blood mononuclear cells (PBMCs) and control samples to quantitatively measure the expression level of key transcripts of interest. A total of 1367 differentially expressed genes (DEGs) were observed in response to a 6 h exposure to PPE57, with 685 being up-regulated and 682 down-regulated. Immune-related gene functions and pathways associated with these genes were evaluated, revealing that the type I IFN signaling pathway was the most significantly enriched pathway in our RNA-seq dataset, with 14 DEGs identified therein including ISG15, MX2, IRF9, IFIT3, IFIT2, OAS3, IFIT1, IFI6, OAS2, OASL, RSAD2, OAS1, IRF7, and MX1. These PPE57-related transcriptomic profiles have implications for a better understanding of host global immune mechanisms underlying MTB infection outcomes. However, more studies regarding these DEGs and type I IFN signaling in this infectious context are necessary to more fully clarify the underlying mechanisms that arise in response to PPE57 during MTB infection.

Characterization of the whole transcriptome of spleens from Chinese indigenous breed Ningxiang pig reveals diverse coding and non-coding RNAs for immunity regulation

The spatio-temporal expression patterns of RNA and comparisons between different developmental stages have been one of the useful techniques for studying animal physiology and functional gene regulations. A Chinese indigenous breed Ningxiang pig is known for its quality meat production, disease resistance and slow growth performances in pig industry. To gain a better understanding of pig immunity and disease resistance, we comprehensively analyzed the whole transcriptome of the spleens from three important developmental nodes of Ningxiang pig at 30, 90 and 210 days of age. By three ways of comparisons (30vs 90 days, 30 vs 210 days and 90 vs 210 days), a total of 364to 865 differentially expressed mRNAs, 37 to 98 differentially expressed miRNAs,220 to 278 lncRNAs, and 96 to 113 circRNAs were identified. Further analysis of expression patterns, potential function and interactions with miRNAs identified the potential non-coding RNAs related to immunomodulation such as ssc-miRNA-150, ssc-miRNA-497, MSTRG24160, MSTRG18646. The results revealed that miRNAs and circRNAs may have evolved to regulate a large set of biological processes of spleen function in Ningxiang pigs, and circRNAs play a role of miRNA sponges. The results from study is the first report of whole transcriptome analysis of Ningxiang pig spleen and provide new insights into the expression changes of RNAs during the spleen development, which contribute to the phenotypic formation of immunity and disease resistancesin Chinese indigenous pig breeds.

Selection and validation of differentially expressed metabolic and immune genes in weaned Ghurrah versus crossbred piglets

In the present investigation, differentially expressed genes (DEGs) were studied using RNA sequencing (RNA-seq) technique in porcine peripheral blood mononuclear cells (PBMC) of weaned Ghurrah and crossbred piglets at 3-month age. Transcriptomic analysis was done using three different packages, namely, EBSeq, DESeq2, and edgeR, to identify the DEGs between Ghurrah and crossbred piglets. Total 7717 DEGs were commonly identified by all three packages, out of which 4151 genes found to be up-regulated, and 3566 genes were down-regulated. Functional annotation of these DEGs indicated metabolism as the most commonly enriched category followed by the immune response. Genes related to metabolism and growth were up-regulated in crossbred piglets as compared with Ghurrah piglets, whereas immunity-related genes were up-regulated in Ghurrah piglets elucidating the disease resistance nature of this indigenous breed over crossbred counterparts. Further, eight DEGs, namely, LRP-1, ADCY4, ERRFI1, LDHD, ARG1, OASL, MGARP, and S100A8, were validated by qRT-PCR in a separate set of biological samples and found to be in concordance with RNA-seq results. Finding in the present study provides insight into genes and their molecular mechanisms governing difference in growth performance between Ghurrah and crossbred pigs.

Characterization of immune pleiotropy of ESR1 gene in pigs

It is well known that the estrogen receptor alpha gene (ESR1) affects the reproductive traits of pigs; however, the immune role of ESR1 gene has not yet been resolved. Here, we characterized the pleiotropic aspects of ESR1 gene in immunity using the pig model. Tissue expression profile showed that the ESR1 gene had a broad ectopic expression in multiple reproductive and immune-related tissues/organs, which provided the tissue-level spatial fundamental of ESR1 gene that might function as a pleiotropic immune regulator. Using the peripheral blood cell model, a coupling transcriptome analytical strategy was proposed and verified that there existed strong positive or negative correlations of ESR1 gene with hundreds of differentially expressed genes that were involved in the immune regulation, indicating that the ESR1 gene might affect or be affected by, directly or indirectly, dozens of immune-related genes in the peripheral blood cells. Furthermore, the results of genetic association analysis showed that the SmaI-polymorphism of ESR1 gene had significant or highly significant associations with multiple immune traits, including platelet (PLT), hematocrit (HCT), the number of CD4-CD8-CD3- cells, plateletcrit (PCT), mean corpuscular volume (MCV), and mean corpuscular hemoglobin concentration (MCHC). Multiple evidences supported the immune pleiotropic roles of ESR1 gene in pigs. The study advances our understanding of the cross-species immune pleiotropic landscape of ESR1 gene and also provides a potential pleiotropic molecular marker for disease-resistant breeding in pigs.

The muscle development transcriptome landscape of ovariectomized goat

In practical production, almost all rams and about 50% of ewes are used to fatten. Researchers have proved that ewe ovariectomy could improve the productivity significantly, but the specific molecular mechanism is still unknown. In this study, five independent cDNA libraries (three and two from ovariectomized and normal ewe longissimus dorsi samples, respectively) were constructed to thoroughly explore the global transcriptome, further to reveal how the ovariectomized ewes influence muscle development by Illumina2000 sequencing technology. As a result, 205 358 transcripts and 118 264 unigenes were generated. 15 490 simple sequence repeats (SSRs) were revealed and divided into six types, and the short repeat sequence SSR (monomers, dimers, trimers) was the domain type. Single nucleotide polymorphism analysis found that the number of transition was greater than the number of transversion among the five libraries. Furthermore, 1612 differently expressed genes (DEGs) (Log2fold_change > 1 and p < 0.05) were revealed between ovariectomized and normal ewe groups, in which 903 genes were expressed commonly in the two groups, and 288 and 421 genes were uniquely expressed in normal and ovariectomized ewe groups, respectively. Gene Ontology (GO) analysis categorized all unigenes into 555 GO terms and 56 DEGs were significantly categorized into 43 GO terms (p < 0.05). KEGG enrichment analysis annotated 12 976 genes (containing 137 DEGs) to 86 pathways, among them 24 and 11 DEGs involved in development and reproduction associated pathways, respectively. To validate the reliability of the RNA-seq analysis, 22 candidate DEGs were randomly selected to perform quantitative real-time polymerase chain reaction. The result showed that 9 and 1 genes were significantly and approximately significantly expressed in control and treatment group, respectively, and the results of RNA-seq are believable in this study. Overall, these results were helpful for elucidating the molecular mechanism of muscle development of ovariectomized animals and the application of female ovariectomy in fattening.

Whole Blood Transcriptome Sequencing Reveals Gene Expression Differences between Dapulian and Landrace Piglets

There is little genomic information regarding gene expression differences at the whole blood transcriptome level of different pig breeds at the neonatal stage. To solve this, we characterized differentially expressed genes (DEGs) in the whole blood of Dapulian (DPL) and Landrace piglets using RNA-seq (RNA-sequencing) technology. In this study, 83 DEGs were identified between the two breeds. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses identified immune response and metabolism as the most commonly enriched terms and pathways in the DEGs. Genes related to immunity and lipid metabolism were more highly expressed in the DPL piglets, while genes related to body growth were more highly expressed in the Landrace piglets. Additionally, the DPL piglets had twofold more single nucleotide polymorphisms (SNPs) and alternative splicing (AS) than the Landrace piglets. These results expand our knowledge of the genes transcribed in the piglet whole blood of two breeds and provide a basis for future research of the molecular mechanisms underlying the piglet differences.

MicroRNA Transcriptome of Poly I:C-Stimulated Peripheral Blood Mononuclear Cells Reveals Evidence for MicroRNAs in Regulating Host Response to RNA Viruses in Pigs

MicroRNAs (miRNAs) are one family of small noncoding RNAs that function to modulate the activity of specific mRNA targets in animals. To understand the role of miRNAs in regulating genes involved in the host immune response to RNA viruses, we profiled and characterized the miRNAs of swine peripheral blood mononuclear cells (PBMC) stimulated with poly I:C, a synthetic dsRNA analog, by miRNA-sequencing (miRNA-seq). We identified a total of 905 miRNAs, of which 503 miRNAs were firstly exploited herein with no annotation in the latest miRBase 21.0. Expression analysis demonstrated that poly I:C stimulation can elicit significantly differentially expressed (DE) miRNAs in Dapulian (n = 20), one Chinese indigenous breed, as well as Landrace (n = 23). By integrating the mRNA expression profiles of the same sample with miRNA profiles, we carried out function analyses of the target genes of these DE miRNAs, with the results indicating that target genes were most enriched in some immune-related pathways and gene ontology (GO) terms, suggesting that DE miRNAs play an important role in the regulation of host to poly I:C stimulation. Furthermore, we also detected 43 and 61 significantly DE miRNAs between the two breeds in the control sample groups and poly I:C stimulation groups, respectively, which may be involved in regulation of the different characteristics of the two breeds. This study describes for the first time the PBMC miRNA transcriptomic response to poly I:C stimulation in pigs, which not only contributes to a broad view of the pig miRNAome but improves our understanding of miRNA function in regulating host immune response to RNA viruses.