Novel bioinformatics method for identification of genome-wide non-canonical spliced regions using RNA-Seq data

Maresin 1 attenuates NAFLD by suppression of endoplasmic reticulum stress via AMPK-SERCA2b pathway

Maresin 1 (MAR1), which is derived from docosahexaenoic acid biosynthesized by macrophages, has been reported to improve insulin resistance. Recently, it has been documented that MAR1 could ameliorate inflammation and insulin resistance in obese mice. These findings led us to investigate the effects of MAR1 on hepatic lipid metabolism. We found that MAR1 could stimulate AMP-activated protein kinase (AMPK), thereby augmenting sarcoendoplasmic reticulum Ca²⁺-ATPase 2b (SERCA2b) expression. This stimulation suppressed lipid accumulation by attenuating the endoplasmic reticulum (ER) stress in hepatocytes under hyperlipidemic conditions. Attenuation was mitigated by knockdown of AMPK or thapsigargin, a SERCA2b inhibitor. We also demonstrated that MAR1 administration resulted in increased hepatic AMPK phosphorylation and Serca2b mRNA expression, whereas hepatic ER stress was reduced in high-fat diet (HFD)-fed mice. Moreover, MAR1 treatment suppressed hepatic lipid synthesis, thereby attenuating hepatic steatosis in HFD-fed mice. In conclusion, our results suggest that MAR1 ameliorates hepatic steatosis via AMPK/SERCA2b-mediated suppression of ER stress. Therefore, MAR1 may be an effective therapeutic strategy for treating non-alcoholic fatty liver disease (NAFLD) via regulation of ER stress-induced hepatic lipogenesis.

VIH from the mud crab is specifically expressed in the eyestalk and potentially regulated by transactivator of Sox9/Oct4/Oct1

Vitellogenesis-inhibiting hormone (VIH) is known to regulate ovarian maturation by suppressing the synthesis of vitellogenin (Vtg) in crustaceans, which belongs to a member of crustacean hyperglycemic hormone (CHH) family synthesized and secreted from the X-organ/sinus gland complex of eyestalks. In this study, the cDNA, genomic DNA (gDNA) and the 5'-upstream regulatory (promoter region) sequences of VIH gene were obtained by conventional PCR, genome walker and tail-PCR techniques according to our transcriptomic database of Scylla paramamosain. The full-length cDNA of SpVIH is 634bp including 105bp 5'UTR, 151bp 3'UTR and 378bp ORF that encodes a peptide of 125 amino acids. The full length gDNA of SpVIH is 790bp containing two exons and one intron. The 5'-flanking promoter regions of SpVIH we isolated are 3070bp from the translation initiation (ATG) and 2398bp from the predicted transcription initiation (A), which consists of putative core promoter region and multiple potential transcription factor binding sites. SpVIH was only expressed in eyestalk. The expression level of SpVIH in eyestalk of female crab decreased gradually along with the development of ovary. As there is not cell line of crabs available, we chose the mature transfection system HEK293FT cell lines to explore the mechanism of transcription regulation of SpVIH in crabs. Sequential deletion assays using luciferase reporter gene in HEK293FT cells revealed that the possible promoter activity regions (including positive and negative transcription factors binding sites simultaneously) presented between pSpVIH-4 and pSpVIH-6. In order to further identify the crucial transcription factors binding site in this region, the site-directed mutagenesis of Sox9/Oct4/Oct1 binding site of pSpVIH-4 was created. The results demonstrated that the transcriptional activity of pSpVIH-4△ decreased significantly (p<0.05). Thus, it is reasonable to deduce that the Sox9/Oct4/Oct1 may be the essential positive transcription factors which regulate the expression of SpVIH.

Identification of genome-wide non-canonical spliced regions and analysis of biological functions for spliced sequences using Read-Split-Fly

Background

It is generally thought that most canonical or non-canonical splicing events involving U2- and U12 spliceosomes occur within nuclear pre-mRNAs. However, the question of whether at least some U12-type splicing occurs in the cytoplasm is still unclear. In recent years next-generation sequencing technologies have revolutionized the field. The “Read-Split-Walk” (RSW) and “Read-Split-Run” (RSR) methods were developed to identify genome-wide non-canonical spliced regions including special events occurring in cytoplasm. As the significant amount of genome/transcriptome data such as, Encyclopedia of DNA Elements (ENCODE) project, have been generated, we have advanced a newer more memory-efficient version of the algorithm, “Read-Split-Fly” (RSF), which can detect non-canonical spliced regions with higher sensitivity and improved speed. The RSF algorithm also outputs the spliced sequences for further downstream biological function analysis.

