Alternative Splicing Regulated by Butyrate in Bovine Epithelial Cells

doi:10.1371/journal.pone.0039182

PLoS One. 2012; 7(6): e39182.

Published online 2012 June 14. doi: 10.1371/journal.pone.0039182

PMCID: PMC3375255

Copyright This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

Alternative Splicing Regulated by Butyrate in Bovine Epithelial Cells

Sitao Wu,¹ Congjun Li,² Wen Huang,³ Weizhong Li,¹ and Robert W. Li²^*

¹Center for Research in Biological Systems, University of California San Diego, San Diego, California, United States of America

²USDA-ARS, Bovine Functional Genomics Laboratory, Beltsville, Maryland, United States of America

³Department of Genetics, North Carolina State University, Raleigh, North Carolina, United States of America

I. King Jordan, Editor

Georgia Institute of Technology, United States of America

* E-mail: robert.li/at/ars.usda.gov

Conceived and designed the experiments: RWL CL. Performed the experiments: CL RWL. Analyzed the data: SW WH WL RWL. Wrote the paper: RWL.

Received March 30, 2012; Accepted May 21, 2012.

Abstract

As a signaling molecule and an inhibitor of histone deacetylases (HDACs), butyrate exerts its impact on a broad range of biological processes, such as apoptosis and cell proliferation, in addition to its critical role in energy metabolism in ruminants. This study examined the effect of butyrate on alternative splicing in bovine epithelial cells using RNA-seq technology. Junction reads account for 11.28 and 12.32% of total mapped reads between the butyrate-treated (BT) and control (CT) groups. 201,326 potential splicing junctions detected were supported by ≥3 junction reads. Approximately 94% of these junctions conformed to the consensus sequence (GT/AG) while ~3% were GC/AG junctions. No AT/AC junctions were observed. A total of 2,834 exon skipping events, supported by a minimum of 3 junction reads, were detected. At least 7 genes, their mRNA expression significantly affected by butyrate, also had exon skipping events differentially regulated by butyrate. Furthermore, COL5A3, which was induced 310-fold by butyrate (FDR <0.001) at the gene level, had a significantly higher number of junction reads mapped to Exon#8 (Donor) and Exon#11 (Acceptor) in BT. This event had the potential to result in the formation of a COL5A3 mRNA isoform with 2 of the 69 exons missing. In addition, 216 differentially expressed transcript isoforms regulated by butyrate were detected. For example, Isoform 1 of ORC1 was strongly repressed by butyrate while Isoform 2 remained unchanged. Butyrate physically binds to and inhibits all zinc-dependent HDACs except HDAC6 and HDAC10. Our results provided evidence that butyrate also regulated deacetylase activities of classical HDACs via its transcriptional control. Moreover, thirteen gene fusion events differentially affected by butyrate were identified. Our results provided a snapshot into complex transcriptome dynamics regulated by butyrate, which will facilitate our understanding of the biological effects of butyrate and other HDAC inhibitors.

Introduction

Alternative splicing (AS) is an important driving force behind vast protein diversity and the evolution of phenotypic complexity in mammals. AS generates various protein isoforms from single genes with different biological properties [1] and plays a critical role in development and disease [2]–[4]. Furthermore, AS is under strict regulatory control, including an accurate recognition of the splice junction that defines intron and exon boundary [5]. AS is mediated by elaborate molecular machinery, the spliceosome, which consists of at least 5 small nuclear RNAs (snRNAs) and over 100 accessory proteins [6]. Mutations in the core element of the spliceosome cause disease [7]. Up to 94% of human genes are estimated to undergo AS [4]; hence, aberrant AS is frequently associated with numerous human diseases [1]. Therefore, splicing modulation has been touted as a therapeutic means for treating genetic diseases caused by splicing mutations [8].

