I am looking to compare concordance in gene expression by whole transcriptome RNAseq in primary tumor specimens and patient tumor derived organoids. My plan is as follows:

RNAseq counts will be normalized. Differential gene expression analysis will be performed relative to a reference normal specimen to determine a tumor enriched transcriptional profile.

For each organoid/tumor pair, we will plot normalized counts of each differentially expressed gene and calculate pearson correlation between organoid and tumor gene expression.

This will generate a correlation coefficient and p-value for each pair.

Beyond calculating summary statistics for the correlation output, is there any other statistical analysis that should be performed?

Any other approaches recommended?