Skip to main content
  • Research article
  • Open access
  • Published:

Meta-analysis of Chicken – Salmonella infection experiments

Abstract

Background

Chicken meat and eggs can be a source of human zoonotic pathogens, especially Salmonella species. These food items contain a potential hazard for humans. Chickens lines differ in susceptibility for Salmonella and can harbor Salmonella pathogens without showing clinical signs of illness. Many investigations including genomic studies have examined the mechanisms how chickens react to infection. Apart from the innate immune response, many physiological mechanisms and pathways are reported to be involved in the chicken host response to Salmonella infection. The objective of this study was to perform a meta-analysis of diverse experiments to identify general and host specific mechanisms to the Salmonella challenge.

Results

Diverse chicken lines differing in susceptibility to Salmonella infection were challenged with different Salmonella serovars at several time points. Various tissues were sampled at different time points post-infection, and resulting host transcriptional differences investigated using different microarray platforms. The meta-analysis was performed with the R-package metaMA to create lists of differentially regulated genes. These gene lists showed many similarities for different chicken breeds and tissues, and also for different Salmonella serovars measured at different times post infection. Functional biological analysis of these differentially expressed gene lists revealed several common mechanisms for the chicken host response to Salmonella infection. The meta-analysis-specific genes (i.e. genes found differentially expressed only in the meta-analysis) confirmed and expanded the biological functional mechanisms.

Conclusions

The meta-analysis combination of heterogeneous expression profiling data provided useful insights into the common metabolic pathways and functions of different chicken lines infected with different Salmonella serovars.

Background

Chicken meat and eggs for human consumption can be contaminated with several Salmonella species, and therefore chicken-derived food products can be regarded as a source of human zoonotic pathogens. Although proper food preparation should kill the pathogens, the food items contain a potential hazard for humans. In chicken both acute fatal and chronic Salmonellosis occurs depending upon the infecting Salmonella serovar [14]. Broad host range Salmonella serovars used most often in studies – including the studies used for this meta-analysis, S. Typhimurium and S. Enteritidis, do not cause fatal infections when chickens older than one day post hatch are orally challenged. Chickens can harbor Salmonella pathogen without showing clinical signs of illness [3, 5]. Many investigations have examined the mechanisms how chickens react to infection, the mechanism of transfer to humans and host immunity to infection [3, 6].

Diverse host species may react differently to Salmonella infection [7]. While one-day old chickens may succumb to broad host range Salmonella infection, older chickens often show no clinical signs. Furthermore, specific chicken lines have been shown to differ in their susceptibility for Salmonella [811]. These clear genetic differences in susceptibility may be due to pleiotropic effects, or to unknown selection-related mechanisms. In the last decade, gene expression profiling studies using microarrays have been widespread in animal genomics and have enabled researchers to monitor the effects of pathogens on host cells and tissues with the aim of gaining insights into the molecular mechanisms that are involved in the host-pathogen interactions. Several genes involved in Salmonella susceptibility in chicken have been determined [1217]. Apart from the innate immune response, many physiological mechanisms and pathways were reported to be involved in the chicken host response to Salmonella infection which are also active in uninfected cells, including energy metabolism, cell shape, and others [1820].

Each of these independent experiments showed how individual hosts within the specific experimental conditions reacted to Salmonella infection. Meta-analysis of these experiment may reveal a common genetic background for the chicken host reaction to the Salmonella infection. Furthermore, the age-related differences in the mechanisms and the outcome of the host immune-response to Salmonella infection suggests that different immune-reactions are possible, and are likely to be age related [21, 22]. Taken together this indicates a complex interplay between chicken host genetics and Salmonella serovars [3, 10, 11, 18, 23].