Butyrate is a preferred substrate for gut epithelial cells. In ruminants, butyrate contributes to 70% of energy metabolism. In monogastric species, butyrate also plays an important role in energy metabolism in the hindgut [9]. Moreover, butyrate acts as a histone deacetylase inhibitor and can induce apoptosis and inhibit cell proliferation in vitro [10]. As a signaling molecule, butyrate induces a profound change in gene expression and results in a significant change in transcript abundance of ~50% genes transcribed in the transcriptome of epithelial cells [11]. The effect of butyrate on growth rate of ovarian carcinoma cells in vitro is associated with a significant decrease of proto-oncogene MYC. Down-regulation of MYC is accomplished at least partially by accelerated MYC mRNA degradation and by inhibiting splicing [12]. Furthermore, butyrate is known to have an impact on AS of vascular endothelial growth factor (VEGF) gene [13]. Butyrate significantly increases the expression levels of both mRNA and protein of transcript variants that are anti-angiogenic in human lung endothelial cells while the pro-angiogenic isoform is not detectable [13]. However, the effect of butyrate on AS has not been systematically investigated in cattle. In this study, we presented evidence that butyrate can influence AS using high-throughput transcriptome technology and bioinformatics.

Results

1. Splice Junction Inference

The total number of reads mapped to the bovine genome (UMD3.1) using the Genomic Short-read Nucleotide Alignment Program (GSNAP) [14] were 52.25 and 47.73 million for the butyrate-treated (BT) and control (CT) groups, respectively (Table 1). For the sake of simplicity, only reads uniquely mapped to the genome (i.e., hitting only one locus) were counted. Among them, junction reads–those spanning exon-exon junctions with at least 8 bp of overhangs in each exon–accounted for 11.28 and 12.32% for the butyrate and control groups, respectively (P<0.001). These junction coverage results were in agreement with a published report [15], which shows 10–16% reads containing splice junctions in breast tumor samples. A total of 799,589 potential splicing junction sites were detected. Of them, 201,326 splice junction sites were supported by a minimum of 3 junction reads (mean; N=8). Approximately 94% of all splicing sites conformed to the consensus sequence, GT/AG while ~3% of junction sites belonged to the most common non-consensus splice junction GC/AG. The percentage of GC/AG junction observed in bovine epithelial cells was higher than that observed in humans and other species [16] where the GT-GC conversation represents only a small fraction at the donor sites (<1%). No AT/AC splice junction sites were observed in the bovine epithelial cell in this study (Table 1).

Table 1

Summary of junction reads.

2. Exon Skipping

A total of 2,834 potentially consistent exon skipping events were detected using GSNAP. A minimum of three junction reads (mean across all eight samples tested) supported each of these events. These junction reads were mapped to the same chromosome and the same strand in the correct order within a close proximity of 20,000 bp as defined (consistent), as compared to other common events, such as scrambles, where the two halves of junction reads were mapped to the same chromosome and on the same strand but in the wrong order. The 2,504 events were involved in 2,113 genes. Table 2 lists select exon skipping events with the difference in mean junction read counts greater than 20.0 between butyrate-treated and control groups. The gene expression level of all genes involved in these exon skipping events was significantly regulated by butyrate treatment at a stringent false discovery rate (FDR <0.0001) [11]. For example, collagen, type V, alpha 3 (COL5A3) was significantly up-regulated (~310 fold) by butyrate (P<0.0001 and FDR <0.0001). The mean number of junction reads mapped to exon #8 (donor) and exon#11 (acceptor) in butyrate-treated cells was 16.00 (±5.10, sd) while no junction reads spanning this region (0.00±0.00) were detected in control cells. This event had the potential to result in the formation of an mRNA isoform with 2 of the COL5A3 69 exons missing (exon#9 and exon#10). Quantitative RT-PCR confirmed this event (data not shown). Similarly, Exon#2 of S100A1 gene may be skipped, supported by multiple junction reads in BT, and the resultant new isoform was also differentially regulated by butyrate (Table 2).

Table 2

Exon skipping events supported by multiple junction reads detected using GSNAP.

3. Gene Fusion Events Significantly Affected by Butyrate

Only consistent splice junctions–both halves of junction reads mapped to the same chromosome strand in the correct order within a default distance, 20,000 bp–were considered in this study. Translocation (two halves mapped to different chromosomes), inversion (two halves mapped to the same chromosome but on opposite strands), scramble (on the same chromosome strand but in an incorrect order) as well as the genomic deletion [17] were excluded. The local gene fusion events or read-through in the bovine epithelial cells were readily detected by GSNAP. Of these, 13 events were significantly affected by butyrate (P<0.05; Table 3). The genes involved in these events were within the genomic interval set by the default parameter, -w =20kb. As the table shows, about half of the events are involved in the members of gene families. For example, the gene fusion event involving keratin 8 (donor gene) and keratin 4 (acceptor gene) was supported by 236 and 165 junction reads for BT and CT groups, respectively (P<0.0001). This event was verified independently using RT-PCR. Similarly, a gene fusion between homeobox A6 (HOXA6) and homeobox A5 (HOXA5) was supported by 3.47 and 7.19 junction reads for BT and CT, respectively (P<0.05). Interestingly, this gene fusion also started from an alternative 5′ splice junction (donor site), changing the 3′ boundary of the upstream exon.

