RNA-Seq Blog

Small RNAs derived from structural non-coding RNAs

RNA-Seq Blog Administrator — Thu, 23 May 2013 12:07:45 +0000

It has been shown in small RNA sequencing-based studies that some small RNA fragments are specifically processed from known structural non-coding RNAs, either through Dicer-dependent or Dicer-independent pathways. Although these small RNAs are often less abundant compared to microRNAs in normal mammalian tissues, they are always present in all sequenced libraries. In this paper, researchers from the Institut Curie, France use the ncPRO-seq pipeline, to describe different profiles of these small RNA fragments, and to discuss their potential processing pathways and functions. To assess whether more small RNA fragments can be detected in small RNA sequencing datasets, they decided to focus on small nuclear RNAs, abbreviated as snRNAs, which are associated with Sm ribonucleoproteins to form functional RNA-protein complexes. Here, they describe a group of small RNA fragments derived from snRNAs, which are typically highly enriched in regions bound by Sm proteins. Based on this, they propose the existence of a potential novel small RNA family associated with Sm proteins.

Chen CJ, Heard E. (2013) Small RNAs derived from structural non-coding RNAs. Methods [Epub ahead of print]. [astract]

Small RNAs derived from structural non-coding RNAs is a post from: RNA-Seq Blog

Incoming search terms:

rna-seq adaptation

Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells

RNA-Seq Blog Administrator — Thu, 23 May 2013 11:51:32 +0000

Recent molecular studies have shown that, even when derived from a seemingly homogenous population, individual cells can exhibit substantial differences in gene expression, protein levels and phenotypic output, with important functional consequences. Existing studies of cellular heterogeneity, however, have typically measured only a few pre-selected RNAs or proteins simultaneously, because genomic profiling methods could not be applied to single cells until very recently.

Here, a team led by researchers at Harvard University used single-cell RNA sequencing to investigate heterogeneity in the response of mouse bone-marrow-derived dendritic cells (BMDCs) to lipopolysaccharide. They found extensive, and previously unobserved, bimodal variation in messenger RNA abundance and splicing patterns, which they validated by RNA-fluorescence in situ hybridization for select transcripts. In particular, hundreds of key immune genes are bimodally expressed across cells, surprisingly even for genes that are very highly expressed at the population average. Moreover, splicing patterns demonstrate previously unobserved levels of heterogeneity between cells. Some of the observed bimodality can be attributed to closely related, yet distinct, known maturity states of BMDCs; other portions reflect differences in the usage of key regulatory circuits. For example, they identified a module of 137 highly variable, yet co-regulated, antiviral response genes. Using cells from knockout mice, the researchers show that variability in this module may be propagated through an interferon feedback circuit, involving the transcriptional regulators Stat2 and Irf7. This study demonstrates the power and promise of single-cell genomics in uncovering functional diversity between cells and in deciphering cell states and circuits.

Shalek AK, Satija R, Adiconis X, Gertner RS, Gaublomme JT, Raychowdhury R, Schwartz S, Yosef N, Malboeuf C, Lu D, Trombetta JT, Gennert D, Gnirke A, Goren A, Hacohen N, Levin JZ, Park H, Regev A. (2013) Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells. Nature [Epub ahead of print]. [abstract]

Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells is a post from: RNA-Seq Blog

Incoming search terms:

quantification of splicejunctionjs
splice index NGS

Graphite Web: web tool for gene set analysis exploiting pathway topology

RNA-Seq Blog Administrator — Wed, 22 May 2013 12:08:28 +0000

Graphite web is a novel web tool for pathway analyses and network visualization for gene expression data of both microarray and RNA-seq experiments. Several pathway analyses have been proposed either in the univariate or in the global and multivariate context to tackle the complexity and the interpretation of expression results. These methods can be further divided into ‘topological’ and ‘non-topological’ methods according to their ability to gain power from pathway topology. Biological pathways are, in fact, not only gene lists but can be represented through a network where genes and connections are, respectively, nodes and edges. To this day, the most used approaches are non-topological and univariate although they miss the relationship among genes. On the contrary, topological and multivariate approaches are more powerful, but difficult to be used by researchers without bioinformatic skills.

