Application of RNA subcellular fraction estimation method to explore RNA localization regulation

RNA localization is involved in multiple biological processes. Recent advances in subcellular fractionation-based sequencing approaches uncovered localization pattern on a global scale. Most of existing methods adopt relative localization ratios (such as ratios of separately normalized transcripts per millions of different subcellular fractions without considering the difference in total RNA abundances in different fractions), however, absolute ratios may yield different results on the preference to different cellular compartment. Experimentally, adding external Spike-in RNAs to different fractionation can be used to obtain absolute ratios. In addition, a spike-in independent computational approach based on multiple linear regression model can also be used. However, currently, no custom tool is available. To solve this problem, we developed a method called subcellular fraction abundance estimator to correctly estimate relative RNA abundances of different subcellular fractionations. The ratios estimated by our method were consistent with existing reports. By applying the estimated ratios for different fractions, we explored the RNA localization pattern in cell lines and also predicted RBP motifs that were associated with different localization patterns. In addition, we showed that different isoforms of same genes could exhibit distinct localization patterns. To conclude, we believed our tool will facilitate future subcellular fractionation-related sequencing study to explore the function of RNA localization in various biological problems.