Table 3

Gene fusion events significantly affected by butyrate in the bovine epithelial cells.

4. Differentially Expressed Isoforms Induced by Butyrate

Mixture-of-Isoforms algorithm (MISO), a probabilistic framework that enables quantification of the expression level of alternatively spliced genes and identification of differentially regulated isoforms from RNA-Seq data [18], was used to detect differentially expressed isoforms induced by butyrate in bovine epithelial cells. In this study, MISO detected a total of 216 isoforms from 98 genes that were differentially regulated by butyrate using a combined cutoff of ΔΨ >0.20 (i.e., the posterior distribution over the change in Psi for each event) and Bayes factor >10. The Bayes factor represents the weight of the evidence in the data in favor of differential expression. For example, Bayes factor =2 means that the isoform is two times more likely to be differentially expressed. Of these events, 50 isoforms from 36 genes were supported by at least 3 pair-wise comparisons (Table 4); 29 of the 36 genes had two mRNA isoforms while the rest had at least 3 isoforms. The mRNA-level expression of defensin, beta 1 gene (DEFB1) was not detected in control cells but significantly up-regulated by butyrate (FDR <0.001). This gene had two isoforms. Isoform #1 had 3 exons and was not detectable in control cells but barely detectable in butyrate-treated cells. However, Isoform#2 had 2 exons and was strongly induced by butyrate (Fig. 1). Similarly, IL-18 (ENSBTAG00000000277) had 2 isoforms with the same number of exons. These 2 isoforms differed in the exon structure, nevertheless. Isoform#1 had a similar number of assigned sequence reads while Isoform#2 had approximately 10 times more assigned reads in butyrate-treated cells (252.08±30.51, mean ± sd) than in control cells (25.94±13.69) (P<0.0001, Fig. 2). Isoform #1 of origin recognition complex, subunit 1 (ORC1, transcript# ENSBTAT00000052505) had mean read counts of 168.26±52.75 (± sd) and 1239.75±48.02 in butyrate and control groups, respectively, and was strongly repressed by butyrate. Isoform#2 of ORC1 remained unchanged. Histone deacetylase 10 (HDAC10) had 2 isoforms with the read count ratio between dominant (ENSBTAT00000014602) and minor (ENSBTAT00000042583) isoforms greater than 150 in the control cells. Similarly, the dominant isoform was significantly repressed (P<0.00001), while the minor isoform was not affected by butyrate. However, the expression of both isoforms of HDAC5 was strongly enhanced by butyrate treatment (data not shown). For genes with ≥3 isoforms, one or 2 isoforms were often differentially regulated by butyrate. As Figure 3 shows, the expression level of Isoforms# 1 and 3 of aminoacylase 1gene (ACY1) was unchanged but that of Isoform#2 was strongly repressed by butyrate (Fig. 3). On the other hand, Isoform#2 of coiled-coil domain containing 24 (CCDC24) was little changed in response to butyrate treatment while Isoforms#1 and 3 were significantly induced by butyrate (Fig. 4). Moreover, differential expression of 6 selected transcript isoforms induced by butyrate, including transcripts such as ENSBTAT00000009005 (LGALS9), ENSBTAT00000010687 (CLSTN3), ENSBTAT00000012035 (MX1), ENSBTAT00000020401 (NASP), ENSBTAT00000036545 (GUK1), and ENSBTAT00000061158 (TACC2), were confirmed using RT-PCR.

Table 4

Differentially expressed isoforms detected using MISO in the bovine epithelial cells.

Figure 1

Relative abundance of transcript isoforms of β-defensin 1 gene induced by butyrate in bovine epithelial cells.

Figure 2

Relative abundance of transcript isoforms of IL-18 gene induced by butyrate in bovine epithelial cells.

Figure 3

Relative abundance of transcript isoforms of aminoacylase 1 gene.