Here, researchers from the University of Padova, Italy present Graphite web, the first public web server for pathway analysis on gene expression data that combines topological and multivariate pathway analyses with an efficient system of interactive network visualizations for easy results interpretation. Specifically, Graphite web implements five different gene set analyses on three model organisms and two pathway databases.

Availability – Graphite Web is freely available at http://graphiteweb.bio.unipd.it/.

Sales G, Calura E, Martini P, Romualdi C. (2013) Graphite Web: web tool for gene set analysis exploiting pathway topology. Nucleic Acids Res [Epub ahead of print]. [article]

Graphite Web: web tool for gene set analysis exploiting pathway topology is a post from: RNA-Seq Blog

Incoming search terms:

www rna-seqblog com graphite-web-web-tool-for-gene-set-analysis-exploiting-pathway-topology

eQTL Mapping Using RNA-seq Data

RNA-Seq Blog Administrator — Wed, 22 May 2013 12:01:27 +0000

As RNA-seq is replacing gene expression microarrays to assess genome-wide transcription abundance, gene expression Quantitative Trait Locus (eQTL) studies using RNA-seq have emerged. RNA-seq delivers two novel features that are important for eQTL studies. First, it provides information on allele-specific expression (ASE), which is not available from gene expression microarrays. Second, it generates unprecedentedly rich data to study RNA-isoform expression. In this paper, the authors review current methods for eQTL mapping using ASE and discuss some future directions. They also review existing works that use RNA-seq data to study RNA-isoform expression and we discuss the gaps between these works and isoform-specific eQTL mapping.

Sun W, Hu Y. (2013) eQTL Mapping Using RNA-seq Data. Stat Biosci 5(1), 198-219. [article]

eQTL Mapping Using RNA-seq Data is a post from: RNA-Seq Blog

Incoming search terms:

www rna-seqblog com eqtl-mapping-using-rna-seq-data

Voom! Precision weights unlock linear model analysis tools for RNA-Seq read counts

RNA-Seq Blog Administrator — Tue, 21 May 2013 12:02:39 +0000

Voom: variance modelling at the observation-level

In the past few years, RNA-seq has emerged as a revolutionary new technology for expression profiling. RNA-seq expression data consists of read counts, and many recent publications have argued therefore that RNA-seq data should be analysed by statistical methods designed specifically for counts. Yet all the statistical methods developed for RNA-seq counts rely on approximations of various kinds.

This article revisits the idea of applying normal-based microarray-like statistical methods to RNA-seq read counts, with the idea that it is more important to model the mean-variance relationship correctly than it is to specify the exact probabilistic distribution of the counts. Log-counts per million are used as expression values. The voom method estimates the mean-variance relationship robustly and generates a precision weight for each individual normalized observation. The normalized log-counts per million and associated precision weights are then entered into the limma analysis pipeline, or indeed into any statistical pipeline for microarray data that is precision weight aware. This opens access for RNA-seq analysts to a large body of methodology developed for microarrays, allowing RNA-seq and microarray data to be analysed in closely comparable ways. The performance of voom and related limma-based pipelines is compared to that of edgeR, DESeq, baySeq, TSPM, PoissonSeq, and DSS. Simulation studies show that voom out-performs previous RNA-seq methods even when the data is generated according to the assumptions of the earlier methods. This is especially true when the sequence depths vary between RNA samples. Several data sets are also analysed to demonstrate how voom can handle heterogeneous data and complex experiments as well as facilitating pathway analysis and gene set testing methods.

Incoming search terms:

The RNA-seq Tuxedo pipeline
www rna-seqblog com voom-precision-weights-unlock-linear-model-analysis-tools-for-rna-seq-read-counts

