Short Chain Fatty Acids (SCFA) Reprogram Gene Expression in Human Malignant Epithelial and Lymphoid Cells.

The effect of short chain fatty acids (SCFAs) on gene expression in human, malignant cell lines was investigated, with a focus on signaling pathways. The commensal microbial flora produce high levels of SCFAs with established physiologic effects in humans. The most abundant SCFA metabolite in the human microflora is n-butyric acid. It is well known to activate endogenous latent Epstein-Barr virus (EBV), that was used as a reference read out system and extended to EBV+ epithelial cancer cell lines. N-butyric acid and its salt induced inflammatory and apoptotic responses in tumor cells of epithelial and lymphoid origin. Epithelial cell migration was inhibited. The n-butyric gene activation was reduced by knock-down of the cell membrane transporters MCT-1 and -4 by siRNA. N-butyric acid show biologically significant effects on several important cellular functions, also with relevance for tumor cell phenotype.

Introduction

The metabolism of the human microbiota is intimately linked with that of the host, especially in mucosal tissues like the gut or the nasopharynx. A feature of the colonic microbiota metabolism is the fermentation of complex carbohydrates [1-3]. One important product of this metabolism is the production of short-chain fatty acids (SCFAs), which can have local effects at the site of production as well as systemic ones, through blood circulation [4-5].

SCFAs refer to free fatty acids with short (less than 6 carbons) aliphatic chains. They include formic acid, valeric acid, caproic acid and butyric acid and its structural isomers [6]. The SCFAs are taken up by blood and affect nutrition and the immune system [7]. N-Butyric acid is a 4-carbon straight chain SCFA, most interesting due to its high production by the microbiota. It reaches a concentration of 20mM in the colon. The metabolism of butyrate (salt of butyric acid) has been estimated to provide about 50% of the daily energy requirements of the gastrointestinal mucosa [8-9].

Although the establishment of a healthy gut microbiota, where bifido- and lactobacteria are prevailing, often coincides with an increase in butyrate concentration, neither lactobacilli nor bifidobacteria produce butyrate [10]. The majority of isolates producing high levels of butyrate (more than 10mM) are related to the Coccoides-Eubacterium phylae, which are other dominant members of the gut microbiota [11-12]. SCFAs are naturally found in foods as well. Thus, by modulation of a diet in favor of the proper microbiota one can modulate butyric acid levels locally and systemically [13].

Cells can be affected by SCFAs in three different ways. SCFA bind cell receptors that regulate cell proliferation and differentiation. SCFAs can enter cells through specific transporters and involve directly in the cellular metabolism, thus affecting cell energy status and signaling processes [14]. SCFAs can inhibit HDAC activity in the nuclei. All major SCFAs have HDAC inhibitory activity at large enough concentrations as shown in in vitro studies [15]. Inhibition of HDAC activity will promote access of transcription factors to promoters and activate gene expression. This, in turn can affect inflammatory and even carcinogenic processes at the gene-expression level [16-17].

We employed an Epstein-Barr virus (EBV) model system as a positive control in our study of effects of SCFAs on cells. More than 95% of adult human population carry EBV virus. It is well established that butyric acid can induce lytic EBV production and switch latency programs in EBV infected B cell lines [18]. Butyrate acts via histone deacetylation to induce lytic EBV replication and lysis of cells [19-21]. The first step of the switch from latency to the lytic virus cycle is the expression of immediate early transactivator genes, BZLF1 and BRLF1, which in concert, activate the subsequent viral lytic cascade [22-23]. The role of the nasopharyngeal microbiome and its metabolites for NPC-risk and–progression is will be of future major interest.

A panel of SCFAs focusing on butyric acid was tested. The expression of the specific transporters for butyric acid entry, MCT1 and 4 were investigated. Further genome-wide expression profiling of cells exposed to butyric acid was analyzed. Thus we could demonstrate a multifaceted effect of butyric acid involving several important host cell signaling pathways.

Materials and Methods

Chemicals

