Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Why does the same transcription factor have many regulons? #535

Open
M-jk22 opened this issue Mar 22, 2024 · 4 comments
Open

Why does the same transcription factor have many regulons? #535

M-jk22 opened this issue Mar 22, 2024 · 4 comments
Labels
question Further information is requested

Comments

@M-jk22
Copy link

M-jk22 commented Mar 22, 2024

c417a7a2f527b0657af819afe6f0b74
1711077120810
In the analysis results, I found that the same transcription factor might correspond to many target gene sets, each with its own AUC. So, can I understand that even the same TF can have many regulons, and that the TF activity in each cell is the result of processing these regulon scores?

@M-jk22 M-jk22 added the question Further information is requested label Mar 22, 2024
@VincentGardeux
Copy link

Hi @M-jk22,

From what I understood these are not the regulons yet, but the modules that were computed on the adjacencies graph.
But then you have another step that summarizes this information at the TF level, thus creating one unique regulon per module.
I think you miss a step of the pipeline? It's like you've run the pruning step, but not the df2regulon step?

@M-jk22
Copy link
Author

M-jk22 commented Mar 22, 2024

Hi @VincentGardeux,
Thank you for your reply. Yes, this is the result after running pyscenic ctx. I thought that during the analysis process, one transcription factor would correspond to a group of target genes. However, from the current analysis results, it seems that after computing with the adjacency graph, one transcription factor is associated with many modules. Then, these modules are integrated using pyscenic aucell to obtain the AUCell score for the transcription factor regulon. I'm not sure if my description is correct. What does the AUC value of each module represent then? Can we use the AUC to filter these modules to select target genes?

@VincentGardeux
Copy link

VincentGardeux commented Mar 22, 2024

Hi @M-jk22,
I am not running the CLI tool, but the notebook version. I did not know that the pyscenic ctx was not returning the summarized regulons but instead the modules (I am no author of this tool, just a user).

There is a preprocessing step, before running the AUCell, where the modules are integrated into regulons, by concatenating all target genes of each module associated with a certain TF, and taking the top motif (higher NES). Then, the AUCell is computed on this summarized TF. So I assume that in the CLI version, the preprocessing is generated while running the AUCell step.

@M-jk22
Copy link
Author

M-jk22 commented Mar 22, 2024

Hi @VincentGardeux,
Thanks for your reply! I've learned a lot from chatting with you. Wishing you a happy and joyful life!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
question Further information is requested
Projects
None yet
Development

No branches or pull requests

2 participants