Figure 4

Relative abundance of transcript isoforms of coiled-coil domain containing 24 (CCDC24) gene in bovine epithelial cells.

Discussion

Butyrate is one of important short-chain fatty acids (SCFAs). Produced in the rumen and hindgut by gut microorganisms, butyrate is rapidly absorbed and utilized by the rumen and colon epithelium, contributing to 75% of the total energy requirement in ruminants and ~10% for humans [9], [19]. Additionally, butyrate is shown to reinforce intestinal barriers and modulate motility and visceral sensitivity of the intestine [20]. As a signaling molecule, butyrate induces apoptosis and inhibits cell proliferation, and therefore, it has antitumorigenic properties. Butyrate is an inhibitor of histone deacetylases (HDACs), one of three classes of enzymes epigenetically modifying chromatin histones. Mounting evidence suggests histone acetylation plays a major role in controlling transcriptional activities of genes [21], [22]. HDAC inhibitors such as butyrate induce hyperacetylation of histones, and therefore increase transcriptional activities [2]. This was supported by our observation that the total number of genes transcribed in butyrate-treated cells (mean ± sd =19,322±155) was significantly higher than in untreated control cells (17,626±125) (P<0.00001). Transcriptional effects of butyrate have been extensively investigated in various cell types, in both in vitro and in vivo models [19], [23]–[25]. Indeed, previous studies have demonstrated that butyrate down-regulates genes controlling cell proliferation in vitro, resulting in the inhibition of the proliferation of epithelial cells [25], [26]. Microarray results show that butyrate has a profound effect in global gene expression, including up-regulation of genes related to apoptosis and differentiation [10], [26]. Our recent transcriptome studies using RNA-seq technology demonstrate that a 24-h butyrate treatment significantly affects the transcript abundance of 11,408 of the 17,625 genes detected in the bovine epithelial cell, representing ~65% of the entire transcriptome [11], [19]. However, these studies focus on the transcription at a gene level. The effect of butyrate on individual transcript isoforms and alternative splicing has been systematically studied only recently in human cells [27]. In this study, we examined the regulation of alternative splicing by butyrate in bovine epithelial cells. Our results should facilitate a better understanding of alternative splicing in the development of epithelial cells-derived diseases.

Of four classes of histone deacetylases, butyrate inhibits enzymatic activities of most HDACs in Class I, II, and IV, which are zinc-dependent, except HDAC6 and HDAC10 [28]. Class III HDACs (also called sirtuins or SIRTs) depend on nicotinamide adenine dinucleotide for their catalytic activity [29]. SIRTs are associated with chromtain regulation and affect genome stability in yeast and may represent pivotal regulators of lifespan and aging [30]. SIRTs catalyze two major biochemical reactions: deacetylation on lysine residues of target proteins by altering cellular [NAD⁺]/[NADH] ratios (SIRT1, SIRT2, SIRT3, SIRT5, and SIRT7) and ADP-ribosylation (SIRT4 and SIRT6) [31]. In neuronal cells, SIRT1, SIRT5, and SIRT6 are down-regulated, whereas SIRT2, SIRT4, and SIRT7 up-regulated by butyrate [32]. Our RNA-seq data suggest that butyrate regulated the transcript abundance or gene expression of the majority of HDACs (Table 5). Butyrate significantly increased the expression of HDAC3, HDAC5, and HDAC11. On the other hand, the expression level of HDACs7-10 was significantly down-regulated. The mRNA levels of SIRT4 and SIRT6 were strongly up-regulated while SIRT1 was significantly down-regulated by butyrate (Table 5). However, the relative abundance of HDAC1, HDAC2, HDAC4, and HDAC6 remained unchanged by butyrate. In addition to its effect on the expression at the whole gene level, butyrate selectively regulated the transcript abundance of different mRNA isoforms. While the abundance of both short and long isoforms of HDAC5 was significantly enhanced by butyrate, only the long and dominant isoform of HDAC10 was significantly down-regulated. Intriguingly, butyrate is unable to inhibit the catalytic activity of HDAC10 [28]. Nevertheless, butyrate may still exert its control on the deacetylase activity of HDAC10 via transcriptional regulation at the mRNA level. Our future work will focus on the biological relevance of various HDAC isoforms induced by butyrate, especially various SIRTs and their roles in cell senescence and aging.