The SCFAs caproic, 2-ethylbutyric, n-butyric, isobutyric, isovaleric, n-valeric acids, were obtained from Merck (Darmstadt, Germany). Formic acid, sodium butyrate, hydrochloric acid were from Sigma-Aldrich (St. Louis, Missouri, USA). 12-O-Tetradecanoylphorbol-13-acetate (TPA) from Calbiochem (San Diego, California, USA) was used at 20 ng/ml.

Cell lines and culture

C666-1[24], AGS[25], HONE1[26] and CNE2[27] are cell lines of human epithelial origin. C666-1 has been established from undifferentiated nasopharyngeal carcinoma (NPC). EBV latency II program is expressed in C666-1. AGS was established from gastric carcinoma and carry EBV in latency stage I. AGS cell line was kindly provided by professor Takada K., at University of Pittsburgh, USA. HONE1 and CNE2 are EBV negative cell lines. Nasopharyngeal carcinoma cell lines were kindly provided by professor Maria Li Lung’s laboratory at the University of Hong Kong, China. Raji and Rael are cell lines derived from Burkitt’s lymphomas and carry EBV genomes in latency type III and I, respectively [28,29]. BL cell lines we received from George Klein, MTC, KI. Raji and Rael cells were cultured in RPMI 1640 (HyClone ™, Thermo Scientific ™) containing 10% fetal bovine serum (HyClone ™, Thermo Scientific ™), penicillin (50 U/ml), streptomycin (50 U/ml) [29]. C666-1, HONE1, CNE2 and AGS cells were cultured in Iscove's Modified Dulbecco's Medium (HyClone ™, Thermo Scientific ™) prepared the same way but with15% FBS for the C666-1 line. All cells were grown at 37°C in an atmosphere with high humidity and 5% CO2. B-cells were subcultured at a concentration of 106 cells/ml a night before treatments. 1,5 x 106/ml cells were treated for 24 hrs by reagents as indicated. NPC cells were seeded on a 12 well plate 24hr before treatment to achieve 90% confluency (ca 0,76x106/well) at the day of treatment. Cells were treated by 10 mM SCFA solutions for 24hr. Cell viability was measured by counting live cells with the red coloured exclusion dye, erythrosin B [30].

RNA isolation, cDNA synthesis, and PCR

RNA was purified with RNeasy Mini Kit (Qiagen, USA). Reverse transcription was performed using RevertAid First Strand cDNA Synthesis Kit (Thermo scientific, EU). Taq DNA Polymerase, Recombinant (Thermo scientific, Sweden) wss employed for amplification. The results were analyzed by agarose gel electrophoresis and recorded by system Gel Doc™ XR ChemiDoc™ XRS—Bio-Rad using Image Lab™ Software 170–82657. Signal intensities was quantitated with ImageJ software and normalized to the housekeeping gene density as indicated.

EBV infection of Raji B cells

EBV-containing conditioned media was collected from a gastric carcinoma AGS cell line, cultured as a monolayer, carrying recombinant EBV-GFP. Viral production in AGS cells was induced by 40h incubation with 10 mM Butyric acid. Media was centrifuged before application to Raji B cells. 1x10 6 Raji cells, cultured in suspension, were pelleted by centrifugation at 1000 prm for 5 min at RT, resuspended in 2 ml virus–containing supernatant with 8 mkg/ml polybrene and incubated in 6-well plates in a humidified incubator at 37°C with 5% CO2 for 90 h with shaking every 24hr. Cells were harvested by centrifugation at 1000 rpm for 5 min at RT and analysed using FACSCalibur flow cytometer (BD Biosciences) using the CellQuest Pro program.

Western blot analysis

Cells were lysed in the presence of 1% NP-40 detergent and protein samples were prepared as described previously [31]. Protein contents were equalized after measurement with Bradford [32]. Total cell extracts were separated on 4–10% polyacrylamide gels (Bio-Rad) with sodium dodecyl sulfate (SDS) and were transferred to nitrocellulose membranes (Bio-Rad). Protein transfer was visualized by PonceauS staining [33]. Western blot signal intensities were quantitated with ImageJ software and normalized to the actin band density. Images of membranes were recorded by the system Gel Doc™ XR ChemiDoc™ XRS—Bio-Rad using Image Lab™ Software 170–82657. The following antibodies were used at 1:1000 dilutions if not otherwise noted: anti-PARP (AB3565, Millipore), anti-Phospho-IkappaBalpha Ser 32/36 (5A5, Cell Signaling), anti-IkappaB-alpha (C-21), anti-MCT4 (H-90), anti-beta-actin (I-19), Santa Cruz Biotechnology.

