Read-Split-Run: an improved bioinformatics pipeline for identification of genome-wide non-canonical spliced regions using RNA-Seq data.

BACKGROUND: Most existing tools for detecting next-generation sequencing-based splicing events focus on generic splicing events. Consequently, special types of non-canonical splicing events of short mRNA regions (IRE1α targeted) have not yet been thoroughly addressed at a genome-wide level using bioinformatics approaches in conjunction with next-generation technologies. During endoplasmic reticulum (ER) stress, the gene encoding the RNase Ire1α is known to splice out a short 26 nt region from the mRNA of the transcription factor Xbp1 non-canonically within the cytosol. This causes an open reading frame-shift that induces expression of many downstream genes in reaction to ER stress as part of the unfolded protein response (UPR). We previously published an algorithm termed "Read-Split-Walk" (RSW) to identify non-canonical splicing regions using RNA-Seq data and applied it to ER stress-induced Ire1α heterozygote and knockout mouse embryonic fibroblast cell lines. In this study, we have developed an improved algorithm "Read-Split-Run" (RSR) for detecting genome-wide Ire1α-targeted genes with non-canonical spliced regions at a faster speed. We applied the RSR algorithm using different combinations of several parameters to the previously RSW tested mouse embryonic fibroblast cells (MEF) and the human Encyclopedia of DNA Elements (ENCODE) RNA-Seq data. We also compared the performance of RSR with two other alternative splicing events identification tools (TopHat (Trapnell et al., Bioinformatics 25:1105-1111, 2009) and Alt Event Finder (Zhou et al., BMC Genomics 13:S10, 2012)) utilizing the context of the spliced Xbp1 mRNA as a positive control in the data sets we identified it to be the top cleavage target present in Ire1α (+/-) but absent in Ire1α (-/-) MEF samples and this comparison was also extended to human ENCODE RNA-Seq data. RESULTS: Proof of principle came in our results by the fact that the 26 nt non-conventional splice site in Xbp1 was detected as the top hit by our new RSR algorithm in heterozygote (Het) samples from both Thapsigargin (Tg) and Dithiothreitol (Dtt) treated experiments but absent in the negative control Ire1α knock-out (KO) samples. Applying different combinations of parameters to the mouse MEF RNA-Seq data, we suggest a General Linear Model (GLM) for both Tg and Dtt treated experiments. We also ran RSR for a human ENCODE RNA-Seq dataset and identified 32,597 spliced regions for regular chromosomes. TopHat (Trapnell et al., Bioinformatics 25:1105-1111, 2009) and Alt Event Finder (Zhou et al., BMC Genomics 13:S10, 2012) identified 237,155 spliced junctions and 9,129 exon skipping events (excluding chr14), respectively. Our Read-Split-Run algorithm also outperformed others in the context of ranking Xbp1 gene as the top cleavage target present in Ire1α (+/-) but absent in Ire1α (-/-) MEF samples. The RSR package including source codes is available at http://bioinf1.indstate.edu/RSR and its pipeline source codes are also freely available at https://github.com/xuric/read-split-run for academic use. CONCLUSIONS: Our new RSR algorithm has the capability of processing massive amounts of human ENCODE RNA-Seq data for identifying novel splice junction sites at a genome-wide level in a much more efficient manner when compared to the previous RSW algorithm. Our proposed model can also predict the number of spliced regions under any combinations of parameters. Our pipeline can detect novel spliced sites for other species using RNA-Seq data generated under similar conditions.

The IRE1α/XBP1s Pathway Is Essential for the Glucose Response and Protection of ß Cells.

