Messenger RNA Sequencing and Pathway Analysis Provide Novel Insights Into the Susceptibility to Salmonella enteritidis Infection in Chickens

Differential responses to avian pathogenic E. coli and the regulatory role of splenic miRNAs in APEC infection in Silkie chickens

Avian pathogenic Escherichia coli (APEC) is a bacterial disease that harms the poultry industry worldwide, but its effect on Chinese Silkie has not been reported. Studies on whether there are differences in Silkie individual resistance to APEC and the regulatory role of spleen miRNAs lay the foundation for strategies against APEC. Therefore, 270 Silkie chickens were infected with the median lethal dose of an E. coli O1, O2, and O78 mixture. These chickens were divided into a susceptible group (Group S) and a recovery group (Group R) according to whether they survived 15 days postinfection (dpi). Moreover, 90 uninfected APEC Silkie served as controls (Group C). The splenic miRNA expression profile was examined to evaluate the role of miRNAs in the APEC infection response. Of the 270 Silkies infected with APEC, 144 were alive at 15 dpi. Cluster analysis and principal component analysis (PCA) of splenic miRNAs revealed that the four Group R replicates were clustered with the three Group C replicates and were far from the three Group S replicates. Differentially expressed (DE) miRNAs, especially gga-miR-146b-5p, play essential roles in immune and inflammatory responses to APEC. Functional enrichment analyses of DEmiRNAs suggested that suppression of immune system processes (biological processes) might contribute to susceptibility to APEC and that FoxO signaling pathways might be closely associated with the APEC infection response and postinfection repair. This study paves the way for screening anti-APEC Silkies and provides novel insights into the regulatory role of miRNAs in APEC infection.

Identification of key transcription factors and their functional role involved in Salmonella typhimurium infection in chicken using integrated transcriptome analysis and bioinformatics approach

Salmonella enterica serovar typhimurium is the cause of significant morbidity and mortality worldwide that causes economic losses to poultry and is able to cause infection in humans. Indigenous chicken breeds are a potential source of animal protein and have the added advantage of being disease resistant. An indigenous chicken, Kashmir favorella and commercial broiler were selected for understanding the mechanism of disease resistance. Following infection in Kashmir favorella, three differentially expressed genes Nuclear Factor Kappa B (NF-κB1), Forkhead Box Protein O3 (FOXO3) and Paired box 5 (Pax5) were identified. FOXO3, a transcriptional activator, is the potential marker of host resistance in Salmonella infection. NF-κB1 is an inducible transcription factor which lays the foundation for studying gene network of the innate immune response of Salmonella infection in chicken. Pax5 is essential for differentiation of pre-B cells into mature B cell. The real time PCR analysis showed that in response to Salmonella Typhimurium infection a remarkable increase of NF-κB1 (P˂0.01)_, FOXO3 (P˂0.01) gene expression in liver and Pax5 (P˂0.01) gene expression in spleen of Kashmir favorella was observed. The protein-protein interaction (PPI) and protein-TF interaction network by STRINGDB analysis suggests that FOXO3 is a hub gene in the network and is closely related to Salmonella infection along with NF-κB1. All the three differentially expressed genes (NF-κB1, FOXO3 and PaX5) showed their influence on 12 interacting proteins and 16 TFs, where cyclic adenosine monophosphate Response Element Binding protein (CREBBP), erythroblast transformation-specific (ETSI), Tumour-protein 53(TP53I), IKKBK, lymphoid enhancer-binding factor-1 (LEF1), and interferon regulatory factor-4 (IRF4) play role in immune responses. This study shall pave the way for newer strategies for treatment and prevention of Salmonella infection and may help in increasing the innate disease resistance.

Phosphoproteomics Profile of Chicken Cecum in the Response to Salmonella enterica Serovar Enteritidis Inoculation