Human IL-8/CXCL8 protein secretion

This was analyzed with a DuoSet ELISA kit, catalog# DY208-05, R&D Systems, accordingly to a manufacturer protocol.

Annexin V and Propidium Iodine staining

This was analyzed with Annexin V-FITC Apoptosis Detection Kit, catalog#:K101-25, Biovision, accordingly to the manufacturer´s protocol.

Both, for ELISA and Annexin V staining, Raji or C666-1 cells were plated on a 6 well plate at 1x106 /well overnight before treatment. Cells were left untreated or exposed to either 10mM HCl or 10mM butyric acid, for 18hr, if not pointed otherwise.

Scratch wound assay

EBV positive and negative epithelial cell lines, CNE2 and HONE1, were seeded in a 6 well plate for 24hr to each 90% confluency at the day of treatment. Cells were then treated with SCFAs for 24hr. A scratch with a P10 pipette tip was made at the time of treatment [34,35]. Images of scratches were recorded by Olympus CKX 41 microscope at 4x magnification and processed by a DpxView Pro [1.14.3] software at 16hr. The remaining width of the scratch gap at 16hr was normalized to the gap measured at 0 hours.

Knock-down of butyrate transporters

siRNAs from Santa-Cruz was employed, against MCT1(sc-37235) and MCT4(sc-45892), with Lipofectamine RNAiMAX reagent, catalog#1377–100, Invitrogen, accordingly to the manufacturer´s protocol. In brief, 10nM siMCT1 or siMCT4, together with 1 ug plasmid, coding for a red fluorescent protein pTomato (cat#632532, Clontech, US), to monitor transfection efficiency, was introduced into C666-1 cells seeded one day earlier in a 6 well plate and treated 24 hrs later with either 10mM HCl or n-butyric acid.

Library construction, sequencing and analysis

The libraries were produced from the expressed mRNAs in parental HONE1 and stably LMP2A transfected NPC cell line before and after exposure to 10mM butyric acid for 24hr. cDNA libraries were prepared by Smartseq2 [36] protocol and sequenced on Illumina HiSeq 2000 platform. QC analysis was performed on the resulting sample reads using FastQC tool [37]. The reads were mapped to the UCSC reference human genome (version hg19) using STAR mapper [38]. Gene read counts and expression level were expressed in RPKMs (reads per million mapped read per kb of transcript length) and calculated based on uniquely mapped reads using an in-house pipeline [39].

Network enrichment analysis

The procedure of network enrichment analysis, NEA, [40] is similar to that of Gene Set Enrichment Analysis, GSEA [41]. In comparison to GSEA, NEA has the advantage of analyzing any DE genes in the global network, regardless of having a pathway/GO annotation. Another difference is that genes from the latter are analyzed in terms of network connectivity to the DE genes and do not have to be differentially expressed themselves. The software and the global interaction network for NEA were those presented by Merid et al., 2014 [42].

Results

Short chain fatty acids induced EBV lytic cycle in B-cell lymphoma and epithelial cancer derived cells

