Oral administration of Lactobacillus plantarum 299v modulates gene expression in the ileum of pigs: prediction of crosstalk between intestinal immune cells and sub-mucosal adipocytes
- Marcel Hulst†1Email author,
- Gabriele Gross†1, 4,
- Yaping Liu1, 5,
- Arjan Hoekman1,
- Theo Niewold2,
- Jan van der Meulen1, 3 and
- Mari Smits1, 3
© The Author(s) 2015
Received: 19 November 2014
Accepted: 28 March 2015
Published: 11 April 2015
To study host–probiotic interactions in parts of the intestine only accessible in humans by surgery (jejunum, ileum and colon), pigs were used as model for humans. Groups of eight 6-week-old pigs were repeatedly orally administered with 5 × 1012 CFU Lactobacillus plantarum 299v (L. plantarum 299v) or PBS, starting with a single dose followed by three consecutive daily dosings 10 days later. Gene expression was assessed with pooled RNA samples isolated from jejunum, ileum and colon scrapings of the eight pigs per group using Affymetrix porcine microarrays. Comparison of gene expression profiles recorded from L. plantarum 299v-treated pigs with PBS-treated pigs indicated that L. plantarum 299v affected metabolic and immunological processes, particularly in the ileum. A higher expression level of several B cell-specific transcription factors/regulators was observed, suggesting that an influx of B cells from the periphery to the ileum and/or the proliferation of progenitor B cells to IgA-committed plasma cells in the Peyer’s patches of the ileum was stimulated. Genes coding for enzymes that metabolize leukotriene B4, 1,25-dihydroxyvitamin D3 and steroids were regulated in the ileum. Bioinformatics analysis predicted that these metabolites may play a role in the crosstalk between intestinal immune cells and sub-mucosal adipocytes. Together with regulation of genes that repress NFKB- and PPARG-mediated transcription, this crosstalk may contribute to tempering of inflammatory reactions. Furthermore, the enzyme adenosine deaminase, responsible for the breakdown of the anti-inflammatory mediator adenosine, was strongly down-regulated in response to L. plantarum 299v. This suggested that L. plantarum 299v-regulated production of adenosine by immune cells like regulatory T cells may also be a mechanism that tempers inflammation in the ileum, and perhaps also in other parts of the pig’s body.
KeywordsLactobacillus plantarum 299v Pigs Intestine Gene expression Crosstalk
Probiotic bacteria may exert a beneficial effect on the host’s health by steering immunological reactions in the gastrointestinal (GI) tract. One example for a bacterial strain that is applied in fermented probiotic drinks is Lactobacillus plantarum (L. plantarum: De Vries et al. 2006; Molin 2001). L. plantarum originates from healthy human colonic mucosa and belongs to the predominant mucosal Lactobacillus flora in the human intestine and survives in gastrointestinal passage (Johansson et al. 1993; Molin 2001). However, knowledge about the mechanisms underlying protective effects of probiotic bacteria like Lactobacilli in the complex environment of the GI tract is still limited (reviewed in reference, van der Meulen et al. 2010).
In gut mucosa of mammals, specialized immunological cells (dendritic cells, M cells, macrophages, etc.) constantly survey the lumen for the presence of potential pathogenic microorganisms. Conserved non-self molecular structures/motifs common for groups of microbes (microbe-associated molecular patterns, MAMP) are recognized by toll-like receptors (TLRs) and pattern recognition receptors (PRRs) of these immune cells. Recognition by TLR–PRRs discriminates between MAMP of pathogenic microbes (pathogen-associated molecular patterns, PAMP) and of commensal bacteria that assist in the fermentation and absorption of food. While PAMP activates TLR-mediated innate defence mechanism and a pro-inflammatory epithelial response to protect the host from infection, recognition of MAMP of commensal bacteria by TLR’s is ignored and prevents a pro-inflammatory response. For probiotic bacterial species, several surface structures and molecules (e.g. flagellins, polysaccharides and lipoteichoic acids) are identified as MAMP (Lebeer et al. 2010). Similarly, as was observed for commensal and pathogenic bacteria, small structural variations between MAMP of different probiotic bacterial species determine the type of response mediated by TLR–PRR activation (reviewed Lebeer et al. 2010). To understand the complex mechanisms how these MAMP of probiotic bacteria interact with cells of the gut, in vivo studies are needed in which heterogeneous multi-cellular systems interact mutually with microorganisms. Such in vivo studies will provide knowledge about how probiotic bacteria exert their beneficial effects in the gut, and how induced (immune) responses in the mucosa are transmitted to other parts of the host body to improve health.
In humans, transcriptional changes in duodenal mucosa upon continuous intraduodenal infusion of L. plantarum strain WCFS1 for 6 h have been studied by collecting duodenal biopsies (van Baarlen et al. 2009). Due to its safety and survival in the human gastrointestinal tract and the availability of the complete genome sequence, strain WCFS1 is particularly suitable to explore probiotic modes of action on the molecular level in humans (De Vries et al. 2006; van Baarlen et al. 2009). However, such studies are invasive for volunteers and only allow sampling in the duodenum and distal colon and not in parts of the intestine like the ileum and jejunum. With regard to digestion of food and other physiological functions, among all mammals, the intestine of pigs probably best resembles that of humans (van der Meulen et al. 2010). Therefore, in the current study, we investigated host–microbe interactions in the jejunum, ileum and colon of pigs after repeated oral administrations of L. plantarum 299v, starting with a single dose followed by three consecutive daily dosings 10 days later. For assessment of the transcriptional response in the mucosa of these intestinal tissues, gene expression profiles were generated using porcine whole-genome microarrays and compared with those obtained from pigs administered with PBS. The results of this study showed that L. plantarum 299v modulates in vivo transcriptional profiles mainly in the ileum. Bioinformatics analysis of response genes enabled us to generate a predictive overview of processes/pathways activated on a transcriptional level in the ileum after repeated consumption of L. plantarum 299v.
The animal experiment was approved by the institutional animal experiment committee “Dier Experimenten Commissie Lelystad” (DEC Protocol 2006110.c) in accordance with the Dutch regulations on animal experiments. Weaned male 4-week-old pigs (Duroc x Topigs 20; eight per group) were obtained from a commercial Dutch piggery. Starting directly after weaning until the end of the experiment, animals were fed ad libitum with commercial dry pig weaning feed. After 14 days of adaptation to housing conditions, faecal samples were collected and the experiment started on the next day (day 0) with a single oral administration of bacterial suspensions to the animals: L. plantarum 299v wild-type strain (Johansson et al. 1993), approximately 1 × 1012 colony-forming units (CFU) per ml in phosphate-buffered saline (PBS), or PBS as a negative control (each 5 ml). Using a syringe, the suspension was injected in the oral cavity of the piglets. After a period of 10 days during which faecal samples of the animals were collected for microbiological analysis, the treatment was repeated for three consecutive days. Twelve hours after the last oral administration, the animals were euthanized by barbiturate overdose. Mucosal scrapings were collected at the mid of the jejunum, ileum and colon (not rinsed), immediately frozen in liquid nitrogen, and stored at −80 °C until total RNA was extracted.
Bacterial strains and culture media
For the animal experiment, a rifampicin-resistant derivative of the wild-type L. plantarum 299v strain was selected as described earlier (Johansson et al. 1993; Bron et al. 2004; Gross et al. 2008). L. plantarum 299v was grown overnight at 37 °C in Man Rogosa Sharp (MRS) broth (Merck, Darmstadt, Germany) containing 100 μg/ml rifampicin. Bacterial suspensions in PBS were prepared by 2 rounds of centrifugation of the bacterial culture and washing/re-suspending of the bacterial pellet in PBS. Bacterial counts of the L. plantarum 299v suspensions in PBS administered to the pigs were determined by plating of serial dilutions on appropriate agar plates. Faecal samples of 1 cm2 were homogenized in physiological salt solution and stored at 4 °C before serial dilutions of these samples were plated on appropriate agar plates. Samples of 1 cm2 frozen mucosal scrapings collected from the jejunum, ileum and colon were homogenized in PBS using the TisuPrep Homogenizer Omni TP TH220P and stored at 4 °C before serial dilutions of these samples were prepared and plated on appropriate agar plates. Serial dilutions of homogenized faecal and mucosal scrapings were plated on MRS agar containing 100 μg/ml rifampicin and incubated anaerobically at 37 °C for 48 h before CFU were counted. Results from bacterial counts were expressed as mean log10 CFU/g faecal sample ± standard error of the mean (SEM).
Isolation of total RNA and microarray analysis
Total RNA was isolated from approximately 0.5 g of frozen mucosal scrapings collected from jejunum, ileum and colon of individual pigs. The tissue was homogenized in 2 ml TRIzol reagent (Invitrogen, Breda, the Netherlands), and the extractions were essentially performed according to the instructions of the manufacturer with additional steps to remove proteoglycan and polysaccharide contaminations and DNAse treatment as described earlier (Niewold et al. 2005). Finally, RNA was dissolved in RNAse-free water and stored at −80 °C until further use. All samples were checked for RNA integrity by agarose gel electrophoresis and UV-spectrophotometry.
After additional RNA clean-up (NucleoSpin RNA II, Macherey–Nagel, Düren, Germany), the quality and integrity of the samples were analysed using Agilent Lab-on-a-Chip and Bioanalyzer (Agilent Technologies, Amstelveen, the Netherlands). All samples scored a RNA integrity number (RIN value) of ≥9. Samples from the same intestinal location of individual animals of the same treatment group were pooled by mixing equal amounts of RNA (n = 8 for all groups, except for the jejunum PBS control: n = 7). Quality of the RNA pools was assessed again, all scoring a RIN value of ≥9. For gene expression analysis, Affymetrix Porcine Genome Microarray Chips (Affymetrix, Santa Clara, CA) were used, following Affymetrix protocols for hybridization, washing, staining and scanning of the arrays. Briefly, biotin-labelled cRNA was synthesized from mRNA by Oligo(dT)-priming with the Affymetrix One-Cycle Target Labeling and Control Kit using 2 μg of pooled total RNA as input. 12.5–20 μg cRNA was further fragmented, and finally, 10 μg was hybridized to the arrays. Quality control and pooling of the RNA samples, Affymetrix microarray analysis and quality control of the arrays were performed externally by ServiceXS (Leiden, the Netherlands). Also data processing, array normalization and analysis of the array data were performed at Service XS using GeneSpring Software (Agilent Technologies). Briefly, spots with a raw intensity ≤50 were filtered out, and a 50th percentile normalization was performed for each chip. In addition, for each individual probe/gene, the intensity was normalized using the median intensity calculated from all chips/samples. The expression profiles induced by the treatment with L. plantarum 299v were compared with the PBS (control) treatment separately for the jejunum, ileum and colon. A threshold level of >2× and <0.5× fold change (FC) of differential expressed probes was used. Normalized array data were posted in the EMBL-EBI ArrayExpress database under deposit number E-MEXP-2198. Note that the results presented in this manuscript were part of a larger study. In E-MEXP-2198 sources, MSA 09, 13, and 17 represent the array data of PBS RNA pools of the jejunum, ileum and colon, respectively, and sources MSA 10, 14, and 18 represent the L. plantarum 299v wild-type RNA pools of the jejunum, ileum and colon, respectively.
The relative concentration of interleukin-1-beta (IL1B) and carbonic anhydrase 2 (CA2) mRNA in all individual RNA samples extracted from mucosal scrapings from the ileum was determined by quantitative real-time PCR (QRT-PCR). RT reactions were performed with Superscript III (Invitrogen) and random hexamer primers (pdN6) according to the manufacturer’s instructions using 250 ng of RNA template. The gene-specific primers and cycle specifications for quantification of IL1B cDNA were applied as described recently (Schreur et al. 2011). Gene-specific primers for CA2 were developed from the pig mRNA sequence XM_001927805.1. A CA2-specific cDNA fragment of 127 bp was amplified using a 20-mer forward (CAACGGCCACTCTTTCAACG) and reverse primer (TTGCCCATCAGATGAACCCC), using the same cycle protocol as described for IL1B (Schreur et al. 2011). The quantity of 18S rRNA in each RNA sample was determined by QRT-PCR using the above described RT reaction products (Durand et al. 2009) and used to normalize IL1B and CA2 data. The quantity of 18S ribosomal RNA showed no essential differences among all individual RNA samples extracted from scrapings.
Bioinformatics and functional analysis
In addition to the annotation provided by Affymetrix, oligonucleotide sequences of differentially expressed probes not annotated yet, or annotated as Unigene, tentative consensus sequences (TC) or mRNA accession number, were compared with the NCBI non-redundant nucleotide databases using blastn to assign a gene name to these probes. Probes that did not produce a significant match with any other eukaryotic mRNA/gene were excluded from gene lists used for functional analysis. Throughout this manuscript, official human gene symbols (HUGO Gene Nomenclature Committee: http://www.genenames.org) were used in the text and in all (supplementary) figures and tables. In all results, paragraphs beneath, information about the biological functions of genes was retrieved by consulting the “GeneCards” (Weizmann Institute of Science) and the NCBI Gene reports, and literature linked to these reports. For genes with multiple functions, mostly only information is provided about their function related to processes/pathways identified in this study, and we refer for information about additional functions to these reports.
The Database for Annotation, Visualization and Integrated Discovery (DAVID version 6.7) website (da Huang et al. 2009) and the “Set Distiller” module of GeneDecks (Stelzer et al. 2009) were used to assign genes to a specific pathway. Because far more human genes are annotated and more information in databases is available for humans than for pigs, the human background was used for this functional analysis. From DAVID, pathways (KEGG) with a p value of <0.1 (EASE score) were retrieved. In GeneDecks, pathways (KEGG, MLPR, CST, GeneGlobe Pathway Central, Invitrogen and Ingenuity) were called significant with a p value <0.05 using the Set Distiller algorithm. KEGG pathways retrieved from DAVID were only listed when not called significant by GeneDecks, or in case more genes were listed by DAVID than by GeneDecks. Pathways that contained <5 regulated genes were not retrieved. Associations of genes with chemical compounds with immune-modulating properties were retrieved from GeneDecks (p value of <0.05). Using the “tissue expression” module of DAVID, regulated genes were explored for expression in specific (subsets of) immune cells and organs by enrichment analysis using the GNF_U133A_QUARTILE cDNA libraries (http://biogps.org/#goto=welcome). Cell-specific libraries enriched with a p value of <0.05 were considered significant. From DAVID, “Functional Annotation charts”, transcription factors or genes involved in regulation of transcription were identified by gene-ontology analysis and uploaded as sub-list in GNCPro (free online software developed and maintained by SABiosciences Inc. a Qiagen company) to establish associations between these genes (i.e. to build a network of transcription factors and regulators). Using the protein interaction tool of DAVID, regulated gene sets were enriched for specific transcription factor binding sequences (UCSC_TFBS module). Seven transcription factors for which the enrichment analysis of their corresponding transcription recognition sequence was most significant (p value of <0.05) were added to the DAVID sub-list loaded in GNCPro. Non-interacting genes were omitted from the displayed network.
Functional associations between immune-modulating chemicals, proteins encoded by differentially expressed genes and enzyme substrates/products linked to these proteins were established using the (protein)–protein–chemical interaction web tool STITCH 4.0 beta (Kuhn et al. 2009). Relevant chemicals were uploaded together with list of genes in STITCH. Associations with a confidence score of ≥0.4 (medium level) were selected from output files and displayed in a network. In supplementary file Table S3, the type and confidence level of each association are listed.
When appropriate, results are indicated ± SE of the mean (SEM) and the significance of the difference between results of the treatment and control conditions was calculated using Student’s t test (two sided; considered statistically significant when p is <0.05).
L. plantarum persistence characteristics in the pig intestine
Gene expression patterns induced by L. plantarum
Compared with the ileum, a significantly lower number of differentially expressed genes affected by L. plantarum 299v were detected in the jejunum and colon, i.e. a total of 54 and 125 (up and down), respectively. To gain insight in the processes induced by L. plantarum 299v, sets of differentially expressed genes were analysed in DAVID and GeneDecks to retrieve significant pathways for each part of the intestine. In contrast to the ileum (see below), only one significant immune-related pathway (PPAR signalling) could be retrieved from DAVID (EASE score <0.1) and GeneDecks (p value of <0.05) for the colon, and none for the jejunum. Only 12 out of the 407 genes regulated in the ileum were also found regulated in the colon. Two of these genes were part of the PPAR signalling pathway (APOC3 and FABP2), and only 2 genes showed a defined function associated with immune cells (CCL25 and CR2). Regulation of several apolipoproteins in the colon (APOA4, APOB, APOA1 and APOC3) suggested that absorption of lipid vitamins and/or metabolism of lipids and lipoproteins was regulated in the colon. Clearly, these parts of the intestine showed a lower responsiveness towards L. plantarum 299v than the ileum. Therefore, further bioinformatics and functional analysis were conducted only for differentially expressed genes in the ileum.
Pathway analysis of genes induced by L. plantarum 299v
Pathway analysis (summarized)
Cell cycle/checkpoint control
ATR, BUB1B, CDC45, CDK2, KIF11, MCM2, MCM3, MSH2, MSH6, ORC1, PCNA, PRIM2, RCC2, TMPO, VRK1, WEE1
ACTG2, DES, GNB2L1, LMNB1, MST4, MYH10, PFN1, PPP1R12A, RCC2, ROCK1, STMN1, TMPO, VAV1
BCL6, HCK, INPP5D, LCP1, MECOM, NQO1, PLCG2, PSIP1, SPI1
Cell cycle_Role of APC in cell cycle regulation
BUB1, BUB1B, CCNB2, CDC27, CDK2, FBXO5, ORC1
Cell cycle_Transition and termination of DNA replication
CDK2, LIG1, MCM2, PCNA, PRIM2, RFC3
Immune response_BCR pathway
BCL6, CD79B, CR2, INPP5D, PIK3R1, PLCG2, VAV1
Chemokine signalling pathway (GeneDecks-KEGG)
CCL28, CXCL13, CXCR4, HCK, PIK3R1, ROCK1, TIAM1, VAV1
Leucocyte transendothelial migration (GeneDecks-KEGG)
CDH5, CXCR4, PIK3R1, PLCG2, ROCK1, VAV1
Focal adhesion (GeneDecks-KEGG)
CAPN2, COL1A1, COL3A1, PIK3R1, PPP1R12A, ROCK1, VAV1
NFKB signalling Targets
AGT, CD40, CFB, LTB, NQO1
RFC3, LIG1, PRIM2, PCNA, MCM2, MCM3
MSH6, RFC3, MSH2, LIG1, PCNA
hsa04662:B cell receptor signalling pathway
CR2, PLCG2, PIK3AP1, CD79B, INPP5D, VAV1, PIK3R1
ENPP3, RRM2, AK1, PRIM2, GUCY1B3, RRM2B, PAPSS2, ADA
hsa04666:Fc gamma R-mediated phagocytosis
ACTG2, HCK, PLCG2, PIKFYVE, ARPC4, INPP5D, VAV1, PIK3R1
hsa04672:Intestinal immune network for IgA production
CXCR4, TNFRSF17, PIGR, CD40, CCL28
hsa04115:p53 signalling pathway
CCNB2, RRM2, RRM2B, ATR, CDK2
hsa03320:PPAR signalling pathway
HMGCS2, APOC3, PPARG, FABP2, MMP1
Genes mapped to the PPAR signalling pathway were all down-regulated by L. plantarum 299v; among them, the nuclear receptor PPARG which forms a receptor complex with the retinoic X receptor (RXR). Binding of ligands like 9-cis-retinoic acid (RA), fatty acids or eicosanoids to this complex activates adipocyte differentiation, lipid metabolism and negatively regulates macrophage activation. In conjunction with this, genes associated with steroids, leukotriene B4 and RA were called significant by GeneDecks (Table S2).
Several NFKB targets were up- or down-regulated in animals treated with L. plantarum 299v. CFB, AGT and NQO1 were down-regulated, and CD40 and LTB were up-regulated. AGT is a potent activator of cortisol by stimulating transcription factor nuclear receptor subfamily 4, group A (Romero et al. 2010). The Ba fragment of complement factor CFB (down-regulated here) inhibits the proliferation of pre-activated B-lymphocytes. NQO1 null-mice suffer from myeloid hyperplasia, suggesting that down-regulation of this quinone reductase supports lymphocyte proliferation. Binding of the LTB–LTA complex to the LT-beta receptor and binding of CD40 to its receptor on the surface of lymphocytes mediate TRAF-dependent NFKB signalling, which may initiate stimulation of lymphocyte attachment, B cell development/survival and T cell co-stimulation. The CD40 receptor also plays a role in the production of IgA in the intestine, a pathway called significant in this study. CCL28, a chemokine that may stimulate migration of IgA-committed plasma cells (B cells) from the Peyer’s patches to the lamina propia, and the Fc receptor PIGR, responsible for translocation of IgA complexes over the epithelial layer, were both down-regulated in response to L. plantarum 299v. In addition to the IgA pathway, B cell receptor signalling pathways were called significant, suggesting that an influx and/or proliferation of subsets of B cells was indeed regulated. Furthermore, pathway analysis suggested that phagocytosis (engulfment of antigens by macrophages, neutrophils, etc.) and activation/proliferation of NK cells were stimulated. The two most dominant pathways, cytoskeleton rearrangement and cell division were activated, most likely, to support the transformation of progenitor cells to active immune cells or to support their migration and homing.
The tissue expression module of DAVID predicted that a high percentage of regulated genes in the ileum are expressed in a lymph node library, CD4+ cells and CD14+ monocytes, suggesting that cell division in both lymphoid- and myeloid-derived progenitor cells was stimulated (see Table S2). The regulation of several B cell-specific markers/genes (e.g. CR2, CD40, EBF1, CXCR4, TNFRSF17, CXCL13, BCL6, BCL7A, CD79B, PIK3AP1, PIK3R1, UBE2V2), along with the calling of the IgA production pathway, suggests that this also applies for B cells. In addition, a high percentage of genes regulated by L. plantarum 299v were suggested to be expressed in adipocytes (Table S2). Except for FABP2, all down-regulated genes in the PPAR pathway were predicted to be expressed in adipocytes. PPARG itself and 9 out of the 20 genes associated with steroid metabolism were predicted to be expressed in CD14+ monocytes (shown in red in Table S2), whereas only 4 of them were predicted to be expressed in CD4+ T cells (shown in blue in Table S2).
Transcriptional regulation of L. plantarum 299v-induced gene expression
Immune-modulatory compounds predicted to be associated with genes regulated by L. plantarum 299v
In most clinical trials with human cohorts, beneficial effect of Lactobacilli intake is recorded after several weeks of intake of a fixed dose before each meal. To ensure that all piglets received a similar volume of the L. plantarum suspension, we delivered it directly in the oral cavity using a syringe. However, administration in the oral cavity of young pigs in a comparable frequency as used in human trials (e.g. 3 times a day for several weeks) would be too stressful for young pigs. Therefore, we started with a single high dose of L. plantarum 299v followed by three consecutive daily dosings 10 days later. Despite this different scheme of intake, information gained in this study in pigs might be valuable for understanding of the mechanism of action of this probiotic bacterium in parts of the human intestine (jejunum and ileum) that can only be sampled in using surgery.
Modulation of gene expression was observed mainly in the ileum, the immunologically most active part of the intestine of pigs and of clusters of genes involved in specific immunological (or immune-related) pathways/processes characteristic for the ileum and also described in previous studies in which the transcriptional response in the intestine was measured after consumption of L. plantarum and other probiotic bacteria. QRT-PCR analysis indicated that for particular genes, a considerable natural variation in expression levels of genes in the ileum exists between individual pigs. This difference in biological/genetic background might have influenced the response to L. plantarum 299v. However, identification of meaningful pathways retrieved in this study suggests that 8 pigs per group/pool were sufficient to allow for the detection of significant effects despite this natural variation.
Inflammatory mediators like IL8, IL1B, FOS, NFKB and acute-phase response proteins like REG3A, for which up-regulation of expression was observed in several animal studies in which the pigs intestine was challenged with bacterial pathogens like ETEC and Salmonella (Niewold et al. 2010; Hulst et al. 2013, and reference herein), did not responded at all to L. plantarum 299v (see also Fig. 2; IL1B expression). This clearly showed that the intestinal mucosa of pigs responded differently, and less vigorous, to exposure to L. plantarum 299v bacteria than to hostile enteric bacteria. The fact that L. plantarum 299v did not colonize in the pig’s intestine may be related to this mild response.
In agreement with other in vivo studies conducted with L. plantarum strain WCFS1 in humans (Bujalance et al. 2007; Troost et al. 2008; van Baarlen et al. 2009), our results also indicate that L. plantarum regulates the influx and/or enhances the proliferation and differentiation of lymphocytes in the intestine without inducing NFKB-mediated inflammatory processes. According to our predictions, and indicated by up-regulation of cell-specific marker genes, gene expression in CD4+ lymphocytes, CD14+ monocytes, adipocytes and B cells were regulated in response to L. plantarum 299v. The latter cells can be progenitor B cells that proliferate to IgA-plasma cell progenitors in the Peyer’s patches. With respect to this B cell proliferation, similar observations have been described in studies in which mice and humans were administrated with other probiotic bacteria (Fukushima et al. 1998, 1999; Tsai et al. 2010). However, L. plantarum 299v could also stimulate an influx of CD14+ monocytes and B cells from the periphery to the ileum.
Cell wall component teichoic acid (TA) was identified as a dominant MAMP of gram-positive bacteria including several probiotic Lactobacillus bacteria (Lebeer et al. 2010). Recognition of L. rhamnosus GG TA’s by TLR2/6 dimers results in NFKB-mediated cytokine production (Claes et al. 2012). In this study, only up-regulation of TLR9 expression was observed, which recognizes bacterial CpG DNA, but no regulation of TLR2/6 gene expression. A recent study, in which healthy non-compromised mice were administered with a mutant of L. plantarum strain WCFS1 containing much less d-alanine incorporated in its cell wand TA, indicated that this dominant MAMP not only influenced a pro-inflammatory response, but also a systemic anti-inflammatory response characterized by an increase in regulatory T cells (Treg) in the spleen (Smelt et al. 2013). Interestingly, the gene adenosine deaminase (ADA) was found down-regulated more than tenfold in our study. Extracellular adenosine produced by Tregs acts as anti-inflammatory mediator on the immune system. A recent study indicated that accumulation of adenosine produced by Tregs in ADA-deficient SCID patients contributes to autoimmune manifestations (Sauer et al. 2012). Therefore, L. plantarum-induced down-regulation of ADA in the ileum of pigs and increase in splenic Tregs in mice (Smelt et al. 2013) suggest that adenosine-mediated immune suppression by Tregs may be initiated in the intestines and transmitted to the periphery. Further studies in non-compromised animals with probiotic bacteria with altered TA structures on their surface focusing on this ADA-Treg immune mechanism could reveal valuable information about how probiotic bacteria communicate their health beneficial effects to other parts of the host body.
By building a network of transcription factors, we identified a set of B cell-specific transcription factors (BCL6, BCOR, SPIB, SPI1, TCF4 and IKZ3). Up-regulation of BCOR may suppress BCL6-mediated transcription in order to prevent excessive B cell maturation. In addition, down-regulation of CCL28 and IgA-transporter PIGR may also counteract excessive B cell migration and unnecessary secretion of produced IgA into the lumen, respectively.
NFKB-mediated transcription can be repressed by up-regulation of CLU, a complement cytolysis inhibitor (Takase et al. 2008) and several ubiquitin-specific proteases, like USP48 (Tzimas et al. 2006) and USP21 (Xu et al. 2010). In addition, genes like DMBT1 (a target gene for the intracellular pathogen recognition receptor NOD2), AGT, CFB and BPI (bactericidal/permeability-increasing protein), normally found up-regulated in intestine infected by bacterial pathogens, were all down-regulated after L. plantarum administration. Furthermore, ubiquitin-conjugating enzymes UBE2E2, UBE2V2 and UBE2T were up-regulated. In conjunction with anti-apoptotic c-IAPs (like BIRC5; up-regulated in this study), human UBE2s promote polyubiquitination of RIPK1 and suppresses DDX58 (alias RIG-I) ubiquitination by RNF125 (Dynek et al. 2010). Both cytosolic receptors, RIPK1 and DDX58, are essential for recognition of viral and bacterial antigens and subsequent activation of interferon regulatory factors (IRF’s) or NFKB-mediated transcription of chemotactic and pro-inflammatory factors like IL1B.
Down-regulation of the enzyme HSD11B2, which catalyses the inter-conversion of the glucocorticoids (GCs) cortisol and cortisone, and the regulation of several other corticosteroid-related genes (SGK2, OSBPL3, OSBPL8, SOAT1 and AKR1C1) suggest that L. plantarum 299v also affected the concentration of endocrine GCs in the ileum. GCs are potent inhibitors of inflammatory cytokine production by macrophages. Interestingly, in the promoter of porcine MIF (macrophage migration inhibitory factor), two potential PPARG binding sites were identified recently (Chen et al. 2010). MIF antagonizes the inhibition of GC-induced cytokine production in macrophages by restoring the transcriptional activity in particular of NFKB. Therefore, down-regulation of PPARG and/or modulation of its transcriptional activity by NRIP1, in conjunction with regulation of GC levels, could be one of the mechanisms of L. plantarum 299v to control the inflammatory activity of macrophages in the ileum.
NRIP1-regulated PPARG signalling has also been shown to affect lipid metabolism and differentiation of adipocytes (Debevec et al. 2007). Our bioinformatics analysis predicted that most of the response genes related to PPARG signalling and lipid/steroid metabolism are expressed by adipocytes and CD14+ monocytes. White adipose tissue is found in the sub-mucosa, a layer just underneath the lamina propia that surrounds the Peyer’s patches embedded in the mucosal layer. Crosstalk between cultured macrophages and adipocytes after exposure to heat-killed L. gasseri and L. rhamnosus GG resulted in down-regulation of PPARG, and subsequently, to an altered cytokine response in macrophages (Miyazawa et al. 2012). A similar mechanism of crosstalk between macrophages and adipocytes in the sub-mucosal layer of the ileum of pigs can also direct polarization of macrophages towards an anti-inflammatory phenotype (Odegaard et al. 2007; Lumeng et al. 2007).
According to Fig. 4, 1,25-dihydroxyvitamin D3 metabolism may play a key role in crosstalk between the different types of cells in the ileum that respond to L. plantarum 299v. Endogenous 1,25-dihydroxyvitamin D3 can block expression of the transcription factor/regulators C/EBP-alpha, sterol regulatory element-binding protein-1 and PPARG in adipocytes, suppressing cell differentiation in these cells (Wood 2008). Production of 1,25-dihydroxyvitamin D3 by human B cells is stimulated by activation of the CD40 receptor and regulated by a complex feedback mechanism in which 1,25-dihydroxyvitamin D3 itself, transient-receptor-potential-cation-channels (like TRPV6; down-regulated here) and the enzyme that degrades 1,25-dihydroxyvitamin D3 (CYP24A; down-regulated here) are involved (Heine et al. 2008).
The inflammatory mediator leukotriene B4 is produced by polymorphonuclear leucocytes and macrophages. A previous studies indicated that leukotriene B4 can activate T cells that inhibit Epstein-Barr virus-induced proliferation of cord blood-derived B-lymphocytes (Liu et al. 2008) and that living and heat-killed preparations of lactic acid bacterial strains were able to inhibit leukotriene B4 production in murine macrophages (Kimoto-Nira et al. 2009). In the current study, we observed a moderate down-regulation (twofold) of PTGR1 (alias LTB4DH), the enzyme that catalyses the conversion of inflammatory mediator leukotriene B4 to biologically less active 12-oxo-leukotriene B4, suggesting that the activity of leukotriene B4 was raised in response to L. plantarum 299v. This is in contrast with results of two in vitro studies cited above and may suggest that in vivo regulation of leukotriene B4 metabolism in the ileum may be dependent on additional factors and/or specific types of cells not present in the above used in vitro cell systems.
Calves treated with corticosteroids produced low levels of IgA in the intestine (Husband et al. 1973). Also, human asthma patients treated with inhaled corticosteroids produce less salivary IgA (Fukushima et al. 2005). Furthermore, studies in ferrets indicated that development of Peyer’s patches and/or delay in the migration of IgA-plasma cell progenitors from this lymphoid organ into the lamina propia was depressed by corticosteroids (Carlile and Beck 1983). This suggests that besides 1,25-dihydroxyvitamin D3 and leukotriene B4, steroids may also play a role in the crosstalk between cells that respond to L. plantarum 299v.
In conclusion, our data suggested that the probiotic bacterium L. plantarum 299v regulates the activity of adipocytes and/or different subsets of B cell in the ileum of pigs. Either infiltration of the ileum by these cells or their maturation/differentiation to functional cells is stimulated. Progenitor B cells may mature to IgA-committed plasma cells and travel from the Peyer’s patches to the lamina propia in response to L. plantarum 299v. In case external challenges are detected, these IgA-producing cells can respond directly by secretion of IgA antibodies into the lumen at the site of damage. In addition, in this study, induced repression of NFKB-mediated transcription and PPARG signalling may implicate that L. plantarum 299v has a prophylactic effect on development of inflammation in the ileum of pigs. Crosstalk between specific immune cells present in the ileum and sub-mucosal adipocytes may play an important role in this repression. The role of leukotriene B4, 1,25-dihydroxyvitamin D3, steroids and adenosine in this crosstalk is an interesting issue for further research.
This work was supported by the Centre for Human Nutrigenomics, the Netherlands. The L. plantarum 299v strain was kindly provided by Johannes Snel and Roelof van der Meer of NIZO Food Research, Ede, the Netherlands. The authors thank Peter van Rossum, Klaas Jan Visser en Johan Meijer for excellent handling of the animal experiment, Dirkjan Schokker for advice on microarray data analysis, and Johannes Snel and Roelof van der Meer for their advice and comments.
Conflict of interest
The author(s) declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
- Biskobing DM, Fan D, Fan X, Rubin J (1997) Induction of carbonic anhydrase II expression in osteoclast progenitors requires physical contact with stromal cells. Endocrinology 138:4852–4857PubMedGoogle Scholar
- Biswas P, Mengozzi M, Mantelli B, Delfanti F, Brambilla A, Vicenzi E, Poli G (1998) 1,25-Dihydroxyvitamin D3 upregulates functional CXCR4 human immunodeficiency virus type 1 coreceptors in U937 minus clones: NF-kappaB-independent enhancement of viral replication. J Virol 72:8380–8383PubMed CentralPubMedGoogle Scholar
- Bron P, Grangette C, Mercenier A, de Vos W, Kleerebezem M (2004) Identification of Lactobacillus plantarum genes that are induced in the gastrointestinal tract of mice. J Bacteriol 186:5721–5729View ArticlePubMed CentralPubMedGoogle Scholar
- Bujalance C, Moreno E, Jimenez-Valera M, Ruiz-Bravo A (2007) A probiotic strain of Lactobacillus plantarum stimulates lymphocyte responses in immunologically intact and immunocompromised mice. Int J Food Microbiol 113:28–34View ArticlePubMedGoogle Scholar
- Carlile AE, Beck F (1983) Maturation in the ferret ileal epithelium and the effect of cortisone acetate. J Anat 136:401–405PubMed CentralPubMedGoogle Scholar
- Chen XD, Chen YL, Wang SB, Lei T, Gan L, Yang ZQ (2010) Molecular characterization and functional analysis of porcine macrophage migration inhibitory factor (MIF) gene. Cytokine 50:84–90View ArticlePubMedGoogle Scholar
- Claes I, Segers M, Verhoeven T, Dusselier M, Sels B, Keersmaecker S, Vanderleyden J, Lebeer S (2012) Lipoteichoic acid is an important microbe-associated molecular pattern of Lactobacillus rhamnosus GG. Microb Cell Fact 11:161View ArticlePubMed CentralPubMedGoogle Scholar
- da Huang W, Sherman B, Lempicki R (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4:44–57View ArticleGoogle Scholar
- De Vries M, Vaughan E, Kleerebezem M, de Vos W (2006) Lactobacillus plantarum—survival, functional and potential probiotic properties in the human intestinal tract. Int Dairy J 16:1018–1028View ArticleGoogle Scholar
- Debevec D, Christian M, Morganstein D, Seth A, Herzog B, Parker M, White R (2007) Receptor interacting protein 140 regulates expression of uncoupling protein 1 in adipocytes through specific peroxisome proliferator activated receptor isoforms and estrogen-related receptor alpha. Mol Endocrinol 21:1581–1592View ArticlePubMed CentralPubMedGoogle Scholar
- Durand SV, Hulst MM, de Wit AA, Mastebroek L, Loeffen WL (2009) Activation and modulation of antiviral and apoptotic genes in pigs infected with classical swine fever viruses of high, moderate or low virulence. Arch Virol 154:1417–1431View ArticlePubMed CentralPubMedGoogle Scholar
- Dynek JN, Goncharov T, Dueber EC, Fedorova AV, Izrael-Tomasevic A, Phu L, Helgason E, Fairbrother WJ, Deshayes K, Kirkpatrick DS, Vucic D (2010) c-IAP1 and UbcH5 promote K11-linked polyubiquitination of RIP1 in TNF signalling. EMBO J 29:4198–4209View ArticlePubMed CentralPubMedGoogle Scholar
- Fukushima Y, Kawata Y, Hara H, Terada A, Mitsuoka T (1998) Effect of a probiotic formula on intestinal immunoglobulin A production in healthy children. Int J Food Microbiol 42:39–44View ArticlePubMedGoogle Scholar
- Fukushima Y, Kawata Y, Mizumachi K, Kurisaki J, Mitsuoka T (1999) Effect of bifidobacteria feeding on fecal flora and production of immunoglobulins in lactating mouse. Int J Food Microbiol 46:193–197View ArticlePubMedGoogle Scholar
- Fukushima C, Matsuse H, Saeki S, Kawano T, Machida I, Kondo Y, Kohno S (2005) Salivary IgA and oral candidiasis in asthmatic patients treated with inhaled corticosteroid. J Asthma 42:601–604View ArticlePubMedGoogle Scholar
- Gross G, van der Meulen J, Snel J, van der Meer R, Kleerebezem M, Niewold TA, Hulst MM, Smits MA (2008) Mannose-specific interaction of Lactobacillus plantarum with porcine jejunal epithelium. FEMS Immunol Med Microbiol 54:215–223View ArticlePubMedGoogle Scholar
- Heine G, Niesner U, Chang HD, Steinmeyer A, Zügel U, Zuberbier T, Radbruch A, Worm M (2008) 1,25-dihydroxyvitamin D(3) promotes IL-10 production in human B cells. Eur J Immunol 38:2210–2218View ArticlePubMedGoogle Scholar
- 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:18View ArticleGoogle Scholar
- Husband AJ, Brandon MR, Lascelles AK (1973) The effect of corticosteroid on absorption and endogenous production of immunoglobulins in calves. Aust J Exp Biol Med Sci 51:707–710View ArticlePubMedGoogle Scholar
- Johansson M, Molin G, Jeppsson B, Nobaek S, Ahrne S, Bengmark S (1993) Administration of different Lactobacillus strains in fermented oatmeal soup: in vivo colonization of human intestinal mucosa and effect on the indigenous flora. Appl Environ Microbiol 59:15–20PubMed CentralPubMedGoogle Scholar
- Kaetzel CS, Robinson JK, Chintalacharuvu KR, Vaerman JP, Lamm ME (1991) The polymeric immunoglobulin receptor (secretory component) mediates transport of immune complexes across epithelial cells: a local defense function for IgA. Proc Natl Acad Sci USA 88:8796–8800View ArticlePubMed CentralPubMedGoogle Scholar
- Kimoto-Nira H, Suzuki C, Kobayashi M, Sasaki K, Mizumachi K (2009) Inhibition of leukotriene B4 production in murine macrophages by lactic acid bacteria. Int J Food Microbiol 129:321–324View ArticlePubMedGoogle Scholar
- Kong J, Li YC (2006) Molecular mechanism of 1,25-dihydroxyvitamin D3 inhibition of adipogenesis in 3T3-L1 cells. Am J Physiol Endocrinol Metab 290:E916–E924View ArticlePubMedGoogle Scholar
- Kuhn M, Szklarczyk D, Franceschini A, Campillos M, von Mering C, Jensen LJ, Beyer A, Bork P (2009) STITCH 2: an interaction network database for small molecules and proteins. Nucleic Acids Res 38(Database issue):D552–D556PubMed CentralPubMedGoogle Scholar
- Lebeer S, Vanderleyden J, Keersmaecker S (2010) Host interactions of probiotic bacterial surface molecules: comparison with commensals and pathogens. Nat Rev Microbiol 8:171–184. doi:https://doi.org/10.1038/nrmicro2297 View ArticlePubMedGoogle Scholar
- Liu A, Claesson HE, Mahshid Y, Klein G, Klein E (2008) Leukotriene B4 activates T cells that inhibit B-cell proliferation in EBV-infected cord blood-derived mononuclear cell cultures. Blood 111:2693–2703View ArticlePubMedGoogle Scholar
- Lumeng CN, Bodzin JL, Saltiel AR (2007) Obesity induces a phenotypic switch in adipose tissue macrophage polarization. J Clin Invest 117:175–184View ArticlePubMed CentralPubMedGoogle Scholar
- Mangell P, Thorlacius H, Syk I, Ahrné S, Molin G, Olsson C, Jeppsson B (2012) Lactobacillus plantarum 299v does not reduce enteric bacteria or bacterial translocation in patients undergoing colon resection. Dig Dis Sci 57:1915–1924View ArticlePubMedGoogle Scholar
- Miyazawa K, He F, Yoda K, Hiramatsu M (2012) Potent effects of, and mechanisms for, modification of crosstalk between macrophages and adipocytes by lactobacilli. Microbiol Immunol 56:847–854View ArticlePubMedGoogle Scholar
- Molin G (2001) Probiotics in foods not containing milk or milk constituents, with special reference to Lactobacillus plantarum 299v. Am J Clin Nutr 73:380S–385SPubMedGoogle Scholar
- Moriwaki Y, Yamamoto T, Higashino K (1999) Enzymes involved in purine metabolism—a review of histochemical localization and functional implications. Histol Histopathol 14:1321–1340PubMedGoogle Scholar
- Niewold T, Kerstens H, van der Meulen J, Smits M, Hulst M (2005) Development of a porcine small intestinal cDNA micro-array: characterization and functional analysis of the response to enterotoxigenic E. coli. Vet Immunol Immunopathol 105:317–329View ArticlePubMedGoogle Scholar
- Niewold TA, van der Meulen J, Kerstens HH, Smits MA, Hulst MM (2010) Transcriptomics of enterotoxigenic Escherichia coli infection. Individual variation in intestinal gene expression correlates with intestinal function. Vet Microbiol 141:110–114View ArticlePubMedGoogle Scholar
- Odegaard JI, Ricardo-Gonzalez RR, Goforth MH, Morel CR, Subramanian V, Mukundan L, Red Eagle A, Vats D, Brombacher F, Ferrante AW, Chawla A (2007) Macrophage-specific PPARgamma controls alternative activation and improves insulin resistance. Nature 447:1116–1120View ArticlePubMed CentralPubMedGoogle Scholar
- Romero DG, Gomez-Sanchez EP, Gomez-Sanchez CE (2010) Angiotensin II-regulated transcription regulatory genes in adrenal steroidogenesis. Physiol Genomics 42A:259–266View ArticlePubMed CentralPubMedGoogle Scholar
- Sauer A, Brigida I, Carriglio N, Hernandez R, Scaramuzza S, Clavenna D, Sanvito F, Poliani P, Gagliani N, Carlucci F, Tabucchi A, Roncarolo M, Traggiai E, Villa A, Aiuti A (2012) Alterations in the adenosine metabolism and CD39/CD73 adenosinergic machinery cause loss of Treg cell function and autoimmunity in ADA-deficient SCID. Blood 119:1428–1439View ArticlePubMed CentralPubMedGoogle Scholar
- Schreur PJW, Rebel JM, Smits MA, van Putten JP, Smith HE (2011) Lgt processing is an essential step in Streptococcus suis lipoprotein mediated innate immune activation. PLoS One 6:e22299View ArticleGoogle Scholar
- Smelt M, de Haan B, Bron P, van Swam I, Meijerink M, Wells J, Kleerebezem M, Faas M, de Vos P (2013) The impact of Lactobacillus plantarum WCFS1 teichoic acid d-alanylation on the generation of effector and regulatory T-cells in healthy mice. PLoS One 8(4):e63099. doi:https://doi.org/10.1371/journal.pone.0063099 View ArticlePubMed CentralPubMedGoogle Scholar
- Stelzer G, Inger A, Olender T, Iny-Stein T, Dalah I, Harel A, Safran M, Lancet D (2009) GeneDecks: paralog hunting and gene-set distillation with GeneCards annotation. OMICS 13:477–487View ArticlePubMedGoogle Scholar
- Takase O, Minto A, Puri T, Cunningham P, Jacob A, Hayashi M, Quigg R (2008) Inhibition of NF-kappaB-dependent Bcl-xL expression by clusterin promotes albumin-induced tubular cell apoptosis. Kidney Int 73:567–577View ArticlePubMedGoogle Scholar
- Troost F, van Baarlen P, Lindsey P, Kodde A, de Vos W, Kleerebezem M, Brummer R (2008) Identification of the transcriptional response of human intestinal mucosa to Lactobacillus plantarum WCFS1 in vivo. BMC Genomics 9:374View ArticlePubMed CentralPubMedGoogle Scholar
- Tsai YT, Cheng PC, Liao JW, Pan TM (2010) Effect of the administration of Lactobacillus paracasei subsp. paracasei NTU 101 on Peyer’s patch-mediated mucosal immunity. Int Immunopharmacol 10:791–798View ArticlePubMedGoogle Scholar
- Tzimas C, Michailidou G, Arsenakis M, Kieff E, Mosialos G, Hatzivassiliou EG (2006) Human ubiquitin specific protease 31 is a deubiquitinating enzyme implicated in activation of nuclear factor-kappaB. Cell Signal 18:83–92View ArticlePubMedGoogle Scholar
- van Baarlen P, Troost F, van Hemert S, van der Meer C, de Vos W, de Groot P, Hooiveld G, Brummer R, Kleerebezem M (2009) Differential NF-kappaB pathways induction by Lactobacillus plantarum in the duodenum of healthy humans correlating with immune tolerance. Proc Natl Acad Sci USA 106:2371–2376View ArticlePubMed CentralPubMedGoogle Scholar
- van der Meulen J, Hulst M, Smits M, Schuurman T (2010) Small intestinal segment perfusion test in piglets: future applications in studying probiotics-gut crosstalk in infectious diarrhoea? Benef Microbes 1:439–445View ArticlePubMedGoogle Scholar
- Wood RJ (2008) Vitamin D and adipogenesis: new molecular insights. Nutr Rev 66:40–46View ArticlePubMedGoogle Scholar
- Xu G, Tan X, Wang H, Sun W, Shi Y, Burlingame S, Gu X, Cao G, Zhang T, Qin J, Yang J (2010) Ubiquitin-specific peptidase 21 inhibits tumor necrosis factor alpha-induced nuclear factor kappaB activation via binding to and deubiquitinating receptor-interacting protein 1. J Biol Chem 285:969–978View ArticlePubMed CentralPubMedGoogle Scholar