Integrative profiling of gene expression and chromatin accessibility elucidates specific transcriptional networks in porcine neutrophils.

Neutrophils are vital components of the immune system for limiting the invasion and proliferation of pathogens in the body. Surprisingly, the functional annotation of porcine neutrophils is still limited. The transcriptomic and epigenetic assessment of porcine neutrophils from healthy pigs was performed by bulk RNA sequencing and transposase accessible chromatin sequencing (ATAC-seq). First, we sequenced and compared the transcriptome of porcine neutrophils with eight other immune cell transcriptomes to identify a neutrophil-enriched gene list within a detected neutrophil co-expression module. Second, we used ATAC-seq analysis to report for the first time the genome-wide chromatin accessible regions of porcine neutrophils. A combined analysis using both transcriptomic and chromatin accessibility data further defined the neutrophil co-expression network controlled by transcription factors likely important for neutrophil lineage commitment and function. We identified chromatin accessible regions around promoters of neutrophil-specific genes that were predicted to be bound by neutrophil-specific transcription factors. Additionally, published DNA methylation data from porcine immune cells including neutrophils were used to link low DNA methylation patterns to accessible chromatin regions and genes with highly enriched expression in porcine neutrophils. In summary, our data provides the first integrative analysis of the accessible chromatin regions and transcriptional status of porcine neutrophils, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the utility of chromatin accessible regions to identify and enrich our understanding of transcriptional networks in a cell type such as neutrophils.

Salmonella Typhimurium induces genome-wide expression and phosphorylation changes that modulate immune response, intracellular survival and vesicle transport in infected neutrophils.

Salmonella Typhimurium is a food-borne pathogen that causes salmonellosis. When in contact with the host, neutrophils are rapidly recruited to act as first line of defense. To better understand the pathogenesis of this infection, we used an in vitro model of neutrophil infection to perform dual RNA-sequencing (both host and pathogen). In addition, and given that many pathogens interfere with kinase-mediated phosphorylation in host signaling, we performed a phosphoproteomic analysis. The immune response was overall diminished in infected neutrophils, mainly JAK/STAT and toll-like receptor signaling pathways. We found decreased expression of proinflammatory cytokine receptor genes and predicted downregulation of the mitogen-activated protein (MAPK) signaling pathway. Also, Salmonella infection inhibited interferons I and II signaling pathways by upregulation of SOCS3 and subsequent downregulation of STAT1 and STAT2. Additionally, phosphorylation of PSMC2 and PSMC4, proteasome regulatory proteins, was decreased in infected neutrophils. Cell viability and survival was increased by p53 signaling, cell cycle arrest and NFkB-proteasome pathways activation. Combined analysis of RNA-seq and phosphoproteomics also revealed inhibited vesicle transport mechanisms mediated by dynein/dynactin and exocyst complexes, involved in ER-to-Golgi transport and centripetal movement of lysosomes and endosomes. Among the overexpressed virulence genes from Salmonella we found potential effectors responsible of these dysregulations, such as spiC, sopD2, sifA or pipB2, all of them involved in intracellular replication. Our results suggest that Salmonella induces (through overexpression of virulence factors) transcriptional and phosphorylation changes that increases neutrophil survival and shuts down immune response to minimize host response, and impairing intracellular vesicle transport likely to keep nutrients for replication and Salmonella-containing vacuole formation and maintenance.

Assessment of DNA methylation in porcine immune cells reveals novel regulatory elements associated with cell-specific gene expression and immune capacity traits.

BACKGROUND: Genetics studies in the porcine immune system have enhanced selection practices for disease resistance phenotypes and increased the efficacy of porcine models in biomedical research; however limited functional annotation of the porcine immunome has hindered progress on both fronts. Among epigenetic mechanisms that regulate gene expression, DNA methylation is the most ubiquitous modification made to the DNA molecule and influences transcription factor binding as well as gene and phenotype expression. Human and mouse DNA methylation studies have improved mapping of regulatory elements in these species, but comparable studies in the pig have been limited in scope. RESULTS: We performed whole-genome bisulfite sequencing to assess DNA methylation patterns in nine pig immune cell populations: CD21+ and CD21- B cells, four T cell fractions (CD4+, CD8+, CD8+CD4+, and SWC6γδ+), natural killer and myeloid cells, and neutrophils. We identified 54,391 cell differentially methylated regions (cDMRs), and clustering by cDMR methylation rate grouped samples by cell lineage. 32,737 cDMRs were classified as cell lowly methylated regions (cLMRs) in at least one cell type, and cLMRs were broadly enriched in genes and regions of intermediate CpG density. We observed strong correlations between differential methylation and expression across immune cell populations, with cell-specific low methylation disproportionately impacting genes exhibiting enriched gene expression in the same cell type. Motif analysis of cLMRs revealed cell type-specific enrichment of transcription factor binding motifs, indicating that cell-specific methylation patterns may influence accessibility by trans-acting factors. Lastly, cDMRs were enriched for immune capacity GWAS SNPs, and many such overlaps occurred within genes known to influence immune cell development and function (CD8B, NDRG1). CONCLUSION: Our DNA methylation data improve functional annotation of the porcine genome through characterization of epigenomic regulatory patterns that contribute to immune cell identity and function, and increase the potential for identifying mechanistic links between genotype and phenotype.