Salmonella enterica serovar Enteritidis (S. Enteritidis) is a foodborne pathogen, which can cause great threats to human health through the consumption of contaminated poultry products. This research combines TMT labeling, HPLC and mass-spectrometry-based phosphoproteomics on cecum of the F1 cross of Guangxi Yao chicken and Jining Bairi chicken. The treated group was inoculated with 0.3 mL inoculum S. Enteritidis, and the control group was inoculated with 0.3 mL phosphate-buffered saline (PBS). A total of 338 differentially phosphorylated modification sites in 243 differentially phosphorylated proteins (DPPs) were chosen for downstream analyses. A total of 213 sites in 146 DPPs were up-regulated and 125 sites in 97 DPPs were down-regulated. Functional analysis was performed for DPPs based on gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and the protein domain. The DPPs were mainly enriched in immune- and metabolic-related GO-BP (biological process) and KEGG pathways. We predicted and classified the subcellular structure and COG/KOG of DPPs. Furthermore, protein-protein interaction network analyses were performed by using multiple algorithms. We identified 71 motifs of the phosphorylated modification sites and selected 18 sites randomly to detect the expression level through parallel reaction monitoring (PRM). S. Enteritidis inoculation caused phosphorylation alteration in immune- and metabolic-related proteins. The invasion of S. Enteritidis may be actualized by inducing cecum cell apoptosis through the endoplasmic reticulum pathway, and chickens could resist the invasion of S. Enteritidis by affecting the function of ECM receptors. The findings herein provide a crucial theoretical foundation to understand the molecular mechanism and epigenetic regulation in response to S. Enteritidis inoculation in chickens.

A transcriptomic analysis of the effects of macrophage polarization and endotoxin tolerance on the response to Salmonella

Salmonella is an intracellular pathogen causing significant morbidity and mortality. Its ability to grow inside macrophages is important to virulence, and is dependent on the activation state of the macrophages. Classically activated M1 macrophages are non-permissive for Salmonella growth, while alternatively activated M2 macrophages are permissive for Salmonella growth. Here we showed that endotoxin-primed macrophages (MEP), such as those associated with sepsis, showed similar levels of Salmonella resistance to M1 macrophages after 2 hr of intracellular infection, but at the 4 hr and 24 hr time points were susceptible like M2 macrophages. To understand this mechanistically, transcriptomic sequencing, RNA-Seq, was performed. This showed that M1 and MEP macrophages that had not been exposed to Salmonella, demonstrated a process termed here as primed activation, in expressing relatively higher levels of particular anti-infective genes and pathways, including the JAK-STAT (Janus kinase-signal transducer and activator of transcription) pathway. In contrast, in M2 macrophages these genes and pathways were largely expressed only in response to infection. Conversely, in response to infection, M1 macrophages, but not MEP macrophages, modulated additional genes known to be associated with susceptibility to Salmonella infection, possibly contributing to the differences in resistance at later time points. Application of the JAK inhibitor Ruxolitinib before infection reduced resistance in M1 macrophages, supporting the importance of early JAK-STAT signalling in M1 resistance to Salmonella.

Comparative RNA-Seq analysis reveals insights in Salmonella disease resistance of chicken; and database development as resource for gene expression in poultry

