- Research article
- Open Access
Transcriptional analysis of porcine intestinal mucosa infected with Salmonella Typhimurium revealed a massive inflammatory response and disruption of bile acid absorption in ileum
Veterinary Research volume 47, Article number: 11 (2016)
Infected pork meat is an important source of non-typhoidal human salmonellosis. Understanding of molecular mechanisms involved in disease pathogenesis is important for the development of therapeutic and preventive strategies. Thus, hereby we study the transcriptional profiles along the porcine intestine during infection with Salmonella Typhimurium, as well as post-transcriptional gene modulation by microRNAs (miRNA). Sixteen piglets were orally challenged with S. Typhimurium. Samples from jejunum, ileum and colon, collected 1, 2 and 6 days post infection (dpi) were hybridized to mRNA and miRNA expression microarrays and analyzed. Jejunum showed a reduced transcriptional response indicating mild inflammation only at 2 dpi. In ileum inflammatory genes were overexpressed (e.g., IL-1B, IL-6, IL-8, IL1RAP, TNFα), indicating a strong immune response at all times of infection. Infection also down-regulated genes of the FXR pathway (e.g., NR1H4, FABP6, APOA1, SLC10A2), indicating disruption of the bile acid absorption in ileum. This result was confirmed by decreased high-density lipoprotein cholesterol in serum of infected pigs. Ileal inflammatory gene expression changes peaked at 2 dpi and tended to resolve at 6 dpi. Furthermore, miRNA analysis of ileum at 2 dpi revealed 62 miRNAs potentially regulating target genes involved in this inflammatory process (e.g., miR-374 and miR-451). In colon, genes involved in epithelial adherence, proliferation and cellular reorganization were down-regulated at 2 and 6 dpi. In summary, here we show the transcriptional changes occurring at the intestine at different time points of the infection, which are mainly related to inflammation and disruption of the bile acid metabolism.
Salmonella enterica subsp. enterica serovar, S. Typhimurium are gram-negative flagellated pathogenic bacteria that cause gastrointestinal disease in animals and humans . Currently, salmonellosis is ranked as the second most common zoonotic disease in the European Union, and most cases of salmonellosis in humans are associated with the consumption of contaminated pork and poultry meat [1–3]. S. Typhimurium the most commonly non-typhoidal serotype isolated from humans worldwide, causes clinical (fever, profuse diarrhea and other gastrointestinal signs) and subclinical disease in pigs [2, 3]. Therefore, prevention and control of salmonellosis in pigs is crucial not only for animal health, decrease in antibiotic use and the reduction of economic losses in the swine industry but also for minimizing the risks to public health [2, 4].
Since the gastrointestinal infection by S. Typhimurium causes similar clinical signs in humans and pigs, and given that the latter have been demonstrated to be a valuable animal model for the study of the human gastrointestinal tract , in vivo experimental infections of pigs with S. Typhimurium will likely reproduce the pathogenesis and the molecular mechanisms underlying this disease in humans. In naturally infected pigs, S. Typhimurium preferentially colonizes ileum, cecum and colon presumably due to pH, not as low in these sections as in the stomach, and a more reduced presence of bile salts than in duodenum or jejunum. Bile salts have antibacterial effects, although it has also been demonstrated that Salmonella shows resistance and tolerance to bile acids [2, 6]. To fight infection, the host defense mechanisms are activated after the adherence of S. Typhimurium to the intestinal epithelial cells. This early pro-inflammatory state can be initiated by the activation of the microbial-specific toll-like receptors (TLRs), which activate nuclear factor kappa B (NF-κB), mitogen-activated protein kinases (MAPK) and caspase-dependent signaling pathways [7, 8]. This induces the expression of inflammatory mediators (e.g., cytokines/chemokines) and antimicrobial peptides (e.g., defensins) . Later on, acquired pathogen-specific responses will be developed with the aim of clearing bacteria. Nevertheless, Salmonella is a very successful enteric pathogen that has developed different virulence strategies to evade detection by the host immune system [10, 11]. Some S. Typhimurium genes responsible of colonizing porcine intestines have been identified and characterized .
Recent studies have demonstrated the importance of certain miRNAs in the modulation of many physiological processes involved in the response to bacterial infections such as signal transduction pathways, membrane trafficking and pro-inflammatory responses [13–15]. miRNAs are small noncoding RNAs that regulate post-transcriptional expression by binding to the 3′ untranslated regions of their target messenger RNAs. It has been reported that a dysregulation of miRNAs occurs in intestinal epithelial cells in response to bacterial pathogens . Also, S. Typhimurium can alter miRNA expression by TLR-independent mechanisms such as secretion of effector proteins . Therefore, a more comprehensive view of the miRNA-mediated regulation of mRNA expression is needed to better understand the gastrointestinal response to invading pathogens .
Although some studies have focused on transcriptional changes of either a reduced number of genes [18, 19] or specific intestinal sections [9, 20], to our knowledge there is limited information about the early transcriptional response to S. Typhimurium infection at the different anatomical portions of the porcine gut. Today, the use of whole-genome approaches such as microarray expression profiling has allowed an unprecedented look about the function of genes and their role in disease . Therefore, in order to perform a comprehensive evaluation of the porcine intestinal response to S. Typhimurium infection, the objective of this work was to investigate the transcriptional profile of different portions of the gut using an in vivo model of Salmonella infection. In addition, the role of miRNAs as post-transcriptional modulators of this immune response was also evaluated.
Materials and methods
Experimental infection and sample processing
Sixteen male and female crossbreed weaned piglets, approximately 4 weeks of age, were used in this study. All piglets were derived from a Salmonella-negative herd and were serologically negative. Pigs were housed in an environmentally controlled isolation facility at 25 °C and under constant light with ad libitum access to feed and water. After an acclimation period of 5 days, four piglets were necropsied (control group), 2 h prior to experimental infection of the other animals. Then, 12 piglets were challenged orally with 108 colony forming units (cfu) of a S. Typhimurium phagetype DT104 strain isolated from a carrier pig . Fever, lethargy and diarrhea were monitored every day. Four randomly chosen infected pigs were necropsied at 1, 2 and 6 days post infection (dpi). Tissue samples were aseptically collected and stored in liquid nitrogen. Fecal samples were collected for bacteriological cultures the day of arrival and the day when the piglets were necropsied. Fecal sample processing and bacteriological analysis was performed following the current EN-ISO standard methodology 6579:2002/Amd 1:2007. Serum samples were obtained from each animal before euthanasia: blood was collected and placed into non-anticoagulated tubes, letting it clot and centrifuging. Obtained serum samples were sent for analysis to an external laboratory (Laboratorio Veterinario Garfia S.L., Cordoba, Spain), from where the following measurements were obtained: total proteins, albumin, blood urea nitrogen (BUN), creatinine, aspartate aminotransferase (AST), alanine transaminase (ALT), alkaline phosphatase, total cholesterol, high density lipoprotein (HDL) cholesterol, low density lipoprotein (LDL) cholesterol, triglycerides, glucose, immunoglobulins (G, A and M), complement component C3, and Salmonella spp. antibodies. Additional information regarding detection methods is included in Additional file 1. Descriptive statistics and normality tests were followed by ANOVA using Dunnet’s post-test.
Sections from jejunum, ileum and colon were collected, sectioned into pieces of around 10 cm and immediately frozen in liquid nitrogen for mucosa isolation and RNA purification. All procedures involving animals were approved by the institutional bioethical committee, and performed according to European regulations regarding animal welfare and protection of animals used for experimental and other scientific purposes.
For mucosa isolation and RNA purification, intestinal tissue samples stored at −80 °C were treated with RNAlater ®-ICE (Ambion, Inc, Austin, TX, USA) and cut into 2 cm pieces, according to manufacturer's instructions. Ileum mucosa was scraped from the intestinal luminal surface with a razor, and was immediately disrupted and homogenized in RLT buffer (RNeasy Mini Kit, QIAGEN, Valencia, CA, USA) using a rotor–stator homogenizer. Further RNA extraction was done using the RNeasy Mini Kit according to manufacturer instructions. For miRNA studies, RNA from ileum (0 and 2 dpi) was isolated using mirVana miRNA isolation kit (Ambion, Inc, Austin, TX, USA). Eluted RNA was treated with DNase using TURBO DNA-free™ Kit (Ambion, Inc, Austin, TX, USA). RNA integrity was assessed in the Agilent Bioanalyzer 2100 (Agilent Technologies, Palo Alto, CA, USA). Only samples with RNA integrity numbers (RIN) ≥7 were used for further analysis.
Microarray hybridization and analysis
Gene expression analysis was carried out using the GeneChip Porcine Genome Array (Affymetrix, Inc., Santa Clara, CA, USA) at the Unidad Científico-Técnica de Apoyo (UCTS) of the Institut de Recerca del Hospital Universitario Vall d’Hebron (Barcelona, Spain). The One-Cycle Eukaryotic Target Labeling Assay (Expression Analysis Technical Manual, Affymetrix, Inc., Santa Clara, CA, USA) was used to obtain biotinylated cRNA from individual mucosal mRNA samples. Then, they were hybridized to the GeneChip Porcine Genome Array and processed using manufacturer’s instructions. Data analysis was conducted using in-house algorithms in R (v. 2.7.0). Quality control analysis of the mRNA array was performed using the robust multi-array analysis (RMA)  included in the affy library of Bioconductor package . Differentially expressed (DE) genes were obtained by paired comparisons using the limma package (moderated t test of linear models after an empirical Bayes correction ); three independent comparisons were carried out for each combination of control vs. infected samples (i.e., 0 dpi 1 dpi, 0 dpi 2 dpi and 0 dpi 6 dpi). Only genes with a fold-change (FC) >1.5 or <1.5 were considered for further investigation. Due to the lack of complete annotation of the GeneChip Porcine Genome Array, the entire dataset was re-annotated using Blast2GO .
Samples for miRNA analysis were hybridized to the Human miRNA Microarray (V3) 8 × 15 K (Agilent Technologies, Inc., Santa Clara, CA, USA) and processed at the Andalusian Center for Molecular Biology and Regenerative Medicine (CABIMER) Genomics Core Facility. For the miRNA microarray data, background correction  and quantile normalization [28, 29] were performed. Differential expression of miRNAs was calculated using the RankProd method , which is a non-parametric method based on the estimated percentage of false predictions (PFP). P values were adjusted for multiple testing using the Benjamini and Hochberg method for false discovery rate , and adjusted P < 0.05 were considered to be statistically significant. As in the mRNA array, miRNA FC threshold was set at 1.5.
Salmonella detection in tissues
Immunohistochemical analysis of intestinal tissue samples of jejunum, ileum and colon at 0, 1, 2 and 6 dpi was performed as previously described, using a S. Typhimurium specific antibody .
Systems biology analysis
Functional analysis of DE genes was carried out using Ingenuity Pathway Analysis (IPA, Ingenuity Systems® Inc, Redwood City, CA, USA). Genes differentially expressed in each intestine section/time point were uploaded into IPA and analyzed separately. Obtained results included biological functions and canonical pathways, which were filtered by setting a threshold of P < 0.05.
Quantitative real-time PCR (qPCR)
For mRNA expression analysis, a panel of 16 selected genes involved in bile acid absorption was assayed by qPCR as previously reported  using gene-specific primers (Additional file 2). PCR conditions were: 5 min at 95 °C followed by 35 cycles of 30 s at 94 °C, 30 s at 57 °C and 45 s at 72 °C. Melting curve analyses were performed to ensure specificity of each assay. Two reference genes, beta-actin and cyclophilin A, were used to normalize mRNA expression values.
For the miRNA expression analysis, qPCR of six differentially expressed miRNAs was performed. Briefly, 100 ng of total RNA per animal was reverse transcribed to cDNA as previously reported [33, 34] and diluted 1:8 times. The 10 µL final PCR reaction mix contained 1 µL of cDNA, 5 µL of PerfeCTa™ SYBR® Green Supermix for iQ™ (Quanta BioSciences, Inc.), and 10 µM of each primer. Cycling conditions were 10 min at 95 °C followed by 40 cycles of 5 s at 95 °C, and 60 s at 60 °C; a final melting curve analysis was performed (60–99 °C). The miRNA-specific primers were designed according to guidelines set by Balcells et al. , and using publicly available software miRprimer  (Additional file 2). After evaluation of its stability by geNorm v3.5 algorithm , five miRNAs (ssc-miR-26a, ssc-let-7a, ssc-miR-103, ssc-miR-17-5p and ssc-miR-16-5p) were used as reference to normalize expression . Relative gene expression was measured in controls and infected pigs at the three times post infection, and expression ratios were calculated according to the 2−ΔΔCt method . Also, miRNA relative expression at 2 dpi (versus control, 0 dpi) was calculated using the same method.
Statistical differences in expression values among groups were assessed using a Kruskal–Wallis test (mRNA) or Student’s t test (miRNA) (Graphpad Prism 6, Graphpad Software Inc, La Jolla, CA, USA). Statistical significance was set at P < 0.05.
All infected animals tested positive for S. Typhimurium in feces, and developed clinical signs characteristic of the disease such as fever (a peak of fever at 2 dpi returning to normal values at 6 dpi), lethargy and diarrhea . Control animals tested negative to S. Typhimurium in feces prior to their necropsy. The sera biochemistry profile along the time course of infection (Additional file 1) showed decreased albumin (at 1 dpi, P = 0.039), total proteins (at 1 dpi, P = 0.009), glucose (at 2 dpi, P = 0.024), high density lipoprotein cholesterol (HDL, at 1 dpi, P = 0.046; 2 dpi, P = 0.003 and 6 dpi, P = 0.041), alanine transaminase (ALT, at 1 dpi, P = 0.030 and 6 dpi, P = 0.001) and porcine IgM (at 2 dpi, P = 0.012 and 6 dpi, P = 0.008).
Complete microarray results can be found in Additional file 3. In jejunum, expression changes were found only in day 2 after infection (2 dpi), and most genes were up-regulated (71%). Expression changes in colon occurred at 2 and 6 dpi, and most genes were found to be down-regulated (70 and 73%, respectively). Ileum was the gut portion where the vast majority of DE genes were found (over 2300 different genes), and the proportion of up- and down-regulated genes was similar at all points of the infection time course (Figure 1).
Ileal transcriptomic response during Salmonella infection
Salmonella colonization observed in ileum at 1 and 2 dpi (Figure 2) is in agreement with previous reports [32, 39]. Functional analysis of DE genes revealed that affected biological pathways included inflammation and immune response, lipid metabolism and cell death and survival (Figure 3; Additional file 4). Immune-related biological pathways were altered at all times of the infection time course, showing a chronological progress of events. The infection triggered an antimicrobial response at 1 dpi, and later the immune response was replaced by an alteration of the immune cell trafficking and inflammatory response (more pronounced at 1 and 2 dpi), cell-to-cell signaling and interaction, infectious disease and cell-mediated immune response. Genes affecting lymphoid tissue structure were impaired at all times, but changes were more evident at 6 dpi, where humoral response was also observed. This differential expression implicated dysregulation of inflammatory/immune processes such as IL-6 signaling pathway, LPS/IL-1 mediated inhibition of RXR function, granulocyte adhesion and diapedesis, IL-10 signaling, differential regulation of cytokine production in intestinal epithelial cells by IL-17A and IL-17F and IL-12 signaling and production in macrophages (Figure 3; Additional file 4). At 6 dpi, we observed a dysregulation of pathways related to proliferation and reposition of damaged tissue (Figure 3).
Along with inflammation and immune response, lipid metabolism and its related functions (e.g., free radical scavenging, endocrine system development and function) were highly impaired at 1 and 2 dpi after experimental infection, dysregulating retinoid X receptor (RXR) related pathways such as LPS/IL-1 mediated inhibition of RXR function, LXR/RXR activation, FXR/RXR activation and hepatic cholestasis. Genes showing the largest transcriptional changes at 1 and 2 dpi (Additional file 3) are involved in these pathways. DE genes described in Table 1 are involved in bile acid metabolism and intestinal absorption of bile through the ileum mucosa (e.g., FABP2, a fatty acid transporter, and FABP6, an ileum-specific bile acid transporter). In order to confirm the disruption of the normal lipid/bile acid absorption pathway in ileum due to S. Typhimurium infection and its impact in overall homeostasis, expression of genes involved in this pathway were further studied using qPCR (Table 1). At 2 dpi, we confirmed up-regulation of IL-1B, IL-6, TLR2, TLR4, TNFα and PPARG, along with down-regulation of ASBT, FABP6, FABP2, FXR, RXRG and APOA1. The overall gene dysregulation found in this pathway tended to resolve at 6 dpi (Table 1), and some regulatory genes (FXR and RXRG) even changed direction of expression (from down to up-regulated).
Transcriptomic response in jejunum and colon during Salmonella infection
Jejunum showed a reduced transcriptional response, with absence of response at 1 and 6 dpi (Figure 1). At 2 dpi we observed scarce changes in genes involved in inflammatory response (e.g., CXCL2 overexpression) but with few consequences in inflammatory/immune response signaling pathways (Additional file 4). In colon, we found no transcriptional changes at 1 dpi, while at 2 and 6 dpi a general down-regulation (70% of genes) of cellular proliferation pathways was found, especially actin-based Rho signaling pathways (Additional file 3).
miRNA regulation of immune response against Salmonella in ileum at 2 dpi
Based on the highest misregulation of gene expression (Additional file 3) and bacterial colonization found in ileum at 2 dpi (Figure 2), the miRNA expression profile was investigated at this time point. A total of 62 miRNAs (FC ≥ 1.5, P < 0.05) were found DE in ileum after S. Typhimurium infection, from which 37 were up-regulated and 25 down-regulated (Additional file 5). The prediction of their target genes revealed that these 62 miRNAs are potential regulators of 880 genes (Additional file 6); these genes are involved in many biological functions such as cellular growth and proliferation, cell death and survival, inflammatory response, immune cell trafficking and gastrointestinal disease (Additional file 7). Although we validated the array results by qPCR (Table 2), we could observe that in general miRNA expression values were very moderate compared to mRNA results, and only miR-451 was found to be statistically significant and biologically meaningful (FC > 2, P < 0.001).
The present experimental time course of infection demonstrated that S. Typhimurium infection in pigs lead to an inflammatory response and activation of immune mechanisms, as shown by massive transcriptional dysregulation. These changes were more pronounced at day 1 and day 2 after experimental infection, and tended to disappear by day 6 after bacterial challenge. Most gene expression changes occurred in the ileum, indicating an early immune response to infection (1 and 2 dpi) followed by cell proliferation (6 dpi).
Activation of TLR2 and TLR4 along the time course of the experimental infection triggered early intestinal (innate) immune response at 1 and 2 dpi. In parallel, initial exposure to pathogens also induced secretion of antimicrobial substances (at 1 and 2 dpi) such as antimicrobial peptides (AMP) and production of reactive oxygen species (ROS) and reactive nitrogen species (RNS) at the site of infection. Antimicrobial elements differentially expressed in this study included S100A9 and S100A8 (highly over-expressed), as well as DEFB1, EDN, S100A12 and LTF, showing the intense local response occurring in the porcine gut after S. Typhimurium infection (in agreement with previous reports ). Stimulation of TLR by bacterial recognition activated cytokine cascades such as IL-6 and IL-10 signaling pathways predominantly at 1 and 2 dpi, which lead to recruitment of other immune cells (e.g., neutrophils, dendritic cells) to the infection site in order to clear infection. This early inflammatory response disappeared at 6 dpi, and it was replaced by tissue regeneration and proliferation. Interestingly, and even though TLR2 and TLR4 expression remained elevated, the host inflammatory response weakened at 6 dpi. The role of TLR2, TLR4 and TLR9 in S. Typhimurium infection has been recently characterized by Arpaia et al. , who found that a failure in host cell TLR recognition of the bacteria impairs formation of the Salmonella-containing vacuole (SCV) and activation of the Salmonella Pathogenicity Island-2 genes .
In the present study, we observed that intestinal inflammation at 1 and 2 dpi was associated with a down-regulation of the farnesoid X receptor (FXR, NR1H4), increasing expression of NF-κB dependent genes (IL-1A, IL-1B, ILRN, IL1R1, IL1RAP, CD14, IL-33, TNFα, TNFRSF1B, IL-6) in agreement with previous reports [40, 42]. FXR, also called bile acid (BA) receptor, is a nuclear receptor that locally modulates intestinal immune response via regulation of cholesterol and BA metabolism. It has been described that FXR forms a heterodimer with the retinoid X factor (RXR), maintaining BA homeostasis in gut and liver . Under physiological conditions, bile acids are synthetized in the liver and secreted into the intestine (i.e., duodenum) for digestion and absorption of dietary fat. Since over 90% of secreted BA are reabsorbed through the ileum to be recycled , FXR has an important role in regulation of the intestinal absorption for protecting the cells from biliary damage. Most inflammatory NF-kB dependent genes were not differentially expressed at 6 dpi, which concurs with a FXR change to overexpression. Therefore, our findings support previous studies , claiming that modulation of the inflammatory process mediated by FXR prevents further tissue damage and avoids disease progression. Intracellular FXR activation by presence of BA in ileocytes has been described to inhibit further absorption by up-regulation of SHP (NR0B2), which represses the apical enterocyte BA transporter ASBT (SLC10A2) limiting additional entrance, and inhibits CYP7A1, the hepatic enzyme that synthetizes BA from cholesterol in the liver . We report here that S. Typhimurium infection in pigs shut down the FXR pathway, subsequently down-regulating FXR target genes; these results are in agreement with previous murine studies . We observed that FXR repression at 1 and 2 dpi (along with NR2B3 down-regulation, its heterodimer partner) during the early inflammatory response lead to a decreased expression of ASBT and the fatty and bile acid transporters FABP2 and FABP6, impairing normal bile absorption in ileum. Repression of ileal BA transporters in the S. Typhimurium infected gut, contrary to the physiological regulatory mechanism, has been associated with decreased BA production in the liver (i.e., down-regulated CYP7A1) . To our knowledge, changes in FXR expression had not been previously reported in pigs after S. Typhimurium experimental infection.
Down-regulation of the RXR-LXR pathway in ileum was evident at 1 and 2 dpi, with repression of its target genes ABCG8, APOA1, APOC3, and LPL. These findings indicate deficiency in cholesterol absorption, and are supported by serum concentrations of HDL cholesterol in the studied pigs, which was significantly lower during the infection time course, especially at 2 dpi, demonstrating a disruption in cholesterol carriage from tissues to the liver. Alteration of lipid metabolism (BA and cholesterol) is tightly linked to the inflammatory response triggered by the S. Typhimurium infection, since activation of the IL-1 signaling cascade down-regulates ASBT , limiting BA absorption and therefore FXR-mediated transcription of the target genes previously mentioned.
After the initial period of non-specific immune response, the acquired pathogen-specific response is activated in order to clear bacteria (although this mechanism can fail, leading to persistence). Antigen presenting cells (e.g., dendritic cells, macrophages) will then stimulate T cells into different types, each of them pathogen- or toxin-specific. It has been shown that Salmonella infection decreases MHC class II molecules’ expression (swine leucocyte antigen, SLA) by inducing polyubiquitination of SLA-DR in infected cells, limiting pathogen recognition . Our data are in agreement with that, indicating a late overstimulation of the antigen presenting function in ileum at 6 dpi, demonstrated by up-regulation of SLA-DRA, SLA-DRB2, SLA-DRB4, SLA-B, SLA-DQA1, and SLA-DQB2. In parallel, we found a down-regulation of B-cell differentiation and function suggesting impairment of a proper humoral immune response, as indicated by down-regulation of molecules involved in B-cell pathways such as FOXO1, CD19, BLNK and EBF1 [49–52].
Transcriptional changes found in jejunum (e.g., overexpression of CXCL2, S100A9), which occurred only at 2 dpi, are product of the early inflammatory/acute phase response. This innate inflammatory response to Salmonella has been previously reported using experimental jejunal loop infection model in pigs, where overexpression of inflammatory molecules was detected from 2 to 8 h after Salmonella perfusion . Colonic response to S. Typhimurium infection started at 2 dpi, and was still present at 6 dpi. Our results indicate a general repression of pathways related to cellular proliferation and disruption of cell junctions, which has also been reported in mouse colon 4 days post Salmonella infection . Gastrointestinal epithelial turnover is a needed physiological process that helps maintaining gut homeostasis; although it can be accelerated due to cell injury , it has been demonstrated that bacterial infections counteract the renewal of epithelial cells due to the effect of bacterial proteins on different host genes . Epithelial cell turnover is mainly triggered by cell oxidative burst that occurs in intestinal epithelial cells during infections (mainly in intestinal stem cells, located in the crypts and in charge of turnover), due to JAK-STAT and JNK pathways’ stimulation . In our study, we detected a down-regulation of LIMS2 at 2 and 6 dpi. This gene encodes an integrin-linked kinase (ILK) binding protein. S. Typhimurium EspO1 protein deregulates cell shedding by acting on host ILK , therefore the down-regulation of LIMS2 could be explained by this phenomenon. Also, we found down-regulation of tissue remodeling-related genes (especially those specific of tight junctions, such as ACTC1, ACTG2, and ACTA1), which can be explained by the action of Salmonella pathogenicity island 1 effector proteins (e.g., SopB, SopE, SopE2 and SipA) that disrupt cell–cell junctions . AvrA has been shown to have an anti-inflammatory effect, acting on mitogen-activated protein kinase kinases (MAPKKs) , explaining the down-regulation of MAPKAP1 (MAPK associated protein 1) found at 2 dpi.
Our microarray results on miRNAs revealed a total of 62 miRNAs differentially expressed in ileum 2 days after S. Typhimurium infection. Assessment of some of the best candidates by qPCR confirmed the microarray results but the fold changes obtained by qPCR were very modest compared with the fold changes found in the microarray data, with the exception of miR-451 which showed highly significant differential expression. In general, changes in expression levels are much more moderate in miRNAs compared to mRNAs. Nevertheless, several miRNAs can target the same mRNA, working cooperatively and making these fold changes accumulative. Moreover, target prediction analysis for the 62 miRNAs found DE in the microarray study revealed 880 genes. These target genes are mainly involved in biological functions such as cellular growth and proliferation, cell death, inflammatory response, immune cell trafficking and gastrointestinal disease, confirming the findings at the mRNA level. Agreeing with our present study, it has been previously described that let-7b overexpression is induced by NF-kB activation . Similarly, we observed miR-374a up-regulation, which has been described as CEBPB regulator , controlling the expression of IL-6, IL-8 and other acute phase inflammatory genes during S. Typhimurium infection . We also found miR-451 overexpression in ileum, in agreement with a recent study of porcine blood miRNA profile after S. Typhimurium infection . Target prediction analysis indicated that this miRNA controls the expression of ATP-binding cassette B1 (ABCB1) [63, 64], highly expressed in the apical surface of epithelial cell in ileum, where it contributes to the luminal efflux of cholesterol . Repression of miR-451 has been associated to ABCB1 up-regulation, impairing S. Typhimurium ability to invade host cells by reducing adhesion to epithelial cells [66, 67]; on the other hand, ABCB1 down-regulation (as in the present study) is associated with inflammatory reaction (TNF activation) in the gut in response to bacterial infections [66, 68]. Additionally, it has been shown that overexpression of ABCB1 is frequently associated with an increase in intracellular pH, therefore down-regulation of this gene could cause a decrease in intracellular pH, which is beneficial to intracellular survival of Salmonella .
To conclude, our biological approximation showed a global deregulation of a number of genes upon S. Typhimurium infection along the porcine gut. Our results reveal that jejunum, ileum and colon respond differently to infection. We observed a slight transcriptional response in jejunum, affected only at 2 dpi. The majority of immune response was observed in ileum, where we found a high inflammatory reaction and a repression of the genes involved in bile acid absorption and metabolism, although this tended to resolve at 6 dpi. In colon, we found down-regulation at 2 and 6 dpi of some genes involved in Salmonella invasion pathways, actin filaments organization and signaling by Rho family GTPases. Our results provide also information about the role of a number of miRNAs in the regulation of bacteria induced immune and inflammatory responses in porcine ileum.
European Food Safety Authority (EFSA) (2013) The European Union summary report on trends and sources of zoonoses, zoonotic agents and food-borne outbreaks in 2011. EFSA J 11:3129–3379
Boyen F, Haesebrouck F, Maes D, Van Immerseel F, Ducatelle R, Pasmans F (2008) Non-typhoidal Salmonella infections in pigs: a closer look at epidemiology, pathogenesis and control. Vet Microbiol 130:1–19
Foley SL, Lynne AM, Nayak R (2008) Salmonella challenges: prevalence in swine and poultry and potential pathogenicity of such isolates. J Anim Sci 86:E149–162
Barton MD (2014) Impact of antibiotic use in the swine industry. Curr Opin Microbiol 19:9–15
Zhang Q, Widmer G, Tzipori S (2013) A pig model of the human gastrointestinal tract. Gut Microbes 4:193–200
van Velkinburgh JC, Gunn JS (1999) PhoP-PhoQ-regulated loci are required for enhanced bile resistance in Salmonella spp. Infect Immun 67:1614–1622
Akira S, Takeda K (2004) Toll-like receptor signalling. Nat Rev Immunol 4:499–511
Hulst M, Smits M, Vastenhouw S, de Wit A, Niewold T, van der Meulen J (2013) Transcription networks responsible for early regulation of Salmonella-induced inflammation in the jejunum of pigs. J Inflamm (Lond) 10:18
Zhou R, O’Hara SP, Chen XM (2011) MicroRNA regulation of innate immune responses in epithelial cells. Cell Mol Immunol 8:371–379
Grassl GA, Finlay BB (2008) Pathogenesis of enteric Salmonella infections. Curr Opin Gastroenterol 24:22–26
Ravindran R, McSorley SJ (2005) Tracking the dynamics of T-cell activation in response to Salmonella infection. Immunology 114:450–458
Carnell SC, Bowen A, Morgan E, Maskell DJ, Wallis TS, Stevens MP (2007) Role in virulence and protective efficacy in pigs of Salmonella enterica serovar typhimurium secreted components identified by signature-tagged mutagenesis. Microbiology 153:1940–1952
Eulalio A, Schulte L, Vogel J (2012) The mammalian microRNA response to bacterial infections. RNA Biol 9:742–750
Witwer KW, Sisk JM, Gama L, Clements JE (2010) MicroRNA regulation of IFN-beta protein expression: rapid and sensitive modulation of the innate immune response. J Immunol 184:2369–2376
Gantier MP, Stunden HJ, McCoy CE, Behlke MA, Wang D, Kaparakis-Liaskos M, Sarvestani ST, Yang YH, Xu D, Corr SC, Morand EF, Williams BR (2012) A miR-19 regulon that controls NF-κB signaling. Nucleic Acids Res 40:8048–8058
Maudet C, Mano M, Sunkavalli U, Sharan M, Giacca M, Forstner KU, Eulalio A (2014) Functional high-throughput screening identifies the miR-15 microRNA family as cellular restriction factors for Salmonella infection. Nat Commun 5:4718
Bao H, Kommadath A, Plastow GS, Tuggle CK, Guan le L, Stothard P (2014) MicroRNA buffering and altered variance of gene expression in response to Salmonella infection. PLoS One 9:e94352
Collado-Romero M, Arce C, Ramirez-Boo M, Carvajal A, Garrido JJ (2010) Quantitative analysis of the immune response upon Salmonella typhimurium infection along the porcine intestinal gut. Vet Res 41:23
Meurens F, Berri M, Auray G, Melo S, Levast B, Virlogeux-Payant I, Chevaleyre C, Gerdts V, Salmon H (2008) Early immune response following Salmonella entericasubspecies entericaserovar typhimurium infection in porcine jejunal gut loops. Vet Res 40:5
Hoeke L, Sharbati J, Pawar K, Keller A, Einspanier R, Sharbati S (2013) Intestinal Salmonella typhimurium infection leads to miR-29a induced caveolin 2 regulation. PLoS One 8:e67300
Hossain H, Tchatalbachev S, Chakraborty T (2006) Host gene expression profiling in pathogen-host interactions. Curr Opin Immunol 18:422–429
Garcia-Feliz C, Collazos JA, Carvajal A, Vidal AB, Aladuena A, Ramiro R, de la Fuente M, Echeita MA, Rubio P (2007) Salmonella enterica infections in Spanish swine fattening units. Zoonoses Public Health 54:294–300
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP (2003) Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 4:249–264
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J (2004) Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 5:R80
Smyth GK (2004) Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 3:Article3
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M (2005) Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21:3674–3676
Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, Smyth GK (2007) A comparison of background correction methods for two-colour microarrays. Bioinformatics 23:2700–2707
Pradervand S, Weber J, Thomas J, Bueno M, Wirapati P, Lefort K, Dotto GP, Harshman K (2009) Impact of normalization on miRNA microarray expression profiling. RNA 15:493–501
Zhao Y, Wang E, Liu H, Rotunno M, Koshiol J, Marincola FM, Landi MT, McShane LM (2010) Evaluation of normalization methods for two-channel microRNA microarrays. J Transl Med 8:69
Hong F, Breitling R, McEntee CW, Wittner BS, Nemhauser JL, Chory J (2006) RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis. Bioinformatics 22:2825–2827
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I (2001) Controlling the false discovery rate in behavior genetics research. Behav Brain Res 125:279–284
Collado-Romero M, Martins RP, Arce C, Moreno A, Lucena C, Carvajal A, Garrido JJ (2012) An in vivo proteomic study of the interaction between Salmonella typhimurium and porcine ileum mucosa. J Proteomics 75:2015–2026
Balcells I, Cirera S, Busk PK (2011) Specific and sensitive quantitative RT-PCR of miRNAs with DNA primers. BMC Biotechnol 11:70
Podolska A, Anthon C, Bak M, Tommerup N, Skovgaard K, Heegaard PM, Gorodkin J, Cirera S, Fredholm M (2012) Profiling microRNAs in lung tissue from pigs infected with Actinobacillus pleuropneumoniae. BMC Genomics 13:459
Busk PK (2014) A tool for design of primers for microRNA-specific quantitative RT-qPCR. BMC Bioinformatics 15:29
Mestdagh P, Van Vlierberghe P, De Weer A, Muth D, Westermann F, Speleman F, Vandesompele J (2009) A novel and universal method for microRNA RT-qPCR data normalization. Genome Biol 10:R64
Timoneda O, Balcells I, Cordoba S, Castello A, Sanchez A (2012) Determination of reference microRNAs for relative quantification in porcine tissues. PLoS One 7:e44413
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25:402–408
Meyerholz DK, Stabel TJ, Ackermann MR, Carlson SA, Jones BD, Pohlenz J (2002) Early epithelial invasion by Salmonella enterica serovar typhimurium DT104 in the swine ileum. Vet Pathol 39:712–720
Bao H, Kommadath A, Liang G, Sun X, Arantes AS, Tuggle CK, Bearson SM, Plastow GS, Stothard P, le Guan L (2015) Genome-wide whole blood microRNAome and transcriptome analyses reveal miRNA-mRNA regulated host response to foodborne pathogen Salmonella infection in swine. Sci Rep 5:12620
Arpaia N, Godec J, Lau L, Sivick KE, McLaughlin LM, Jones MB, Dracheva T, Peterson SN, Monack DM, Barton GM (2011) TLR signaling is required for Salmonella typhimurium virulence. Cell 144:675–688
Vavassori P, Mencarelli A, Renga B, Distrutti E, Fiorucci S (2009) The bile acid receptor FXR is a modulator of intestinal innate immunity. J Immunol 183:6251–6261
Eloranta JJ, Kullak-Ublick GA (2008) The role of FXR in disorders of bile acid homeostasis. Physiology (Bethesda) 23:286–295
Guyton AC, Hall JE (2006) Secretory functions of the alimentary tract. In: Guyton AC, Hall JE (eds) Textbook of medical physiology. Elsevier Saunders, Philadelphia, pp 791–807
Matsubara T, Li F, Gonzalez FJ (2013) FXR signaling in the enterohepatic system. Mol Cell Endocrinol 368:17–29
Romain G, Tremblay S, Arena ET, Antunes LCM, Covey S, Chow MT, Finlay BB, Menendez A (2013) Enterohepatic bacterial infections dysregulate the FGF15-FGFR4 endocrine axis. BMC Microbiol 13:238
Neimark E, Chen F, Li X, Magid MS, Alasio TM, Frankenberg T, Sinha J, Dawson PA, Shneider BL (2006) c-Fos is a critical mediator of inflammatory-mediated repression of the apical sodium-dependent bile acid transporter. Gastroenterology 131:554–567
Van Parys A, Boyen F, Verbrugghe E, Leyman B, Bram F, Haesebrouck F, Pasmans F (2012) Salmonella typhimurium induces SPI-1 and SPI-2 regulated and strain dependent downregulation of MHC II expression on porcine alveolar macrophages. Vet Res 43:52
Szydlowski M, Jablonska E, Juszczynski P (2014) FOXO1 transcription factor: a critical effector of the PI3 K-AKT axis in B-cell development. Int Rev Immunol 33:146–157
Hasegawa M, Fujimoto M, Poe JC, Steeber DA, Lowell CA, Tedder TF (2001) A CD19-dependent signaling pathway regulates autoimmunity in Lyn-deficient mice. J Immunol 167:2469–2478
Minegishi Y, Rohrer J, Coustan-Smith E, Lederman HM, Pappu R, Campana D, Chan AC, Conley ME (1999) An essential role for BLNK in human B cell development. Science 286:1954–1957
Vilagos B, Hoffmann M, Souabni A, Sun Q, Werner B, Medvedovic J, Bilic I, Minnich M, Axelsson E, Jaritz M, Busslinger M (2012) Essential role of EBF1 in the generation and function of distinct mature B cell types. J Exp Med 209:775–792
Liu X, Lu R, Xia Y, Sun J (2010) Global analysis of the eukaryotic pathways and networks regulated by Salmonella typhimurium in mouse intestinal infection in vivo. BMC Genomics 11:722
Pitsouli C, Apidianakis Y, Perrimon N (2009) Homeostasis in infected epithelia: stem cells take the lead. Cell Host Microbe 6:301–307
Kim M, Ashida H, Ogawa M, Yoshikawa Y, Mimuro H, Sasakawa C (2010) Bacterial interactions with the host epithelium. Cell Host Microbe 8:20–35
Buchon N, Broderick NA, Poidevin M, Pradervand S, Lemaitre B (2009) Drosophila intestinal response to bacterial infection: activation of host defense and stem cell proliferation. Cell Host Microbe 5:200–211
Kim M, Ogawa M, Fujita Y, Yoshikawa Y, Nagai T, Koyama T, Nagai S, Lange A, Fassler R, Sasakawa C (2009) Bacteria hijack integrin-linked kinase to stabilize focal adhesions and block cell detachment. Nature 459:578–582
Boyle EC, Brown NF, Finlay BB (2006) Salmonella enterica serovar typhimurium effectors SopB, SopE, SopE2 and SipA disrupt tight junction structure and function. Cell Microbiol 8:1946–1957
Jones RM, Wu H, Wentworth C, Luo L, Collier-Hyams L, Neish AS (2008) Salmonella AvrA coordinates suppression of host immune and apoptotic defenses via JNK pathway blockade. Cell Host Microbe 3:233–244
Teng GG, Wang WH, Dai Y, Wang SJ, Chu YX, Li J (2013) Let-7b is involved in the inflammation and immune responses associated with Helicobacter pylori infection by targeting Toll-like receptor 4. PLoS One 8:e56709
Pan S, Zheng Y, Zhao R, Yang X (2013) miRNA-374 regulates dexamethasone-induced differentiation of primary cultures of porcine adipocytes. Horm Metab Res 45:518–525
Huang TH, Uthe JJ, Bearson SM, Demirkale CY, Nettleton D, Knetter S, Christian C, Ramer-Tait AE, Wannemuehler MJ, Tuggle CK (2011) Distinct peripheral blood RNA responses to Salmonella in pigs differing in Salmonella shedding levels: intersection of IFNG TLR and miRNA pathways. PLoS One 6:e28768
Zhu H, Wu H, Liu X, Evans BR, Medina DJ, Liu CG, Yang JM (2008) Role of MicroRNA miR-27a and miR-451 in the regulation of MDR1/P-glycoprotein expression in human cancer cells. Biochem Pharmacol 76:582–588
Kovalchuk O, Filkowski J, Meservy J, Ilnytskyy Y, Tryndyak VP, Chekhun VF, Pogribny IP (2008) Involvement of microRNA-451 in resistance of the MCF-7 breast cancer cells to chemotherapeutic drug doxorubicin. Mol Cancer Ther 7:2152–2159
Le May C, Berger JM, Lespine A, Pillot B, Prieur X, Letessier E, Hussain MM, Collet X, Cariou B, Costet P (2013) Transintestinal cholesterol excretion is an active metabolic process modulated by PCSK9 and statin involving ABCB1. Arterioscler Thromb Vasc Biol 33:1484–1493
Siccardi D, Mumy KL, Wall DM, Bien JD, McCormick BA (2008) Salmonella enterica serovar typhimurium modulates P-glycoprotein in the intestinal epithelium. Am J Physiol Gastrointest Liver Physiol 294:G1392–1400
Crowe A (2011) The role of P-glycoprotein and breast cancer resistance protein (BCRP) in bacterial attachment to human gastrointestinal cells. J Crohns Colitis 5:531–542
Belliard AM, Lacour B, Farinotti R, Leroy C (2004) Effect of tumor necrosis factor-alpha and interferon-gamma on intestinal P-glycoprotein expression, activity, and localization in Caco-2 cells. J Pharm Sci 93:1524–1536
Chakravortty D, Rohde M, Jäger L, Deiwick J, Hensel M (2005) Formation of a novel surface structure encoded by Salmonella Pathogenicity Island 2. EMBO J 24:2043–2052
JHU and MCR carried out the molecular genetic studies. JHU, MCR and SZL performed data analysis and interpretation of results, and drafted the manuscript. AC performed the experimental infection, and AC, MCR and CA collected and processed the tissue samples. RB and MGC designed and carried out the statistical analyses of both microarrays. SC designed the miRNA validation assays and participated in the interpretation of results. JJG conceived and designed the project, and participated in the interpretation and discussion of the results, as well as in the writing of the manuscript. All authors read and approved the final manuscript.
We thank Erena Ruiz Mora, Juana Molina and Reyes Alvarez for their technical assistance, and Eloisa Andújar and Mónica Pérez from the Genomic Unit of CABIMER for their excellent array technical assistance. This work was supported by the Spanish Ministry of Economy and Competitiveness (AGL2011-28904 and AGL2014-54089-R). JHU is a predoctoral researcher supported by the FPI Research Program of the Spanish Ministry of Economy and Competitiveness (BES-2012-058642). SZL is a postdoctoral researcher supported by the Postdoctoral Trainee Program of the Spanish Ministry of Economy and Competitiveness (FPDI-2013-15619).
The authors declare that they have no competing interests.
Juber Herrera Uribe and Melania Collado-Romero contributed equally