RNA-seq analysis of allele-specific expression, hybrid effects, and regulatory divergence in hybrids

RNA-Seq Blog Administrator — Tue, 21 May 2013 11:51:08 +0000

Hybridization is a prominent process among natural plant populations that can result in phenotypic novelty, heterosis, and changes in gene expression. The effects of intraspecific hybridization on F1 hybrid gene expression were investigated using parents from divergent, natural populations of Cirsium arvense, an invasive Compositae weed.

Using an RNA-seq approach, the expression of 68,746 unigenes was quantified in parents and hybrids. The expression levels of 51% of transcripts differed between parents, a majority of which had <1.25x fold-changes. More unigenes had higher expression in the invasive parent (P1) than the non-invasive parent (P2). Of those that were divergently expressed between parents, 10% showed additive and 81% showed non-additive (transgressive or dominant) modes of gene action in the hybrids. A majority of the dominant cases had P2-like expression patterns in the hybrids. Comparisons of allele-specific expression also enabled a survey of cis- and trans-regulatory effects. Cis- and trans-regulatory divergence was found at 70% and 68% of 62,281 informative SNP sites, respectively. Of the 17% of sites exhibiting both cis- and trans- effects, a majority (70%) had antagonistic regulatory interactions (cis x trans); trans-divergence tended to drive higher expression of the P1 allele whereas cis-divergence tended to increase P2 transcript abundance. Trans-effects correlated more highly than cis- with parental expression divergence and accounted for a greater proportion of the regulatory divergence at sites with additive compared to non-additive inheritance patterns. This study explores the nature of, and types of mechanisms underlying, expression changes that occur in upon intraspecific hybridization in natural populations.

Bell GD, Kane NC, Rieseberg LH, Adams KL. (2013) RNA-seq analysis of allele-specific expression, hybrid effects, and regulatory divergence in hybrids compared with their parents from natural populations. Genome Biol Evol [Epub ahead of print]. [abstract]

RNA-seq analysis of allele-specific expression, hybrid effects, and regulatory divergence in hybrids is a post from: RNA-Seq Blog

Incoming search terms:

remove rrna biotin probe plos one
zggwcnpdyzw

A bi-Poisson model for clustering gene expression profiles by RNA-seq

RNA-Seq Blog Administrator — Mon, 20 May 2013 12:17:59 +0000

With the availability of gene expression data by RNA-seq, powerful statistical approaches for grouping similar gene expression profiles across different environments have become increasingly important. A team led by researchers at Penn State University describe and assess a computational model for clustering genes into distinct groups based on the pattern of gene expression in response to changing environment. The model capitalizes on the Poisson distribution to capture the count property of RNA-seq data. A two-stage hierarchical expectation-maximization (EM) algorithm is implemented to estimate an optimal number of groups and mean expression amounts of each group across two environments. A procedure is formulated to test whether and how a given group shows a plastic response to environmental changes. The impact of gene-environment interactions on the phenotypic plasticity of the organism can also be visualized and characterized. The model was used to analyse an RNA-seq dataset measured from two cell lines of breast cancer that respond differently to an anti-cancer drug, from which genes associated with the resistance and sensitivity of the cell lines are identified. They performed simulation studies to validate the statistical behaviour of the model. The model provides a useful tool for clustering gene expression data by RNA-seq, facilitating understanding of gene functions and networks.

Wang N, Wang Y, Hao H, Wang L, Wang Z, Wang J, Wu R. (2013) A bi-Poisson model for clustering gene expression profiles by RNA-seq. Brief Bioinform [Epub ahead of print]. [abstract]

A bi-Poisson model for clustering gene expression profiles by RNA-seq is a post from: RNA-Seq Blog

Incoming search terms:

rna-seq identified a super-long intergenic transcript functioning in adipoge
gene expression heart vertebrate
illumina sequencing scriptseq
star alignment r rnaseq
www rna-seqblog com a-bi-poisson-model-for-clustering-gene-expression-profiles-by-rna-seq

