Quartz-Seq: a highly reproducible and sensitive single-cell RNA sequencing method, reveals non-genetic gene-expression heterogeneity

Development of a highly reproducible and sensitive single-cell RNA sequencing (RNA-seq) method would facilitate the understanding of the biological roles and underlying mechanisms of non-genetic cellular heterogeneity. In this study, we report a novel single-cell RNA-seq method called Quartz-Seq that has a simpler protocol and higher reproducibility and sensitivity than existing methods. We show that single-cell Quartz-Seq can quantitatively detect various kinds of non-genetic cellular heterogeneity, and can detect different cell types and different cell-cycle phases of a single cell type. Moreover, this method can comprehensively reveal gene-expression heterogeneity between single cells of the same cell type in the same cell-cycle phase.

Quantitative expression profile of distinct functional regions in the adult mouse brain

The adult mammalian brain is composed of distinct regions with specialized roles including regulation of circadian clocks, feeding, sleep/awake, and seasonal rhythms. To find quantitative differences of expression among such various brain regions, we conducted the BrainStars (B*) project, in which we profiled the genome-wide expression of ∼50 small brain regions, including sensory centers, and centers for motion, time, memory, fear, and feeding. To avoid confounds from temporal differences in gene expression, we sampled each region every 4 hours for 24 hours, and pooled the samples for DNA-microarray assays. Therefore, we focused on spatial differences in gene expression. We used informatics to identify candidate genes with expression changes showing high or low expression in specific regions. We also identified candidate genes with stable expression across brain regions that can be used as new internal control genes, and ligand-receptor interactions of neurohormones and neurotransmitters. Through these analyses, we found 8,159 multi-state genes, 2,212 regional marker gene candidates for 44 small brain regions, 915 internal control gene candidates, and 23,864 inferred ligand-receptor interactions. We also found that these sets include well-known genes as well as novel candidate genes that might be related to specific functions in brain regions. We used our findings to develop an integrated database (http://brainstars.org/) for exploring genome-wide expression in the adult mouse brain, and have made this database openly accessible. These new resources will help accelerate the functional analysis of the mammalian brain and the elucidation of its regulatory network systems.

A functional genomics strategy reveals clockwork orange as a transcriptional regulator in the Drosophila circadian clock

The Drosophila circadian clock consists of integrated autoregulatory feedback loops, making the clock difficult to elucidate without comprehensively identifying the network components in vivo. Previous studies have adopted genome-wide screening for clock-controlled genes using high-density oligonucleotide arrays that identified hundreds of clock-controlled genes. In an attempt to identify the core clock genes among these candidates, we applied genome-wide functional screening using an RNA interference (RNAi) system in vivo. Here we report the identification of novel clock gene candidates including clockwork orange (cwo), a transcriptional repressor belonging to the basic helix-loop-helix ORANGE family. cwo is rhythmically expressed and directly regulated by CLK-CYC through canonical E-box sequences. A genome-wide search for its target genes using the Drosophila genome tiling array revealed that cwo forms its own negative feedback loop and directly suppresses the expression of other clock genes through the E-box sequence. Furthermore, this negative transcriptional feedback loop contributes to sustaining a high-amplitude circadian oscillation in vivo. Based on these results, we propose that the competition between cyclic CLK-CYC activity and the adjustable threshold imposed by CWO keeps E-box-mediated transcription within the controllable range of its activity, thereby rendering a Drosophila circadian clock capable of generating high-amplitude oscillation.

An improved single-cell cDNA amplification method for efficient high-density oligonucleotide microarray analysis

A systems-level understanding of a small but essential population of cells in development or adulthood (e.g. somatic stem cells) requires accurate quantitative monitoring of genome-wide gene expression, ideally from single cells. We report here a strategy to globally amplify mRNAs from single cells for highly quantitative high-density oligonucleotide microarray analysis that combines a small number of directional PCR cycles with subsequent linear amplification. Using this strategy, both the representation of gene expression profiles and reproducibility between individual experiments are unambiguously improved from the original method, along with high coverage and accuracy. The immediate application of this method to single cells in the undifferentiated inner cell masses of mouse blastocysts at embryonic day (E) 3.5 revealed the presence of two populations of cells, one with primitive endoderm (PE) expression and the other with pluripotent epiblast-like gene expression. The genes expressed differentially between these two populations were well preserved in morphologically differentiated PE and epiblast in the embryos one day later (E4.5), demonstrating that the method successfully detects subtle but essential differences in gene expression at the single-cell level among seemingly homogeneous cell populations. This study provides a strategy to analyze biophysical events in medicine as well as in neural, stem cell and developmental biology, where small numbers of distinctive or diseased cells play critical roles.