We tested expression of EBV immediate early genes in Raji cells treated with n-butyric acid. In comparison to untreated Raji cells (Fig 1, lanes 7–8) all treated cells showed significant increase in expression of BZLF1 and BRLF1, the EBV immediate early genes (two upper panels, Fig 1) in relation to expression of the house-keeping GAPDH gene and an internal control for EBV latent gene expression, the EBV transcript BART A (Fig 1, the lowest panel). TPA, an established inducer of EBV lytic replication, was used as a control of induction. The observed pattern of EBV lytic genes was not due to decrease of pH due to the exposure to acids, as in cells treated with low molarity of hydrochloric acid (HCl), with comparable pK, we did not observe induction of EBV IE genes expression (Fig 1, lanes 3–4). Thus, we confirmed that 10 mM n-butyric acid induced the early phase of the lytic cycle. We then used BZLF1 expression as a read-out for testing a whole panel of SCFAs on EBV-carrying C666-1 cells. In comparison to non-treated C666-1 cells (Fig 2, lane 9) the most significant upregulation of BZLF1, was obtained with 10 mM n-butyric acid (Fig 2, lane 1). The BZLF1 response to SCFAs in C666-1 cells was highly reproducible (Fig 3). We employed one more EBV model cell line to confirm the effect of n-butyric acid on EBV lytic cycle replication. We took the advantage of a recombinant EBV virus, which also carries the gene for green fluorescent protein (GFP). We exposed this AGS cell line carrying EBV-GFP to 10mM n-butyric acid, collected supernatants and applied these to the Raji B cell line, accordingly to the protocol, described in Material and Methods. After exposure to n-butyric acid the AGS cell supernatants resulted in infection of Raji, as indicated by appearance of GFP positive Raji cells, while no infection was seen with supernatants of mock or HCl treated AGS cells (Fig 4).

Fig 1

Expression of EBV lytic IE genes in Raji cells.

RT-PCR analysis upon overnight treatments with 10mM n-butyric acid, 10mM HCl or 20ng/ml TPA.

Fig 2

Analysis of short chain fatty acids induced BZLF1 gene expression in C666-1 NPC cells.

RT-PCR analysis upon overnight treatment with 10mM of short chain fatty acids as indicated. A. BZLF1 expression. B. BART A expression. C. GAPDH expression. Cells were treated with 10mM of SCFAs as indicated or were left untreated. D. Digital analysis of data from panels A-C, demonstrating BZLF1 expression in relation to GAPDH.

Fig 3

Analysis of n-butyric acid induced BZLF1 gene expression in C666-1 NPC cells.

RT-PCR analysis upon treatment with 10mM of n-butyric acid or HCl. A. BZLF1 expression. B. GAPDH expression. C. BART A expression. D. Digital analysis of relative BZLF1 expression calibrated to GAPDH, calculated from data in panels A-C.

Fig 4

N-butyric acid induces EBV virus production.

Flow cytometric analysis of GFP fluorescence in Raji cells upon superinfection with EBV-GFP recombinant virus produced by EBV-GFP carrying AGS epithelial cells, upon treatment with n-butyric acid as stated in the M&M.

N-Butyric acid induced inflammatory response in B-cell lymphoma and epithelial cancer cells

N-Butyric acid significantly induced expression of proinflammatory cytokines IL-6 and IL-8, which we demonstrated by qRT-PCR (Fig 5A) and ELISA analysis. We performed ELISA for the strongest induced cytokine, IL-8 (Fig 5B). While upregulation of IL-6 and IL-8 genes is detectable already after overnight exposure to butyric acid in both epithelial and BL cell lines, secretion of IL-8 cytokine, is increased only in C666-1 epithelial cells after 18hr with a peak at 48hrs. To substantiate our observation of expression of the inflammatory cytokines, we analysed activation of NFkB, a crucial transcription factor responsible for the cytokine expression. We detected a decrease in the expression of the NFkB inhibitory factor, IkB, in our cell lines of epithelial and B cell origin upon treatment with n-butyric acid (Fig 6) and an increase of NFkB expression to the same level as after TPA-exposure, an established inducer of NFkB activation [19].

Fig 5

N-butyric acid induces proinflammatory cytokine expression.

A. qRT-PCR analysis of Il-6 and IL-8 expression after overnight exposure to 10mM n-butyric acid. Values are adjusted to the EF1α gene expression. B. ELISA analysis of IL-8 secretion from C666-1 cells after 48hr of exposure to 10mM n-butyric acid.

Fig 6

N-butyric acid induces NFkappaB signaling.

Western blot analysis of phospho-IkappaB and NFkappaB expression in treated C666-1 (A) and Raji cells (B), adjusted to IkappaB and GAPDH expression, respectively.

N-Butyric acid induced apoptosis in epithelial and B cell lines

