Wei Wang

A fundamental assumption, common to the vast majority of high-throughput transcriptome analyses, is that the expression of most genes is unchanged among samples and that total cellular RNA remains constant. As the number of analyzed experimental systems increases however, different independent studies demonstrate that this assumption is often violated. We present a calibration method using RNA spike-ins that allows for the measurement of absolute cellular abundance of RNA molecules. We apply the method to pooled RNA from cell populations of known sizes. For each transcript, we compute a nominal abundance that can be converted to absolute by dividing by a scale factor determined in separate experiments: the yield coefficient of the transcript relative to that of a reference spike-in measured with the same protocol. The method is derived by maximum likelihood theory in the context is a complete statistical model for sequencing counts contributed by cellular RNA and spike-ins. The counts are based on a sample from a fixed number of cells to which a fixed population of spike-in molecules has been added. We illustrate and evaluate the method with applications to two global expression data sets, one from the model eukaryote Saccharomyces cerevisiae, proliferating at different growth rates, and differentiating cardiopharyngeal cell lineages in the chordate Ciona robusta. We tested the method in a technical replicate dilution study, and in a k-fold validation study. We present a complete statistical model for the analysis of RNA-seq data from a population of cells using external RNA spike-ins and a maximum-likelihood method for genome-wide estimation of transcripts per cell. The model includes biological variability of cellular transcript number and sampling noise. We derive an unbiased estimator of transcripts per cell for every transcript, given by simply multiplying the count by a library-dependent, but transcript-independent, scale factor. This is a nominal estimate that can be converted to an absolute estimate by dividing by the transcript’s relative yield coefficient, measured in a separate experiment. A negative binomial probability mass function with novel normalization (size) factors allows for parametric testing of hypotheses concerning dependence of the absolute abundance of each transcript on experimental condition. Our method integrates information from every RNA-seq experiment across all replicates and experimental conditions to determine the calibration constants. We test the method with a dilution study and a k-fold cross-validation study. We illustrate our method with applications to two independent data sets from yeast and the sea squirt that were derived by different library preparation protocols. We show that our methods detect genome-wide amplification of expression, and we compare our method to others.

In embryos, multipotent progenitors divide to produce distinct progeny and express their full potential. In vertebrates, multipotent cardiopharyngeal progenitors produce second-heart-field-derived cardiomyocytes, and branchiomeric skeletal head muscles. However, the mechanisms underlying these early fate choices remain largely elusive. The tunicate Ciona emerged as an attractive model to study early cardiopharyngeal development at high resolution: through two asymmetric and oriented divisions, defined cardiopharyngeal progenitors produce distinct first and second heart precursors, and pharyngeal muscle (aka atrial siphon muscle, ASM) precursors. Here, we demonstrate that differential FGF-MAPK signaling distinguishes between heart and ASM precursors. We characterize a feed-forward circuit that promotes the successive activations of essential ASM determinants, Hand-related, Tbx1/10 and Ebf. Finally, we show that coupling FGF-MAPK restriction and cardiopharyngeal network deployment with cell divisions defines the timing of gene expression and permits the emergence of diverse cell types from multipotent progenitors.