https://orcid.org/0000-0002-5834-2869
Figures
Abstract
Managing Salmonella enterica Enteritidis (SE) carriage in chicken is necessary to ensure human food safety and enhance the economic, social and environmental sustainability of chicken breeding. Salmonella can contaminate poultry products, causing human foodborne disease and economic losses for farmers. Both genetic selection for a decreased carriage and gut microbiota modulation strategies could reduce Salmonella propagation in farms. Two-hundred and twenty animals from the White Leghorn inbred lines N and 61 were raised together on floor, infected by SE at 7 days of age, transferred into isolators to prevent oro-fecal recontamination and euthanized at 12 days post-infection. Caecal content DNA was used to measure individual Salmonella counts (ISC) by droplet digital PCR. A RNA sequencing approach was used to measure gene expression levels in caecal tonsils after infection of 48 chicks with low or high ISC. The analysis between lines identified 7516 differentially expressed genes (DEGs) corresponding to 62 enriched Gene Ontology (GO) Biological Processes (BP) terms. A comparison between low and high carriers allowed us to identify 97 DEGs and 23 enriched GO BP terms within line 61, and 1034 DEGs and 288 enriched GO BP terms within line N. Among these genes, we identified several candidate genes based on their putative functions, including FUT2 or MUC4, which could be involved in the control of SE infection, maybe through interactions with commensal bacteria. Altogether, we were able to identify several genes and pathways associated with differences in SE carriage level. These results are discussed in relation to individual caecal microbiota compositions, obtained for the same animals in a previous study, which may interact with host gene expression levels for the control of the caecal SE load.
Citation: Cazals A, Rau A, Estellé J, Bruneau N, Coville J-L, Menanteau P, et al. (2022) Comparative analysis of the caecal tonsil transcriptome in two chicken lines experimentally infected with Salmonella Enteritidis. PLoS ONE 17(8): e0270012. https://doi.org/10.1371/journal.pone.0270012
Editor: Michael H. Kogut, USDA-Agricultural Research Service, UNITED STATES
Received: May 31, 2022; Accepted: August 3, 2022; Published: August 17, 2022
Copyright: © 2022 Cazals et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All sequencing data files are available from the NCBI Sequence Read Archive (SRA) (Bioproject accession number PRJNA649900).
Funding: Joint program of the Institutes Carnot Santé Animale (ICSA) and Pasteur Maladies Infectieuses (PMI). The PhD of Anais Cazals was funded by Région Ile de France (50%) and the Animal Genetics Division of INRAE (50%). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Salmonella is a zoonotic pathogen that can cause human foodborne disease. In 2019, more than 80,000 human salmonellosis cases were confirmed in Europe with 140 reported deaths [1]. Among these cases, more than 9000 were associated with 926 food-borne outbreaks (FBOs), with a large majority (72.4%) caused by the serovar S. Enterica Enteritidis (SE). Eggs produced by infected layer hens seem to be the major food vehicle, representing more than 37% of the FBOs. In parallel, in spite of strict hygiene control in farms, the systematic detection of Salmonella serovars, and the use of vaccination, the prevalence of Salmonella in laying hen flocks increased from 2.07% in 2014 to 3.44% in 2019 [1]. In chicken, the carriage is asymptomatic. The bacteria can persist a long time in the gut and can quickly spread within a contaminated farm via oro-fecal recontaminations between birds [2]. Understanding the impact of factors such as host genetics or gut microbiota on Salmonella carriage, and even more their combined impact, could lead to innovative strategies to reduce Salmonella transmission and ensure human food safety.
The caecal tonsil is a major barrier controlling the entry of bacteria in the organism [3, 4] and is therefore a tissue particularly relevant for identifying host factors potentially involved in the control of SE. Several studies have been conducted on the caecal tonsil transcriptome. In particular, they have helped to identify biological processes associated with resistance to S. Enteritidis [5–7], S. Typhimurium [8, 9] and S. Pullorum [10]. Nevertheless, in these studies, the impact of host genetics on gene expression was not examined. The expression of specific immune genes has been compared between the two experimental inbred chicken lines 61 and 15I, but not whole transcriptome [11]. The impact of host genetic variations was considered in a recent study of the caecal tissue transcriptome after Campylobacter colonisation. Comparisons between the experimental White Leghorn inbred chicken lines 61 and N led to the identification of a large number of differentially expressed genes, which may underlie variation in heritable resistance to the pathogen [12].
The host genetic background is an important factor for the outcome of Salmonella infection in chicken. A number of quantitative trait loci (QTL) and candidate genes associated with Salmonella resistance have been identified [13]. In the inbred lines N (resistant) and 61 (susceptible) in particular, several QTLs with low to moderate effects were identified [14–16]. However, no causal gene could be pinpointed due to the large size of the QTL genomic regions. More generally, only a few genes have been identified for their direct implication in the control of SE load in chicken, and knowledge is lacking about the mechanisms leading to genetic resistance. For the studies conducted on Salmonella carriage in the N and 61 lines, birds were reared together on floor after infection, thus allowing Salmonella oro-fecal recontamination between birds. In the present study, we used another infection model, making use of isolators. Previously tested on the experimental White Leghorn line PA12, this model showed a strong reduction of oro-fecal recontaminations, leading to much increased Salmonella individual variation among birds [2]. It is therefore an interesting model to identify birds with highly contrasted carriage levels, in order to facilitate the identification of host genes involved in these differences.
In the current study, we performed an analysis using information about caecal tonsil gene expression in two distinct genetic lines: the inbred chicken lines N and 61, respectively resistant and susceptible to SE infection. The objectives of this study were to:
- identify differentially expressed genes between genetic lines in the caecal tonsils after SE infection, in order to identify potential pathways involved in the genetic resistance to SE;
- identify genes and pathways associated with SE resistance within line (low vs high carriers);
Results
Two-hundred and forty animals from the two experimental White Leghorn inbred lines N and 61 were raised together on floor until 7 days of age. Then, chicks were challenged with Salmonella enterica Enteritidis (SE) LA5 by oral infection and separated into four isolators. Two independent replicates (n = 120) were conducted with a total number of 240 chicks. No clinical signs of disease were observed on the animals. Caecal contents and caecal tonsils were collected at 12 days post infection.
The abundance of SE in caecal contents was measured by Droplet Digital PCR (ddPCR) and, as described previously, significant differences of Salmonella abundance were observed between lines for the two experiments [17]. The observed variability of carriage allowed us to identify extreme low and high carriers within each line, and 48 extreme animals were selected for the caecal tonsil RNA extraction, balancing the “experiment”, “isolator” and “sex” factors. Different groups were defined as described in Fig 1 according to the line, class and experiment factors. Means, standard deviations and p-values according to these groups are also given in Fig 1.
Mean (red points) and standard deviation (SD; red bars) of the Salmonella Enteritidis abundance at 12 dpi in caecal contents (log10/g of caceal contents) of chicken groups infected with SE according to the line and experiment. Difference in mean between low and high carriers according to the line and experiment and t-test p-value. Difference in mean between lines and t-test p-value. Results of the DEG and BP GO term enrichment analyses. Results of the 16S analysis from the study [17].
Differentially expressed genes (DEGs) in caecal tonsils between lines
On average, more than 40M reads were sequenced for each of the 48 samples. After quality control, a total of 24,356 expressed genes were identified and used for the following analyses. Using all 48 samples, a principal component analysis (PCA) showed a distinct clustering between lines (Fig 2). Two ANOVA analyses on the PCA1 and PCA2 components showed that gene expression is significantly affected by line and sex (S1 File). A differential analysis with DESeq2 allowed the identification of 7,516 DEGs between the two lines (p-adj<0.05) among which 3,944 were up- and 3,572 were down-regulated (S1 Fig and S1 Table).
The first principal component contains 29% of the variance and may be attributed to the line effect and the second principal component contains 25% of the variance and may be attributed to the sex effect.
DEGs in caecal tonsils between low and high carriers within lines and experiments
Gene expression levels between low and high Salmonella carriers within each line and each experiment (L6L1/L6H1, L6L2/L6H2 and LNL1/LNH1, LNL2/LNH2) were compared. A PCA showed a clustering between the high and lower carriers within each of the N and 61 lines in experiment 1 (Fig 3; LNH1 vs LNL1 and L6H1 vs L6L1), but not in the equivalent groups of experiment 2. ANOVA analysis on the PCA1 and PCA2 components showed that gene expression is significantly (P<0.05) affected by low/high classes and by sex in experiment 1, but not in experiment 2 (S1 File).
(A) Principal component analysis (PCA) of the gene expression in caecal tonsils within line 61 in experiment 1. (B) PCA of the gene expression in caecal tonsils within line N in experiment 1.
A differential analysis with DESeq2 between low and high carriers within line 61 in experiment 1 (L6L1/L6H1) allowed the identification of 97 DEGs (p-adj < 0.05), among which 42 were up- and 55 were down-regulated (S2 Fig and S2 Table). A similar analysis performed between low and high carriers within line N in experiment 1 (LNL1/LNH1) allowed the identification of 1,034 DEGs (p-adj< 0.05) with 794 up- and 240 down-regulated genes (S3 Fig and S3 Table). Only 1 DEG was shared between these two comparisons (Fig 4). In experiment 2, no significant DEGs were found between low and high carriers regardless of the line (S4 and S5 Tables). The results are summarized in Fig 1. Merging both experiments in a single analysis, including a fixed effect for the experiment did not provide conclusive results.
Differentially expressed genes identified by comparing different groups of animals: (A) L6L1/L6H1 and LNL1/LNH1 (low vs high carriers of lines 61 and N, respectively, in experiment 1), (B) L6L1/L6H1 and LN/L6 (low vs high carriers of line 61 in experiment 1, and line N vs line 61), and (C) LNL1/LNH1 and LN/L6 (low vs high carriers of line N in experiment 1, and line N vs line 61).
DEGs common to intra-line and between-line analyses
When comparing DEGs identified between lines 61 and N, and within line 61 between low and high carriers, 38 genes were shared. Among these shared genes, 9 genes appeared to be regulated in the same direction: 4 DEGs (CEMIP, DMXL1, FUT2, NOS2) were upregulated in both the resistant line N (low carriage) and in low carriers in line 6, and 5 DEGs (PLEKHS1, RPS6KB2, CYP4B7, CYP2D6, CYP2AC1) were downregulated in both line N and in low carriers in line 61 (Fig 4 and Tables 1 and 2).
When comparing DEGs identified between lines 61 and N and within line N between low and high carriers, 560 genes were shared. Among these shared genes, 58 appeared to be regulated in the same direction: 43 DEGs were upregulated in both the resistant line N (low carriage) and in low carriers in line N, and 15 DEGs were downregulated in both line N and in low carriers in line N (Fig 4 and Tables 3 and 4).
Functional enrichment analysis
Enrichment analyses were performed with topGO on DEGs between all animals from lines N and 61 and between low and high carriers within each line in experiment 1. The comparison between lines N and 61 led to the identification of 62 significantly (p-value < 0.05) enriched Biological Processes (BP) Gene Ontology (GO) terms (S6 Table). The comparison between low and high carriers of lines 61 (L6L1/L6H1) led to the identification of 23 BP GO terms (S7 Table). The comparison between low and high carriers of line N (LNH1/ LNL1) led to the identification of 288 BP GO terms (S8 Table).
The results are summarized in Fig 1. Three BP GO terms were shared between the L6L1/L6H1 and N/6 analyses (Table 5), fourteen between the LNL1/LNH1 and LN/L6 analyses (Table 6) and four between the L6L1/L6H1 and LNL1/LNH1 analyses (Table 7). Interestingly, the 3 GO BP terms enriched and shared between L6L1/L6H1 and N/6 were related to the response to biotic stimulus and other organisms.
Discussion
We explored the caecal tonsil transcriptome after Salmonella Enteritidis infection by comparing samples from two genetic lines displaying contrasted levels of Salmonella carriage after infection: lines N and 61. These inbred lines of chicken have been used in many infection studies, with results in apparent contradiction with ours. While N is more resistant to Salmonella Enteritidis in this study and in our previous experiments [14] in comparison to line 61, it was more susceptible with Salmonella Typhimurium [18, 19] or Campylobacter [12]. However, many factors may explain this difference: a potential genetic drift of the line (maintained for years at different experimental stations), environmental differences (different experimental farms), the infection model (the route of infection, the phenotype measured, the organ targeted, etc) [13]. However, to our opinion the most important factor is the pathogen itself: ST, Campylobacter or SE. Animals can display different resistance mechanisms to distinct pathogens, and even to distinct serotypes of the same pathogen. This was evidenced for instance in a previous study comparing genetic locations of QTLs for resistance to SE and to ST: while two QTLs were probably common to both serotypes, several QTLs were specific of the serotype, thus revealing the probable existence of partially distinct mechanisms for those two serotypes of Salmonella [20].
A large difference of gene expression between lines
We showed a strong impact of the genetic background on gene expression in caecal tonsils, with 7,516 significant DEGs and 62 GO BP terms identified between the two lines. This strong difference in gene expression between lines is likely the result of genetic differences between the two lines. Some of these genes could explain the differences in susceptibility to Salmonella between lines. Thus, several genes display functions which may explain the higher resistance of line N, in which they were expressed to a greater extent. The latter genes code for the major histocompatibility complex I or II (MHCIBF2, MHCIYF5, MHCIIBLB1, MHCIIBLB2), antimicrobial peptides such as granzyme A and K, or the avian beta-defensins 10, 13 and 14. In the same way, the genes TLR1,2,4,7,15 [21–24], NOS2, Gal 13, PSAP and IGL, which have already been associated with SE resistance in a genetic study in chicken [13], were more expressed in the resistant line N.
Interestingly, we showed previously that the caecal microbiota composition of these animals was highly different between individuals from lines N and 61 [17]. Do some of these DEGs in caecal tonsils indirectly impact Salmonella carriage through the modulation of the caecal microbiota composition? Conversely, do differences of microbiota composition indirectly impact Salmonella carriage through the modulation of gene expression in caecal tonsils? Further elements are needed to answer these questions.
A difference of gene expression between low and high carriers within line
The comparisons between low and high carriers within each line in experiment 1 (L6L1/L6H1 and LNL1/LNH1) revealed only one DEG and four BP GO terms in common, leading us to the conclusion that host pathways leading to a higher resistance could differ between lines. It may be explained by their genetic differences. However, one of the four common BP GO terms may be directly related to the Salmonella infection process: the transmembrane receptor protein tyrosine kinase signalling pathway (GO:0007169), which is enriched in low carriers in both lines. Several studies have shown that the activation of protein tyrosine kinases and other specific transcription factors directly affects the innate immune response during SE infection [5, 25]. Therefore, the transcriptional up-regulation of this pathway may be one of the mechanisms by which the animals from both lines control the infection.
Interestingly, in spite of the smaller difference between low and high carriers observed in line N compared to line 61, we identified many more DEGs and BP GO terms in line N, compared to line 61 (1034 DEGs and 288 BP GO terms in line N vs. 97 DEGs and 23 BP GO terms in line 61). Previously, we showed that the microbiota composition was significantly different between low and high carriers in line 61 but not between groups in line N [17]. One tentative hypothesis could be that the largest variability in Salmonella carriage in line 61 could be explained by differences in gut microbiota composition, as demonstrated with the PA12 chicken line [26]. The resistance against Salmonella in line 61 may thus be driven by mechanisms involving the intestinal microbiota, whereas in line N resistance mechanisms could be triggered by other factors, such as intra-line genetic variations.
Even in controlled environmental conditions with inbred lines, variability in the level of Salmonella carriage was identified within each line, as expected with the use of isolators preventing inter-individual recontamination [2]. This variability may derive from the residual genetic variability remaining within each line, or it may be caused by variations in the caecal microbiota composition, which is highly variable between lines and within line 61 [17, 26]. However, this variability was much higher in experiment 1, allowing the identification of more contrasted low and high classes than in experiment 2. This higher variability, and thus greater contrast in host response to Salmonella infection between low and high carriers, is probably the reason why significant DEGs and BP GO terms could be identified in experiment 1 but not experiment 2. Experimental conditions were highly controlled and similar between experiments but some unchecked factors, such as the environmental microbial exposure at hatch, may differ and explain the observed differences between the two. Results in experiment 2 could be used as a “tendency validation” of the experiment 1. For example, in experiment 2, MUC4 was identified as more expressed in resistant animals, as in experiment 1 but with a higher pvalue, closer to the significance threshold (pvalue = 0.08 and padj = 0.9).
Identification of genes associated with the resistance to Salmonella infection
To highlight host genes that may be consistently associated with resistance to SE infection, we decided to focus on common DEGs between the intra-line and inter-line analyses. Indeed, genes more highly expressed in both the resistant line N and in low carriers in the susceptible line 61, may be more likely to be associated with the resistance to SE infection. We further investigated if some of these genes have already been associated with the control of Salmonella infection in previous studies, supporting the reliability of our analysis.
Four genes were identified as more highly expressed in both the resistant line N and in low carriers of the susceptible line 61 and are therefore associated with Salmonella resistance. Three of these genes have functions that are related to Salmonella: NOS2, DMXL1 and FUT2. First, NOS2 is a gene coding for the inducible nitric oxide synthase (iNOS) protein involved in macrophage inflammatory response [27]. Its function in the innate immunity against bacteria, viruses, fungi and parasites is well established, especially against Salmonella Typhimurium in mice models [28–33] and in in vitro models [34, 35]. In chicken, a transcriptome analysis of caeca showed an increase of NOS2 expression in the SE infected group [6], while genetic studies showed an association between NOS2 gene alleles and spleen SE bacterial load [36, 37]. An in vitro test showed the implication of iNOS between pathogen and macrophage cells during SE infection [38]. Finally, CNP (chitosan-nanoparticle) vaccination seems to increase NOS2 expression and protect against SE infection in chicken [39]. Second, DMXL1 is a gene involved in the phagosome acidification in macrophages. It seems to have an impact on innate immunity and macrophage bacteria killing after an activation by the TPL-2 kinase [40]. Finally, FUT2 is a gene coding for the α-1,2-fucosyltransferase enzyme involved in the glycosylation profile of the gastrointestinal tract. In mice and human, it has been shown that a “non-secretor” individual heterozygous for a loss-of-function mutation is more susceptible to chronic intestinal diseases such as Crohn’s disease and to pathogen infection [41]. More specifically, it has been shown that “non secretor” mice (Fut-/-) show an increase of Salmonella Typhimurium in caecal tissue compared to wild-type mice [42, 43]. Moreover, a FUT2 polymorphism was associated with the human faecal microbiota composition and diversity [44], which could explain host-microbe interactions and susceptibility to infection [45]. Finally, it has been shown that FUT2 was associated with the abundance of Christensenellaceae [46], a bacteria family we identified as associated with Salmonella Enteritidis resistance in the same animals [17]. Thus, FUT2 may be a gene indirectly associated with SE resistance through the modulation of the microbiota composition, modulating the abundance of competitive bacteria against Salmonella.
Forty-two genes were identified as more expressed both in the resistant line N and in low carriers of the resistant line N and may therefore be associated with Salmonella resistance. Eight of these genes have functions of interest: MHCIBF2, CLDN7, SIRT5, ENTPD1, SYT7, SLC22A23, S100B and MUC4. The MHCIBF2 (Major histocompatibility complex class I antigen BF2) gene is the predominant ligand of cytotoxic T lymphocytes [47]. It has been shown that a particular MHC I haplotype may contribute to control the response to SE infection in chicken [48, 49]. The CLDN7 (claudin 7) gene is involved in the formation of tight junctions between epithelial cells. It seems that a downregulation of CLDN7 by pathogen could facilitate translocation of invasive bacteria across the epithelium [50]. The SIRT5 (sirtuin 5) gene could have large impact on cellular homeostasis and is more expressed in colorectal cancer [51]. The ENTPD1 (ectonucleoside triphosphate diphosphohydrolase 1 or CD39) gene has an impact on inflammatory bowel disease (IBD) and regulation of pro-inflammatory responses and pathogen colonization [52–54], as does the SLC22A23 (Solute Carriers family) gene, which is associated with intestinal inflammation in human [55]. The SYT7 (synaptotagmin VII) gene is associated with the control of cytotoxic granule fusion in lymphocytes, and mice lacking syt7 have reduced ability to clear an infection [56]. The following two genes could be indirectly associated with SE resistance through the modulation of the microbiota composition. The S100B (S100 calcium binding protein B) gene codes for a signalling molecule which could be implicated in the communication mechanisms between microbiota and gut, and could explain differences between healthy and pathological microbiota in human [57]. Finally, it has been established that MUC4 (mucin 4, cell surface associated) and more generally MUC genes have an important role in preventing pathogens infections [58]. Many studies in chicken showed the association of the expression of mucin genes, especially MUC2, with SE infection [59–61]. Mucins favour the establishment and the maintenance of a commensal microbiota, and form a protection barrier against pathogens. In pigs, a decrease of MUC4 gene expression has been associated with ST infection [62], and genetic variants in this gene have been associated with the increase of gene expression relative to immune function and gut homoeostasis [63, 64].
Identification of genes associated with the susceptibility to Salmonella infection
Five genes were identified as less expressed both in the resistant line N and in low carriers of the susceptible line 61 and may therefore be associated with Salmonella susceptibility. Among these genes, three belong to the cytochrome P450 family: CYP4B7, CYP2D6 and CYP2AC1. CYP enzymes play a key role in metabolic processes in the intestine as the metabolism of xenobiotic substances. Metabolism of CYP enzymes is closely connected with infection, inflammation and intestinal microbiota in human [65, 66].
Fifteen genes were identified as less expressed both in the resistant line N and in low carriers of the resistant line N and may therefore be associated with Salmonella susceptibility. Four of these genes have immune functions: LECT2, DEFB4A, LOC121107850 and EXP. The LECT2 (leukocyte cell derived chemotaxin 2) gene has an antibacterial function, is expressed in chicken heterophils, and increases in abundance in macrophage after SE infection [53]. It is also more expressed in vaccinated chicken again SE. The DEFB4A (defensin beta 4A) gene has an important role in innate immunity in mucosal tissues through its antimicrobial activity against various microorganisms. It has been shown that beta-defensin plays a role in immunoprotection against Salmonella Enteritidis in in vitro embryonic chicken cell model [54] and in the development of innate immunity in gastrointestinal tract of newly hatched chicks [55]. The LOC121107850 (T-cell-interacting, activating receptor on myeloid cells protein 1-like) gene codes for an activating receptor on myeloid cells protein 1-like, and the EPX (eosinophil peroxidase) gene codes for an enzyme in myeloid, which has a function in bacterial destruction [67].
The increase in the expression of these genes is probably related to the high level of Salmonella infection. However, with this type of experimental design it is not possible to decipher whether host gene expression is responsible for high levels of Salmonella infection, or whether it is a consequence of the high infection.
Comparing the results obtained in the analyses between genetic lines on the one hand, and between high and low carriers within each line on the other hand, seems interesting to highlight genes having a meaningful impact in the response to SE infection. All of the genes discussed in these two last sections appear to be associated with SE infection and a modulation of the expression of these genes could prevent the colonization of the pathogen.
Identification of BP GO terms associated with the response to Salmonella infection
Some of the BP GO terms identified are implicated in immune response: notably immune system process (GO:0002376), defense response (GO:0006952) or cell surface receptor signalling pathway (GO:0007166). These BP GO terms have been already identified in another study working on the jejunum of a commercial genetic line and comparing high and low Salmonella infection using a kinome peptide array [68]. Thus, despite differences in experimental designs and methods used, similar results were found.
The three BP GO terms identified both between lines N and 61 and between low and high carriers in line 61 (LN/L6 and L6L1/L6H1) have interesting links to immunity: response to biotic stimulus, response to external biotic stimulus and response to other organisms (Table 5). Seven of the fourteen BP GO terms identified both between lines N and 61 and between low and high carriers in line N (LN/L6 and LNL1/LNH1) are related to the regulation of the polymerisation or depolymerisation of proteins as actin (Table 6). It has been shown that a cytoskeletal actin rearrangement is induced by SE invasion in host cells [69]. Indeed, the actin cytoskeleton is targeted by Salmonella to promote its invasion, survival and growth in cells [70]. Genes implicated in these pathways may contribute to response to Salmonella infection and could be under genetic control.
Conclusion
The two experimental genetic lines N and 61, displaying contrasted levels of Salmonella carriage after infection, showed a large difference of caecal tonsil gene expression associated with the outcome to SE infection. The comparison of resistant chicks (from line N or from low carriers within both lines) and susceptible chicks (from line 61 or from high carriers within both lines) allowed us to identify several genes and pathways associated with Salmonella resistance. Different mechanisms seem to be involved in the response to SE between these two experimental lines. A lower number of DEGs is associated with a larger inter-individual variability in line 61 compared to line N.
Materials and methods
Experimental design
All animal procedures were authorised by the Ethic committee: APAFIS#5833-2016062416362298v3. Animals from the two experimental White Leghorn inbred lines N and 61 were provided by the experimental unit PEAT (Pole d’Expérimentation Avicole de Tours, Nouzilly, France). As described previously [17], animals were raised together on the floor until infection at the PFIE unit (Plateforme d’Infectiologie Expérimentale, INRA, Nouzilly, France) with free access to food and water. At 7 days of age, chicks were orally infected with Salmonella enterica Enteritidis (Strain 775 [LA5 Nal20Sm500], 5.104 cfu/0.2 mL/chick) and immediately separated into isolators. Four isolators were used for each experiment: two isolators for chicks from line N and two others for chicks from line 61, with 30 birds per isolator. Caecal contents and caecal tonsils were collected at 12 days post infection (19 days of age) after the animal sacrifice and were immediately frozen in liquid nitrogen and stored at -80°C until use. Two experiments were conducted with a total number of 240 chicks.
DNA extraction, Salmonella count by Droplet Digital PCR and choice of low and high carriers
As described previously [17], individual caecal DNA was extracted from an average of 200 mg of frozen caecal contents and DNA samples were stored at -20°C. Individual abundances of Salmonella Enteritidis in caecal contents were obtained by Droplet Digital PCR (ddPCR) using the QX200 Droplet Digital PCR system (Bio-Rad) at the @bridge platform (INRAE, Jouy-en-Josas, France). The ddPCR method has been proven to reliably quantify the amount of salmonella spp. in a sample [71]. We targeted and amplified a region of the InvA gene specific to SE [72]. Data were analysed with a log transformation of the copies of Salmonella. Analyses of variance (ANOVA) were performed to test the significance of differences of the copies of Salmonella according to different factors (line, sex, experiment or isolator) using the anova function from base R (Type I sum of squares). We selected 48 chicks, either low or high Salmonella carriers based on data from ddPCR, balancing the experiment, sex, isolator and line factors.
RNA extraction and sequencing
Caecal tonsils from the 48 low/high carrier chicks were first grinded using ULTRA-TURRAX T25 (IKA). RNAs were then extracted using the NucleoSpin RNA Kit (MACHEREY-NAGEL) according to the manufacturer’s protocol. The quantity of RNA was measured using a Nanodrop spectrophotometer (Thermo Scientific) and its quality was assessed using a 2100 Bioanalyzer Expert system using a total RNA nano Kit (Agilent), all samples displayed an RNA integrity number (RIN) > 7. RNAs were sent to the genomic platform GeT-Plage (Toulouse, France) for the cDNA library preparation (TruSeq Stranded mRNA kit, Illumina) and the sequencing (NovaSeq 6000 Sequencing System, Illumina). The sequencing data analysed during the current study are available in the NCBI Sequence Read Archive (SRA) database under the Bioproject accession number PRJNA649900.
Transcriptome analysis
The quality control of raw reads was performed using the FastQC program. Trimming was performed using Sickle (version 1.33). Reads were then mapped on the Gallus reference genome (Gallus_gallus.GRCg6a.98.gtf) with STAR (version 2.5.3a) [73]. Genes were counted with htseq-count (version 0.12.4) [74].
The metadata and table of gene counts have been included as Additional files 6 and 7, respectively. PCA and ANOVA tests on PCA components were performed with the stats package (3.6.1). DESeq2 (version 1.26.0) was used with R version 3.6.1 to perform differential gene expression analyses [75]. In our model, the sex and the isolator effects were fixed. Genes with an adjusted p-value < 0.05 were considered to be differentially expressed genes (DEGs). The ggplot2 packages (3.2.1) were used for visualisation. The functional enrichment analyses were performed with topGO package version 2.38.1 [76].
Supporting information
S1 File. ANOVA analysis results on PCA1 and PCA2 components.
https://doi.org/10.1371/journal.pone.0270012.s001
(TXT)
S2 File. Metadata associated with all samples.
https://doi.org/10.1371/journal.pone.0270012.s002
(TXT)
S1 Fig. DESeq2 results for the comparison of the N and 61 lines.
Distribution of raw p-values (A), dispersion plot (B), and MA plot highlighting 1492 and 1285 up- and down-regulated genes with a p-adj < 0.05 (C).
https://doi.org/10.1371/journal.pone.0270012.s004
(TIFF)
S2 Fig. DESeq2 results for the comparison of the low and high carriers within line 61 in experiment 1.
Distribution of raw p-values (A), dispersion plot (B), and MA plot highlighting 42 and 55 up- and down-regulated genes with a p-adj < 0.05 (C).
https://doi.org/10.1371/journal.pone.0270012.s005
(TIFF)
S3 Fig. DESeq2 results for the comparison of the low and high carriers within line N.
Distribution of raw p-values (A), dispersion plot (B), and MA plot highlighting 794 and 240 up- and down-regulated genes with a p-adj < 0.05 (C).
https://doi.org/10.1371/journal.pone.0270012.s006
(TIFF)
S1 Table. Differentially expressed genes between LN and L6 groups.
https://doi.org/10.1371/journal.pone.0270012.s007
(XLSX)
S2 Table. Differentially expressed genes between L6H1 and L6L1 groups in experiment 1.
https://doi.org/10.1371/journal.pone.0270012.s008
(XLSX)
S3 Table. Differentially expressed genes between LNL1 and LNH1 groups in experiment 1.
https://doi.org/10.1371/journal.pone.0270012.s009
(XLSX)
S4 Table. Differentially expressed genes between L6H2 and L6L2 groups in experiment 2.
https://doi.org/10.1371/journal.pone.0270012.s010
(TXT)
S5 Table. Differentially expressed genes between LNL2 and LNH2 groups in experiment 2.
https://doi.org/10.1371/journal.pone.0270012.s011
(TXT)
S6 Table. Enrichment analyse results performed on DEGs between line N and 61.
https://doi.org/10.1371/journal.pone.0270012.s012
(XLSX)
S7 Table. Enrichment analyse results performed on DEGs between low and high carriers in line 61—experiment 1.
https://doi.org/10.1371/journal.pone.0270012.s013
(XLSX)
S8 Table. Enrichment analyse results performed on DEGs between low and high carriers in line N—experiment 1.
https://doi.org/10.1371/journal.pone.0270012.s014
(XLSX)
Acknowledgments
We thank our colleagues from the experimental unit PEAT who raised the parents of the animals studied and provided the chicks used in our study. We also thank colleagues from the experimental unit PFIE, who efficiently monitored the experiments and collected the samples, and colleagues from the SPVB research team, who helped collect the samples.
References
- 1. EFSA. The European Union One Health 2019 Zoonoses Report. EFSA J 2021. 2021;19(2):6406. pmid:33680134
- 2. Menanteau P, Kempf F, Trotereau J, Virlogeux-Payant I, Gitton E, Dalifard J, et al. Role of systemic infection, cross contaminations and super-shedders in Salmonella carrier state in chicken. Environ Microbiol. 2018;20(9):3246–60. pmid:29921019
- 3. Setta AM, Barrow PA, Kaiser P, Jones MA. Early immune dynamics following infection with Salmonella enterica serovars Enteritidis, Infantis, Pullorum and Gallinarum: cytokine and chemokine gene expression profile and cellular changes of chicken cecal tonsils. Comp Immunol Microbiol Infect Dis. 2012;35(5):397–410. pmid:22512820
- 4. Berndt A, Wilhelm A, Jugert C, Pieper J, Sachse K, Methner U. Chicken cecum immune response to Salmonella enterica serovars of different levels of invasiveness. Infect Immun. 2007;75(12):5993–6007. pmid:17709416
- 5. Kogut MH, Swaggerty CL, Byrd JA, Selvaraj R, Arsenault RJ. Chicken-Specific Kinome Array Reveals that Salmonella enterica Serovar Enteritidis Modulates Host Immune Signaling Pathways in the Cecum to Establish a Persistence Infection. Int J Mol Sci. 2016;17(8):1207. pmid:27472318
- 6. Matulova M, Varmuzova K, Sisak F, Havlickova H, Babak V, Stejskal K, et al. Chicken innate immune response to oral infection with Salmonella enterica serovar Enteritidis. Vet Res. 2013;44(1):37. pmid:23687968
- 7. Luan DQ, Chang GB, Sheng ZW, Zhang Y, Zhou W, Li ZZ, et al. Analysis of Gene Expression Responses to a Salmonella Infection in Rugao Chicken Intestine Using GeneChips. Asian-Australas J Anim Sci. 2012;25(2):278–85. pmid:25049563
- 8. Wang F, Zhang J, Zhu B, Wang J, Wang Q, Zheng M, et al. Transcriptome Analysis of the Cecal Tonsil of Jingxing Yellow Chickens Revealed the Mechanism of Differential Resistance to Salmonella. Genes. 2019;10(12). pmid:31795199
- 9. Khan S, Chousalkar KK. Transcriptome profiling analysis of caeca in chicks challenged with Salmonella Typhimurium reveals differential expression of genes involved in host mucosal immune response. Appl Microbiol Biotechnol. 2020;104(21):9327–42. pmid:32960293
- 10. Ma T, Chang G, Chen R, Sheng Z, Dai A, Zhai F, et al. Identification of key genes in the response to Salmonella enterica Enteritidis, Salmonella enterica Pullorum, and poly(I:C) in chicken spleen and caecum. BioMed Res Int. 2014;2014:154946. pmid:24707473
- 11. Sadeyen JR, Trotereau J, Velge P, Marly J, Beaumont C, Barrow PA, et al. Salmonella carrier state in chicken: comparison of expression of immune response genes between susceptible and resistant animals. Microbes Infect. 2004;6(14):1278–86. pmid:15555534
- 12. Russell KM, Smith J, Bremner A, Chintoan-Uta C, Vervelde L, Psifidi A, et al. Transcriptomic analysis of caecal tissue in inbred chicken lines that exhibit heritable differences in resistance to Campylobacter jejuni. BMC Genomics. 2021;22(1):411. pmid:34082718
- 13. Calenge F, Kaiser P, Vignal A, Beaumont C. Genetic control of resistance to salmonellosis and to Salmonella carrier-state in fowl: a review. Genet Sel Evol. 2010;42(1):11. pmid:20429884
- 14. Tilquin P, Barrow PA, Marly J, Pitel F, Plisson-Petit F, Velge P, et al. A genome scan for quantitative trait loci affecting the Salmonella carrier-state in the chicken. Genet Sel Evol GSE. 2005;37(5):539–61. pmid:16093014
- 15. Calenge F, Lecerf F, Demars J, Feve K, Vignoles F, Pitel F, et al. QTL for resistance to Salmonella carrier state confirmed in both experimental and commercial chicken lines. Anim Genet. 2009;40(5):590–7. pmid:19422366
- 16. Calenge F, Vignal A, Demars J, Fève K, Menanteau P, Velge P, et al. New QTL for resistance to Salmonella carrier-state identified on fowl microchromosomes. Mol Genet Genomics. 2011;285(3):237–43. pmid:21279652
- 17. Cazals A, Estellé J, Bruneau N, Coville JL, Menanteau P, Rossignol MN, et al. Differences in caecal microbiota composition and Salmonella carriage between experimentally infected inbred lines of chickens. Genet Sel Evol GSE. 2022;54(1):7. pmid:35093028
- 18. 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. 2004;132(1):117–26. pmid:14979597
- 19. Fife MS, Howell JS, Salmon N, Hocking PM, van Diemen PM, Jones MA, et al. Genome-wide SNP analysis identifies major QTL for Salmonella colonization in the chicken: Major QTL for Salmonella colonization in the chicken. Anim Genet. 2011;42(2):134–40. pmid:20579012
- 20. Thanh-Son T, Catherine B, Nigel S, Mark F, Pete K, Elisabeth LBD, et al. A maximum likelihood QTL analysis reveals common genome regions controlling resistance to Salmonella colonization and carrier-state. BMC Genomics. 2012;13(1):198. pmid:22613937
- 21. Nerren JR, Swaggerty CL, MacKinnon KM, Genovese KJ, He H, Pevzner I, et al. Differential mRNA expression of the avian-specific toll-like receptor 15 between heterophils from Salmonella-susceptible and -resistant chickens. Immunogenetics. 2009;61(1):71–7. pmid:19002681
- 22. Abasht B, Kaiser MG, Lamont SJ. Toll-like receptor gene expression in cecum and spleen of advanced intercross line chicks infected with Salmonella enterica serovar Enteritidis. Vet Immunol Immunopathol. 2008;123(3‑4):314–23. pmid:18394716
- 23. MacKinnon KM, He H, Nerren JR, Swaggerty CL, Genovese KJ, Kogut MH. Expression profile of toll-like receptors within the gastrointestinal tract of 2-day-old Salmonella enteriditis-infected broiler chickens. Vet Microbiol. 2009;137(3‑4):313–9. pmid:19201111
- 24. Swaggerty CL, He H, Genovese KJ, Duke SE, Kogut MH. Loxoribine pretreatment reduces Salmonella Enteritidis organ invasion in 1-day-old chickens. Poult Sci. 2012;91(4):1038–42. pmid:22399745
- 25. Swaggerty CL, He H, Genovese KJ, Pevzner IY, Kogut MH. Protein tyrosine kinase and mitogen-activated protein kinase signalling pathways contribute to differences in heterophil-mediated innate immune responsiveness between two lines of broilers. Avian Pathol J WVPA. 2011;40(3):289–97. pmid:21711188
- 26. Kempf F, Menanteau P, Rychlik I, Kubasová T, Trotereau J, Virlogeux-Payant I, et al. Gut microbiota composition before infection determines the Salmonella super- and low-shedder phenotypes in chicken. Microb Biotechnol. 2020;13(5):1611–30. pmid:32639676
- 27. McNeill E, Crabtree MJ, Sahgal N, Patel J, Chuaiphichai S, Iqbal AJ, et al. Regulation of iNOS function and cellular redox state by macrophage Gch1 reveals specific requirements for tetrahydrobiopterin in NRF2 activation. Free Radic Biol Med. 2015;79:206–16. pmid:25451639
- 28. Mastroeni P, Vazquez-Torres A, Fang FC, Xu Y, Khan S, Hormaeche CE, et al. Antimicrobial actions of the NADPH phagocyte oxidase and inducible nitric oxide synthase in experimental salmonellosis. II. Effects on microbial proliferation and host survival in vivo. J Exp Med. 2000;192(2):237–48. pmid:10899910
- 29. Cherayil BJ, Antos D. Inducible nitric oxide synthase and Salmonella infection. Microbes Infect. 2001;3(9):771–6. pmid:11489426
- 30. Sebastiani G, Blais V, Sancho V, Vogel SN, Stevenson MM, Gros P, et al. Host Immune Response to Salmonella enterica Serovar Typhimurium Infection in Mice Derived from Wild Strains. Infect Immun. 2002;70(4):1997–2009. pmid:11895964
- 31. Roy MF, Malo D. Genetic regulation of host responses to Salmonella infection in mice. Genes Immun. 2002;3(7):381–93. pmid:12424619
- 32. Jennewein J, Matuszak J, Walter S, Felmy B, Gendera K, Schatz V, et al. Low-oxygen tensions found in Salmonella-infected gut tissue boost Salmonella replication in macrophages by impairing antimicrobial activity and augmenting Salmonella virulence. Cell Microbiol. 2015;17(12):1833–47. pmid:26104016
- 33. Yadav S, Pathak S, Sarikhani M, Majumdar S, Ray S, Chandrasekar BS, et al. Nitric oxide synthase 2 enhances the survival of mice during Salmonella Typhimurium infection-induced sepsis by increasing reactive oxygen species, inflammatory cytokines and recruitment of neutrophils to the peritoneal cavity. Free Radic Biol Med. 2018;116:73–87. pmid:29309892
- 34. Cristina Cerquetti M, Hovsepian E, Sarnacki SH, Goren NB. Salmonella enterica serovar Enteritidis dam mutant induces low NOS-2 and COX-2 expression in macrophages via attenuation of MAPK and NF-κB pathways. Microbes Infect. 2008;10(14):1431–9. pmid:18801455
- 35. Nairz M, Schleicher U, Schroll A, Sonnweber T, Theurl I, Ludwiczek S, et al. Nitric oxide–mediated regulation of ferroportin-1 controls macrophage iron homeostasis and immune function in Salmonella infection. J Exp Med. 2013;210(5):855–73. pmid:23630227
- 36. Malek M, Lamont SJ. Association of INOS, TRAIL, TGF-beta2, TGF-beta3, and IgL genes with response to Salmonella enteritidis in poultry. Genet Sel Evol GSE. 2003;35 Suppl 1:S99-111.
- 37. Tohidi R, Idris IB, Panandam JM, Bejo MH. The effects of polymorphisms in IL-2, IFN-γ, TGF-β2, IgL, TLR-4, MD-2, and iNOS genes on resistance to Salmonella enteritidis in indigenous chickens. Avian Pathol J WVPA. 2012;41(6):605–12.
- 38. He H, Arsenault RJ, Genovese KJ, Johnson C, Kogut MH. Chicken macrophages infected with Salmonella (S.) Enteritidis or S. Heidelberg produce differential responses in immune and metabolic signaling pathways. Vet Immunol Immunopathol. 2018;195:46–55. pmid:29249317
- 39. Acevedo-Villanueva K, Renu S, Gourapura R, Selvaraj R. Efficacy of a nanoparticle vaccine administered in-ovo against Salmonella in broilers. PloS One. 2021;16(4):e0247938. pmid:33822791
- 40. Breyer F, Härtlova A, Thurston T, Flynn HR, Chakravarty P, Janzen J, et al. TPL-2 kinase induces phagosome acidification to promote macrophage killing of bacteria. EMBO J. 2021;40(10):e106188. pmid:33881780
- 41. Galeev A, Suwandi A, Cepic A, Basu M, Baines JF, Grassl GA. The role of the blood group-related glycosyltransferases FUT2 and B4GALNT2 in susceptibility to infectious disease. Int J Med Microbiol IJMM. 2021;311(3):151487. pmid:33662872
- 42. Goto Y, Obata T, Kunisawa J, Sato S, Ivanov II, Lamichhane A, et al. Innate lymphoid cells regulate intestinal epithelial cell glycosylation. Science. 2014;345(6202):1254009. pmid:25214634
- 43. Suwandi A, Galeev A, Riedel R, Sharma S, Seeger K, Sterzenbach T, et al. Std fimbriae-fucose interaction increases Salmonella-induced intestinal inflammation and prolongs colonization. PLoS Pathog. 2019;15(7):e1007915. pmid:31329635
- 44. Wacklin P, Tuimala J, Nikkilä J, Sebastian Tims null, Mäkivuokko H, Alakulppi N, et al. Faecal microbiota composition in adults is associated with the FUT2 gene determining the secretor status. PloS One. 2014;9(4):e94863. pmid:24733310
- 45. Hall AB, Tolonen AC, Xavier RJ. Human genetic variation and the gut microbiome in disease. Nat Rev Genet. 2017;18(11):690–9. pmid:28824167
- 46. Waters JL, Ley RE. The human gut bacteria Christensenellaceae are widespread, heritable, and associated with health. BMC Biol. 2019;17(1):83. pmid:31660948
- 47. Kaufman J. From Chickens to Humans: The Importance of Peptide Repertoires for MHC Class I Alleles. Front Immunol. 2020;11:601089. pmid:33381122
- 48. Zhou H, Lamont SJ. Chicken MHC class I and II gene effects on antibody response kinetics in adult chickens. Immunogenetics. 2003;55(3):133–40. pmid:12743657
- 49. Liu W, Miller MM, Lamont SJ. Association of MHC class I and class II gene polymorphisms with vaccine or challenge response to Salmonella enteritidis in young chicks. Immunogenetics. 2002;54(8):582–90. pmid:12439621
- 50. Clarke TB, Francella N, Huegel A, Weiser JN. Invasive Bacterial Pathogens Exploit TLR-Mediated Downregulation of Tight Junction Components to Facilitate Translocation across the Epithelium. Cell Host Microbe. 2011;9(5):404–14. pmid:21575911
- 51. Kumar S, Lombard DB. Functions of the sirtuin deacylase SIRT5 in normal physiology and pathobiology. Crit Rev Biochem Mol Biol. 2018;53(3):311–34. pmid:29637793
- 52. Vuerich M, Mukherjee S, Robson SC, Longhi MS. Control of Gut Inflammation by Modulation of Purinergic Signaling. Front Immunol. 2020;11:1882. pmid:33072065
- 53. Zeng J, Ning Z, Wang Y, Xiong H. Implications of CD39 in immune-related diseases. Int Immunopharmacol. 2020;89(Pt A):107055. pmid:33045579
- 54. Alam MS, Costales MG, Cavanaugh C, Williams K. Extracellular adenosine generation in the regulation of pro-inflammatory responses and pathogen colonization. Biomolecules. 2015;5(2):775–92. pmid:25950510
- 55. Serrano León A, Amir Shaghaghi M, Yurkova N, Bernstein CN, El-Gabalawy H, Eck P. Single-nucleotide polymorphisms in SLC22A23 are associated with ulcerative colitis in a Canadian white cohort. Am J Clin Nutr. 2014;100(1):289–94. pmid:24740203
- 56. Sleiman M, Stevens DR, Chitirala P, Rettig J. Cytotoxic Granule Trafficking and Fusion in Synaptotagmin7-Deficient Cytotoxic T Lymphocytes. Front Immunol. 2020;11:1080. pmid:32547563
- 57. Orsini M, Di Liddo R, Valeriani F, Mancin M, D’Incà R, Castagnetti A, et al. In Silico Evaluation of Putative S100B Interacting Proteins in Healthy and IBD Gut Microbiota. Cells. 2020;9(7):E1697. pmid:32679810
- 58. Ijaz A, Veldhuizen EJA, Broere F, Rutten VPMG, Jansen CA. The Interplay between Salmonella and Intestinal Innate Immune Cells in Chickens. Pathogens. 2021;10(11):1512. pmid:34832668
- 59. Xie S, Zhang H, Matjeke RS, Zhao J, Yu Q. Bacillus coagulans protect against Salmonella enteritidis-induced intestinal mucosal damage in young chickens by inducing the differentiation of goblet cells. Poult Sci. 2022;101(3):101639. pmid:35016049
- 60. Song J, Li Q, Everaert N, Liu R, Zheng M, Zhao G, et al. Effects of inulin supplementation on intestinal barrier function and immunity in specific pathogen-free chickens with Salmonella infection. J Anim Sci. 2020;98(1):skz396. pmid:31894241
- 61. de Barros Moreira Filho AL, de Oliveira CJB, de Oliveira HB, Campos DB, Guerra RR, Costa FGP, et al. High Incubation Temperature and Threonine Dietary Level Improve Ileum Response Against Post-Hatch Salmonella Enteritidis Inoculation in Broiler Chicks. PLoS ONE. 2015;10(7):e0131474. pmid:26131553
- 62. Kim CH, Kim D, Ha Y, Cho KD, Lee BH, Seo IW, et al. Expression of mucins and trefoil factor family protein-1 in the colon of pigs naturally infected with Salmonella typhimurium. J Comp Pathol. 2009;140(1):38–42. pmid:19064270
- 63. Luise D, Motta V, Bertocchi M, Salvarani C, Clavenzani P, Fanelli F, et al. Effect of Mucine 4 and Fucosyltransferase 1 genetic variants on gut homoeostasis of growing healthy pigs. J Anim Physiol Anim Nutr. 2019;103(3):801–12. pmid:30734380
- 64. Jørgensen CB, Cirera S, Anderson SI, Archibald AL, Raudsepp T, Chowdhary B, et al. Linkage and comparative mapping of the locus controlling susceptibility towards E. COLI F4ab/ac diarrhoea in pigs. Cytogenet Genome Res. 2003;102(1‑4):157–62. pmid:14970696
- 65. Stavropoulou E, Pircalabioru GG, Bezirtzoglou E. The Role of Cytochromes P450 in Infection. Front Immunol. 2018;9:89. pmid:29445375
- 66. Bezirtzoglou EEV. Intestinal cytochromes P450 regulating the intestinal microbiota and its probiotic profile. Microb Ecol Health Dis. 2012;23. pmid:23990816
- 67. Sakamaki K, Tomonaga M, Tsukui K, Nagata S. Molecular cloning and characterization of a chromosomal gene for human eosinophil peroxidase*. J Biol Chem. 1989;264(28):16828–36. pmid:2550461
- 68. Swaggerty CL, Kogut MH, He H, Genovese KJ, Johnson C, Arsenault RJ. Differential Levels of Cecal Colonization by Salmonella Enteritidis in Chickens Triggers Distinct Immune Kinome Profiles. Front Vet Sci. 2017;4:214. pmid:29322049
- 69. La Ragione RM, Cooley WA, Velge P, Jepson MA, Woodward MJ. Membrane ruffling and invasion of human and avian cell lines is reduced for aflagellate mutants of Salmonella enterica serotype Enteritidis. Int J Med Microbiol. 2003;293(4):261–72. pmid:14503791
- 70. Guiney DG, Lesnick M. Targeting of the actin cytoskeleton during infection by Salmonella strains. Clin Immunol Orlando Fla. 2005;114(3):248–55. pmid:15721835
- 71. Villamil C, Calderon MN, Arias MM, Leguizamon JE. Validation of Droplet Digital Polymerase Chain Reaction for Salmonella spp. Quantification. Front Microbiol. 2020;11:1512. pmid:32733415
- 72. Daum LT, Barnes WJ, McAvin JC, Neidert MS, Cooper LA, Huff WB, et al. Real-Time PCR Detection of Salmonella in Suspect Foods from a Gastroenteritis Outbreak in Kerr County, Texas. J Clin Microbiol. 2002;40(8):3050–2. pmid:12149377
- 73. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinforma Oxf Engl. 2013;29(1):15–21. pmid:23104886
- 74. Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9. pmid:25260700
- 75. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. pmid:25516281
- 76. Alexa A, Rahnenfuhrer J. topGO: Enrichment Analysis for Gene Ontology [Internet]. Bioconductor version: Release (3.12); 2021. https://bioconductor.org/packages/topGO/