I have been following the section of Gene set enrichment and pathway analysis in single-cell best practices. Specifically, I am trying to use decoupler.run_gsea to perform pathway enrichment with scRNA-seq data. The results look very good and consistent with our understanding of the data. But I wish I can include the leading edge genes in the output somehow? I searched around but could not really figure out how to achieve that, at least using decoupler. Maybe GSEApy is better route for this? Any suggestions would be very much appreciated. Sorry if I missed anything obvious.