Salmonella, one of the major infectious diseases in poultry, causes considerable economic losses in terms of mortality and morbidity, especially in countries that lack effective vaccination programs. Besides being resistant to diseases, indigenous chicken breeds are also a potential source of animal protein in developing countries. For understanding the disease resistance, an indigenous chicken line Kashmir faverolla, and commercial broiler were selected. RNA-seq was performed after challenging the chicken with Salmonella Typhimurium. Comparative differential expression results showed that following infection, a total of 3153 genes and 1787 genes were differentially expressed in the liver and spleen, respectively. The genes that were differentially expressed included interleukins, cytokines, NOS2, Avβ-defensins, toll-like receptors, and other immune-related gene families. Most of the genes and signaling pathways involved in the innate and adaptive immune responses against bacterial infection were significantly enriched in the Kashmir faverolla. Pathway analysis revealed that most of the enriched pathways were MAPK signaling pathway, NOD-like receptor signaling pathway, TLR signaling pathway, PPAR signaling pathway, endocytosis, etc. Surprisingly some immune-related genes like TLRs were upregulated in the susceptible chicken breed. On postmortem examination, the resistant birds showed small lesions in the liver compared to large necrotic lesions in susceptible birds. The pathological manifestations and RNA sequencing results suggest a balancing link between resistance and infection tolerance in Kashmir faverolla. Here we also developed an online Poultry Infection Database (https://skuastk.org/pif/index.html), the first publicly available gene expression resource for disease resistance in chickens. The available database not only shows the data for gene expression in chicken tissues but also provides quick search, visualization and download capacity.

Differential immunological response detected in mRNA expression profiles among diverse chicken lines in response to Salmonella challenge

Salmonella enterica serovar Enteritidis is a bacterial pathogen that contributes to poultry production losses and human foodborne illness. The bacterium elicits a broad immune response involving both the innate and adaptive components of the immune system. Coordination of the immune response is largely directed by cytokines. The objective of the current study was to characterize the expression of a select set of cytokines and regulatory immune genes in three genetically diverse chicken lines after infection with S. Enteritidis. Leghorn, Fayoumi and broiler day-old chicks were orally infected with pathogenic S. Enteritidis or culture medium. At 2 and 18 h postinfection, spleens and ceca were collected and mRNA expression levels for 7 genes (GM-CSF, IL2, IL15, TGF-β1, SOCS3, P20K, and MHC class IIβ) were evaluated by real-time quantitative PCR. Genetic line had a significant effect on mRNA expression levels of IL15, TGF-β1, SOCS3 and P20K in the spleen and on P20K and MHC class IIβ in the cecum. Comparing challenged vs. unchallenged birds, the expression of SOCS3 and P20K mRNA were significantly higher in the spleen and cecum, while MHC class IIβ mRNA was significantly lower in spleen. Combining the current RNA expression results with those of previously reported studies on the same samples reveals distinct RNA expression profiles among the three genetic chicken lines and the 2 tissues. This study illustrates that these diverse genetic lines have distinctively different immune response to S. Enteritidis challenge within the spleen and the cecum.

Quantitative trait loci and transcriptome signatures associated with avian heritable resistance to Campylobacter

Campylobacter is the leading cause of bacterial foodborne gastroenteritis worldwide. Handling or consumption of contaminated poultry meat is a key risk factor for human campylobacteriosis. One potential control strategy is to select poultry with increased resistance to Campylobacter. We associated high-density genome-wide genotypes (600K single nucleotide polymorphisms) of 3000 commercial broilers with Campylobacter load in their caeca. Trait heritability was modest but significant (h² = 0.11 ± 0.03). Results confirmed quantitative trait loci (QTL) on chromosomes 14 and 16 previously identified in inbred chicken lines, and detected two additional QTLs on chromosomes 19 and 26. RNA-Seq analysis of broilers at the extremes of colonisation phenotype identified differentially transcribed genes within the QTL on chromosome 16 and proximal to the major histocompatibility complex (MHC) locus. We identified strong cis-QTLs located within MHC suggesting the presence of cis-acting variation in MHC class I and II and BG genes. Pathway and network analyses implicated cooperative functional pathways and networks in colonisation, including those related to antigen presentation, innate and adaptive immune responses, calcium, and renin–angiotensin signalling. While co-selection for enhanced resistance and other breeding goals is feasible, the frequency of resistance-associated alleles was high in the population studied and non-genetic factors significantly influenced Campylobacter colonisation.

Transcriptional response of blood leukocytes from turkeys challenged with Salmonella enterica serovar Typhimurium UK1

Non-typhoidal Salmonella is one of the most common causes of bacterial foodborne disease and consumption of contaminated poultry products, including turkey, is one source of exposure. Minimizing Salmonella colonization of commercial turkeys could decrease the incidence of Salmonella-associated human foodborne illness. Understanding host responses to these bacteria is critical in developing strategies to minimize colonization and reduce food safety risk. In this study, we evaluated bacterial load and blood leukocyte transcriptomic responses of 3-week-old turkeys challenged with the Salmonella enterica serovar Typhimurium (S. Typhimurium) UK1 strain. Turkeys (n = 8/dose) were inoculated by oral gavage with 10⁸ or 10¹⁰ colony forming units (CFU) of S. Typhimurium UK1, and fecal shedding and tissue colonization were measured across multiple days post-inoculation (dpi). Fecal shedding was 1-2 log₁₀ higher in the 10¹⁰ CFU group than the 10⁸ CFU group, but both doses effectively colonized the crop, spleen, ileum, cecum, colon, bursa of Fabricius and cloaca without causing any detectable clinical signs in either group of birds. Blood leukocytes were isolated from a subset of the birds (n = 3-4/dpi) both pre-inoculation (0 dpi) and 2 dpi with 10¹⁰ CFU and their transcriptomic responses assayed by RNA-sequencing (RNA-seq). At 2 dpi, 647 genes had significant differential expression (DE), including large increases in expression of immune genes such as CCAH221, IL4I1, LYZ, IL13RA2, IL22RA2, and ACOD1. IL1β was predicted as a major regulator of DE in the leukocytes, which was predicted to activate cell migration, phagocytosis and proliferation, and to impact the STAT3 and toll-like receptor pathways. These analyses revealed genes and pathways by which turkey blood leukocytes responded to the pathogen and can provide potential targets for developing intervention strategies or diagnostic assays to mitigate S. Typhimurium colonization in turkeys.

Transcriptome profiling analysis of caeca in chicks challenged with Salmonella Typhimurium reveals differential expression of genes involved in host mucosal immune response

Temporal regulation of global gene expression in the caeca of chickens infected with Salmonella Typhimurium has not been investigated previously. In this study, we performed the transcriptome analysis of the caeca of Salmonella Typhimurium challenged chicks to understand the regulation of the mucosal immune system in a temporal manner. The Salmonella infection resulted in the activation of the caecal immune system by the upregulation of the differentially expressed genes (DEGs; false discovery rate (FDR) < 0.05; log₂ fold change > 1) involved in biological pathways such as Toll-like receptor signaling pathway, Salmonella infection, cytokine-cytokine receptor interaction, phagosome, apoptosis and intestinal immune network for IgA production. The activation of biological pathways such as RIG-I-like receptor signaling pathway, ErbB signaling pathway and cellular senescence showed a time-dependent response of the host immune system. A 49% increase in the DEGs on day 7 compared with day 3 post-infection (p.i.) suggested a time-dependent role of multiple genes such as AvBD1, AvBD2, AvBD7, IL2, IL10, IL21, SIVA1, CD5, CD14 and GPR142 in the regulation of the immune system. Nested network analysis of the individual biological pathways showed that IL6 played a significant role in the immune system regulation by activating the pathways, including Toll-like receptor signaling pathway, Salmonella infection, intestinal immune network for IgA production and C-type lectin receptor signaling pathway. The downregulated DEGs (FDR < 0.05; log₂ fold change < -1) showed that Salmonella challenge affected the functions of pathways, such as tryptophan metabolism, retinol metabolism, folate biosynthesis and pentose and glucoronate interconversions, suggesting the disruption of cellular mechanisms involved in nutrient synthesis, absorption and metabolism. Overall, the immune response was temporally regulated through the activation of Toll-like signaling receptor pathway, cytokine-cytokine interactions and Salmonella infection, where IL6 played a significant role in the modulation of caecal immune system against Salmonella Typhimurium. KEY POINTS: • The immune response to Salmonella Typhimurium challenge was temporally regulated in the caeca of chickens. • Many newly identified genes have been shown to be involved in the activation of the immune system. • Toll-like receptors and interleukins played a key role in immune system regulation.

Genome-Wide Epistatic Interaction Networks Affecting Feed Efficiency in Duroc and Landrace Pigs

Interactions among genomic loci have often been overlooked in genome-wide association studies, revealing the combinatorial effects of variants on phenotype or disease manifestation. Unexplained genetic variance, interactions among causal genes of small effects, and biological pathways could be identified using a network biology approach. The main objective of this study was to determine the genome-wide epistatic variants affecting feed efficiency traits [feed conversion ratio (FCR) and residual feed intake (RFI)] based on weighted interaction SNP hub (WISH-R) method. Herein, we detected highly interconnected epistatic SNP modules, pathways, and potential biomarkers for the FCR and RFI in Duroc and Landrace purebreds considering the whole population, and separately for low and high feed efficient groups. Highly interacting SNP modules in Duroc (1,247 SNPs) and Landrace (1,215 SNPs) across the population and for low feed efficient (Duroc-80 SNPs, Landrace-146 SNPs) and high feed efficient group (Duroc-198 SNPs, Landrace-232 SNPs) for FCR and RFI were identified. Gene and pathway analyses identified ABL1, MAP3K4, MAP3K5, SEMA6A, KITLG, and KAT2B from chromosomes 1, 2, 5, and 13 underlying ErbB, Ras, Rap1, thyroid hormone, axon guidance pathways in Duroc. GABBR2, GNA12, and PRKCG genes from chromosomes 1, 3, and 6 pointed towards thyroid hormone, cGMP-PKG and cAMP pathways in Landrace. From Duroc low feed efficient group, the TPK1 gene was found involved with thiamine metabolism, whereas PARD6G, DLG2, CRB1 were involved with the hippo signaling pathway in high feed efficient group. PLOD1 and SETD7 genes were involved with lysine degradation in low feed efficient group in Landrace, while high feed efficient group pointed to genes underpinning valine, leucine, isoleucine degradation, and fatty acid elongation. Some SNPs and genes identified are known for their association with feed efficiency, others are novel and potentially provide new avenues for further research. Further validation of epistatic SNPs and genes identified here in a larger cohort would help to establish a framework for modelling epistatic variance in future methods of genomic prediction, increasing the accuracy of estimated genetic merit for FE and helping the pig breeding industry.

Association of SPOP Expression with the Immune Response to Salmonella Infection in Chickens

Simple Summary

Salmonella frequently causes human illness via the consumption of contaminated meat or eggs. At present, studies about how the host immune response against Salmonella is regulated are limited. Speckle-type POZ (poxvirus and zinc finger) protein (SPOP) is a specific adaptor of Cul3-based ubiquitin ligase, which catalyzes the ubiquitination and degrades the substrates. However, its role in the immune response is still unknown. Therefore, this study measured expression of SPOP and the proinflammatory cytokines interleukin-1β and interleukin-8 in chicken macrophage cells stimulated with a bacterial substitute and assessed their relationship using the quantitative polymerase chain reaction. We then validated the results in chickens infected with Salmonella. Notably, SPOP expression gradually decreased and then gradually increased in cells after challenging the bacterial substitute, indicating its potential involvement in the regulation of the immune response. Additionally, SPOP expression was negatively correlated with expression of interleukin 1β and interleukin-8 both in vivo and in vitro. More importantly, SPOP expression was related to immunoglobulin (Ig) A production and bacterial loads in chickens infected with Salmonella. These results indicate that SPOP could be a potential marker of the immune response in chickens.

Abstract

Salmonellosis is a zoonosis that is not only harmful to the health of poultry but also poses a threat to human health. Although many measures have been put in place to reduce morbidity, they have not provided satisfactory results. Therefore, it is necessary to clarify the immune mechanisms involved in improving the resistance of chickens against Salmonella. BTB (Broad-complex Tramtrack and Bric-a-brac) Speckle-type POZ (poxvirus and zinc finger) protein (SPOP) regulates protein expression by promoting substrate ubiquitination and degradation. The correlation between SPOP expression and the immune response has not been fully described. Therefore, the aim of this study was to clarify this relationship. In vitro, we stimulated chicken macrophage cells (HD11) with lipopolysaccharide, then analyzed the correlation between SPOP and IL1β or IL8 expression using quantitative real-time polymerase chain reaction (qRT-PCR). In vivo, we infected 7-days-old chickens with Salmonella Typhimurium, then analyzed the association between SPOP expression and the immune response, including IL1β and IL8 expression, IgA production, and bacterial loads. We found that SPOP may participate in the regulation of the immune response in macrophage cells. SPOP expression was negatively correlated with IL-1β and IL-8 expression both in vivo and in vitro. SPOP expression was also negatively related to bacterial loads and immunoglobulin (Ig) A production. These results indicate that SPOP may have important functions in the response to Salmonella infection.

Translating 'big data': better understanding of host-pathogen interactions to control bacterial foodborne pathogens in poultry

Recent technological advances has led to the generation, storage, and sharing of colossal sets of information ('big data'), and the expansion of 'omics' in science. To date, genomics/metagenomics, transcriptomics, proteomics, and metabolomics are arguably the most ground breaking approaches in food and public safety. Here we review some of the recent studies of foodborne pathogens (Campylobacter spp., Salmonella spp., and Escherichia coli) in poultry using big data. Genomic/metagenomic approaches have reveal the importance of the gut microbiota in health and disease. They have also been used to identify, monitor, and understand the epidemiology of antibiotic-resistance mechanisms and provide concrete evidence about the role of poultry in human infections. Transcriptomics studies have increased our understanding of the pathophysiology and immunopathology of foodborne pathogens in poultry and have led to the identification of host-resistance mechanisms. Proteomic/metabolomic approaches have aided in identifying biomarkers and the rapid detection of low levels of foodborne pathogens. Overall, 'omics' approaches complement each other and may provide, at least in part, a solution to our current food-safety issues by facilitating the development of new rapid diagnostics, therapeutic drugs, and vaccines to control foodborne pathogens in poultry. However, at this time most 'omics' approaches still remain underutilized due to their high cost and the high level of technical skills required.

Transcriptome Analysis of the Cecal Tonsil of Jingxing Yellow Chickens Revealed the Mechanism of Differential Resistance to Salmonella

Salmonella is one of the most common food-borne pathogens. It can be transmitted between chickens, as well as to people by contaminated poultry products. In our study, we distinguished chickens with different resistances mainly based on bacterial loads. We compared the cecal tonsil transcriptomes between the susceptible and resistant chickens after Salmonella infection, aiming to identify the crucial genes participating in the antibacterial activity in the cecal tonsil. A total of 3214 differentially expressed genes (DEGs), including 2092 upregulated and 1122 downregulated genes, were identified between the two groups (fold change ≥ 2.0, padj < 0.05). Many DEGs were mainly involved in the regulation of two biological processes: crosstalk between the cecal tonsil epithelium and pathogenic bacteria, such as focal adhesion, extracellular-matrix-receptor interaction, and regulation of the actin cytoskeleton and host immune response including the cytokine-receptor interaction. In particular, the challenged resistant birds exhibited strong activation of the intestinal immune network for IgA production, which perhaps contributed to the resistance to Salmonella infection. These findings give insight into the mRNA profile of the cecal tonsil between the two groups after initial Salmonella stimulation, which may extend the known complexity of molecular mechanisms in chicken immune response to Salmonella.

A new chicken 55K SNP genotyping array

BACKGROUND

China has the richest local chicken breeding resources in the world and is the world's second largest producer of meat-type chickens. Development of a moderate-density SNP array for genetic analysis of chickens and breeding of meat-type chickens taking utility of those resources is urgently needed for conventional farms, breeding industry, and research areas.

RESULTS

Eight representative local breeds or commercial broiler lines with 3 pools of 48 individuals within each breed/line were sequenced and supplied the major SNPs resource. There were 7.09 million - 9.41 million SNPs detected in each breed/line. After filtering using multiple criteria such as preferred incorporation of trait-related SNPs and uniformity of distribution across the genome, 52.18 K SNPs were selected in the final array. It consists of: (i) 19.22 K SNPs from the genomes of yellow-feathered, cyan-shank partridge and white-feathered chickens; (ii) 5.98 K SNPs related to economic traits from the Illumina 60 K SNP Bead Chip, which were found as significant associated SNPs with 15 traits in a Beijing-You crossed Cobb F2 resource population by genome-wide association study analysis; (iii) 7.63 K SNPs from 861 candidate genes of economic traits; (iv) the 0.94 K SNPs related to residual feed intake; and (v) 18.41 K from chicken SNPdb. The polymorphisms of 9 extra local breeds and 3 commercial lines were examined with this array, and 40 K - 47 K SNPs were polymorphic (with minor allele frequency > 0.05) in those breeds. The MDS result showed that those breeds can be clearly distinguished by this newly developed genotyping array.

CONCLUSIONS

We successfully developed a 55K genotyping array by using SNPs segregated from typical local breeds and commercial lines. Compared to the existing Affy 600 K and Illumina 60 K arrays, there were 21,41 K new SNPs included on our Affy 55K array. The results of the 55K genotyping data can therefore be imputed to high-density SNPs genotyping data. The array offers a wide range of potential applications such as genomic selection breeding, GWAS of interested traits, and investigation of diversity of different chicken breeds.