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

Comparative Transcriptome Analysis of Deep-Rooting and Shallow-Rooting Potato (Solanum tuberosum L.) Genotypes under Drought Stress

The selection and breeding of deep rooting and drought-tolerant varieties has become a promising approach for improving the yield and adaptability of potato (Solanum tuberosum L.) in arid and semiarid areas. Therefore, the discovery of root-development-related genes and drought tolerance signaling pathways in potato is important. In this study, we used deep-rooting (C119) and shallow-rooting (C16) potato genotypes, with different levels of drought tolerance, to achieve this objective. Both genotypes were treated with 150 mM mannitol for 0 h (T0), 2 h (T2), 6 h (T6), 12 h (T12), and 24 h (T24), and their root tissues were subjected to comparative transcriptome analysis. A total of 531, 1571, 1247, and 3540 differentially expressed genes (DEGs) in C16 and 1531, 1108, 674, and 4850 DEGs in C119 were identified in T2 vs. T0, T6 vs. T2, T12 vs. T6, and T24 vs. T12 comparisons, respectively. Gene expression analysis indicated that a delay in the onset of drought-induced transcriptional changes in C16 compared with C119. Functional enrichment analysis revealed genotype-specific biological processes involved in drought stress tolerance. The metabolic pathways of plant hormone transduction and MAPK signaling were heavily involved in the resistance of C16 and C119 to drought, while abscisic acid (ABA), ethylene, and salicylic acid signal transduction pathways likely played more important roles in C119 stress responses. Furthermore, genes involved in root cell elongation and division showed differential expression between the two genotypes under drought stress. Overall, this study provides important information for the marker-assisted selection and breeding of drought-tolerant potato genotypes.

miRNA-150-5p promotes hepatic stellate cell proliferation and sensitizes hepatocyte apoptosis during liver fibrosis

Aim: To explore the role of miRNA-150-5p (miR-150-5p) in liver fibrosis. Materials & methods: miRNA expression profiles, CCl₄-induced liver fibrosis progression and regression rodent models, quantitative real-time PCR, miR-150-5p mimics and inhibitors, cell proliferation and apoptosis detection, RNA sequencing and bioinformatics analysis were employed. Results: Liver tissue miR-150-5p expression was positively associated with liver fibrosis progression and regression; however, miR-150-5p exhibited a cell-specific expression pattern, namely, it was enhanced in hepatocytes but reduced in hepatic stellate cells (HSCs) during liver fibrosis; miR-150-5p overexpression promoted HSC apoptosis and sensitized hepatocyte apoptosis; miR-150-5p mimic had a larger influence on the transcriptomic stability of HSCs than that of hepatocytes; miR-150-5p mediated activation of interferon signaling pathways might be responsible for HSC apoptosis. Conclusion: miR-150-5p exhibited an opposite regulation and function pattern between HSCs and hepatocytes during liver fibrosis.

A Brief Overview of lncRNAs in Endothelial Dysfunction-Associated Diseases: From Discovery to Characterization

Long non-coding RNAs (lncRNAs) are a novel class of regulatory RNA molecules and they are involved in many biological processes and disease developments. Several unique features of lncRNAs have been identified, such as tissue-and/or cell-specific expression pattern, which suggest that they could be potential candidates for therapeutic and diagnostic applications. More recently, the scope of lncRNA studies has been extended to endothelial biology research. Many of lncRNAs were found to be critically involved in the regulation of endothelial function and its associated disease progression. An improved understanding of endothelial biology can thus facilitate the discovery of novel biomarkers and therapeutic targets for endothelial dysfunction-associated diseases, such as abnormal angiogenesis, hypertension, diabetes, and atherosclerosis. Nevertheless, the underlying mechanism of lncRNA remains undefined in previous published studies. Therefore, in this review, we aimed to discuss the current methodologies for discovering and investigating the functions of lncRNAs and, in particular, to address the functions of selected lncRNAs in endothelial dysfunction-associated diseases.

ARHGEF4-mediates the actin cytoskeleton reorganization of hepatic stellate cells in 3-dimensional collagen matrices

The actin cytoskeleton of hepatic stellate cells (HSCs) is reorganized when they are cultured in 3D collagen matrices. Here, we investigated the molecular mechanism of actin cytoskeleton reorganization in HSCs cultured in 3D floating collagen matrices (FCM) compared to those on 2D polystyrene surfaces (PS). First, we found that the generation of dendritic cellular processes was controlled by Rac1. Next, we examined the differential gene expression of HSCs cultured on 2D PS and in 3D FCM by RNA-Seq and focused on the changes of actin cytoskeleton reorganization-related molecular components and guanine nucleotide exchange factors (GEFs). The results showed that the expression of genes associated with actin cytoskeleton reorganization-related cellular components, filopodia and lamellipodia, were significantly decreased, but podosome-related genes was significantly increased in 3D FCM. Furthermore, we found that a Rac1-specific GEF, ARHGEF4, played roles in morphological changes, migration and podosome-related gene expression in HSCs cultured in 3D FCM.

Abbreviations: 2D PS: 2-dimensional polystyrene surface; 3D FCM: 3-dimensional floating collagen matrices; ARHGEF4: Rho guanine nucleotide exchange factor 4; ARHGEF6: Rho guanine nucleotide exchange factor 6; GEF: guanine nucleotide exchange factor; HSC: hepatic stellate cell

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.

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.

The impact of the RBM4-initiated splicing cascade on modulating the carcinogenic signature of colorectal cancer cells

A growing body of studies has demonstrated that dysregulated splicing profiles constitute pivotal mechanisms for carcinogenesis. In this study, we identified discriminative splicing profiles of colorectal cancer (CRC) cells compared to adjacent normal tissues using deep RNA-sequencing (RNA-seq). The RNA-seq results and cohort studies indicated a relatively high ratio of exon 4-excluded neuro-oncological ventral antigen 1 (Nova1⁻⁴) and intron 2-retained SRSF6 (SRSF6^{+intron 2}) transcripts in CRC tissues and cell lines. Nova1 variants exhibited differential effects on eliminating SRSF6 expression in CRC cells by inducing SRSF6^{+intron 2} transcripts which were considered to be the putative target of alternative splicing-coupled nonsense-mediated decay mechanism. Moreover, the splicing profile of vascular endothelial growth factor (VEGF)165/VEGF165b transcripts was relevant to SRSF6 expression, which manipulates the progression of CRC calls. These results highlight the novel and hierarchical role of an alternative splicing cascade that is involved in the development of CRC.

Intelligent biology and medicine in 2015: advancing interdisciplinary education, collaboration, and data science

We summarize the 2015 International Conference on Intelligent Biology and Medicine (ICIBM 2015) and the editorial report of the supplement to BMC Genomics. The supplement includes 20 research articles selected from the manuscripts submitted to ICIBM 2015. The conference was held on November 13-15, 2015 at Indianapolis, Indiana, USA. It included eight scientific sessions, three tutorials, four keynote presentations, three highlight talks, and a poster session that covered current research in bioinformatics, systems biology, computational biology, biotechnologies, and computational medicine.