The Effect of Age, Gender and Comorbidities Upon SARS-CoV-2 Spike Antibody Induction After Two Doses of Sinopharm Vaccine and the Effect of a Pfizer/BioNtech Booster Vaccine.

Severe acute respiratory syndrome coronavirus (SARS-CoV)-2 emerged in China in 2019 and has since travelled the world infecting millions. SARS-CoV-2 causes Corona Virus Disease (COVID-19), that has to date taken over 4 million lives. The Kingdom of Bahrain's vaccine roll-out has consisted of Sinopharm's BBIBP-CorV (Sinopharm) and Pfizer/BioNtech's BNT162b2 (Pfizer/BioNtech). Testing for SARS-CoV-2 anti-Spike (S) antibodies is a useful technique in estimating an individual's immune protection against the infection. In this study we evaluated S antibody levels by electro-chemiluminescence immunoassay in 379 individuals double vaccinated with Sinopharm and 15 of whom were given a booster with the Pfizer/BioNtech vaccine. Among our double vaccinated cohort, we found a spectrum of S antibody levels. Indeed, we found that a significant proportion of individuals with low S antibody levels had clinical conditions, which were mainly immune-related disorders. Furthermore, a significant proportion of individuals with low S antibody levels were above 50 years of age. Finally, we observed a significant increase in S antibody levels after the Pfizer/BioNtech booster was administered. These findings reveal that while a large proportion of Sinopharm vaccinated individuals did not develop high levels of antibodies against the S protein, a booster dose of the Pfizer/BioNtech vaccine significantly enhances S antibody levels, revealing this "triple dose" vaccination strategy as a useful method of ensuring protective immunity against SARS-CoV-2.

Study of microRNA expression in Salmonella Typhimurium-infected porcine ileum reveals miR-194a-5p as an important regulator of the TLR4-mediated inflammatory response.

Infection with Salmonella Typhimurium (S. Typhimurium) is a common cause of food-borne zoonosis leading to acute gastroenteritis in humans and pigs, causing economic losses to producers and farmers, and generating a food security risk. In a previous study, we demonstrated that S. Typhimurium infection produces a severe transcriptional activation of inflammatory processes in ileum. However, little is known regarding how microRNAs regulate this response during infection. Here, small RNA sequencing was used to identify 28 miRNAs differentially expressed (DE) in ileum of S. Typhimurium-infected pigs, which potentially regulate 14 target genes involved in immune system processes such as regulation of cytokine production, monocyte chemotaxis, or cellular response to interferon gamma. Using in vitro functional and gain/loss of function (mimics/CRISPR-Cas system) approaches, we show that porcine miR-194a-5p (homologous to human miR-194-5p) regulates TLR4 gene expression, an important molecule involved in pathogen virulence, recognition and activation of innate immunity in Salmonella infection.

Reference Transcriptomes of Porcine Peripheral Immune Cells Created Through Bulk and Single-Cell RNA Sequencing.

