Identification of positive regulators per condition

[DaianeH]

Hi,

I have a scATAC-seq dataset from samples with TP53 Mut and TP53 WT cells. I am interested in finding the difference between the two groups of samples. When making the archR object, I used a scRNA-seq dataset to transfer cell-type labels and create a GeneIntegrationMatrix. The scRNA-seq and the scATAC-seq datasets were generated on the same samples, but not on the exact same cells (different barcodes).

I would like to identify positive regulators of the TP53 Mut group of samples, and positive regulators of the TP53 WT group of samples, to see if the different conditions have distinct positive regulators.

How can I do that? Do I need to re-run the entire workflow twice, creating an archR project for each condition, separately, or can I do that on my current archR object that contains all samples (both Mut and WT), and somehow define that I'm interested in the positive regulators of each condition?

Thank you so much in advance!