- Research article
- Open Access
Butyrate induces profound changes in gene expression related to multiple signal pathways in bovine kidney epithelial cells
BMC Genomics volume 7, Article number: 234 (2006)
Global gene expression profiles of bovine kidney epithelial cells regulated by sodium butyrate were investigated with high-density oligonucleotide microarrays. The bovine microarray with 86,191 distinct 60mer oligonucleotides, each with 4 replicates, was designed and produced with Maskless Array Synthesizer technology. These oligonucleotides represent approximately 45,383 unique cattle sequences.
450 genes significantly regulated by butyrate with a median False Discovery Rate (FDR) = 0 % were identified. The majority of these genes were repressed by butyrate and associated with cell cycle control. The expression levels of 30 selected genes identified by the microarray were confirmed using real-time PCR. The results from real-time PCR positively correlated (R = 0.867) with the results from the microarray.
This study presented the genes related to multiple signal pathways such as cell cycle control and apoptosis. The profound changes in gene expression elucidate the molecular basis for the pleiotropic effects of butyrate on biological processes. These findings enable better recognition of the full range of beneficial roles butyrate may play during cattle energy metabolism, cell growth and proliferation, and possibly in fighting gastrointestinal pathogens.
The most common short-chain fatty acids (SCFA) are natural microbial fermentation products in the gastrointestinal tract. SCFA, including propionic, butyric and valeric acids, each with three, four and five carbons, respectively, contribute to the energy balance of all mammalian species [1, 2]. The major sources of these carbohydrates are hemicelluloses and fiber, which consists of plant cell wall polysaccharides such as cellulose and pectins . In ruminants, SCFA are a major energy source and contribute up to 70% of their energy requirements . Beyond their nutritional impact, SCFA, especially butyrate, have a multitude of cellular regulatory effects that modulate cell differentiation, proliferation, and motility. All three major components of SCFA induce apoptosis and inhibit cell proliferation, however, butyrate has the most potent effect [3, 4]. Roles for butyrate have been established in cell differentiation, proliferation, motility and in particular induction of cell cycle arrest and apoptosis [5–7]. Apoptosis is a genetically regulated cellular suicide mechanism that plays a crucial role in development and in the defense of homeostasis of animals . The mechanism(s) by which butyrate induces cellular differentiation and suppresses growth has not been elucidated. Similarly, the mechanism through which butyrate and other short-chain fatty acids induce the cell cycle regulatory and apoptotic effects and the mechanism by which the decision between cell death and survival is arbitrated are poorly understood [9–12].
In a previous study , potential biological roles of butyrate were investigated using the established Madin-Darby bovine kidney epithelial cell line (MDBK). The study focused on determining whether normal bovine cells in a standard cell culture condition were sensitive to the growth inhibitory effects of butyrate. The data suggested that sodium butyrate (SB) could induce apoptosis and cell cycle arrest in MDBK cells. Up to 38% of cells became apoptotic after 24 hours of treatment with 10 mM of butyrate. Butyrate also blocked the surviving cells at two distinct stages, G1 and M/G2. However, more studies are needed to better understand the relationship between butyrate and alterations in the expression of genes involved in cell cycle, apoptosis, and transcriptional regulation. Recent advances in high-throughput genomic tools such as microarray technology allowed us to examine the genome-wide effects of sodium butyrate on MDBK cells.
Butyrate Induces Cell Cycle Arrest and Hyperacetylation of Histone 3 in MDBK Cells
We previously reported that butyrate induced cell cycle arrest in MDBK cells . Prior to microarray analysis, the butyrate induced cell cycle arrest was reconfirmed. As shown in Figure 1, after butyrate treatment for 24 h, cell population profiles changed significantly. In the surviving cell population there was a significant increase in the number of cells in G1 (cells with 2C DNA contents) whereas those in S phase (cells falling between the 2C and 4C DNA contents) were decreased. This result confirmed our prior observation that cells were arrested at the G1/S boundary and DNA replication was blocked by the butyrate treatment. We also confirmed accumulation of acetylated histone 3 (H3) due to the butyrate treatment (Fig. 2). H3 acetylation was selected as the marker for the accumulation of acetylated histones because H3 is one of the core histones (H2A, H2B, H3, and H4), and is highly conserved across a wide range of organisms. The antibody against the acetyl H3(lys18) is also readily available commercially. To determine whether these biochemical attributes of butyrate were also observed in MDBK cells, two specific antibodies, the monoclonal antibody against the acetyl H3 and the monoclonal antibody against acetyl-phospho H3, were used to evaluate the histone deacetylase inhibitory activity of butyrate. Histone deacetylase (HDAC) inhibitors induce the hyperacetylation of nucleosomal histones. Evidence suggested that H3 phosphorylation is restricted to a small fraction of highly acetylated H3 histones. H3 phosphorylation is cell cycle dependent and may be associated with induced FOS and MYC oncogenes. No evidence indicates H3 phosphorylation directly results from HADC inhibitory activities. As shown in Fig. 2, butyrate treatment induced not only accumulation of hyperacetylation of H3 (Lys18) but also phosphor (Ser10)-acetyl (Lys9) H3. The quantified relative densities from the Western Blot show approximately a 2-fold increase of acetylated H3 protein and 2-fold increase of phosphor-acetyl H3 (Fig. 2B). Therefore it is reasonable to assume the increase detected by the antibody against phosph(Ser10)-Acetyl (lys9)-H3 is mostly due to the acetylation of H3 induced by butyrate.
Butyrate Induces Profound Changes in Gene Expression in MDBK Cells
In our previous study  and prior to the current microarray experiment, we monitored butyrate induced cell death and cell cycle arrest in MDBK cells in a time/dose-dependent manner using flow cytometry and Western blotting. For the purpose of this investigation and for determining the differential gene expression induced by butyrate, we selected a single dose proven to be able to generate maximum biological impacts in cell cycle arrest. Due to research budget constraints, we selected only a single time point as the initial global expression screening. We intend to verify the genes identified with the microarray in a time/dose-dependent manner using real-time RT-PCR in subsequent studies. We identified 450 genes significantly regulated by sodium butyrate at a very stringent false discovery rate (FDR) = 0% (see Additional file 1). However, when relaxing the stringency threshold to FDR = 10%, which should still be acceptable in most cases, there were 3662 genes significantly regulated (3662/45383 = 8%, data not shown). This percentage was consistent with a previous report  in which the authors used cDNA microarrays consisting of ~8000 sequences to demonstrate that approximately 7% of sequences assayed exhibited alteration by butyrate in human colon carcinoma cells.
Cell Cycle Control
The single largest category of genes regulated by butyrate (103) was cell cycle control related (Table 1). Butyrate repressed the vast majority of these genes including cyclins, cyclin-dependent kinases, histone deacetylases, helicases, chromosomal structure proteins as well as kinesins. Ubiquitin conjugating enzyme E2 C was also down-regulated. However, Max interacting protein 1 (MXI1) was up-regulated 12.6-fold by butyrate.
Apoptosis and Extracellular Matrix
Fourteen genes related to apoptosis and extracellular matrix (ECM) were significantly regulated by butyrate (Table 2). Generally, butyrate induced pro-apoptotic genes and repressed anti-apoptotic genes. Inhibin, beta A (INHBA) and adrenomedullin were among the induced genes while apoptosis inhibitors such as survivin and FAIM were down-regulated. Interestingly, Insulin-like growth factor 2 (IGF2) was up-regulated while IGF binding proteins 4 and 6 were down-regulated. Collagens, such as COL5A2 and COL3A1, and extracellular matrix protein Spondin 1 were repressed by butyrate while tissue inhibitor of metalloproteinase 2 (TIMP2) was significantly up-regulated.
Real-Time RT PCR
Thirty genes that represent different expression levels and functional classes were selected for real-time RT-PCR confirmation (Table 3). The real-time PCR data generally confirmed the microarray analysis. Linear regression analysis demonstrated a strong positive correlation between the two technological platforms with R = 0.867.
Cell cycle regulatory and apoptotic effects of butyrate and other short-chain fatty acids at the cellular and molecular levels in normal bovine cells have not been studied thus far but would serve as a principle launch point to validate the need for further study of these phenomena in cattle. Furthermore, utilization of particular signaling pathways in regulating cellular function or inducing gene expression appears to be dependent on, among other factors, the type of stimulus and cell examined. In our previous study, an important question asked and answered was whether normal bovine cells in a standard cell culture condition are sensitive to the growth inhibitory effects of butyrate. In this follow-up study, we utilize microarray technique to examine the genome-wide effects of sodium butyrate on MDBK cells as an important component of our efforts to understand the mechanisms of this phenomenon. Our data presented in this report indicate that the effects of butyrate are mediated through coordinated changes in gene expression that are the outcome of interactions between transduction pathways. We identified 450 genes significantly regulated by sodium butyrate at a very stringent false discovery rate (FDR) = 0%. However, since many of the genes have no direct links to the cell cycle arrest or apoptosis, their involvement in these biological effects certainly warrants more investigation. Therefore, instead of speculating, our discussion is concentrated primarily in cell cycle and apoptosis.
Sodium butyrate (SB) exerts a very broad range of effects on many biological pathways via its inhibitory ability on HDAC. SB is a potent inducer of a G1 cell cycle arrest. It also provokes apoptosis by activating both the death receptor and intrinsic apoptotic pathway. It regulates the cell cycle via down-regulation of cyclins and activation of CDK . In addition, it was reported that butyrate modulates host immune responses by activating neutrophils  and enhancing IL-4-dependent IgE production . Butyrate has immune suppression  and anti-inflammatory properties, in part by suppressing nuclear factor NF-κB activity [18–20]. Regulation of enzymes involved in cytoskeleton and cell membranes by butyrate has also been reported . The majority of its effects directly result from HDAC inhibition. However, butyrate appears to be involved in signal transduction via its own receptor GPR41, resulting in inosityl 1,4,5-triphosphate generation, intracellular Ca2+ release, ERK1/2 activation, and inhibition of cAMP accumulation . Because of their ability to inhibit cell proliferation, HDAC inhibitors (HDI) including butyrate have been extensively exploited as anti-tumor agents . Butyrate appears to have pleiotropic effects on various biological processes. The genes we identified (see Additional file 1), which fall within a broad range of functional categories, appeared to provide the molecular basis for its pleiotropic effects.
Butyrate seems able to inhibit all class I HDACs . It would be safe to assume that, like trichostatin A (TSA), butyrate can indeed inhibit the activities of HDAC8 and HDAC10, even though butyrate may have a different mechanism of action. Our studies suggested that butyrate indeed repressed histone deacetylase 8 (HDAC8) mRNA expression. The missing link is why this inhibition of enzymatic activities in turn down-regulates their own expression in mRNA levels. In mouse neural cells, it was observed that HDAC inhibitors affect the expression of HDACs themselves . In these cells, both TSA and SB indeed elevated the expression of HDAC1, HDAC3, HDAC5 and HDAC6 whereas mRNA levels of HDAC 2 and HDAC7 did not change. The mRNA levels of HDAC8 and HDAC10 were not detectable in these cells. While the mechanism and biological relevance of HDI regulation of HDAC expression remains unclear, it appears that there indeed exists an auto-regulatory feedback loop to the expression of several HDACs after their activities are inhibited.
It appears that the effects of HDI such as TSA and SB on MMP (matrix metalloproteinase) expression are specific to cell types. In mouse 3T3 fibroblasts, TSA represses MMP2 expression , while in human colonic cells DHD/K12, MMP production is inhibited by butyrate . In HT1080 tumor cells, both protein and mRNA levels of TIMP1, TIMP2, MMP2 and MMP9 are increased by butyrate treatment . Based on their data of limited modulation of MMP by butyrate in human SW1116 colon cancer cells , Emenaker et al suggested that SCFA, such as those derived from dietary fiber, may protect against invasive colon cancer through stimulation of TIMP and inhibition of uPA (plasminogen activator, urokinase) activities rather than their effects on MMP activities. In the present study, we found that both TIMP2 and MMPs, such as MMP1, MMP9 and MMP13, were induced by butyrate. The expression levels of TIMP2 induced by butyrate are similar (10.82, 11.28, and 14.05 folds, respectively) as detected by the three sequences that represent this gene on the microarray [TIGR: TC260960, TIGR: TC289409 and TIGR: TC289410]. However, changes in MMP levels were not considered significant due to the high stringency cutoff used in this study (FDR = 0%). The net effect of this induction on both MMPs and their inhibitor (TIMP2) is still unclear. TIMP2 can promote apoptosis in an in vivo colorectal cancer model yet can protect B16 melanoma cells from apoptosis . The elucidation of the mechanisms involved in controlling these distinctly opposing phenotypic effects of TIMP2 is of paramount importance.
Insulin-like growth factor binding proteins (IGFBP) modulate IGF action and regulate cell growth and apoptosis by preventing IGF from interacting with their own receptors. In our study, insulin-like growth factor (IGF-2) was upregulated by SB, which is consistent with other published data . Our microarray and real-time RT-PCR results confirmed that IGFBP6, which has a 100-fold higher affinity for IGF2 than IGF1 , was down-regulated by butyrate. It seems paradoxical that while IGF2 is up-regulated significantly by SB, its highest affinity binding protein is down-regulated. Our results may suggest different functions of various IGFBP members in regulating apoptosis and cell cycle progression. It would be intriguing to see how IGFBP exert their actions in cell growth and apoptosis via an IGF-independent fashion.
Despite the fact that the effect of HDAC inhibitors such as SB on the expression of cell cycle regulatory genes, such as cyclins and cyclin-related kinase (CDK), were documented and a few attempts were made to use high-throughput approaches, such as microarrays , differential display  and SAGE, to study the effects of HDAC inhibitors on cell cycle control, the extent of the effect of these inhibitors on cell proliferation and cell cycle has not been fully realized in part due to limited gene representation on these microarrays [19, 30, 33, 34]. Due to significant differences in gene representation, species, microarray platforms, and data analysis tools used in these studies, a direct comparison between our results and those published is seemingly difficult. In general, our results are in good agreement with published reports. For example, we detected ~8% of all genes were significantly regulated by butyrate in bovine MDBK cells, which is consistent with a previous report  in which the authors used cDNA microarray consisting of ~8000 sequences to demonstrate that approximately 7% of sequences assayed exhibited alteration by butyrate in human colon carcinoma cells. Down-regulation of cyclins, PCNA, CDKs, and upregulation of IGF2, MMPs, and TIMP2 were also confirmed by previous reports [19, 34]. However, many genes, such as Aurora kinases, BUB1 and BUB1B, centromere proteins, kinesins, Max interacting protein 1 (MXI1), minichromosomal maintenance deficient proteins (MCMs), and spindle pole body components were not previously recognized to be regulated by SB. Our efforts in this study are among the first to systematically categorize the butyrate-regulated genes related with cell cycle control with a genome-wide approach in farm animals. While the vast majority of these genes are down-regulated, MXI1 is up-regulated. As a key component of the mitotic checkpoint, MXI1 binds with MAX to form a sequence-specific DNA-binding protein complex and acts as a transcriptional repressor. Up-regulation of MXI1 by SB could result in down-regulation of cyclins, which in turn negatively regulates centromere proteins.
Accumulation of cells with 2C and 4C DNA contents suggests inhibition by butyrate of cell cycle at G1 and M/G2 phases and suggests that a common responding element in genes responsive to the treatment of butyrate is required for progression of both phases G1 and M/G2. This observation is also consistent with the previous report that the inhibition of G1 progression by butyrate is not restricted to a specific mitogenic signaling pathway , but may also include the inhibitory effect on initiation of DNA replication. In this report, multiple genes such as minichromosome maintenance (MCM) proteins 2, 3, 4, 5, and 6, as well as Orc1 (Origin Recognition Complex largest subunit) are significantly down-regulated. The products of these genes has been shown to be rate limiting for initiation of DNA replication in eukaryotic cell lines and are essential for the assembly of the pre-replication complex (pre-RC) [36–38]. This finding indicates that in someway, butyrate treatment directly targets these genes and down regulates the genes that are essential for initiation of DNA replication. CDC2/Cdk1 and related cyclins are also significantly down-regulated. At least four roles have now been recognized for cyclin A activated Cdk1 protein kinases in regulating cell cycle events. First, Cdk2/cyclin A is responsible for activating pre-replication complexes at the beginning of S-phase in order to begin DNA synthesis . Second, Cdk1/cyclin A inhibits assembly of new pre-replication complexes during S-phase  by inactivating Cdc6 . Third, Cdk2/cyclin A is required for the G2 to M-phase transition . Finally, Li and DePamphilis  recently revealed that Cdk1/cyclin A is required for preventing Orc1 in mammals or ORC in Xenopus from binding to chromatin during mitosis. Our results show that targeted destruction of cdc6 and cdc2/cdk1 may be involved in the apoptosis and cell cycle arrest induced by butyrate. They are consistent with the growing body of evidence suggesting that disruption of the coordination between regulation of DNA synthesis and cyclin-dependent kinase activity is an important feature of apoptosis. It is of importance and interest for us to understand the mechanism(s) of how butyrate targets these genes and causes these changes in gene expressions. However, it will require a great deal of efforts and certainly is out of the scope of this report.
Understanding the mechanism of butyrate in affecting apoptosis, cell proliferation and differentiation and especially the difference of its mode of action compared to other HDAC inhibitors will facilitate designing novel and more potent HDAC inhibitors that target specific cellular processes that are dysregulated in neoplastic cells. Dissection of the pathways regulated by butyrate will provide a basis for better utilization of its anti-tumor, anti-metastatic, immune-mediating, and anti-inflammatory properties. In addition, butyrate has been shown to up-regulate transcription levels of muslins, which are the major components of gastrointestinal mucosa considered to be the first line of defense against pathogens . Because of its ability to affect functions of monocyte-derived dendritic cells and macrophage  and to activate neutrophils, butyrate and its related biological pathways could be manipulated to fight against gastrointestinal pathogens. Up-regulation of glutathione S-transferases [43, 44], genes known to be involved in defense against oxidative stress, by butyrate provides evidence of a favorable modulation of toxicological defense systems. The results presented in this paper not only help to better dissect HDAC inhibition and mechanism in apoptosis and cell cycle control but also help to understand ruminant metabolism and physiology, which in turn could lead to improvement in energy efficiency of cattle.
The present research identified 450 genes significantly regulated by sodium butyrate in bovine kidney epithelial cells. The genes related to multiple signal pathways such as cell cycle control and apoptosis were presented. The profound changes in gene expression elucidate the molecular basis for the pleiotropic effects of butyrate on biological processes. The results presented in this paper can provide clues on the mechanism of histone deacetylase inhibition by butyrate and resulting alterations in the expression of genes involved in cell cycle, apoptosis, and transcriptional regulation. Since butyrate functions as both a nutrient and signaling molecule regulating the cell growth and proliferation, these findings enable better recognition of the full range of roles butyrate may play during cattle energy metabolism, cell growth and proliferation.
Cell Culture and Treatments
The Madin-Darby bovine kidney epithelial cells (MDBK, American Type Culture Collection, Manassas, VA., Catalog No. CCL-22) were cultured in Eagle's minimal essential medium supplemented with 5% fetal bovine serum (Invitrogen, Carlsbad, CA) in 25 cm2 flasks with medium renewal twice per week. Cell cultures were maintained in a water-jacked incubator with 5% CO2 at 37°C. Sub-cultivations were performed when cells attained 80 to 90 % confluence, according to the product information supplied by American Type Culture Collection. At approximately 50% confluence (during the exponential phase), the cells were treated with 10 mM of sodium butyrate for 24 h (Calbiochem, San Diego, CA. Three replicate flasks of cells for both treatment and control groups (a total of 6 samples) were used for the flow cytometry and microarray experiments. The harvested cells were snap frozen in liquid N2 and stored at -80°C until RNA extraction.
Flow Cytometric Analysis of Cells
The detailed procedures were described in a previous publication . Briefly, cells collected by trypsinization were washed and resuspended in PBS buffer. Two volumes of ice-cold 100% ethanol were added drop wise into tubes and mixed with cells in suspension by slow vertexing. After incubation with RNase I, cells were then stained with propidium iodide (PI). Measuring the fluorescence by flow cytometry provided a measure of the amount of PI taken up by the cells and, indirectly, the amount of DNA content. Cell DNA content was analyzed using a flow cytometer (FC500, Beckman Coulter, Palatine, IL) and collected data were analyzed using Cytomics RXP (Beckman Coulter). At least 10,000 cells per sample were analyzed.
Preparation of Cell Extracts and Western Blot Analysis
Preparations of cells and cell extracts, SDS-PAGE and Western Blot analysis were described previously . Briefly, the protein from different samples was separated by SDS PAGE on two identical 4 to 20% polyacrylamide gradient gels. One gel was stained with SimpleBlue (Invitrogen) and one was transferred to a membrane and probed with monoclonal anti- acetyl-phospho H3 and anti- acetyl H3 antibodies. The target bands on the Western Blots from three experiments were quantified with a NIH Image software. The relative densities were measured and corrected with the stained protein density.
Isolation of Total RNA
Total RNA was extracted using Trizol following the manufacturer's recommendations (Invitrogen). Trace genomic DNA in the crude total RNA samples was removed by incubation with 4–10 units DNase I per 100 μg total RNA (Ambion, Austin, TX) at 37°C for 30 min. Total RNA was further purified using an RNeasy Mini kit (Qiagen, Valenica, CA). The concentration of the total RNA was determined using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Rockland, DE) and RNA integrity was verified using a Bioanalyzer 1000 (Agilent, Palo Alto, CA).
Generation of Biotin-labeled cRNA
Biotin-labeled cRNA was generated with a modified procedure of the Superscript Choice System (Invitrogen) for double-strand (ds) cDNA synthesis followed by in vitro transcription. Briefly, the 1st strand cDNA was synthesized from 4.0 μg total RNA by 1.0 unit SuperScript II reverse transcriptase (Invitrogen) in the presence of 100 pmoles T7 promoter Oligo dT primer. After 2nd strand synthesis, the DNA was purified with a DNA Clean & Concentrator-5 kit (Zymo Research, Orange, CA) and eluted with 8 to 16 μl of deionized (dd) H2O. The recovered ds cDNA was further concentrated down to 3 μl by a speed vacuum device. cRNA was synthesized with a MEGAscript in vitro Transcription kit (Ambion). The in vitro transcription reaction was carried out in a total volume of 23.0 μl consisting of 3.0 μl of ds cDNA, 2.3 μl 10X Ambion reaction buffer, 2.3 μl 10X Ambion T7 enzyme mix, and 15.4 μl NTP labeling mix (7.5 mM ATP, 7.5 mM GTP, 5.625 mM UTP, 5.625 mM CTP and 1.875 mM biotin-16-UTP and 1.875 mM biotin-11 CTP). The in vitro transcription reaction was incubated at 37°C for ~16 hours in a thermocycler. The cRNA was purified with an RNeasy mini kit (Qiagen). Generally, 40 to 60 μg of cRNA can be obtained from 4.0 μg of input total RNA. The size range of the cRNA, expected to be between 300 to 3000 bp with the maximum intensity centered at least 1000 bp, was verified using a Bioanalyzer 1000. The biotinylated cRNA was fragmented to 50 to 200 bp by heating cRNA in a buffer consisting of 40 mM Tris-acetate, pH 8.0, 100 mM potassium acetate, and 30 mM magnesium acetate at 95°C for 35 min.
Oligonucleotide Microarray, Hybridization, Image Acquisition and Data Analysis
The bovine microarray platform used was described previously . Briefly, a total of 86,191 unique 60mer oligonucleotides were designed and synthesized it in situ using photo deprotection chemistry . Each unique oligonucleotide was repeated 4 times on the array (a total of ~340,000 features). These oligonucleotides represented 45,383 unique bovine sequences/genes, including 40,808 Tentative Consensus sequences (TCs) from TIGR Bos taurus gene index  and 4,575 singletons.
The microarrays were pre-hybridized with 1X MES hybridization buffer (100 mM MES, 1.0 M Na+, 20 mM EDTA, 0.01% Tween20), 40 μg herring sperm DNA and 200 μg acetylated BSA at 45°C for 15 min followed by hybridization with 10 μg denatured and fragmented cRNA per microarray at 45°C for 16 – 20 h with constant rotation. After hybridization, the microarrays were immediately washed extensively under non-stringent conditions (6× SSPE, 0.01% Tween20) at room temperature (RT) followed by a stringent wash (100 mM MES salt and free acid solution, 0.1 M Na+, 0.01% Tween20) at 45°C. After the final rinse with the non-stringent wash buffer, the microarrays were stained with 1× Stain buffer (100 mM MES, 1 M Na+, 0.05% Tween20, 50 mg/ml of BSA, and 1 mg/ml of Cy3-streptavidin) at RT for 25 min. The stain buffer was removed and the microarrays were rinsed once more with non-stringent wash buffer. The microarrays were immediately dried under a stream of argon gas and scanned using an Axon GenePix 4000B scanner (Molecular Devices Corp., Union City, CA) at 5 μM resolution. The data were extracted from the raw images using NimbleScan software (NimbleGen, Madison, WI). The control and butyrate treatment each had 3 replicates and a total of 6 microarrays were used in the experiment (GEO Accession GSE3970). The microarray data are available as accession GSE3970 in the Gene Expression Omnibus repository at the National Center for Biotechnology Information .
Relative signal intensities (log2) for each feature were generated using the Robust Multi-Array Average (RMA) algorithm [49, 50]. The data were processed based on quantile normalization method  using the R package . This normalization method aims to make the distribution of intensities for each array in a set of arrays the same. The method assumes that a quantile-quantile plot of two data vectors with the same distribution will have a straight diagonal line. The method performed better in dealing with bias and reducing variability across arrays compared to other methods . The background-adjusted, normalized, and log transformed intensity values were then analyzed using the Significance Analysis of Microarrays method  with two-class unpaired design (SAM version 2.20 ). SAM is the most popular method for microarray analysis with 635 citations of the original publication as of October 2004 . SAM ranks genes based on a modified t-test statistic. The unique features of SAM include implementing permutation testing, and the ability to estimate a global false discovery rate (FDR, an expected percentage of false positives among the claimed positives) and a gene error chance (q-value). A sequence was declared to be significant when it met a stringent median false discovery rate (FDR) cutoff at 0 % (see Additional file 1). A BLAST search was conducted for all sequences that met the threshold to remove possible redundancy. When a gene was represented by multiple sequences, the fold change with q value of only one sequence was selected to represent this gene.
Real-time RT-PCR analysis was carried out with the iQ SYBR Green Supermix kit (Biorad) using 200 nM of each amplification primer (see Additional file 2) and the 1st-strand cDNA (80 ng of the input total RNA equivalents) in a 25 μl reaction volume. The amplification was carried out on an iCycler iQ™ Real Time PCR Detection System (BioRad) with the following profile: 95°C – 60s; 40 cycles of 94°C-15s, 60°C -30s, and 72°C -30s. The melting curve analysis was performed for each primer pair. Expression levels of β-actin remained constant (within 0.5 Ct between samples) and were used as endogenous controls. Relative gene expression data were calculated using the 2ΔΔC Tmethod .
Miller SJ: Cellular and physiological effects of short-chain fatty acids. Mini Rev Med Chem. 2004, 4 (8): 839-845.
Sharp WM, Johnson RR, Owens FN: Ruminal VFA production with steers fed whole or ground corn grain. J Anim Sci. 1982, 55 (6): 1505-1514.
Emenaker NJ, Calaf GM, Cox D, Basson MD, Qureshi N: Short-chain fatty acids inhibit invasive human colon cancer by modulating uPA, TIMP-1, TIMP-2, mutant p53, Bcl-2, Bax, p21 and PCNA protein expression in an in vitro cell culture model. J Nutr. 2001, 131 (11 Suppl): 3041S-6S.
Hague A, Paraskeva C: The short-chain fatty acid butyrate induces apoptosis in colorectal tumour cell lines. Eur J Cancer Prev. 1995, 4 (5): 359-364.
Chen JS, Faller DV, Spanjaard RA: Short-chain fatty acid inhibitors of histone deacetylases: promising anticancer therapeutics?. Curr Cancer Drug Targets. 2003, 3 (3): 219-236. 10.2174/1568009033481994.
Gassull MA, Cabre E: Nutrition in inflammatory bowel disease. Curr Opin Clin Nutr Metab Care. 2001, 4 (6): 561-569. 10.1097/00075197-200111000-00018.
Scheppach W, Bartram HP, Richter F: Role of short-chain fatty acids in the prevention of colorectal cancer. Eur J Cancer. 1995, 31A (7-8): 1077-1080. 10.1016/0959-8049(95)00165-F.
Johnson IT: Anticarcinogenic effects of diet-related apoptosis in the colorectal mucosa. Food Chem Toxicol. 2002, 40 (8): 1171-1178. 10.1016/S0278-6915(02)00051-0.
Dashwood RH, Myzak MC, Ho E: Dietary HDAC inhibitors: time to rethink weak ligands in cancer chemoprevention?. Carcinogenesis. 2006, 27 (2): 344-349. 10.1093/carcin/bgi253.
Saha RN, Pahan K: HATs and HDACs in neurodegeneration: a tale of disconcerted acetylation homeostasis. Cell Death Differ. 2006, 13 (4): 539-550. 10.1038/sj.cdd.4401769.
Scheppach W, Weiler F: The butyrate story: old wine in new bottles?. Curr Opin Clin Nutr Metab Care. 2004, 7 (5): 563-567. 10.1097/00075197-200409000-00009.
Rupnarain C, Dlamini Z, Naicker S, Bhoola K: Colon cancer: genomics and apoptotic events. Biol Chem. 2004, 385 (6): 449-464. 10.1515/BC.2004.053.
Li CJ, Elsasser TH: Butyrate-induced apoptosis and cell cycle arrest in bovine kidney epithelial cells: involvement of caspase and proteasome pathways. J Anim Sci. 2005, 83 (1): 89-97.
Li CJ, Vassilev A, DePamphilis ML: Role for Cdk1 (Cdc2)/cyclin A in preventing the mammalian origin recognition complex's largest subunit (Orc1) from binding to chromatin during mitosis. Mol Cell Biol. 2004, 24 (13): 5875-5886. 10.1128/MCB.24.13.5875-5886.2004.
Kasama T, Miwa Y, Isozaki T, Odai T, Adachi M, Kunkel SL: Neutrophil-derived cytokines: potential therapeutic targets in inflammation. Curr Drug Targets Inflamm Allergy. 2005, 4 (3): 273-279. 10.2174/1568010054022114.
Yamamoto I, Matsunaga T, Sakata K, Nakamura Y, Doi S, Hanmyou F: Histone hyperacetylation plays a role in augmentation of IL-4-induced IgE production in LPS-stimulated murine B-lymphocytes by sodium butyrate. J Biochem (Tokyo). 1996, 119 (6): 1056-1061.
Millard AL, Mertes PM, Ittelet D, Villard F, Jeannesson P, Bernard J: Butyrate affects differentiation, maturation and function of human monocyte-derived dendritic cells and macrophages. Clin Exp Immunol. 2002, 130 (2): 245-255. 10.1046/j.0009-9104.2002.01977.x.
Ahmad MS, Krishnan S, Ramakrishna BS, Mathan M, Pulimood AB, Murthy SN: Butyrate and glucose metabolism by colonocytes in experimental colitis in mice. Gut. 2000, 46 (4): 493-499. 10.1136/gut.46.4.493.
Joseph J, Mudduluru G, Antony S, Vashistha S, Ajitkumar P, Somasundaram K: Expression profiling of sodium butyrate (NaB)-treated cells: identification of regulation of genes related to cytokine signaling and cancer metastasis by NaB. Oncogene. 2004, 23 (37): 6304-6315. 10.1038/sj.onc.1207852.
Segain JP, Raingeard de la Bletiere D, Bourreille A, Leray V, Gervois N, Rosales C, Ferrier L, Bonnet C, Blottiere HM, Galmiche JP: Butyrate inhibits inflammatory responses through NFkappaB inhibition: implications for Crohn's disease. Gut. 2000, 47 (3): 397-403. 10.1136/gut.47.3.397.
Prasad KN, Sinha PK: Effect of sodium butyrate on mammalian cells in culture: a review. In Vitro. 1976, 12 (2): 125-132.
Le Poul E, Loison C, Struyf S, Springael JY, Lannoy V, Decobecq ME, Brezillon S, Dupriez V, Vassart G, Van Damme J, Parmentier M, Detheux M: Functional characterization of human receptors for short chain fatty acids and their role in polymorphonuclear cell activation. J Biol Chem. 2003, 278 (28): 25481-25489. 10.1074/jbc.M301403200.
Piekarz R, Bates S: A review of depsipeptide and other histone deacetylase inhibitors in clinical trials. Curr Pharm Des. 2004, 10 (19): 2289-2298. 10.2174/1381612043383980.
Bertos NR, Wang AH, Yang XJ: Class II histone deacetylases: structure, function, and regulation. Biochem Cell Biol. 2001, 79 (3): 243-252. 10.1139/bcb-79-3-243.
Ailenberg M, Silverman M: Differential effects of trichostatin A on gelatinase A expression in 3T3 fibroblasts and HT-1080 fibrosarcoma cells: implications for use of TSA in cancer therapy. Biochem Biophys Res Commun. 2003, 302 (2): 181-185. 10.1016/S0006-291X(03)00150-5.
Aparicio A, Gardner A, Tu Y, Savage A, Berenson J, Lichtenstein A: In vitro cytoreductive effects on multiple myeloma cells induced by bisphosphonates. Leukemia. 1998, 12 (2): 220-229. 10.1038/sj.leu.2400892.
Zeng H, Briske-Anderson M: Prolonged butyrate treatment inhibits the migration and invasion potential of HT1080 tumor cells. J Nutr. 2005, 135 (2): 291-295.
Emenaker NJ, Basson MD: Short chain fatty acids inhibit human (SW1116) colon cancer cell invasion by reducing urokinase plasminogen activator activity and stimulating TIMP-1 and TIMP-2 activities, rather than via MMP modulation. J Surg Res. 1998, 76 (1): 41-46. 10.1006/jsre.1998.5279.
Valente P, Fassina G, Melchiori A, Masiello L, Cilli M, Vacca A, Onisto M, Santi L, Stetler-Stevenson WG, Albini A: TIMP-2 over-expression reduces invasion and angiogenesis and protects B16F10 melanoma cells from apoptosis. Int J Cancer. 1998, 75 (2): 246-253. 10.1002/(SICI)1097-0215(19980119)75:2<246::AID-IJC13>3.0.CO;2-B.
Chiba T, Yokosuka O, Fukai K, Kojima H, Tada M, Arai M, Imazeki F, Saisho H: Cell growth inhibition and gene expression induced by the histone deacetylase inhibitor, trichostatin A, on human hepatoma cells. Oncology. 2004, 66 (6): 481-491. 10.1159/000079503.
Baxter RC, Saunders H: Radioimmunossay of insulin-like growth factor-binding protein-6 in human serum and other body fluids. J Endocrinol. 1992, 134 (1): 133-139.
Basson MD, Liu YW, Hanly AM, Emenaker NJ, Shenoy SG, Gould Rothberg BE: Identification and comparative analysis of human colonocyte short-chain fatty acid response genes. J Gastrointest Surg. 2000, 4 (5): 501-512. 10.1016/S1091-255X(00)80093-1.
Mariadason JM, Corner GA, Augenlicht LH: Genetic reprogramming in pathways of colonic cell maturation induced by short chain fatty acids: comparison with trichostatin A, sulindac, and curcumin and implications for chemoprevention of colon cancer. Cancer Res. 2000, 60 (16): 4561-4572.
Tabuchi Y, Takasaki I, Doi T, Ishii Y, Sakai H, Kondo T: Genetic networks responsive to sodium butyrate in colonic epithelial cells. FEBS Lett. 2006, 580 (13): 3035-3041. 10.1016/j.febslet.2006.04.048.
Charollais RH, Buquet C, Mester J: Butyrate blocks the accumulation of CDC2 mRNA in late G1 phase but inhibits both the early and late G1 progression in chemically transformed mouse fibroblasts BP-A31. J Cell Physiol. 1990, 145 (1): 46-52. 10.1002/jcp.1041450108.
Sun WH, Coleman TR, DePamphilis ML: Cell cycle-dependent regulation of the association between origin recognition proteins and somatic cell chromatin. Embo J. 2002, 21 (6): 1437-1446. 10.1093/emboj/21.6.1437.
Li CJ, Bogan JA, Natale DA, DePamphilis ML: Selective activation of pre-replication complexes in vitro at specific sites in mammalian nuclei. J Cell Sci. 2000, 113 ( Pt 5): 887-898.
Li CJ, DePamphilis ML: Mammalian Orc1 protein is selectively released from chromatin and ubiquitinated during the S-to-M transition in the cell division cycle. Mol Cell Biol. 2002, 22 (1): 105-116. 10.1128/MCB.22.1.105-116.2002.
Coverley D, Laman H, Laskey RA: Distinct roles for cyclins E and A during DNA replication complex assembly and activation. Nat Cell Biol. 2002, 4 (7): 523-528. 10.1038/ncb813.
Jiang W, Wells NJ, Hunter T: Multistep regulation of DNA replication by Cdk phosphorylation of HsCdc6. Proc Natl Acad Sci U S A. 1999, 96 (11): 6193-6198. 10.1073/pnas.96.11.6193.
Furuno N, den Elzen N, Pines J: Human cyclin A is required for mitosis until mid prophase. J Cell Biol. 1999, 147 (2): 295-306. 10.1083/jcb.147.2.295.
Gaudier E, Forestier L, Gouyer V, Huet G, Julien R, Hoebler C: Butyrate regulation of glycosylation-related gene expression: evidence for galectin-1 upregulation in human intestinal epithelial goblet cells. Biochem Biophys Res Commun. 2004, 325 (3): 1044-1051. 10.1016/j.bbrc.2004.10.141.
Pool-Zobel BL, Selvaraju V, Sauer J, Kautenburger T, Kiefer J, Richter KK, Soom M, Wolfl S: Butyrate may enhance toxicological defence in primary, adenoma and tumor human colon cells by favourably modulating expression of glutathione S-transferases genes, an approach in nutrigenomics. Carcinogenesis. 2005, 26 (6): 1064-1076. 10.1093/carcin/bgi059.
Ranganna K, Yousefipour Z, Yatsu FM, Milton SG, Hayes BE: Gene expression profile of butyrate-inhibited vascular smooth muscle cell proliferation. Mol Cell Biochem. 2003, 254 (1-2): 21-36. 10.1023/A:1027383710582.
Li RW, Meyer MJ, Van Tassel CP, Sonstegard TS, Conner EE, Van Amburgh ME, Boisclair YR, Capuco AV: Identification of estrogen-responsive genes in the parenchyma and fat pad of the bovine mammary gland by microarray analysis. Physiol Genomics. 2006, In press:
Singh-Gasson S, Green RD, Yue Y, Nelson C, Blattner F, Sussman MR, Cerrina F: Maskless fabrication of light-directed oligonucleotide microarrays using a digital micromirror array. Nat Biotechnol. 1999, 17 (10): 974-978. 10.1038/13664.
National Center for Biotechnology Information. [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/geo]
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249.
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31 (4): e15-10.1093/nar/gng015.
Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.
Department of Statistics-Stanford University. [http://www.stat.stanford.edu/ntibs/clickwrap/sam]
Larsson O, Wahlestedt C, Timmons JA: Considerations when using the significance analysis of microarrays (SAM) algorithm. BMC Bioinformatics. 2005, 6 (1): 129-10.1186/1471-2105-6-129.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
The authors thank Joy Castano for her excellent technical assistance. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U. S. Department of Agriculture.
RWL conceived the study, carried out the microarray experiment, provided data analysis, and drafted the manuscript. CJL cultured the cells, carried out the western blotting and flow cytometry analysis, participated in data mining, and helped to draft the manuscript. Both authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Supplementary table 1. This data represents the 466 sequences identified to be significant at FDR = 0%, which represents 450 cattle genes. (XLS 88 KB)
Additional file 2: Supplementary table 2. This data represents the forward/reverse primer sequences for the genes selected for real-time PCR confirmation. (XLS 20 KB)
About this article
Cite this article
Li, R.W., Li, C. Butyrate induces profound changes in gene expression related to multiple signal pathways in bovine kidney epithelial cells. BMC Genomics 7, 234 (2006). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-7-234
- Cell Cycle Control
- Sodium Butyrate
- Sodium Butyrate
- MDBK Cell