Although glucose uniquely stimulates proinsulin biosynthesis in ß cells, surprisingly little is known of the underlying mechanism(s). Here, we demonstrate that glucose activates the unfolded protein response transducer inositol-requiring enzyme 1 alpha (IRE1α) to initiate X-box-binding protein 1 (Xbp1) mRNA splicing in adult primary ß cells. Using mRNA sequencing (mRNA-Seq), we show that unconventional Xbp1 mRNA splicing is required to increase and decrease the expression of several hundred mRNAs encoding functions that expand the protein secretory capacity for increased insulin production and protect from oxidative damage, respectively. At 2 wk after tamoxifen-mediated Ire1α deletion, mice develop hyperglycemia and hypoinsulinemia, due to defective ß cell function that was exacerbated upon feeding and glucose stimulation. Although previous reports suggest IRE1α degrades insulin mRNAs, Ire1α deletion did not alter insulin mRNA expression either in the presence or absence of glucose stimulation. Instead, ß cell failure upon Ire1α deletion was primarily due to reduced proinsulin mRNA translation primarily because of defective glucose-stimulated induction of a dozen genes required for the signal recognition particle (SRP), SRP receptors, the translocon, the signal peptidase complex, and over 100 other genes with many other intracellular functions. In contrast, Ire1α deletion in ß cells increased the expression of over 300 mRNAs encoding functions that cause inflammation and oxidative stress, yet only a few of these accumulated during high glucose. Antioxidant treatment significantly reduced glucose intolerance and markers of inflammation and oxidative stress in mice with ß cell-specific Ire1α deletion. The results demonstrate that glucose activates IRE1α-mediated Xbp1 splicing to expand the secretory capacity of the ß cell for increased proinsulin synthesis and to limit oxidative stress that leads to ß cell failure.

The unfolded protein response transducer IRE1α prevents ER stress-induced hepatic steatosis.

The endoplasmic reticulum (ER) is the cellular organelle responsible for protein folding and assembly, lipid and sterol biosynthesis, and calcium storage. The unfolded protein response (UPR) is an adaptive intracellular stress response to accumulation of unfolded or misfolded proteins in the ER. In this study, we show that the most conserved UPR sensor inositol-requiring enzyme 1 α (IRE1α), an ER transmembrane protein kinase/endoribonuclease, is required to maintain hepatic lipid homeostasis under ER stress conditions through repressing hepatic lipid accumulation and maintaining lipoprotein secretion. To elucidate physiological roles of IRE1α-mediated signalling in the liver, we generated hepatocyte-specific Ire1α-null mice by utilizing an albumin promoter-controlled Cre recombinase-mediated deletion. Deletion of Ire1α caused defective induction of genes encoding functions in ER-to-Golgi protein transport, oxidative protein folding, and ER-associated degradation (ERAD) of misfolded proteins, and led to selective induction of pro-apoptotic UPR trans-activators. We show that IRE1α is required to maintain the secretion efficiency of selective proteins. In the absence of ER stress, mice with hepatocyte-specific Ire1α deletion displayed modest hepatosteatosis that became profound after induction of ER stress. Further investigation revealed that IRE1α represses expression of key metabolic transcriptional regulators, including CCAAT/enhancer-binding protein (C/EBP) ß, C/EBPδ, peroxisome proliferator-activated receptor γ (PPARγ), and enzymes involved in triglyceride biosynthesis. IRE1α was also found to be required for efficient secretion of apolipoproteins upon disruption of ER homeostasis. Consistent with a role for IRE1α in preventing intracellular lipid accumulation, mice with hepatocyte-specific deletion of Ire1α developed severe hepatic steatosis after treatment with an ER stress-inducing anti-cancer drug Bortezomib, upon expression of a misfolding-prone human blood clotting factor VIII, or after partial hepatectomy. The identification of IRE1α as a key regulator to prevent hepatic steatosis provides novel insights into ER stress mechanisms in fatty liver diseases associated with toxic liver injuries.

UPR pathways combine to prevent hepatic steatosis caused by ER stress-mediated suppression of transcriptional master regulators.

The unfolded protein response (UPR) is linked to metabolic dysfunction, yet it is not known how endoplasmic reticulum (ER) disruption might influence metabolic pathways. Using a multilayered genetic approach, we find that mice with genetic ablations of either ER stress-sensing pathways (ATF6alpha, eIF2alpha, IRE1alpha) or of ER quality control (p58(IPK)) share a common dysregulated response to ER stress that includes the development of hepatic microvesicular steatosis. Rescue of ER protein processing capacity by the combined action of UPR pathways during stress prevents the suppression of a subset of metabolic transcription factors that regulate lipid homeostasis. This suppression occurs in part by unresolved ER stress perpetuating expression of the transcriptional repressor CHOP. As a consequence, metabolic gene expression networks are directly responsive to ER homeostasis. These results reveal an unanticipated direct link between ER homeostasis and the transcriptional regulation of metabolism, and suggest mechanisms by which ER stress might underlie fatty liver disease.