For TSSs that were differentially expressed, we performed enrichment tests for FANTOM5 cell specificity of TSSs. We downloaded the table of cell type specificity of CAGE tag clusters from http://fantom.gsc.riken.jp/5/datafiles/phase1.1/extra/Sample_ontology_enrichment_of_CAGE_peaks/. For each TSS, we noted whether it was differentially expressed and if it overlapped (based on genomic coordinates) with a FANTOM5 cell-specific tag cluster. We used the table function from base R to create a 2 × 2 contingency table for differential expression statistics and cell specificity, and tested for independence between rows and columns of this table using a Fisher’s Exact test (via R’s fisher.test). P-values for each cell category were corrected for multiple testing using the Benjamini–Hochberg method. The complete table of statistics for all cell type terms can be found in Supplementary Data 8.
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this article to respond.