Hi, I'm working on scATAC-seq motif analysis using Signac, and I encountered an issue related to motif name redundancy.
When running motif enrichment test by FindMotifs() function, I observe that the motif.name column in the results contains many variations of AP-1 family motifs, such as FOSL1, FOSL1::JUND, FOSL2::JUN, BATF::JUN, FOS::JUN, JUNB, ....
While this reflects the structure of the JASPAR PWM database, it makes the plot difficult to interpret since dozens of highly similar AP-1 motifs dominate the top-ranked results.
I'd like to ask:
- Is there a recommended way to collapse redundant motifs? for example by consolidating them under a single representative gene name?
- Does Signac provide any built-in method to map motif names to gene names?
Thank you in advance for your help! :)