Cell death after SCFA treatment was evaluated counting live Raji cell numbers 24hr after the treatments. Only n-butyric acid and TPA induced EBV lytic cycle, while all SCFAs treatment regimens resulted in considerable decrease in the cell numbers (Fig 7A). Apoptosis was measured by flow cytometry analysis after propidium iodide and Annexin-V-FITC labeling (Fig 7B).

Fig 7

N-butyric acid induces apoptosis in Raji cells.

A.The number of live cells after 24hr of treatment as indicated. B. Flow cytometry dot plot from a representative experiment. A staple diagram in the lower panel represents the apoptotic cells population (Annexin-V-positive cells+Annexin-V/propidium iodide-positive cells) for each condition.

In addition, we evaluated PARP cleavage by western blot (Fig 8). Despite cell death after all treatment regimens, the PARP cleavage product—as a specific indicator of apoptosis—was detected only after 10mM n-butyric acid treatment of Raji and C666-1 cells.

Fig 8

N-butyric acid induces PARP cleavage.

Western blot analysis of cells treated with butyrate or TPA for 24hr, A-C) Raji cells, D-F) C666-1 cells. A,D) PonceauS staining. C,E) Actin. B,F) anti PARP Ab.

N-Butyric acid inhibited NPC cancer cells migration

To evaluate the impact of n-butyric acid on the behavior of highly migratory NPC cancer cells we analyzed C666-1 and HONE1 (data not shown) migration in a scratch wound assay. Reproducibly n-butyric acid inhibited C666-1 NPC cell migration, even at 0.5mM concentration, as judged by degree of gap closing in the scratch wound assay for C666-1 cells (Fig 9).

Fig 9

N-butyric acid inhibits migration of C666-1 NPC cells.

Evaluation of the gap distance in migration assay of C666-1 NPC cells upon treatment as indicated.

Expression of n-butyric acid transporters in epithelial and B cells

We investigated expression of MCT1 and MCT4, cellular SCFA transporters, in epithelial and B cells. We observed lower MCT1 and 4 expression in Raji B cells than in C666-1 epithelial cells (Fig 10A–10C). Treatment with n-butyric acid resulted in increase of the MCT expression in C666-1 cells, which was most evident for MCT4 (Fig 10A). We could not detect changes in MCT expression in Raji cells after treatment with n-butyric acid (Fig 10B). Activation of EBV lytic replication (BZLF1expression) decreased in response to butyric acid in cells, when MCT1 and MCT4-expression were knocked-down by siRNA (Fig 10D, the left panel). The expression of IL-8 cytokine was also decreased in these cells (Fig 10D, the right panel). This shows an important role for these transporters in butyric acid cell entrance.

Fig 10

N-butyric acid enters cells through MCT1 and 4 transporters.

A. RT-PCR analysis of MCT1 and 4 expression in C666-1 NPC cells upon n-butyric acid treatment during 24hr. B. RT-PCR analysis of MCT expression in Raji B cells upon different treatments during 24hr. Lanes represent the same as in A. C. Graph represents MCT1 and 4 gene expression in non treated cells, measured as pixel density and adjusted to GAPDH expression. D. The left panel. Expression of BZLF1, the viral lytic message, and cellular MCT transporters in C666-1 cells, adjusted to BART A and GAPDH, respectively, is decreased after knock-down of MCT, followed by n-butyric acid treatment. The right panel. qRT-PCR analysis of IL-8 expression in these cells. Pixel density values are compared to those after HCl treatment.

Gene expression analysis

In order to quantify changes in gene expression of NPC cells upon n-butyric acid exposure, we compared mRNA sequencing of parental HONE1 and HONE1 LMP2A transfected cells, since the LMP2A-mRNA is commonly detected in EBV associated malignancies [24]. A total of 15716 genes matching mRNA in HONE1 and 15331 genes matching those in HONE1 LMP2A cells passed preprocessing and mapping steps according to [39] and were used for the analysis of differential expression (S1 Table). We counted genes as differentially expressed when they 1) did not have exact '0' values in either column, 2) had a sum of RPKM values in two columns above 25 and 3) the RPKM values changed at least 2-fold in either direction. By using these criteria we obtained two sets of genes with altered expression (altered gene sets, AGS): 309 genes in HONE1 cells and 957 genes in LMP2A expressing HONE1 cells. These lists were analyzed for functional relations to signaling-related terms from Gene Ontology database using the method of Network Enrichment Analysis, NEA[40, 42). At the confidence threshold of false discovery rate (FDR) <0.01, we could identify in HONE1 cells four signaling pathways enriched in network connections to DE genes in HONE1 NPC cells after n-butyric acid treatment—all with highly relevant biological impact.