Meta-analysis methods integrate results of independent studies creating very large datasets with increased statistical power [24, 25]. It allows a more objective appraisal of evidence than individual studies, and has been widely used to interpret contradictory results from diverse studies. Furthermore, this analysis method overcomes the problem of reduced statistical power associated with studies of small sample size (reviewed by [26, 27]. Such methods enable analyses at a higher level than possible on the individual datasets. Host-specific general mechanisms can be determined in addition to mechanisms operating under specific conditions. Thus, using previously published individual datasets we were able to highlight new results that contribute to understanding of common disease mechanisms and physiology. Different experiments were performed under the umbrella of a large EU-funded project called SABRE - Cutting Edge Genomics for Sustainable Animal Breeding [9, 21, 22, 28, 29]. This meta-analysis brings the individual studies together offering the potential to highlight new host-pathogen interaction mechanisms and elucidate possible general host-response mechanisms. The objective of this study was to determine the general chicken host response to Salmonella infection independent of age of the chicken, age at infection and, time post infection and independent of host response time post-infection. The results indicate several common chicken host reaction mechanisms to Salmonella infection.

Methods

Animals and Salmonella challenges

Experiment 1

The original animal experiment was described by Fife et al. [29]. In short, two inbred chicken lines differing in susceptibility to gut pathogens (lines N and 6, with line 6 more resistant than line N, [29]) were at three weeks of age orally infected with 5.1x107-1.97x108 cfu S. Typhimurium according to the method of Barrow et al. [30]. The caecal tonsils and spleens were sampled at 2, 3, and 4 days post infection, (n = 10), and four birds at each time point were used as uninfected controls. Total RNA for these samples was isolated and used for hybridization to the 20.6 K chicken oligo array (ARK genomics; http://www.ark-genomics.org/) microarrays. Infection and infection clearance was determined by ceacal counts of S. Typhimurium (cfu 106l) and differences between the lines investigated. A total of 32 microarrays per line were obtained.

Experiment 2

The original experiment was described by Schokker et al. [31] (GEO data: GSE27069). In short, three commercial chicken lines differing for Salmonella sensitivity were orally infected with 105 cfu S. Enteritidis at the day of hatch. The jejunum was sampled at 8 h and days 1 and 2 post infection, 10 animals each, of which 5 were used for microarray analysis. A reference pool was created from 0.33, 1 and 2 days post infection birds, for all three lines together, as well as control and infected birds. Infection was checked by body weight and liver weight gain and liver clearance, and cloaca swaps. Total RNA was isolated and hybridized to the same microarrays as experiment 1. A total of 45 microarrays were obtained [31].

Experiment 3

The original experiment was described by Schokker et al. [16, 28] (ArrayExpress data: E-MEXP-042). In short, chickens were challenged orally at the day of hatch with 105 cfu S. Enteritidis. The jejunum was sampled at 8 h, and 1, 2, 4, 8, 12, and 21 days post infection, 5 animals each for both control and infected situation. Infection was checked by body weight and liver weight gain and liver clearance. Total RNA was isolated and single color hybridized against Agilent chicken microarrays. A total of 70 microarrays were obtained [16, 21].

Experiment 4

The original experiment was described by van Hemert et al. [21] (GEO data: GSE3702). In short, two chicken lines differing in growth rate and Salmonella sensitivity were orally infected with 105 cfu of S. Enteritidis at one day of age and jejunum samples were taken after 24 h. A non-infected control was used, 5 chicken each. Total RNA was isolated and hybridized against Affymetrix chicken microarrays using group comparison. A total of four microarrays were obtained [21].

Meta-analysis methodology

Pre-processing microarray data

The microarray data pre-processing was carried out using functions from the LIMMA package (version 3.2.1) [32]. The quality of the arrays was evaluated through several diagnostic plots. The “normexp” method [33] was used for background correction, followed by normalization within individual microarrays using the default “print tip loess” method and normalization between arrays using the “quantile” method. The background correction was set to: offset = 50. The offset can be used to add a constant to the intensities before log-transforming, so that the log-ratios are shrunk towards zero at the lower intensities. This may eliminate or reverse the usual 'fanning' of log-ratios at low intensities associated with local background subtraction. Areas with higher than average background were removed from the results. Bad hybridization always removed whole microarrays. Especially in experiment 1 this removed parts of the results. After normalization 20 slides of experiment 1 and one slide of Experiment 2 were deleted due to poor quality hybridization. This will inevitably affect the results, but this procedure ensures that only good quality data were used.

Meta-analysis

Meta-analysis was carried out using the directpvalcombi function from the metaMA package (Meta-analysis for MicroArrays) (version 1.1) in R [24]. The input for the meta-analysis were the individual microarrays of all experiments. The meta-analysis produced lists of gene names with differential expression under specific conditions. The lists of genes were grouped in (1) DE: the list of Differentially Expressed genes (i.e. in the experiments and in the meta-analysis), and (2) IDD (Integration Driven Discoveries): the list of genes that were determined differentially expressed in the meta-analysis that were not identified in any of the individual studies alone (i.e. new differentially expressed genes). For both DE and IDD gene lists (e) experiment (i.e. 4 studies) and (t) time (14 studies) were generated. Subsequently from these (e) and (t) the following groups were also extracted, namely (et): overlap between the (e) and (t) groups, (e-t): genes unique in (e), and (t-e): genes unique in (t) (Figure 1). The groups included results from the different tissues. Since the analyses focus on expression differences related to Salmonella infection no interaction with tissue-specific gene expression can be expected.

Figure 1
figure 1

Overview experiment and time analyses. For both IDD and DE experiment (e) and time (t) gene lists were generated. Thereafter different subsets were extracted, unique time genes (t-e), unique experiment genes (e-t), and overlap between (e) and (t) (et).

Functional bioinformatics analyses

The lists of differentially expressed genes were analyzed for biological functionalities using the DAVID (The Database for Annotation, Visualization and Integrated Discovery) software [3436], version 6.6. The gene lists were analyzed against the gene list of the human genome since the annotation of the human genome, especially with physiological data, is more advanced than the chicken genome. Therefore, all genes were converted to human Entrez identifiers. The false discovery rate, multiple testing correction for statistical significance [37], and the fold enrichment analyses were manually included for all analyses. The tissue-specific profiles and functional annotations and clusterings of the gene lists were investigated.

Results

Meta-analysis

Due to the platform differences used in the individual studies the number of genes available in all studies was reduced as expressed in Figure 2. The Figure shows that 7,643 genes were common to all microarray platforms and thus available for meta-analysis over all studies. Different platforms may use different probes for the same genes, and the probes may differ in hybridization characteristics. However, since differential expression of genes was measured within a platform for each of the experiments this will not affect the meta-analysis.

Figure 2
figure 2

Visualization of the overlap between the microarray platforms. For the three different platforms, ARK-genomics, Affymetrix and Agilent, all probes were mapped to human entrez gene identifiers (EGIDs). Subsequently the overlap between all platforms was calculated.

The results of the meta-analysis DE-group were expressed in several lists of genes with regulated expression in more than one or all experiments (Table 1 – the Table with the gene lists is Additional file 1). The Table shows that approximately 3,000 genes are differently expressed, irrespective of age of infection or Salmonella serovar of infection. Differential expression was between control and infected animals of the same age. Furthermore, most differently expressed genes were the same in all experiments irrespective of sampling time point post-infection (Figure 3). The Figure shows that the large majority of the genes are shared by the (e) and the (t) categories.

Table 1 Number of differentially expressed genes per group
Figure 3
figure 3

Visualization of the differentially expressed genes of the DE-group. The (e), (t), and (e + t) categories show large overlap between the categories. The lists of genes were grouped in (1) DE: the Number of Differentially Expressed genes (i.e. in the experiments and in the meta-analysis), and (2) IDD (Integration Driven Discoveries): the number of genes that were determined differentially expressed in the meta-analysis that were not identified in any of the individual studies alone (i.e. new differentially expressed genes). These two groups were studied in detail in five categories each: (e): experiment - i.e. 4 studies, (t): time after infection – 14 different time point post-infection, (et): overlap between the (e) and (t) groups, (e-t): genes unique in (e), and (t-e): genes unique in (t).

The IDD group (i.e. the list of genes found only differently expressed in the meta-analysis) genes were predominantly in the (t) category (Figure 4). A few genes were differently expressed in the meta-analysis only for the (e) category.

Figure 4
figure 4

Visualization of the differentially expressed genes of the IDD-group. The (e), (t), and (e + t) categories showing overlap between the categories. The lists of genes were grouped in (1) DE: the Number of Differentially Expressed genes (i.e. in the experiments and in the meta-analysis), and (2) IDD (Integration Driven Discoveries): the number of genes that were determined differentially expressed in the meta-analysis that were not identified in any of the individual studies alone (i.e. new differentially expressed genes). These two groups were studied in detail in five categories each: (e): experiment - i.e. 4 studies, (t): time after infection – 14 different time point post-infection, (et): overlap between the (e) and (t) groups, (e-t): genes unique in (e), and (t-e): genes unique in (t).

Functional bioinformatics analyses

The gene lists were analyzed for biological functional groups using the DAVID software. First the expression profile of differentially expressed genes of the DE-group was compared with normal tissue-specific expression profiles from the same time points and the same tissue type (Table 2). Significant results were only obtained for the DE group, (e), (t), and (et) categories. The Table indicates that the sampled intestinal tissue showed expression profiles related to a number of different tissue types. The most significant tissue expression profile is epithelium. This cell type is abundantly present in intestinal tissue. Several other cell types and tissues also showed similarities for tissue-specific expression profiles, some of them not relevant in intestine (data not shown).

Table 2 Tissue specificity of differentially expressed gene profiles
Table 3 Clustering of gene lists using functional annotations

The results of the differently expressed gene lists were then analyzed for functional biological mechanisms. The results are shown in Additional file 2. The results showed lists of biological functions for the (e) and the (t) categories of the DE group. Moreover, the (et) category showed that both lists were largely similar, and the top of the lists were even identical. The top of the lists indicated that phosphorylation of proteins, acetylation in the cytoplasm and lumen of other cellular components, and ATP consuming processes were important biological mechanisms during chicken host reaction to Salmonella infection. The meta-analysis showed additional significant results for both the experiment (e-t) and the time (t-e) categories of the phosphoprotein biological function during chicken host reaction to Salmonella infection. Similarly, in the IDD group the time (t-e) category further indicated additional significant results especially for the phosphoprotein biological function.

Finally, a cluster analysis was performed for the lists of biological functional annotations. The DE-group (e) and (et) categories showed over 600 clusters. Due to the fact that the list of the (t) category in the DE group was longer than 3000 entries, clustering was technically not possible for the DAVID software. Since the lists of functional annotations of the (et) and the (t) categories of the DE group were very similar the (t) category clusters were deduced from the (e) and (et) categories clusters (see below).

Table 3 shows the results for clusters with enrichment scores larger than 2. An enrichment score indicate whether the number of genes in a cluster is equal to the expected number of genes (due to the number of genes of that physiological group in the genome and on the microarray) or higher or lower than expected. A high enrichment score thus indicates that the physiological trait of the cluster may be significant to the trait. The clusters were ordered by enrichment scores. The range of enrichment scores was from over 10 to almost zero. Enrichment scores less than 2 were omitted, leaving between 35 and 40 clusters in the DE-group (e) and (et) categories, respectively. For completeness, the Additional file 3 includes the biological functions per cluster in detail, i.e. larger than 1. The content of each cluster is a group of biological functions taken together from various databases centered on a specific theme. For example (see Table 3), cluster 1 groups biological functions together related to the lumen of cell organelles, especially relating to the nucleus. Cluster 2 groups nucleotide binding and phosphorylation functions (especially serine and threonine phosphorylation), while cluster 3 groups mitochondrial membrane functionalities. The Table contains both cell / tissue morphological clusters and (macro) molecular biogenesis functional clusters. The functions of several protein domains were also relevantly clustered.

Taken together functional clustering analysis showed that these functional annotations can be grouped together in higher order biological morphological structures and biological processes. The clusters can be divided in 21 clusters describing metabolic processes – of which three were related to energy metabolism, eight clusters describing (cell) morphological features, three clusters specifically pointing to protein domains involved, one cluster related to apoptosis and one cluster is a collection of processes, making it difficult to recognize a central theme. Apart from apoptosis, these clusters describe normal cellular physiological processes taking also place in non-infected animals, e.g. during growth and development of the tissues and organs. Nevertheless, these processes also participate in the host reaction to infection with Salmonella.

Apart from small differences in the order of clusters the (e) and (et) categories of the DE group differ only in a few clusters from each other. A specific tyrosine phosphorylation was found in the (et) category but not in the (e) category of the DE group while the (e) category showed a protein domain WD cluster and a cell movement cluster, both not found in the (et) category. Finally, it should be noted that in none of the other categories (DE and IDD groups) a significant cluster with enrichment score of at least 2 could be found.

Discussion

The objective of this study was to determine the general chicken host response to Salmonella infection independent of age of the chicken host at time of bacterial challenge and independent of host response time post-infection, investigating various tissues and using chicken lines differing in susceptibility for Salmonella. The results highlight several biological mechanisms related to energy metabolism, apoptosis, specific protein domains indicating groups of involved proteins, and several cellular morphological structures where the affected processes are taking place. Overall, the reported meta-analysis approach showed successful integration of heterogeneous data sets of limited size by increasing statistical power. Using the results of this study for future biomarker analysis may provide in early diagnosis and warning of potentially hazardous food.

Meta-analysis using data from different sources and different technologies

A meta-analysis is performed using the original raw data from a number of individual experiments. Since the experiments may have different objectives and use different technologies, the experiments or data may not be directly comparable. In our study we compared data from four studies using: (1) genetically different chicken lines differing in Salmonella susceptibility), (2) different Salmonella serovars, (3) different sampling time points, (4) different sampled tissues, (5) different microarray types, and (6) different ages of bacterial challenge. Intuitively, it would be expected that these differences would affect the meta-analysis: (1) Genetically different lines of chicken, differing in Salmonella susceptibility, were expected to differ in reaction mechanism and/or reaction severity. (2) A pathogen specific host reaction was expected to different Salmonella serovars. (3) Sampling at different times post infection was suggested to show different temporal expression patterns related to the stage of infection. (4) Expression patterns are also expected to differ between different tissues or cell types. (5) Finally, different microarrays contained different sets of genes, so results from one study were expected to be missing from another study and vice versa. (6) The age of challenge of the birds would be expected to produce very different responses due to the poorly developed immune system of day old chicks compared to 3 week old birds. Despite of all these differences our meta-analysis indicated that the chicken lines react to Salmonella infection through comparable mechanisms irrespective of Salmonella serovar and tissue type, and therefore it may be concluded that we identified common mechanisms of the host response to the bacterial challenge. However, due to the different experimental ages of the animals used in the diverse studies, this conclusion may be hampered by the developmental differences of tissues and organs in the animals. It can be expected that at least part of the mechanisms found may relate to this. This could have been investigated only if control samples of all experimental ages in the individual datasets would have been available. But often these control samples are only available for the last experimental sampling age. Further experiments are needed to elucidate this point.

Although it is not certain, it can be expected that the results would have been more comprehensive if all experiments were performed under standard procedures. Similarly, the functional annotation analysis to elucidate potential biological mechanisms of the functional reaction of chicken to Salmonella infection would have been more robust.

What does the functional annotation analysis teach us about the chicken host reaction to the infection with a Salmonella bacterium?

The first indication that the chicken host reaction to Salmonella infection was similar between the diverse experiments was obtained from the similarities in the gene lists for the differently conducted experiment (e) and time (t) categories, i.e. the (et) category. One unexpected finding was that the expression profiles related to several different cell types. Intestinal tissue is composed of many different cell types that could be indicated by a mixture of expression profiles. Furthermore, localized infection will change tissue expression profiles, which will be exacerbated by the influx of immune cells, which will further change the overall expression profile. However, the results indicated similarities to the expression profiles of several cell types including many unrelated tissues like liver and brain. The epithelium cell type of the intestine was the highest ranking tissue in all three analysis groups. Also platelet and muscle tissue, and perhaps colon expression profiles were recognizable – these cell types are also included in the intestinal tissue. Other cell types may also be on the list for several reasons. One reason may be that a cell type has a high turnover rate like epithelial cells in intestine tissue. In these cell types the general mechanism for cell division will be activated and therefore all these cell types appear on the list. Especially in developing young-age animals this may be expected. Finally, cell types and tissues may have been included in the list because we used the human physiological information instead of chicken physiological information for the DAVID software to create the list, e.g. lymph tissue may be inserted for that reason (although the chicken intestine contains a limited number of Peyer’s patches as lymphoid tissues [38]).

The functional annotation is the result of the analysis of the DAVID software using the same gene lists to analyze several different databases containing biological function information. Due to the similarities within the gene lists the DE group (e), (t), and (et) categories showed similar functional annotations. Furthermore, the top category functional annotation “phosphoprotein” was also found in the differently expressed genes unique for both the experiments and the time (t) array after infection, and in the time-related genes found specifically in the meta-analysis (IDD-group). These results suggest two mechanisms: (1) the analysis is robust and indicated the same biological functionalities for all experiments despite the experimental differences, and (2) the meta-analysis adds new genes and data to the already existing data, but does indicate new biological mechanisms for the reaction of chicken hosts to the Salmonella infection. On the other hand, protein phosphorylation is an important regulatory mechanism for protein function in normal tissue and changed phosphoprotein content of the cell may have important physiological consequences for cellular metabolism (see below).

The clustering of biological functional annotations showed only in the lower part of the list differences between the DE (e) and DE (t) category. While these differences themselves were statistically significant, the place on the list may suggest that the differences in the reaction of chicken to Salmonella are small. Alternatively, these differences point towards differences in the expression profiles related to time point after infection. However, due to the structure of the dataset these differences may also relate to deviations in the general chicken reaction mechanism caused by different chicken breed/lines, different tissues or different Salmonella serovars (jejunum vs. caecum; S. Enteritidis vs. S. Typhimurium).

While most clusters of differently expressed genes were similar in both the experiment (e) and the time (t) categories, some interesting differences were obtained. It should be noted that clusters found in one category but not in the other may be the result of real missing clusters or clusters failing to reach the enrichment score limit in one of the two categories. Three clusters were found in the differently expressed genes group experiment (e) category, but not in the time (t) category suggesting that these genes were not, or less regulated in time after infection and may be constitutively active during the chicken host reaction to Salmonella infection: (E1) Tyrosine phosphorylation, (E2) Protein domain WD (repeat), and (E3) cell movement. Two clusters were found in the differently expressed genes group time (t) category but not in the experiment (e) category, suggesting that these genes were especially regulated at different moments in time after infection of the chicken: (T1) ubiquitin / proteasome mediated proteolysis, and (T2) Intracellular vesicles.

Tyrosine phosphorylation (E1)

Phosphorylation activates or deactivates many proteins in cellular processes and protein phosphorylation in particular plays a significant role in a wide range of cellular processes [3941]. Tyrosine phosphorylation is considered to be one of the key steps in signal transduction and regulation of enzymatic activity (for a review see [42]. The consequences of the difference between the (e)- and (t) categories (for tyrosine phosphorylation) may be important.

Both signal transduction and enzymatic activity may regulate a variety of important processes in the cell, including immune processes, cellular metabolism, and cell morphology, which may be related to the chicken host reaction to infection with Salmonella, e.g. via changes in the actin cytoskeleton [43, 44].

Protein domain WD (E2)

The WD40 repeat (also known as the WD or beta-transducin repeat) is a short structural motif of approximately 40 amino acids, often terminating in a tryptophan-aspartic acid (W-D) dipeptide [45]. Several of these repeats are combined to form a type of protein domain called the WD domain. WD-containing proteins have 4 to 16 repeating units, all of which are thought to form a circularized beta-propeller structure [46, 47]. WD-repeat proteins are a large protein family found in all eukaryotes and are implicated in a variety of functions ranging from signal transduction and transcription regulation to cell cycle control and apoptosis, which may be directly related to the chicken immune reaction to the Salmonella infection. All these specific functions were also found in other clusters. Thus, the difference between the (e) and (t) categories may induce modulations of the intensities of the processes described in several of the other clusters, thereby representing another mechanism for these proteins to modulate the chicken host response to Salmonella infection. Furthermore, the underlying common function of all WD-repeat proteins is coordinating multi-protein complex assemblies, where the repeating units serve as a rigid scaffold for protein interactions. The specificity of the proteins is determined by the sequences outside the repeats themselves. Several of the clusters relate to macromolecules which may be differently regulated between the (e) and (t) categories. A specific macromolecule includes the E3 ubiquitin ligase suggesting that also proteolysis is regulated [46, 47].

Cell movement (E3)

Cell movement could relate to the influx of immune cells to the site of infection/ tissue. Also in non-infected tissues immune cells move through the tissue, but this process will be enhanced during infection. It may be suggested that regulation of this process may be one of the fundamental mechanisms of the cellular immune response of the chicken host.

Ubiquitin / proteasome mediated proteolysis (T1)

The destination of Ubiquitin tagged proteins is the proteasome for proteolysis. The ubiquination system functions in a wide variety of cellular processes, including the immune response and inflammation, antigen processing, apoptosis and cell cycle. Furthermore, the development and degeneration of several tissues is affected – probably via biogenesis of organelles such as ribosomes and modulation of cell surface receptors, ion channels, and the secretory pathway (for a review see [48]. The ubiquination system is responsive to stress and extracellular modulators such as Salmonella infection [49]. It is clear that the wide variety of cellular metabolic functions regulated by the ubiquitin / proteasome system may affect the chicken host response to Salmonella. Its regulation of expression especially at different time points after infection can modulate the response of the chicken host to Salmonella infection through a variety of mechanisms described in the other clusters.

Intracellular vesicles (T2)

Intracellular vesicles transport material – e.g. (macro) molecules - through the cell – either importing or exporting material, or transporting material to different cellular locations. They deliver molecules both for excretion or to lysosomes for degradation, and may import food components for energy and cellular components synthesis processes. Lotz et al. [50] described that the HSP90 protein is important for the regulation of intracellular vesicle transport. The HSP90 protein is a molecular chaperone regulating the folding and thereby the activity of macromolecules [51]. These functions can be found in several of the other clusters. Thus, the differential expression especially at different time points after infection can modulate the response of the chicken host to Salmonella infection during the cause of the infection through a variety of mechanisms described in the other clusters.

Conclusions

These results shed light on the important biological mechanisms that are active in the chicken gut cells during Salmonella infection – although part of the processes may relate to growth and development of the tissues and organs as discussed above. From our data we conclude that similar host mechanisms apply to S. Enteritidis and S. Typhimurium infection, and that similar biological mechanisms appear underlying the processes regulated during different times after infection.

Abbreviations

DE:

The Number of Differentially Expressed genes

IDD = (Integration Driven Discoveries):

the number of genes that were determined differentially expressed in the meta-analysis that were not identified in any of the individual studies alone (i.e. new differentially expressed genes)

e:

Experiment - i.e. 4 studies

t:

Time after infection – 14 different time point post-infection

et:

Overlap between the (e) and (t) groups

e-t:

Genes unique in (e)

t-e:

Genes unique in (t).

References

  1. Hofstad MS, John BH, Calnek BW, Reid WN, Yoder HW: Diseases of Poultry. 1992, Panima Education Book Agency, New Delhi, Ind, 65-123. 8

    Google Scholar 

  2. Kabir LSM: Avian Colibacillosis and Salmonellosis: A Closer Look at Epidemiology, Pathogenesis, Diagnosis, Control and Public Health Concerns. Int J Environ Res Public Health. 2010, 7: 89-114. 10.3390/ijerph7010089.

    Article  Google Scholar 

  3. Wales AD, Davies RH: A critical review of Salmonella Typhimurium infection in laying Hens. Avian Pathol. 2011, 40: 429-436. 10.1080/03079457.2011.606799.

    Article  CAS  PubMed  Google Scholar 

  4. Chappell L, Kaiser P, Barrow P, Jones MA, Johnston C, Wigley P: The immunobiology of avian systemic salmonellosis. Vet Immunol Immunopathol. 2009, 128: 53-59. 10.1016/j.vetimm.2008.10.295.

    Article  CAS  PubMed  Google Scholar 

  5. Guard-Petter J: The chicken, the egg and Salmonella enteritidis. Environ Microbiol. 2001, 3: 421-430. 10.1046/j.1462-2920.2001.00213.x.

    Article  CAS  PubMed  Google Scholar 

  6. Parker MT: Enteric infections. Topley and Wilson's principles of bacteriology, virology and immunity. Edited by: Parker MT, Collier LH. 1990, Edward Arnold, London, 424-446.

    Google Scholar 

  7. Barrow PA, Huggins MB, Lovell MA: Host Specificity of Salmonella Infection in Chickens and Mice Is Expressed In Vivo Primarily at the Level of the Reticuloendothelial System. Infect Immun. 1994, 62: 4602-4610.

    PubMed Central  CAS  PubMed  Google Scholar 

  8. Bumstead N, Barrow P: Resistance to Salmonella gallinarium, S. pullorium, and S. enteritidis in inbred lines of chicken. Avian Dis. 1993, 37: 189-193. 10.2307/1591473.

    Article  CAS  PubMed  Google Scholar 

  9. Van Hemert S, Hoekman AJW, Smits MA, Rebel JMJ: Gene expression responses to a Salmonella infection in the chicken intestine differ between lines. Vet Immunol Immunopathol. 2006, 114: 247-258. 10.1016/j.vetimm.2006.08.007.

    Article  CAS  PubMed  Google Scholar 

  10. Calenge F, Kaiser P, Vignal A, Beaumont C: Genetic control of resistance to salmonellosis and toSalmonellacarrier-state in fowl: a review. Genet Selec Evol. 2010, 42: 11-10.1186/1297-9686-42-11.

    Article  Google Scholar 

  11. Redmond SB, Chuammitri P, Andreasen CB, Palić D, Lamont SJ: Proportion of circulating chicken heterophils and CXCLi2 expression in response to Salmonella enteritidis are affected by genetic line and immune modulating diet. Vet Immunol Immunopathol. 2011, 140: 323-328. 10.1016/j.vetimm.2011.01.006.

    Article  CAS  PubMed  Google Scholar 

  12. Hu J, Bumstead N, Barrow P, Sebastiani G, Olien L, Morgan K, Malo D: Resistance to Salmonellosis in the Chicken Is Linked to NRAMP1 and TNC. Genome Res. 1997, 7: 693-704.

    CAS  PubMed  Google Scholar 

  13. Leveque G, Forgetta V, Morroll S, Smith AL, Bumstead N, Barrow P, Loredo-Osti JC, Morgan K, Malo D: Allelic Variation in TLR4 Is Linked to Susceptibility to Salmonella enterica Serovar Typhimurium Infection in Chickens. Infect Immun. 2003, 71: 1116-1124. 10.1128/IAI.71.3.1116-1124.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Sadeyen J-R, Trotereau J, Velge P, Marly J, Beaumont C, Barrow PA, Bumstead N, Lalmanach A-C: Salmonella carrier state in chicken: comparison of expression of immune response genes between susceptible and resistant animals. Microbes Infect. 2004, 6: 1278-1286. 10.1016/j.micinf.2004.07.005.

    Article  CAS  PubMed  Google Scholar 

  15. Fife M, Salmon N, Hocking P, Kaiser P: Fine mapping of the chicken salmonellosis resistance locus (SAL1). Anim Genet. 2009, 40: 871-877. 10.1111/j.1365-2052.2009.01930.x.

    Article  CAS  PubMed  Google Scholar 

  16. Pan Z, Fang Q, Geng S, Kang X, Cong Q, Jiao X: Analysis of immune-related gene expression in chicken peripheral blood mononuclear cells following Salmonella enterica serovar Enteritidis infection in vitro. Res Vet Sci. 2012

    Google Scholar 

  17. Coble DJ, Redmond SB, Hale B, Lamont SJ: Distinct lines of chickens express different splenic cytokine profiles in response to Salmonella Enteritidis challenge. Poult Sci. 2011, 90: 1659-1663. 10.3382/ps.2010-01279.

    Article  CAS  PubMed  Google Scholar 

  18. Galan JE, Cossart P: Host-pathogen interactions: a diversity of themes, a variety of molecular machines. Curr Opin Microbiol. 2005, 8: 1-3. 10.1016/j.mib.2004.12.015.

    Article  PubMed  Google Scholar 

  19. Mattoo S, Lee YM, Dixon JE: Interactions of bacterial effector proteins with host proteins. Curr Opin Immunol. 2007, 19: 392-401. 10.1016/j.coi.2007.06.005.

    Article  CAS  PubMed  Google Scholar 

  20. Te Pas MFW, van Hemert S, Hulsegge I, Rebel JMJ, Smits MA: A pathways analysis tool for analyzing microarray data of species with low physiological information. Adv Bioinf. 2008, 10.1155/2008/719468. ID 719468

    Google Scholar 

  21. Van Hemert S, Hoekman AJW, Smits MA, Rebel JMJ: Early host gene expression responses to a Salmonella infection in the intestine of chickens with different genetic background examined with cDNA and oligonucleotide microarrays. Comp Biochem Physiol D. 2006, 1: 292-299.

    Google Scholar 

  22. Schokker D, Hoekman AJW, Smits MA, Rebel JMJ: Gene expression patterns associated with chicken jejunal development. Dev Comp Immunol. 2009, 33: 1156-1164. 10.1016/j.dci.2009.06.002.

    Article  CAS  PubMed  Google Scholar 

  23. Crhanova M, Hradecka H, Faldynova M, Matulova M, Havlickova H, Sisak F, Rychlik I: Immune response of chicken gut to natural colonization by gut microflora and to Salmonella enterica serovar enteritidis infection. Infect Immun. 2011, 79: 2755-2763. 10.1128/IAI.01375-10.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Marot G, Foulley JL, Mayer CD, Jaffrézic F: Moderated effect size and P-value combinations for microarray meta-analyses. Bioinformatics. 2009, 25: 2692-2699. 10.1093/bioinformatics/btp444.

    Article  CAS  PubMed  Google Scholar 

  25. Editorial: Meta-analysis and moderation. Nature Rev Drugs Disc. 2010, 9:747.

  26. Hedges LV, Olkin I: Statistical Methods for Meta-Analysis. 1985, New York, NY, Academic Press

    Google Scholar 

  27. Stangl DK, Berry DA: Meta-Analysis in Medicine and Health Policy. 2000, Marcel Dekker, New York, NY

    Book  Google Scholar 

  28. Schokker D, Smits MA, Hoekman AJW, Parmentier HK, Rebel JMJ: Effects of Salmonella on spatial-temporal processes of jejunal development in chickens. Dev Comp Immunol. 2010, 34: 1090-1100. 10.1016/j.dci.2010.05.013.

    Article  CAS  PubMed  Google Scholar 

  29. Fife MS, Howel JS, Salmon N, Hocking PM, van Diemen PM, Jones MA, Stevens MP, Kaiser P: Genome-wide SNP analysis identifies major QTL for Salmonella colonization in the chicken. Anim Genet. 2010, 42: 134-140.

    Article  Google Scholar 

  30. Barrow PA, Bumstead N, Marston K, Lovell MA, Wigley P: Faecal shedding and intestinal colonization of Salmonella enterica in in-bred chickens: the effect of host genetic background. Epidemiol Infect. 2003, 132: 117-126.

    Article  Google Scholar 

  31. Schokker D, Peters THF, Hoekman AJW, Rebel JMJ, Smits MA: Differences in the early response of hatchlings of different chicken breeding lines to Salmonella enterica serovar Enteritidis infection. Poultry Sci. 2012, in press

    Google Scholar 

  32. Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: No. 1-Article 3

    Google Scholar 

  33. Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, Smyth GK: A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007, 23: 2700-2707. 10.1093/bioinformatics/btm412.

    Article  CAS  PubMed  Google Scholar 

  34. Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nat Protoc. 2009, 4: 44-57.

    Article  CAS  Google Scholar 

  35. Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucl Acid Res. 2009, 37: 1-13. 10.1093/nar/gkn923.

    Article  Google Scholar 

  36. [http://david.abcc.ncifcrf.gov/]

  37. Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc B. 1995, 57: 289-300.

    Google Scholar 

  38. Befus AD, Johnston N, Leslie GA, Bienenstock J: Gut-associated lymphoid tissue in the chicken. I. Morphology, ontogeny, and some functional characteristics of Peyer's patches. J Immunol. 1980, 125: 2626-2632.

    CAS  PubMed  Google Scholar 

  39. Barford D, Das AK, Egloff MP: The structure and mechanism of protein phosphatases: insights into catalysis and regulation. Annu Rev Biophys Biomol Struct. 1998, 27: 133-164. 10.1146/annurev.biophys.27.1.133.

    Article  CAS  PubMed  Google Scholar 

  40. Van Weeren PC, de Bruyn KM, de Vries-Smits AM, van Lint J, Burgering BM: Essential role for protein kinase B (PKB) in insulin-induced glycogen synthase kinase 3 inactivation. Characterization of dominant-negative mutant of PKB. J Biol Chem. 1998, 273: 13150-13156. 10.1074/jbc.273.21.13150.

    Article  CAS  PubMed  Google Scholar 

  41. Blom N, Sicheritz-Pontén T, Gupta R, Gammeltoft S, Brunak S: Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics. 2004, 4: 1633-1649. 10.1002/pmic.200300771.

    Article  CAS  PubMed  Google Scholar 

  42. Daniel JM, Reynolds AB: Tyrosine phosphorylation and cadherin/catenin function. Bioessays. 1997, 19: 883-891. 10.1002/bies.950191008.

    Article  CAS  PubMed  Google Scholar 

  43. Galan JE, Zhou D: Striking a balance: Modulation of the actin cytoskeleton by Salmonella. Proc Natl Acad Sci USA. 2000, 97: 8754-8761. 10.1073/pnas.97.16.8754.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Patel JC, Galan JE: Manipulation of the host actin cytoskeleton by Salmonella–all in the name of entry. Curr Opin Microbiol. 2005, 8: 10-15. 10.1016/j.mib.2004.09.001.

    Article  CAS  PubMed  Google Scholar 

  45. Neer EJ, Schmidt CJ, Nambudripad R, Smith TF: The ancient regulatory-protein family of WD-repeat proteins. Nature. 1994, 371: 297-300. 10.1038/371297a0.

    Article  CAS  PubMed  Google Scholar 

  46. Smith TF, Gaitatzes C, Saxena K, Neer EJ: The WD repeat: a common architecture for diverse functions. Trends Biochem Sci. 1999, 24: 181-185. 10.1016/S0968-0004(99)01384-5.

    Article  CAS  PubMed  Google Scholar 

  47. Li D, Roberts R: WD-repeat proteins: structure characteristics, biological function, and their involvement in human diseases. Cell Mol Life Sci. 2001, 58: 2085-2097. 10.1007/PL00000838.

    Article  CAS  PubMed  Google Scholar 

  48. Pickart CM: Mechanisms underlying ubiquitination. Ann Rev Biochem. 2001, 70: 503-533. 10.1146/annurev.biochem.70.1.503.

    Article  CAS  PubMed  Google Scholar 

  49. Zhang Y, Higashide W, Dai S, Sherman DM, Zhou D: Recognition and ubiquitination of Salmonella type III effector SopA by a ubiquitin E3 ligase, HsRMA1. J Biol Chem. 2005, 280: 38682-38688. 10.1074/jbc.M506309200.

    Article  CAS  PubMed  Google Scholar 

  50. Lotz GP, Brychzy A, Heinz S, Obermann WMJ: A novel HSP90 chaperone complex regulates intracellular vesicle transport. J Cell Sci. 2008, 121: 717-723. 10.1242/jcs.015610.

    Article  CAS  PubMed  Google Scholar 

  51. Csermely P, Schnaider T, Soti C, Prohászka Z, Nardai G: The 90-kDa molecular chaperone family: structure, function, and clinical applications. A comprehensive review. Pharmacol Ther. 1998, 79: 129-168. 10.1016/S0163-7258(98)00013-8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgement

The authors gratefully acknowledge financial participation from the European Community under the Sixth Framework Programme for Research, Technological Development and Demonstration Activities, for the Integrated Project SABRE, Cutting Edge Genomics for Sustainable Animal Breeding, Contract no: FOOD-CT-2006-01625. The views expressed in this publication are the sole responsibility of the author(s) and do not necessarily reflect the views of the European Commission. Neither the European Commission nor any person acting on behalf of the Commission is responsible for the use, which might be made of the information. The information in this document is provided as is and no guarantee or warranty is given that the information is fit for any particular purpose. The user thereof uses the information at its sole risk and liability. Furthermore, additional finances were from the Kennisbasis (Knowledge Base) grants no KB-01-006, KB-04-004-006 and KB-04-004-012 of the Dutch Ministry of Agriculture, Nature and Food security.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Marinus FW te Pas.

Additional information

Competing interests

The authors declare that they have no competing interests

Authors’ contributions

MtP: Coordinated the meta-analysis, participated in the design of the meta-analysis, participated in the analysis of the results of the meta-analysis, and drafted the manuscript. IH: Participated in the design of the meta-analysis, performed the pre-processing of the data, participated in the meta-analysis, and helped in the analysis of the results of the meta-analysis. DS: Carried out the molecular genetic studies, analyzed the molecular results, participated in the design of the meta-analysis, participated in the meta-analysis, and helped in the analysis of the results of the meta-analysis. MS: Participated in the coordination the overall work in the project, participated in the coordination of the meta-analysis, participated in the design and analysis of the meta-analysis, and participated in drafting the manuscript. MF: Carried out the molecular genetic studies, analyzed the molecular results, and participated in drafting the manuscript. RZ: Carried out the animal studies, molecular genetic studies and immunologic reactions. M-LE: Carried out the animal studies, molecular genetic studies and immunologic reactions. AR: Coordinated the overall work in the project, participated in the coordination of the meta-analysis, participated in the design and analysis of the meta-analysis, and participated in drafting the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2011_4058_MOESM1_ESM.xlsx

Additional file 1: The lists of genes with regulated expression in more than one or all experiments of the meta-analysis DE-group. (XLSX 105 KB)

Additional file 2: Functional biological mechanisms of the differently expressed genes. (XLSX 54 KB)

Additional file 3: Details of the biological functions per cluster. (XLSX 26 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

te Pas, M.F., Hulsegge, I., Schokker, D. et al. Meta-analysis of Chicken – Salmonella infection experiments. BMC Genomics 13, 146 (2012). https://doi.org/10.1186/1471-2164-13-146

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-13-146

Keywords