PARRoT- a homology-based strategy to quantify and compare RNA-sequencing from non-model organisms

Rhythms in lipid droplet content driven by a metabolic oscillator are conserved throughout evolution

The biological clock in eukaryotes controls daily rhythms in physiology and behavior. It displays a complex organization that involves the molecular transcriptional clock and the redox oscillator which may coordinately work to control cellular rhythms. The redox oscillator has emerged very early in evolution in adaptation to the environmental changes in O2 levels and has been shown to regulate daily rhythms in glycerolipid (GL) metabolism in different eukaryotic cells. GLs are key components of lipid droplets (LDs), intracellular storage organelles, present in all living organisms, and essential for energy and lipid homeostasis regulation and survival; however, the cell bioenergetics status is not constant across time and depends on energy demands. Thus, the formation and degradation of LDs may reflect a time-dependent process following energy requirements. This work investigated the presence of metabolic rhythms in LD content along evolution by studying prokaryotic and eukaryotic cells and organisms. We found sustained temporal oscillations in LD content in Pseudomonas aeruginosa bacteria and Caenorhabditis elegans synchronized by temperature cycles, in serum-shock synchronized human embryonic kidney cells (HEK 293 cells) and brain tumor cells (T98G and GL26) after a dexamethasone pulse. Moreover, in synchronized T98G cells, LD oscillations were altered by glycogen synthase kinase-3 (GSK-3) inhibition that affects the cytosolic activity of the metabolic oscillator or by knocking down LIPIN-1, a key GL synthesizing enzyme. Overall, our findings reveal the existence of metabolic oscillations in terms of LD content highly conserved across evolutionary scales notwithstanding variations in complexity, regulation, and cell organization.

Rhythmic Component Analysis Tool (RCAT): A Precise, Efficient and User-Friendly Tool for Circadian Clock Genes Analysis

High-throughput next-generation sequencing (NGS) technologies and real-time circadian dynamics reporting systems produce large amounts of experimental data on RNA and protein levels in the field of circadian rhythm and therefore require statistical knowledge and computational skills for quantitative analysis. Although there are many software applications that can process these data, they are often difficult to use and computationally inefficient. Hence, a convenient, user-friendly tool that can accurately acquire rhythmic components (period, amplitude, and phase) of circadian clock genes is necessary. Here, we develop a new analysis tool named rhythmic component analysis tool (RCAT), which has an easily understood interface featuring a one-button operation, that presents all results as tables and images and automatically saves them as CSV files. We use the relative amplitude error (RAE), widely-adopted criteria on the circadian research field to estimate the quality of results. To illustrate the analytical ability of the RCAT under different situations, we generate four groups of time-series data by CircaInSilico (a web server for generating synthetic genome biology data to benchmark statistical methods for studying biological rhythms) with different collection intervals and amplitude ranges and use RCAT to analyze them. To demonstrate the effectiveness of RCAT, we analyze two sets of case studies with time-series data: one set uses microarray and RNA-Seq data from the gene expression omnibus (GEO) repository to identify core clock genes (CCGs) with significant periodicity in the liver, and the other set uses real-time fluorescence reporting data collected by Lumicycle® (a commonly-used luminometer) to calculate the precise period, amplitude and phase. In these examples, most cycling samples are successfully detected by the RCAT within a short collection time, and accurate rhythmic components are also successfully computed. These results indicate that RCAT improves flexibility and convenience in periodic oscillation data analysis. RCAT, is freely available at: https://github.com/lzbbest/Rhythmic-Component-Analysis-Tool/releases . It, as a cross-platform software, can be run not only on Linux, but also on Win10, Win8 and Win7.

Structural variation analysis of 6,500 whole genome sequences in amyotrophic lateral sclerosis

There is a strong genetic contribution to Amyotrophic lateral sclerosis (ALS) risk, with heritability estimates of up to 60%. Both Mendelian and small effect variants have been identified, but in common with other conditions, such variants only explain a little of the heritability. Genomic structural variation might account for some of this otherwise unexplained heritability. We therefore investigated association between structural variation in a set of 25 ALS genes, and ALS risk and phenotype. As expected, the repeat expansion in the C9orf72 gene was identified as associated with ALS. Two other ALS-associated structural variants were identified: inversion in the VCP gene and insertion in the ERBB4 gene. All three variants were associated both with increased risk of ALS and specific phenotypic patterns of disease expression. More than 70% of people with respiratory onset ALS harboured ERBB4 insertion compared with 25% of the general population, suggesting respiratory onset ALS may be a distinct genetic subtype.

