Finding RNA genes associated with X protein expression in CITE-seq Seurat object

Question

[Picasa]

Hi everyone,

I’m working with a CITE-seq dataset using Seurat. In the ADT (protein) assay, I see that some cells express the protein marker X while others don’t.

What I’d like to do is find RNA genes (from the RNA assay) that are associated with the X protein expression levels

Has anyone done something similar or could suggest best practices or tools for correlating continuous ADT vs RNA ?

Thanks a lot for your help !

[ATpoint]

One way would be to make a correlation test between every pair of gene-ADT, and test for significance. Another would be to run something like WGCNA and filter for pairs that end up in the same module. Just thinking aloud, not done exactly this myself so far. Pseudobulk would probably be preferred if possible.