Sensitive high-throughput single-cell RNA-seq reveals within-clonal transcript correlations in yeast populations

Single-cell RNA sequencing has revealed extensive cellular heterogeneity within many organisms, but few methods have been developed for microbial clonal populations. The yeast genome displays unusually dense transcript spacing, with interleaved and overlapping transcription from both strands, resulting in a minuscule but complex pool of RNA that is protected by a resilient cell wall. Here, we have developed a sensitive, scalable and inexpensive yeast single-cell RNA-seq (yscRNA-seq) method that digitally counts transcript start sites in a strand- and isoform-specific manner. YscRNA-seq detects the expression of low-abundance, noncoding RNAs and at least half of the protein-coding genome in each cell. In clonal cells, we observed a negative correlation for the expression of sense-antisense pairs, whereas paralogs and divergent transcripts co-expressed. By combining yscRNA-seq with index sorting, we uncovered a linear relationship between cell size and RNA content. Although we detected an average of ~3.5 molecules per gene, the number of expressed isoforms is restricted at the single-cell level. Remarkably, the expression of metabolic genes is highly variable, whereas their stochastic expression primes cells for increased fitness towards the corresponding environmental challenge. These findings suggest that functional transcript diversity acts as a mechanism that provides a selective advantage to individual cells within otherwise transcriptionally heterogeneous populations.