Results

We used open access ENCODE project RNA-Seq data to search spliced intron sequences against the U12-type spliced intron sequence database to examine whether some events could occur as potential signatures of U12-type splicing. The check was performed by searching spliced sequences against 5’ss and 3’ss sequences from the well-known orthologous U12-type spliceosomal intron database U12DB. Preliminary results of searching 70 ENCODE samples indicated that the presence of 5’ss with U12-type signature is more frequent than U2-type and prevalent in non-canonical junctions reported by RSF. The selected spliced sequences have also been further studied using miRBase to elucidate their functionality. Preliminary results from 70 samples of ENCODE datasets show that several miRNAs are prevalent in studied ENCODE samples. Two of these are associated with many diseases as suggested in the literature. Specifically, hsa-miR-1273 and hsa-miR-548 are associated with many diseases and cancers.

Conclusions

Our RSF pipeline is able to detect many possible junctions (especially those with a high RPKM) with very high overall accuracy and relative high accuracy for novel junctions. We have incorporated useful parameter features into the pipeline such as, handling variable-length read data, and searching spliced sequences for splicing signatures and miRNA events. We suggest RSF, a tool for identifying novel splicing events, is applicable to study a range of diseases across biological systems under different experimental conditions.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1801-y) contains supplementary material, which is available to authorized users.

Physiological/pathological ramifications of transcription factors in the unfolded protein response

Numerous environmental, physiological, and pathological insults disrupt protein-folding homeostasis in the endoplasmic reticulum (ER), referred to as ER stress. Eukaryotic cells evolved a set of intracellular signaling pathways, collectively termed the unfolded protein response (UPR), to maintain a productive ER protein-folding environment through reprogramming gene transcription and mRNA translation. The UPR is largely dependent on transcription factors (TFs) that modulate expression of genes involved in many physiological and pathological conditions, including development, metabolism, inflammation, neurodegenerative diseases, and cancer. Here we summarize the current knowledge about these mechanisms, their impact on physiological/pathological processes, and potential therapeutic applications.

Identification of human short introns

Canonical pre-mRNA splicing requires snRNPs and associated splicing factors to excise conserved intronic sequences, with a minimum intron length required for efficient splicing. Non-canonical splicing-intron excision without the spliceosome-has been documented; most notably, some tRNAs and the XBP1 mRNA contain short introns that are not removed by the spliceosome. There have been some efforts to identify additional short introns, but little is known about how many short introns are processed from mRNAs. Here, we report an approach to identify RNA short introns from RNA-Seq data, discriminating against small genomic deletions. We identify hundreds of short introns conserved among multiple human cell lines. These short introns are often alternatively spliced and are found in a variety of RNAs-both mRNAs and lncRNAs. Short intron splicing efficiency is increased by secondary structure, and we detect both canonical and non-canonical short introns. In many cases, splicing of these short introns from mRNAs is predicted to alter the reading frame and change protein output. Our findings imply that standard gene prediction models which often assume a lower limit for intron size fail to predict short introns effectively. We conclude that short introns are abundant in the human transcriptome, and short intron splicing represents an added layer to mRNA regulation.

Comparison of Alternative Splicing Junction Detection Tools Using RNA-Seq Data

BACKGROUND

Alternative splicing (AS) is a posttranscriptional process that produces differ-ent transcripts from the same gene and is important to produce diverse protein products in response to environmental stimuli. AS occurs at specific sites on the mRNA sequence, some of which have been de-fined. Multiple bioinformatics tools have been developed to detect AS from experimental data.

OBJECTIVES

The goal of this review is to help researchers use specific tools to aid their research and to develop new AS detection tools based on these previously established tools.

METHOD

We selected 15 AS detection tools that were recently published; we classified and delineated them on several aspects. Also, a performance comparison of these tools with the same starting input was conducted.

RESULT

We reviewed the following categorized features of the tools: Publication information, working principles, generic and distinct workflows, running platform, input data requirement, sequencing depth dependency, reads mapped to multiple locations, isoform annotation basis, precise detected AS types, and performance benchmarks.

CONCLUSION

Through comparisons of these tools, we provide a panorama of the advantages and short-comings of each tool and their scopes of application.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2896-7) contains supplementary material, which is available to authorized users.

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.