Modeling expression ranks for noise-tolerant differential expression analysis of scRNA-seq data.

Systematic delineation of complex biological systems is an ever-challenging and resource-intensive process. Single-cell transcriptomics allows us to study cell-to-cell variability in complex tissues at an unprecedented resolution. Accurate modeling of gene expression plays a critical role in the statistical determination of tissue-specific gene expression patterns. In the past few years, considerable efforts have been made to identify appropriate parametric models for single-cell expression data. The zero-inflated version of Poisson/negative binomial and log-normal distributions have emerged as the most popular alternatives owing to their ability to accommodate high dropout rates, as commonly observed in single-cell data. Although the majority of the parametric approaches directly model expression estimates, we explore the potential of modeling expression ranks, as robust surrogates for transcript abundance. Here we examined the performance of the discrete generalized beta distribution (DGBD) on real data and devised a Wald-type test for comparing gene expression across two phenotypically divergent groups of single cells. We performed a comprehensive assessment of the proposed method to understand its advantages compared with some of the existing best-practice approaches. We concluded that besides striking a reasonable balance between Type I and Type II errors, ROSeq, the proposed differential expression test, is exceptionally robust to expression noise and scales rapidly with increasing sample size. For wider dissemination and adoption of the method, we created an R package called ROSeq and made it available on the Bioconductor platform.

Corrigendum: Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis.

[This corrects the article on p. 70 in vol. 4, PMID: 27709111.].

Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis.

The study of single cells has evolved over the past several years to include expression and genomic analysis of an increasing number of single cells. Several studies have demonstrated wide spread variation and heterogeneity within cell populations of similar phenotype. While the characterization of these populations will likely set the foundation for our understanding of genomic- and expression-based diversity, it will not be able to link the functional differences of a single cell to its underlying genomic structure and activity. Currently, it is difficult to perturb single cells in a controlled environment, monitor and measure the response due to perturbation, and link these response measurements to downstream genomic and transcriptomic analysis. In order to address this challenge, we developed a platform to integrate and miniaturize many of the experimental steps required to study single-cell function. The heart of this platform is an elastomer-based integrated fluidic circuit that uses fluidic logic to select and sequester specific single cells based on a phenotypic trait for downstream experimentation. Experiments with sequestered cells that have been performed include on-chip culture, exposure to various stimulants, and post-exposure image-based response analysis, followed by preparation of the mRNA transcriptome for massively parallel sequencing analysis. The flexible system embodies experimental design and execution that enable routine functional studies of single cells.

Isolation and amplification of mRNA within a simple microfluidic lab on a chip.

The major modules for realizing molecular biological assays in a micro-total analysis system (µTAS) were developed for the detection of pathogenic organisms. The specific focus was the isolation and amplification of eukaryotic mRNA within a simple, single-channel device for very low RNA concentrations that could then be integrated with detection modules. The hsp70 mRNA from Cryptosporidium parvum was used as a model analyte. Important points of study were surface chemistries within poly(methyl methacrylate) (PMMA) microfluidic channels that enabled specific and sensitive mRNA isolation and amplification reactions for very low mRNA concentrations. Optimal conditions were achieved when the channel surface was carboxylated via UV/ozone treatment followed by the immobilization of polyamidoamine (PAMAM) dendrimers on the surface, thus increasing the immobilization efficiency of the thymidine oligonucleotide, oligo(dT)25, and providing a reliable surface for the amplification reaction, importantly, without the need for blocking agents. Additional chemical modifications of the remaining active surface groups were studied to avoid nonspecific capturing of nucleic acids and hindering of the mRNA amplification at low RNA concentrations. Amplification of the mRNA was accomplished using nucleic acid sequence-based amplification (NASBA), an isothermal, primer-dependent technique. Positive controls consisting of previously generated NASBA amplicons could be diluted 10(15) fold and still result in successful on-chip reamplification. Finally, the successful isolation and amplification of mRNA from as few as 30 C. parvum oocysts was demonstrated directly on-chip and compared to benchtop devices. This is the first proof of successful mRNA isolation and NASBA-based amplification of mRNA within a simple microfluidic device in relevant analytical volumes.