Thus, 424 gene members of the GO term “apoptosis” had 2834 links in total, i.e. in the global network. These genes had 627 links to the set of differentially expressed genes in our experiment (AGS), which was highly significant (Table 1).

Table 1

Genes expression of which changed in HONE1 cells treated with 10mM Butyric Acid for 24hr, were enriched for four different signaling-related GO terms.

Functional gene set	No. of genes in FGS	No.of links in FGS	No. of links between AGS and FGS	Network enrichment score	False discovery rate
GO_BP:CELL MIGRATION	93	662	167	11.87	5.03–03
GO_BP:APOPTOSIS	424	2834	627	11.50	5.76E-03
GO_BP:IMMUNE RESPONSE	232	1500	345	10.43	8.54E-03
GO_BP:KAPPAB_KINASE_NFKB_CASCADE	107	597	146	8.07	2.45E-02

Foot notes:

FGS–Functional Gene Set, a previously characterized group of genes with a common function.

FDR–False Discovery Rate of the network enrichment analysis

Score–the chi-squared score of network enrichment.

By analyzing enrichment of functional relations to KEGG pathways [42], we could see up-regulation of inflammatory Toll-like receptor signaling (S2 and S3 Tables) and pathways controlling cellular motility, such as axon guidance, actin cytoskeleton, focal adhesion and adherence junction as response to butyric acid (S2 and S3 Tables). To reveal effects of LMP2A expression in HONE1 NPC cells we set the threshold for RPKM fold selection to 50. LMP2A positive cells activate fewer signaling pathways (Fig 11, compare the upper and lower bubble diagrams).

Fig 11

N-butyric acid affects cellular gene expression in LMP2A expressing cells.

Network enrichment analysis of differentially expressed genes in HONE1 and HONE1 LMP2A cells upon n-butyric acid exposure.

Table 2 summarizes the signaling pathways differentially enriched in the parental NPC cells and the LMP2A transfected HONE1 NPC cells in response to n-butyric acid treatment.

Table 2

A comparison of the differentially activated signaling pathways in LMP2A positive cells.

HONE1	HONE1+LMP2A
KEGG_04010_MAPK_SIGNALING_PATHWAY	KEGG_04142_LYSOSOME
KEGG_04360_AXON_GUIDANCE	KEGG_04114_OOCYTE_MEIOSIS
KEGG_04320_DORSO-VENTRAL_AXIS_FORMATION	KEGG_04120_UBIQUITIN_MEDIATED_PROTEOLYSIS
KEGG_04662_B_CELL_RECEPTOR_SIGNALING_PAT	EGG_04930_TYPE_II_DIABETES_MELLITUS
KEGG_04612_ANTIGEN_PROCESSING_AND_PRESEN	KEGG_04960_ALDOSTERONE-REGULATED_SODIUM_
KEGG_04640_HEMATOPOIETIC_CELL_LINEAGE	KEGG_04910_INSULIN_SIGNALING_PATHWAY
KEGG_04650_NATURAL_KILLER_CELL_MEDIATED_
KEGG_04666_FC_GAMMA_R-MEDIATED_PHAGOCYTO
KEGG_04623_CYTOSOLIC_DNA-SENSING_PATHWAY
KEGG_04610_COMPLEMENT_AND_COAGULATION_CA

Discussion

SCFAs are important metabolic products of commensal microbiotas, e.g. in the nasopharynx or the gut. Produced at human mucosa the SCFAs have access to blood circulation and thus can also exert systemic effects. As it is known that some SCFAs have strong physiologic or metabolic effects on human host cells, it was of interest to investigate their effects on cell lines derived from cancers in human mucosal tissues. We employed latently EBV-infected cell lines both as a model and read out systems for SCFAs, as well as to evaluate the effect of SCFAs on EBV positive cells.

