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

Provide a rule to collapse expression of genes with the same key #6

Open
Feat-FeAR opened this issue Jun 17, 2024 · 1 comment
Open
Labels
bug Something isn't working enhancement New feature or request

Comments

@Feat-FeAR
Copy link
Member

In subsetGenes.xSeries, when, for a given key value (e.g., gene symbol), more than one occurrence is found in the annotation table of the xSeries (i.e., the same duplicated SYMBOL for multiple ENSG IDs) no action is taken and each entry is assigned to a different bar in the final chart, even if its name is not unique. It should be better to provide a rule to collapse statistics of genes with the same key value (e.g., keeping only the most expressed or their sum...), provided that user's "subsetting intention" is based on geneset key, rather than on (unique) ENSG IDs...
See the same problem in GOZER

DEGs <- DEGs[order(DEGs$adj_pval), ]
DEGs <- DEGs[!duplicated(DEGs$GeneSymbol), ]
@Feat-FeAR Feat-FeAR changed the title Provide a rule to collapse expression of genes with the same symbol Provide a rule to collapse expression of genes with the same key Jun 17, 2024
@Feat-FeAR Feat-FeAR added bug Something isn't working enhancement New feature or request labels Jun 17, 2024
@Feat-FeAR
Copy link
Member Author

Importantly, this really happens even within the restricted geneset of ion channels, as of org.Hs.eg.db v.3.18.0 (Bioconductor version 3.18 - BiocManager 1.30.23):

KCNMB2

ENSG00000197584 KCNMB2  potassium calcium-activated channel subfamily M regulatory beta subunit 2
ENSG00000275163 KCNMB2  potassium calcium-activated channel subfamily M regulatory beta subunit 2

AQP12

ENSG00000184945 AQP12A  aquaporin 12A   protein-coding
ENSG00000185176 AQP12A,AQP12B   aquaporin 12A,aquaporin 12B     protein-coding

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

1 participant