Table 5

The expression of histone deacetylases affected by butyrate in the bovine epithelial cells.

Distant gene fusion events are well known in tumors, often resulting from genomic abnormalities such as chromosomal translocation. These events, such as BCL-ABL, lead to the formation of a novel chimeric protein with different functions and are one of the common mechanisms for oncogene activation [33]. Recently, a new type of fusion involving two adjacent genes in the same orientation on the same chromosome has been described [34], [35]. Adjacent genes are normally transcripted independently. However, a single transcript can be occasionally formed to include at least part of one exon from each of two or more distinct genes [36]. This phenomenon, Transcription Induced Chimeras (TICs) [35] or Conjoined Genes (CGs) [36], is widespread in mammalian genomes. It is estimated that at least 4%-5% of the tandem gene pairs in the human genome can be transcripted into TICs. Moreover, these TICs may possess novel functions because >70% of them are conserved in other vertebrate genomes [36]. In this study, we detected 13 TICs that were supported by multiple junction reads. Intriguingly, these TICs were also differentially regulated by butyrate in bovine epithelial cells. Approximately 46% of these fusion events were involved in the members of gene families, which is much higher than 11% as previously reported [36]. For example, TICs were formed between 2 homeobox genes, HOXA6 and HOXA5 and between 2 keratin genes, KRT8 and KRT4. In addition, the fusion between zinc finger proteins ZNF865 and ZNF524 was also supported by multiple junction reads, and a significantly higher number of reads was detected in untreated control cells than in butyrate-treated cells. A similar fusion event between ZNF649 and ZNF577 was identified in prostate tumors [37]. A relatively higher percentage of TICs between genes with related functions identified in this study should be further examined. Most importantly, the functional significance of these fusion events, especially their possible role in transcription regulation, should be experimentally determined.

One advantage of GSNAP algorithm is its potential to identify novel splice sites and therefore possibly novel transcript isoforms. The algorithm relies on a maximum entropy model and uses frequencies of nucleotides neighboring a donor and acceptor splice site to discriminate between true and false splicing sites [14]. The power of this approach was exemplified by a case study involving the prohibitin gene (PHB). The prohibitin protein complex, located in mitochondrial inner membrane, is formed by heteromeric binding of both PHB and PHB2 [38] and is involved in transcription regulation and cell cycle progression by blocking the G1/S transition of the cell cycle [39]. Prohibitin induces apoptosis by interacting with the retinoblastoma protein as well as being involved in the p53 pathway. Its 3′ untranslated region (UTR) acts as a novel class of non-coding regulatory RNAs. Additionally, PHB expression is up-regulated in the retina in aging and diabetic models and may serve as an oxidative marker [38]. A recent study using thyroid tumor cell lines demonstrates that butyrate increases PHB mRNA expression. Furthermore, butyrate as well as other HDAC inhibitors, such as trichostatin A, affects PHB splicing [2], leading to the over-expression of the longer isoform with 3′ UTR. Both inhibitors decrease the mRNA levels of the shorter isoform but increase those of the longer isoform, which exerts a growth-suppressive action. Our results showed that butyrate significantly down-regulated the mRNA expression of both PHB and PHB2 in the bovine epithelial cell (FDR <0.0001). No known isoforms in both genes have been annotated in cattle so far. Annotated PHB and PHB2 genes have 7 and 10 exons, respectively. GSNAP correctly identified all normally splicing exon-exon junctions. Moreover, GSNAP detected novel splice sites. For example, several junction reads detected in the untreated control cells suggest a possible exon skipping event that may result in skipping of Exon#2 in PHB2. Such reads were not detectable in the butyrate-treated cells. In the PHB gene, significantly higher numbers of junction reads in the control group than the butyrate-treated group indicated that multiple alternative splicing events involved Exon#1 and Exon#2 and the intron between them. These events occurred in the 5′ UTR and did not seem to alter its primary protein structure. In humans, the 3′ UTR of PHB is attributed to its anti-tumorigenic and anti-proliferative properties [40]. The biological implication of various splicing events in the 5′ UTR of PHB genes in cattle is worthy of further investigation.

Materials and Methods

1. Treatments