EBV-carrying cell lines from nasopharynx and EBV positive B cell lines were exposed to a range of physiological concentrations of a panel of SCFAs. We demonstrated that C666-1, Raji and AGS respond to SCFA by induction inflammatory genes and with induction of early steps of the EBV lytic cycle. It is known that EBV lytic replication is induced by Na-butyrate, and this effect has been linked to the effect of butyrate on histones. 10mM n-butyric acid significantly induces immediate early EBV genes BZLF1 and BRLF1 [43-45]. So did TPA, which we used as a positive control. Interestingly only n-butyric acid, but not TPA, also showed induction of apoptosis, evidenced by a decrease in the cell numbers, morphological changes and PARP cleavage. It has been suggested that EBV lytic induction is mechanistically linked to the triggering of cell death [46]. The concurrent induction of inflammation and apoptotic processes in our model systems points to different responses in two separate cell subpopulations. Although we observed upregulation of inflammatory gene expression in Raji BL cells, we could not detect IL-8 secretion from Raji cells upon exposure to butyric acid. Apparently, IL-8 mRNA produced in B cells without T cell help, would be quickly destroyed. It was more difficult to induce lytic replication in B cell lines than in epithelial cells with n-butyric acid.

In line with the our experimental data and a previous report [47], the effect of 10mM n-butyric acid on the signaling pathways regulating the migratory capacity of NPC cells was confirmed by gene expression analysis. We found that n-butyric acid inhibited migration of NPC cells independent of EBV status.

Butyric acid is actively transported in and out of cells by two specific transporters, MCT1 and MCT4 [48,49]. Epithelial cells express both and to a higher degree than the B cell line, which express mostly MCT1. The less efficient induction of EBV lytic cycle in Raji cells compared to epithelial cells upon n-butyric acid treatment correlated to differential expression of SCFA transporters in these cells. Recently, involvement of butyrate and MCT transporters in cell migration was demonstrated. Integrin b1 cooperates with MCT4 in migration of breast cancer cells [49]. Addition of butyric acid may induce MCT1 and/or MCT4 turnover and thus cause disruption of their cooperation with integrins in epithelial cells. A recent demonstration of the importance of proper spatial organization of cell surface receptors [50] corroborates our suggestion. Integrins a6b4, with which MCT1 and/or MCT4 may interact, were reported to contribute to NPC cell motility [51]. Interaction with butyric acid transporters MCT1 and MCT4 may increase the insight into a6b4 mediated migration.

The analysis of gene expression profiles of n-butyric acid treated and non treated HONE1 NPC cell lines, confirmed the significance of the cell signaling pathways studied, as major effects of the n-butyric acid response.

Virus-targeted lytic induction in EBV-associated malignancies may facilitate antiviral treatment and has been suggested [52]. Recently it was proposed to combine gemcitabine and valproic acid [53]. Valproic acid, a small organic acid is commonly used as an anti-epileptic drug and has several known adverse effects. N-butyric acid is a physiologically occurring SCFA and can be influenced by dietary effects on the composition of the microbiome. To regulate intestinal absorption and to establish an adequate pharmacokinetic approach will be challenges. Utilizing recent improvements in technology, both tasks are within reach and may be tried in a clinical setting.

SCFAs broadly affect host cell gene expression and programs. Several of the effects of butyric acid can be linked also to phenotypes and functions important for tumor cell biology. Possible adverse effects of SCFAs, always found also in normal physiologic conditions, may depend on concentrations and context. The inhibition of specific HDACs, like HDAC1,2,3, is most efficient in inducing apoptosis. The butyric acid inhibits exactly these HDACs, which also confer the activation of EBV [54]. The balance and level of microbial metabolites such as butyric acid may turn out to play a role in disease risks and progression. This can be influenced by dietary effects on the composition of the microflora.

Differential gene expression after n-butyric acid treatment.

(XLSX)

Click here for additional data file.

KEGG pathways in HONE1 NPC cells.

(DOCX)

Click here for additional data file.

KEGG pathways in HONE1 LMP2A NPC cells.

(DOCX)

Click here for additional data file.