Optimizing de novo assembly of short-read RNA-seq data for phylogenomics

RNA-Seq Blog Administrator — Mon, 20 May 2013 11:55:08 +0000

RNA-seq has shown huge potential for phylogenomic inferences in non-model organisms. However, error, incompleteness, and redundant assembled transcripts for each gene in de novo assembly of short reads cause noise in analyses and a large amount of missing data in the aligned matrix. To address these problems, we compare de novo assemblies of paired end 90 bp RNA-seq reads using Oases, Trinity, Trans-ABySS and SOAPdenovo-Trans to transcripts from genome annotation of the model plant Ricinus communis. By doing so we evaluate strategies for optimizing total gene coverage and minimizing assembly chimeras and redundancy.

Researchers at the University of Michigan found that the frequency and structure of chimeras vary dramatically among different software packages. The differences were largely due to the number of trans-self chimeras that contain repeats in the opposite direction. More than half of the total chimeras in Oases and Trinity were trans-self chimeras. Within each package, they found a trade-off between maximizing reference coverage and minimizing redundancy and chimera rate.

In order to reduce redundancy, they investigated three methods:

1) using cap3 and CD-HIT-EST to combine highly similar transcripts,

2) only retaining the transcript with the highest read coverage, or removing the transcript with the lowest read coverage for each subcomponent in Trinity, and

3) filtering Oases single k-mer assemblies by number of transcripts per locus and relative transcript length, and then finding the transcript with the highest read coverage.

The researchers then utilized results from blastx against model protein sequences to effectively remove trans chimeras. After optimization, seven assembly strategies among all four packages successfully assembled 42.9–47.1% of reference genes to more than 200 bp, with a chimera rate of 0.92–2.21%, and on average 1.8–3.1 transcripts per reference gene assembled.

With rapidly improving sequencing and assembly tools, this study provides a framework to benchmark and optimize performance before choosing tools or parameter combinations for analyzing short-read RNA-seq data. The study demonstrates that choice of assembly package, k-mer sizes, post-assembly redundancy-reduction and chimera cleanup, and strand-specific RNA-seq library preparation and assembly dramatically improves gene coverage by non-redundant and non-chimeric transcripts that are optimized for downstream phylogenomic analyses.

Yang Y, Smith SA. (2013) Optimizing de novo assembly of short-read RNA-seq data for phylogenomics. BMC Genomics 14(1), 328. [abstract]

Optimizing de novo assembly of short-read RNA-seq data for phylogenomics is a post from: RNA-Seq Blog

Incoming search terms:

lokus no

Interview – Maarten Leerkes – Genome Analysis Specialist at the (NIAID) Office of Cyber Infrastructure and Computational Biology (OCICB)

RNA-Seq Blog Administrator — Fri, 17 May 2013 14:37:07 +0000

Maarten Leerkes speaks to Izzy Scott Moncrieff in the run up to the upcoming RNA-Seq 2013 Summit, 18th–20^th June 2013, Boston, MA

Maarten Leerkes provides bioinformatics support for various projects at NIAID as an employee of Medical Science and Computing, Inc. His areas of research include the use of bioinformatics to interpret sequencing data and to find patterns that can be extrapolated into diagnostic tools for improving treatment options for patients.

What initially attracted you to RNA-Seq?

During my doctoral training, I worked with open reading frame ESTs and SAGE platforms that aimed towards similar end products as RNA-Seq. The difference being that RNA-Seq has more sequencing depth and its information content is much higher.

A prevailing theme during my Ph.D. was a need for innovative technologies to process extra volumes of data and to answer certain types of research questions that were limited with existing technologies. With the development of RNA-Seq, many challenges to answering some of those questions were overcome. My interest, for example, was to answer questions related to alternative splicing and translocations that lead to fusions between gene products. The increased resolution of RNA Seq really opened up possibilities for me as I pursued my research.

