A study of TP53 RNA splicing illustrates pitfalls of RNA-seq methodology

Abstract

TP53 undergoes multiple RNA-splicing events, resulting in at least nine mRNA transcripts encoding at least 12 functionally different protein isoforms. Antibodies specific to p53 protein isoforms have proven difficult to develop, thus researchers must rely on the transcript information to infer isoform abundance. In this study, we used deep RNA-seq, droplet digital PCR (ddPCR), and real-time quantitative reverse transcriptase PCR (RT-qPCR) from nine human cell lines and RNA-seq data available for tumors in The Cancer Genome Atlas to analyze TP53 splice variant expression. All three methods detected expression of the FL/40TP53α-T1 variant in most human tumors and cell lines. However, other less abundant variants were only detected with PCR-based methods. Using RNA-seq simulation analysis, we determined why RNA-seq is unable to detect less abundant TP53 transcripts and discuss the implications of these findings for the general interpretation of RNA-seq data.