We calculated the tissue specificity of gene abundance reflected by the tau score (τ)102 (ranging from 0 to 1, with 1 for highly tissue-specific genes and 0 for ubiquitously transcribed genes) for each gene with scaled TPM values. For each tissue, we averaged all replicates and then calculated τ to account for unequal numbers of replicates among tissues. We used τ ≥ 0.75 as the cut-off for tissue-specific genes.
We calculated the abundance distribution (i.e., transcriptome complexity) of distinct transcripts across tissues, reflected as the fraction of total RNAs contributed by the most highly expressed genes.
Differential gene expression analysis was performed using edgeR (version 3.22.5)103, with a false discovery rate (FDR) ≤0.05 and log2(fold change) ≥1 as cut-offs for statistical significance.
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.