Pigs are a valuable human biomedical model and an important protein source supporting global food security. The transcriptomes of peripheral blood immune cells in pigs were defined at the bulk cell-type and single cell levels. First, eight cell types were isolated in bulk from peripheral blood mononuclear cells (PBMCs) by cell sorting, representing Myeloid, NK cells and specific populations of T and B-cells. Transcriptomes for each bulk population of cells were generated by RNA-seq with 10,974 expressed genes detected. Pairwise comparisons between cell types revealed specific expression, while enrichment analysis identified 1,885 to 3,591 significantly enriched genes across all 8 cell types. Gene Ontology analysis for the top 25% of significantly enriched genes (SEG) showed high enrichment of biological processes related to the nature of each cell type. Comparison of gene expression indicated highly significant correlations between pig cells and corresponding human PBMC bulk RNA-seq data available in Haemopedia. Second, higher resolution of distinct cell populations was obtained by single-cell RNA-sequencing (scRNA-seq) of PBMC. Seven PBMC samples were partitioned and sequenced that produced 28,810 single cell transcriptomes distributed across 36 clusters and classified into 13 general cell types including plasmacytoid dendritic cells (DC), conventional DCs, monocytes, B-cell, conventional CD4 and CD8 αß T-cells, NK cells, and γδ T-cells. Signature gene sets from the human Haemopedia data were assessed for relative enrichment in genes expressed in pig cells and integration of pig scRNA-seq with a public human scRNA-seq dataset provided further validation for similarity between human and pig data. The sorted porcine bulk RNAseq dataset informed classification of scRNA-seq PBMC populations; specifically, an integration of the datasets showed that the pig bulk RNAseq data helped define the CD4CD8 double-positive T-cell populations in the scRNA-seq data. Overall, the data provides deep and well-validated transcriptomic data from sorted PBMC populations and the first single-cell transcriptomic data for porcine PBMCs. This resource will be invaluable for annotation of pig genes controlling immunogenetic traits as part of the porcine Functional Annotation of Animal Genomes (FAANG) project, as well as further study of, and development of new reagents for, porcine immunology.

Changes in H3K27ac at Gene Regulatory Regions in Porcine Alveolar Macrophages Following LPS or PolyIC Exposure.

Changes in chromatin structure, especially in histone modifications (HMs), linked with chromatin accessibility for transcription machinery, are considered to play significant roles in transcriptional regulation. Alveolar macrophages (AM) are important immune cells for protection against pulmonary pathogens, and must readily respond to bacteria and viruses that enter the airways. Mechanism(s) controlling AM innate response to different pathogen-associated molecular patterns (PAMPs) are not well defined in pigs. By combining RNA sequencing (RNA-seq) with chromatin immunoprecipitation and sequencing (ChIP-seq) for four histone marks (H3K4me3, H3K4me1, H3K27ac and H3K27me3), we established a chromatin state map for AM stimulated with two different PAMPs, lipopolysaccharide (LPS) and Poly(I:C), and investigated the potential effect of identified histone modifications on transcription factor binding motif (TFBM) prediction and RNA abundance changes in these AM. The integrative analysis suggests that the differential gene expression between non-stimulated and stimulated AM is significantly associated with changes in the H3K27ac level at active regulatory regions. Although global changes in chromatin states were minor after stimulation, we detected chromatin state changes for differentially expressed genes involved in the TLR4, TLR3 and RIG-I signaling pathways. We found that regions marked by H3K27ac genome-wide were enriched for TFBMs of TF that are involved in the inflammatory response. We further documented that TF whose expression was induced by these stimuli had TFBMs enriched within H3K27ac-marked regions whose chromatin state changed by these same stimuli. Given that the dramatic transcriptomic changes and minor chromatin state changes occurred in response to both stimuli, we conclude that regulatory elements (i.e. active promoters) that contain transcription factor binding motifs were already active/poised in AM for immediate inflammatory response to PAMPs. In summary, our data provides the first chromatin state map of porcine AM in response to bacterial and viral PAMPs, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the role of HMs, especially H3K27ac, in regulating transcription in AM in response to LPS and Poly(I:C).

Comparative proteomic analysis reveals different responses in porcine lymph nodes to virulent and attenuated homologous African swine fever virus strains.

African swine fever (ASF) is a pathology of pigs against which there is no treatment or vaccine. Understanding the equilibrium between innate and adaptive protective responses and immune pathology might contribute to the development of strategies against ASFV. Here we compare, using a proteomic approach, the course of the in vivo infection caused by two homologous strains: the virulent E75 and the attenuated E75CV1. Our results show a progressive loss of proteins by day 7 post-infection (pi) with E75, reflecting tissue destruction. Many signal pathways were affected by both infections but in different ways and extensions. Cytoskeletal remodelling and clathrin-endocytosis were affected by both isolates, while a greater number of proteins involved on inflammatory and immunological pathways were altered by E75CV1. 14-3-3 mediated signalling, related to immunity and apoptosis, was inhibited by both isolates. The implication of the Rho GTPases by E75CV1 throughout infection is also evident. Early events reflected the lack of E75 recognition by the immune system, an evasion strategy acquired by the virulent strains, and significant changes at 7 days post-infection (dpi), coinciding with the peak of infection and the time of death. The protein signature at day 31 pi with E75CV1 seems to reflect events observed at 1 dpi, including the upregulation of proteosomal subunits and molecules described as autoantigens (vimentin, HSPB1, enolase and lymphocyte cytosolic protein 1), which allow the speculation that auto-antibodies could contribute to chronic ASFV infections. Therefore, the use of proteomics could help understand ASFV pathogenesis and immune protection, opening new avenues for future research.