Long non-coding RNA ZFAS1 sponges miR-486 to promote osteosarcoma cells progression and metastasis in vitro and vivo

Long noncoding RNAs (lncRNAs) have been wildly demonstrated to participate in the osteosarcoma tumorigenesis. ZFAS1 is a novel identified lncRNA, however, its role in osteosarcoma is still unclear. In present study, we utilize lncRNA microarray assay to screen the lncRNA expression profile in osteosarcoma tissue, and investigate the regulatory function of ZFAS1 in osteosarcoma. LncRNA microarray assay revealed that lncRNA ZFAS1 was significantly up-regulated in 3 pairs of osteosarcoma and adjacent non-tumor tissue, which was confirmed by RT-PCR. Furthermore, in 53 pairs of osteosarcoma patient samples, the up-regulated expression of ZFAS1 was closely related to poor prognosis. In vitro, loss-of-function experiments showed that ZFAS1 knockdown significantly suppressed the proliferation, induced cycle arrest at G0/G1 phase and enhance apoptosis. In vivo, ZFAS1 knockdown inhibited the tumor growth. Bioinformatics online programs predicted that ZFAS1 sponge miR-486 at 3’-UTR with complementary binding sites, which was validated using luciferase reporter assay and RNA immunoprecipitation (RIP) assay. Rescue experiments confirmed that miR-486 could reverse the functions of ZFAS1 on osteosarcoma genesis. In conclusion, our results demonstrate that ZFAS1 act as competing endogenous RNA (ceRNA) for miR-486, and act as oncogene in osteosarcoma tumorigenesis, and discover the functional regulatory pathway of ZFAS1 sponging miR-486.

Comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri identifies potential genes related to triterpenoid saponin biosynthesis

Background

Bacopa monnieri commonly known as Brahmi is utilized in Ayurveda to improve memory and many other human health benefits. Bacosides enriched standardized extract of Bacopa monnieri is being marketed as a memory enhancing agent. In spite of its well known pharmacological properties it is not much studied in terms of transcripts involved in biosynthetic pathway and its regulation that controls the secondary metabolic pathway in this plant. The aim of this study was to identify the potential transcripts and provide a framework of identified transcripts involved in bacosides production through transcriptome assembly.

Results

We performed comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri in two independent biological replicate and obtained 22.48 million and 22.0 million high quality processed reads in shoot and root respectively. After de novo assembly and quantitative assessment total 26,412 genes got annotated in root and 18,500 genes annotated in shoot sample. Quality of raw reads was determined by using SeqQC-V2.2. Assembled sequences were annotated using BLASTX against public database such as NR or UniProt. Searching against the KEGG pathway database indicated that 37,918 unigenes from root and 35,130 unigenes from shoot were mapped to 133 KEGG pathways. Based on the DGE data we found that most of the transcript related to CYP450s and UDP-glucosyltransferases were specifically upregulated in shoot tissue as compared to root tissue. Finally, we have selected 43 transcripts related to secondary metabolism including transcription factor families which are differentially expressed in shoot and root tissues were validated by qRT-PCR and their expression level were monitored after MeJA treatment and wounding for 1, 3 and 5 h.

Conclusions

This study not only represents the first de novo transcriptome analysis of Bacopa monnieri but also provides information about the identification, expression and differential tissues specific distribution of transcripts related to triterpenoid sapogenin which is one of the most important pharmacologically active secondary metabolite present in Bacopa monnieri. The identified transcripts in this study will establish a foundation for future studies related to carrying out the metabolic engineering for increasing the bacosides biosynthesis and its regulation for human health benefits.

Electronic supplementary material

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

Bioinformatics and systems biology research update from the 15th International Conference on Bioinformatics (InCoB2016)

The International Conference on Bioinformatics (InCoB) has been publishing peer-reviewed conference papers in BMC Bioinformatics since 2006. Of the 44 articles accepted for publication in supplement issues of BMC Bioinformatics, BMC Genomics, BMC Medical Genomics and BMC Systems Biology, 24 articles with a bioinformatics or systems biology focus are reviewed in this editorial. InCoB2017 is scheduled to be held in Shenzen, China, September 20-22, 2017.