The bovine epithelial cells used in this study were previously described [11]. The culture and sodium butyrate treatment was essentially the same as reported [26]. Briefly, cells were treated with 10 mM sodium butyrate for 24h. As a result of butyrate treatment, cell cycle arrest was notable. The percentage of G1/G0 cells was increased from ~41% in normal cell populations to 79% in butyrate-treated cell populations, in a good agreement with previous results [10]. Cells were then harvested and high-quality total RNA (RNA Integrity number or RIN >9.0) was processed using an Illumina TruSeq RNA sample prep kit following the manufacturer’s instruction (Illumina, San Diego, CA, USA). After quality control procedures, individual RNA-seq libraries were then pooled based on their respective sample-specific 6-bp adaptors and sequenced at 50bp/sequence read using an Illumina HiSeq 2000 sequencer as described previously [41]. The mean number of sequence reads generated per sample was 67,527,111±8,330,388.58 (± sd). A total of eight replicates in two groups, butyrate treated (BT) and untreated (control or CT; N =4 for each group), were used. Raw sequence reads were deposited to the NCBI Sequence Read Archive (accession# SRA051007.1).

2. Data Analysis and Bioinformatics

Raw sequence reads were first checked using our quality control pipeline. Nucleotides of each raw read were scanned for low quality and trimmed using SolexaQA [42]. Input reads were then aligned to the bovine reference genome using GSNAP [14] with default parameters. Two mismatches were tolerated (“-m 2″ flag). The intron size was specified by “-w” =20,000 bp. Read counts, including junction reads, were extracted from GSNAP –N flag output files and further analyzed.

For differentially expressed isoform detection, Mixture-of-Isoforms program (MISO v 0.4.1, released February 1, 2012) [18] was used. Trimmed reads were first aligned to the bovine genome (UMD3.1) with TopHat [43] using default parameters. The TopHat SAM output files were converted into BAM files and input into MISO with the GTF file from ENSEMBL bovine genebuild release 65. Differentially expressed isoforms were detected using the following filtering parameters: the sum of inclusion and exclusion reads is greater than 10 (≥1 inclusion read and ≥1 exclusion read), ΔΨ greater than 0.20, and the Bayes factor ≥10. In addition, each SAM output file from TopHat alignment, along with the GTF file from ENSEMBL bovine genebuild release v65, were used in the Cuffdiff program [44] in the Cufflink package (v1.3.0) to test for differential expression at the gene level, as described previously [19].

3. Quantitative RT-PCR

Ten selected alternative splicing events, including 2 exon-skipping event, 2 gene fusion events, and 6 transcript isoforms, were tested using real-time RT-PCR as described previously [45]. Briefly, the cDNA synthesis was performed using an iScript cDNA Synthesis kit (Bio-Rad, Hercules, CA). Real-time RT-PCR analysis was carried out with an iQ SYBR Green Supermix kit (Biorad) using 200 nM of each primer and the 1^st-strand cDNA (100 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. A melting curve analysis was performed for each primer pair. PCR products were further analyzed on a high-sensitive DNA chip using an Agilent Bioanalyzer 2100 (Agilent, Palo Alto, CA, USA) for product length.

Acknowledgments

We’d like to thank Alicia Beavers for their 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. The USDA is an equal opportunity provider and employer.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: This study was partly supported by Award R01HG005978 to WL from the National Human Genome Research Institute (NHGRI). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NHGRI or the National Institutes of Health. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. No additional external funding was received for this study

References

1. Tazi J, Bakkour N, Stamm S. Alternative splicing and disease. Biochimica et biophysica acta. 2009;1792:14–26. [PubMed]

2. Puppin C, Passon N, Franzoni A, Russo D, Damante G. Histone deacetylase inhibitors control the transcription and alternative splicing of prohibitin in thyroid tumor cells. Oncology reports. 2011;25:393–397. [PubMed]

3. Tazi J, Bakkour N, Marchand V, Ayadi L, Aboufirassi A, et al. Alternative splicing: regulation of HIV-1 multiplication as a target for therapeutic action. The FEBS journal. 2010;277:867–876. [PubMed]

4. Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, et al. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008;456:470–476. [PMC free article] [PubMed]

5. Murray JI, Voelker RB, Henscheid KL, Warf MB, Berglund JA. Identification of motifs that function in the splicing of non-canonical introns. Genome biology. 2008;9:R97. [PMC free article] [PubMed]

6. Will CL, Luhrmann R. Spliceosome structure and function. Cold Spring Harbor perspectives in biology 3. 2011. [PubMed]

7. Boon KL, Grainger RJ, Ehsani P, Barrass JD, Auchynnikava T, et al. prp8 mutations that cause human retinitis pigmentosa lead to a U5 snRNP maturation defect in yeast. Nature structural & molecular biology. 2007;14:1077–1083. [PMC free article] [PubMed]

8. Nissim-Rafinia M, Aviram M, Randell SH, Shushi L, Ozeri E, et al. Restoration of the cystic fibrosis transmembrane conductance regulator function by splicing modulation. EMBO reports. 2004;5:1071–1077. [PubMed]

9. Bergman EN. Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiological reviews. 1990;70:567–590. [PubMed]

10. Li RW, Li C. Butyrate induces profound changes in gene expression related to multiple signal pathways in bovine kidney epithelial cells. BMC genomics. 2006;7:234. [PMC free article] [PubMed]

11. Wu S, Li RW, Li W, Li CJ. Transcriptome characterization by RNA-seq unravels the mechanisms of butyrate-induced epigenomic regulation in bovine cells. submitted. 2012. [PMC free article] [PubMed]

12. Krupitza G, Harant H, Dittrich E, Szekeres T, Huber H, et al. Sodium butyrate inhibits c-myc splicing and interferes with signal transduction in ovarian carcinoma cells. Carcinogenesis. 1995;16:1199–1205. [PubMed]

13. Ciura J, Jagodzinski PP. Butyrate increases the formation of anti-angiogenic vascular endothelial growth factor variants in human lung microvascular endothelial cells. Molecular biology reports. 2010;37:3729–3734. [PubMed]

14. Wu TD, Nacu S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics. 2010;26:873–881. [PMC free article] [PubMed]

15. Wang K, Singh D, Zeng Z, Coleman SJ, Huang Y, et al. MapSplice: accurate mapping of RNA-seq reads for splice junction discovery. Nucleic acids research. 2010;38:e178. [PMC free article] [PubMed]

16. Churbanov A, Winters-Hilt S, Koonin EV, Rogozin IB. Accumulation of GC donor splice signals in mammals. Biology direct. 2008;3:30. [PMC free article] [PubMed]

17. Nacu S, Yuan W, Kan Z, Bhatt D, Rivers CS, et al. Deep RNA sequencing analysis of readthrough gene fusions in human prostate adenocarcinoma and reference samples. BMC medical genomics. 2011;4:11. [PMC free article] [PubMed]

18. Katz Y, Wang ET, Airoldi EM, Burge CB. Analysis and design of RNA sequencing experiments for identifying isoform regulation. Nature methods. 2010;7:1009–1015. [PMC free article] [PubMed]

19. Baldwin RLV, Wu S, Li W, Li CJ, Bequette BJ, et al. Quantification of transcriptome responses of the rumen epithelium to butyrate infusion. submitted. 2012. [PMC free article] [PubMed]

20. Canani RB, Costanzo MD, Leone L, Pedata M, Meli R, et al. Potential beneficial effects of butyrate in intestinal and extraintestinal diseases. World journal of gastroenterology : WJG. 2011;17:1519–1528. [PMC free article] [PubMed]

21. Thiagalingam S, Cheng KH, Lee HJ, Mineva N, Thiagalingam A, et al. Histone deacetylases: unique players in shaping the epigenetic histone code. Annals of the New York Academy of Sciences. 2003;983:84–100. [PubMed]

22. Shin JH, Li RW, Gao Y, Baldwin RT, Li CJ. Genome-wide ChIP-seq mapping and analysis reveal butyrate-induced acetylation of H3K9 and H3K27 correlated with transcription activity in bovine cells. Functional & integrative genomics. 2012;12:119–130. [PubMed]

23. Daly K, Shirazi-Beechey SP. Microarray analysis of butyrate regulated genes in colonic epithelial cells. DNA and cell biology. 2006;25:49–62. [PubMed]

24. Li CJ, Li RW. Butyrate induced cell cycle arrest in bovine cells through targeting gene expression relevant to DNA replication apparatus. Gene regulation and systems biology. 2008;2:113–123. [PMC free article] [PubMed]

25. Ranganna K, Yousefipour Z, Yatsu FM, Milton SG, Hayes BE. Gene expression profile of butyrate-inhibited vascular smooth muscle cell proliferation. Molecular and cellular biochemistry. 2003;254:21–36. [PubMed]

26. Li CJ, Li RW, Wang YH, Elsasser TH. Pathway analysis identifies perturbation of genetic networks induced by butyrate in a bovine kidney epithelial cell line. Functional & integrative genomics. 2007;7:193–205. [PubMed]

27. Hnilicova J, Hozeifi S, Duskova E, Icha J, Tomankova T, et al. Histone deacetylase activity modulates alternative splicing. PloS one. 2011;6:e16727. [PMC free article] [PubMed]

28. Davie JR. Inhibition of histone deacetylase activity by butyrate. The Journal of nutrition. 2003;133:2485S–2493S. [PubMed]

29. Schemies J, Uciechowska U, Sippl W, Jung M. NAD(+) -dependent histone deacetylases (sirtuins) as novel therapeutic targets. Medicinal research reviews. 2010;30:861–889. [PubMed]

30. Corbi G, Conti V, Scapagnini G, Filippelli A, Ferrara N. Role of sirtuins, calorie restriction and physical activity in aging. Frontiers in bioscience. 2012;4:768–778. [PubMed]

31. Carafa V, Nebbioso A, Altucci L. Sirtuins and disease: the road ahead. Frontiers in pharmacology. 2012;3:4. [PMC free article] [PubMed]

32. Kyrylenko S, Kyrylenko O, Suuronen T, Salminen A. Differential regulation of the Sir2 histone deacetylase gene family by inhibitors of class I and II histone deacetylases. Cellular and molecular life sciences : CMLS. 2003;60:1990–1997. [PubMed]

33. Shtivelman E, Lifshitz B, Gale RP, Canaani E. Fused transcript of abl and bcr genes in chronic myelogenous leukaemia. Nature. 1985;315:550–554. [PubMed]

34. Akiva P, Toporik A, Edelheit S, Peretz Y, Diber A, et al. Transcription-mediated gene fusion in the human genome. Genome research. 2006;16:30–36. [PubMed]

35. Parra G, Reymond A, Dabbouseh N, Dermitzakis ET, Castelo R, et al. Tandem chimerism as a means to increase protein complexity in the human genome. Genome research. 2006;16:37–44. [PubMed]

36. Prakash T, Sharma VK, Adati N, Ozawa R, Kumar N, et al. Expression of conjoined genes: another mechanism for gene regulation in eukaryotes. PloS one. 2010;5:e13284. [PMC free article] [PubMed]

37. Kinsella M, Harismendy O, Nakano M, Frazer KA, Bafna V. Sensitive gene fusion detection using ambiguously mapping RNA-Seq read pairs. Bioinformatics. 2011;27:1068–1075. [PMC free article] [PubMed]

38. Lee H, Arnouk H, Sripathi S, Chen P, Zhang R, et al. Prohibitin as an oxidative stress biomarker in the eye. International journal of biological macromolecules. 2010;47:685–690. [PubMed]

39. Gamble SC, Odontiadis M, Waxman J, Westbrook JA, Dunn MJ, et al. Androgens target prohibitin to regulate proliferation of prostate cancer cells. Oncogene. 2004;23:2996–3004. [PubMed]

40. Jupe ER, Liu XT, Kiehlbauch JL, McClung JK, Dell’Orco RT. The 3′ untranslated region of prohibitin and cellular immortalization. Experimental cell research. 1996;224:128–135. [PubMed]

41. Li RW, Schroeder SG. Cytoskeleton remodeling and alterations in smooth muscle contractility in the bovine jejunum during nematode infection. Functional & integrative genomics. 2011. [PubMed]

42. Cox MP, Peterson DA, Biggs PJ. SolexaQA: At-a-glance quality assessment of Illumina second-generation sequencing data. BMC bioinformatics. 2010;11:485. [PMC free article] [PubMed]

43. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–1111. [PMC free article] [PubMed]

44. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature biotechnology. 2010;28:511–515. [PMC free article] [PubMed]

45. Li RW, Rinaldi M, Capuco AV. Characterization of the abomasal transcriptome for mechanisms of resistance to gastrointestinal nematodes in cattle. Veterinary research. 2011;42:114. [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of
Public Library of Science