Regulatory role of microRNA in mesenteric lymph nodes after Salmonella Typhimurium infection.

Salmonellosis is a gastrointestinal disease caused by non-typhoidal Salmonella serovars such as Salmonella Typhimurium. This pathology is a zoonosis, and food animals with subclinical infection constitute a vast reservoir for disease. After intestinal colonization, Salmonella Typhimurium reaches mesenteric lymph nodes (MLN), where infection is controlled avoiding systemic spread. Although the molecular basis of this infection has been extensively studied, little is known about how microRNA (miRNA) regulate the expression of proteins involved in the Salmonella-host interaction. Using small RNA-seq, we examined expression profiles of MLN 2 days after infection with Salmonella Typhimurium, and we found 110 dysregulated miRNA. Among them, we found upregulated miR-21, miR-155, miR-150, and miR-221, as well as downregulated miR-143 and miR-125, all of them previously linked to other bacterial infections. Integration with proteomic data revealed 30 miRNA potentially regulating the expression of 15 proteins involved in biological functions such as cell death and survival, inflammatory response and antigenic presentation. The inflammatory response was found increased via upregulation of miRNA such as miR-21 and miR-155. Downregulation of miR-125a/b, miR-148 and miR-1 were identified as potential regulators of MHC-class I components PSMB8, HSP90B1 and PDIA3, respectively. Furthermore, we confirmed that miR-125a is a direct target of immunoproteasome component PSMB8. Since we also found miR-130 downregulation, which is associated with upregulation of HSPA8, we suggest induction of both MHC-I and MHC-II antigen presentation pathways. In conclusion, our study identifies miRNA that could regulate critical networks for antigenic presentation, inflammatory response and cytoskeletal rearrangements.

Physical training and weight loss in dogs lead to transcriptional changes in genes involved in the glucose-transport pathway in muscle and adipose tissues.

Obesity is a worldwide problem in humans and domestic animals. Interventions, including a combination of dietary management and exercise, have proven to be effective for inducing weight loss in humans. In companion animals, the role of exercise in the management of obesity has received relatively little attention. The aim of the present study was to investigate changes in the transcriptome of key energy metabolism genes in muscle and adipose tissues in response to diet-induced weight loss alone, or combined with exercise in dogs. Overweight pet dogs were enrolled on a weight loss programme, based on calorie restriction and physical training (FD group, n = 5) or calorie restriction alone (DO group, n = 7). mRNA expression of 12 genes and six microRNAs were investigated using quantitative real-time PCR (qPCR). In the FD group, FOXO1 and RAC1 were expressed at lower levels in adipose tissue, whereas ESRRA and AKT2 were more highly expressed in muscle, when compared with the DO group. Comparing expression before and after the intervention, in the DO group, nine genes and three microRNAs showed significant altered expression in adipose tissue (PPARG, ADIPOQ and FOXO1; P < 0.001) and seven genes and two microRNAs were significantly downregulated (NRF2, RAC1, ESRRA, AKT2, PGC1a and mir-23; P < 0.001) in muscle. Thus, calorie restriction causes regulation of several metabolic genes in both tissues. The mild exercise, incorporated into this study design, was sufficient to elicit transcriptional changes in adipose and muscle tissues, suggesting a positive effect on glucose metabolism. The study findings support inclusion of exercise in management of canine obesity.

Wet-lab tested microRNA assays for qPCR studies with SYBR® Green and DNA primers in pig tissues.

MicroRNAs are key post-transcriptional regulators of gene expression that are involved in several biological processes including those that mediate disease pathophysiology. Hence, quantifying microRNA expression levels can provide important and novel insights into disease biology. In recent years, the pig has emerged as an excellent large animal model for studying human diseases and conditions (e.g. obesity) due to similarities in organ size, gastro-intestinal tract, metabolism, immune response, genetics and the availability of relevant tissues that are not normally easily available in humans. We have previously developed two useful tools in the field of microRNA quantitative real time PCR (qPCR): 1) a very specific, sensitive and simple qPCR method based on DNA primers, MiR-specific qPCR; and 2) the free primer-design software miRprimer. The present study integrates in a publicly accessible database all available information on validated porcine microRNA qPCR assays that have utilized these tools. Due to the high phylogenetic conservation in microRNA sequence between pig, humans and other domestic species this database is a very valuable resource for the broader scientist community who are working on microRNAs and want to use readily tested qPCR assays in a simple and cost-effective manner.