Abstract
The precise coordination of important biological processes, such as differentiation and development, relies heavily on the regulation of gene expression. In eukaryotic cells, understanding the distribution of RNA transcripts between the nucleus and cytosol is essential for gaining valuable insights into the process of protein production. The most efficient way to estimate the levels of RNA species genome-wide is through RNA sequencing (RNAseq). While RNAseq can be performed separately in the nucleus and in the cytosol, comparing transcript levels between compartments is challenging since measurements are relative to the unknown total RNA volume. Here, we show theoretically that if, in addition to nuclear and cytosolic RNAseq, whole-cell RNAseq is also performed, then accurate estimations of the localization of transcripts can be obtained. Based on this, we designed a method that estimates, first the fraction of the total RNA volume in the cytosol (nucleus), and then, this fraction for every transcript. We evaluate our methodology on simulated data and nuclear and cytosolic single-cell data available. Finally, we use our method to investigate the subcellular localization of transcripts using bulk RNAseq data from the ENCODE project.