Download a PDF of the entire interview at – http://rna-seqsummit.com/library

Interview – Maarten Leerkes – Genome Analysis Specialist at the (NIAID) Office of Cyber Infrastructure and Computational Biology (OCICB) is a post from: RNA-Seq Blog

Incoming search terms:

Ribo-zero library(LncRNA library)?
foto payudara dan memek hot
kontol dan memek
pofo memek hot

High-Resolution Transcriptome Maps Reveal Strain-Specific Regulatory Features of Multiple Campylobacter jejuni Isolates

RNA-Seq Blog Administrator — Fri, 17 May 2013 14:27:46 +0000

Many species have evolved into diverse strains with phenotypic and genotypic variations that facilitate adaptation to different ecological niches and, in the case of pathogens, to different hosts. Whereas comparison of genome sequences reveals differences and similarities among strains, the consequences of genomic variations can be tracked by studying the functional output from the genome. RNA sequencing has been revolutionizing transcriptome analyses of both pro- and eukaryotes. However, the bioinformatics-based analysis is still lagging behind, and transcriptome features are often manually annotated, which is laborious and time-consuming. This is even more compounded for the analyses of multiple strains.

Here, a team led by researchers at the University of Würzburg and the University of Tübingen, Germany compared the primary transcriptomes of four isolates of Campylobacter jejuni, the leading cause of bacterial gastroenteritis in humans, and provide genome-wide transcriptional start site (TSS) maps using a novel automated annotation method. Their comparative RNA–seq showed that most TSS are conserved in multiple strains, but they also observed SNP–dependent promoter usage. Furthermore, the researchers identified a novel minimal RNA–based CRISPR immune system as well as strain-specific small RNA repertoires. This automated, comparative TSS annotation will facilitate and improve transcriptome annotation for a wider range of organisms and provides insights into the contribution of transcriptome differences to phenotypic variation among closely related species.

Dugar G, Herbig A, Förstner KU, Heidrich N, Reinhardt R, et al. (2013) High-Resolution Transcriptome Maps Reveal Strain-Specific Regulatory Features of Multiple Campylobacter jejuni Isolates. PLoS Genet 9(5), e1003495. [article]

High-Resolution Transcriptome Maps Reveal Strain-Specific Regulatory Features of Multiple Campylobacter jejuni Isolates is a post from: RNA-Seq Blog

Incoming search terms:

memek india
High-Resolution Transcriptome Maps Reveal Strain-Specific Regulatory Features of Multiple Campylobacter jejuni Isolates
MicroRazerS

RNA-Seq Blog

Small RNAs derived from structural non-coding RNAs

Incoming search terms:

Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells

Incoming search terms:

Graphite Web: web tool for gene set analysis exploiting pathway topology

Incoming search terms:

eQTL Mapping Using RNA-seq Data

Incoming search terms:

Voom! Precision weights unlock linear model analysis tools for RNA-Seq read counts

Incoming search terms:

RNA-seq analysis of allele-specific expression, hybrid effects, and regulatory divergence in hybrids

Incoming search terms:

A bi-Poisson model for clustering gene expression profiles by RNA-seq

Incoming search terms:

Optimizing de novo assembly of short-read RNA-seq data for phylogenomics

Incoming search terms:

Interview – Maarten Leerkes – Genome Analysis Specialist at the (NIAID) Office of Cyber Infrastructure and Computational Biology (OCICB)

Maarten Leerkes speaks to Izzy Scott Moncrieff in the run up to the upcoming RNA-Seq 2013 Summit, 18th–20th June 2013, Boston, MA

Incoming search terms:

High-Resolution Transcriptome Maps Reveal Strain-Specific Regulatory Features of Multiple Campylobacter jejuni Isolates

Incoming search terms:

Maarten Leerkes speaks to Izzy Scott Moncrieff in the run up to the upcoming RNA-Seq 2013 Summit, 18th–20^th June 2013, Boston, MA