New normal linear modeling strategies are presented for analyzing read counts from RNA-seq experiments. The voom method estimates the mean-variance relationship of the log-counts, generates a precision weight for each observation and enters these into the limma empirical Bayes analysis pipeline. This opens access for RNA-seq analysts to a large body of methodology developed for microarrays. Simulation studies show that voom performs as well or better than count-based RNA-seq methods even when the data are generated according to the assumptions of the earlier methods. Two case studies illustrate the use of linear modeling and gene set testing methods.

/pdf/voom-precision-weights-unlock-linear-model-analysis-tools-5bgw8lj2b8.pdf

voom: precision weights unlock linear model analysis tools for RNA-seq read counts

Droplet Barcoding for Single-Cell Transcriptomics Applied to Embryonic Stem Cells

Read alignment is an ongoing challenge for the analysis of data from sequencing technologies. This article proposes an elegantly simple multi-seed strategy, called seed-and-vote, for mapping reads to a reference genome. The new strategy chooses the mapped genomic location for the read directly from the seeds. It uses a relatively large number of short seeds (called subreads) extracted from each read and allows all the seeds to vote on the optimal location. When the read length is <160 bp, overlapping subreads are used. More conventional alignment algorithms are then used to fill in detailed mismatch and indel information between the subreads that make up the winning voting block. The strategy is fast because the overall genomic location has already been chosen before the detailed alignment is done. It is sensitive because no individual subread is required to map exactly, nor are individual subreads constrained to map close by other subreads. It is accurate because the final location must be supported by several different subreads. The strategy extends easily to find exon junctions, by locating reads that contain sets of subreads mapping to different exons of the same gene. It scales up efficiently for longer reads.

/pdf/the-subread-aligner-fast-accurate-and-scalable-read-mapping-2r44yfipz6.pdf

The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote

Over the last decade, the introduction of microarray technology has had a profound impact on gene expression research. The publication of studies with dissimilar or altogether contradictory results, obtained using different microarray platforms to analyze identical RNA samples, has raised concerns about the reliability of this technology. The MicroArray Quality Control (MAQC) project was initiated to address these concerns, as well as other performance and data analysis issues. Expression data on four titration pools from two distinct reference RNA samples were generated at multiple test sites using a variety of microarray-based and alternative technology platforms. Here we describe the experimental design and probe mapping efforts behind the MAQC project. We show intraplatform consistency across test sites as well as a high level of interplatform concordance in terms of genes identified as differentially expressed. This study provides a resource that represents an important first step toward establishing a framework for the use of microarrays in clinical and regulatory settings.

/pdf/the-microarray-quality-control-maqc-project-shows-inter-and-284p6wv6is.pdf

The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements

In multicellular organisms, biological function emerges when heterogeneous cell types form complex organs. Nevertheless, dissection of tissues into mixtures of cellular subpopulations is currently challenging. We introduce an automated massively parallel single-cell RNA sequencing (RNA-seq) approach for analyzing in vivo transcriptional states in thousands of single cells. Combined with unsupervised classification algorithms, this facilitates ab initio cell-type characterization of splenic tissues. Modeling single-cell transcriptional states in dendritic cells and additional hematopoietic cell types uncovers rich cell-type heterogeneity and gene-modules activity in steady state and after pathogen activation. Cellular diversity is thereby approached through inference of variable and dynamic pathway activity rather than a fixed preprogrammed cell-type hierarchy. These data demonstrate single-cell RNA-seq as an effective tool for comprehensive cellular decomposition of complex tissues.

/pdf/massively-parallel-single-cell-rna-seq-for-marker-free-uiiq82x1s8.pdf

Massively Parallel Single-Cell RNA-Seq for Marker-Free Decomposition of Tissues into Cell Types

Standard controls and best practice guidelines advance acceptance of data from research, preclinical and clinical laboratories by providing a means for evaluating data quality. The External RNA Controls Consortium (ERCC) is developing commonly agreed-upon and tested controls for use in expression assays, a true industry-wide standard control.

The External RNA Controls Consortium: a progress report.

Despite the widespread use of microarrays, much ambiguity regarding data analysis, interpretation and correlation of the different technologies exists. There is a considerable amount of interest in correlating results obtained between different microarray platforms. To date, only a few cross-platform evaluations have been published and unfortunately, no guidelines have been established on the best methods of making such correlations. To address this issue we conducted a thorough evaluation of two commercial microarray platforms to determine an appropriate methodology for making cross-platform correlations. In this study, expression measurements for 10,763 genes uniquely represented on Affymetrix U133A/B GeneChips® and Amersham CodeLink™ UniSet Human 20 K microarrays were compared. For each microarray platform, five technical replicates, derived from the same total RNA samples, were labeled, hybridized, and quantified according to each manufacturers' standard protocols. The correlation coefficient (r) of differential expression ratios for the entire set of 10,763 overlapping genes was 0.62 between platforms. However, the correlation improved significantly (r = 0.79) when genes within noise were excluded. In addition to levels of inter-platform correlation, we evaluated precision, statistical-significance profiles, power, and noise levels for each microarray platform. Accuracy of differential expression was measured against real-time PCR for 25 genes and both platforms correlated well with r values of 0.92 and 0.79 for CodeLink and GeneChip, respectively. As a result of this study, we recommend using only genes called 'present' in cross-platform correlations. However, as in this study, a large number of genes may be lost from the correlation due to differing levels of noise between platforms. This is an important consideration given the apparent difference in sensitivity of the two platforms. Data from microarray analysis need to be interpreted cautiously and therefore, we provide guidelines for making cross-platform correlations. In all, this study represents the most comprehensive and specifically designed comparison of short-oligonucleotide microarray platforms to date using the largest set of overlapping genes.

/pdf/performance-evaluation-of-commercial-short-oligonucleotide-2me2rqmwzy.pdf

Performance evaluation of commercial short-oligonucleotide microarrays and the impact of noise in making cross-platform correlations.

Parkinson's disease pathogenesis proceeds through several phases, culminating in the loss of dopaminergic neurons of the substantia nigra (SN). Although the 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) model of oxidative SN injury is frequently used to study degeneration of dopaminergic neurons in mice and non-human primates, an understanding of the temporal sequence of molecular events from inhibition of mitochondrial complex 1 to neuronal cell death is limited. Here, microarray analysis and integrative data mining were used to uncover pathways implicated in the progression of changes in dopaminergic neurons after MPTP administration. This approach enabled the identification of small, yet consistently significant, changes in gene expression within the SN of MPTP-treated animals. Such an analysis disclosed dysregulation of genes in three main areas related to neuronal function: cytoskeletal stability and maintenance, synaptic integrity, and cell cycle and apoptosis. The discovery and validation of these alterations provide molecular evidence for an evolving cascade of injury, dysfunction, and cell death.

https://www.jneurosci.org/content/jneuro/24/34/7445.full.pdf

Timothy J Sendera

Papers

The External RNA Controls Consortium: a progress report.

Performance evaluation of commercial short-oligonucleotide microarrays and the impact of noise in making cross-platform correlations.

Dysregulation of Gene Expression in the 1-Methyl-4-Phenyl-1,2,3,6-Tetrahydropyridine-Lesioned Mouse Substantia Nigra