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

Trans_susie mode outputs are cis interaction #179

Open
Hiroaki-S-Ri opened this issue Jan 23, 2025 · 0 comments
Open

Trans_susie mode outputs are cis interaction #179

Hiroaki-S-Ri opened this issue Jan 23, 2025 · 0 comments

Comments

@Hiroaki-S-Ri
Copy link

Hi,
I am trying Trans_susie mode anarysis.
Even if I try to use trans_susie mode, the result of cis interaction will be output.

Looking at the source code, in tensorqtl.py, a function called susie.map_loci is used in tensorqtl.py, and I check map_loci in susie.py, a function called map is used.
In this map function, a function called genotypeio.InputGeneratorCis is used.
I check InputGeneratorCis of genotypeio.py. In the InputGeneratorCis, common between the genotype chromosome and the phenotype chromosome are extracted, and also the snps are extracted with ±window. So I think the outputs become a cis interaction.

Actually, I did such as code.
python3 -m tensorqtl
test.ReplaceSampleID.qc.chr22 pheno.enhancer.chr13.bed.gz chr22SNPs_chr13enhancer.finemap_trans
--mode trans_susie
--susie_loci enhancer.chr13.loci.txt
--covariates cov.txt
--fdr 0.05

Traceback (most recent call last):
File "/opt/local/python/3.6.9/lib/python3.6/runpy.py", line 193, in _run_module_as_main
"main", mod_spec)
File "/opt/local/python/3.6.9/lib/python3.6/runpy.py", line 85, in _run_code
exec(code, run_globals)
File "/opt/local/python/3.6.9/lib/python3.6/site-packages/tensorqtl/main.py", line 2, in
tensorqtl.main()
File "/opt/local/python/3.6.9/lib/python3.6/site-packages/tensorqtl/tensorqtl.py", line 259, in main
maf_threshold=maf_threshold, max_iter=500, window=args.window, summary_only=False)
File "/agora/work2/sg/sasaki/eQTL/202308_enhancerQTL_normalizeValue=promoter+enhancermatrix3/tensorqtl_finemapping/susie.py", line 662, in map
raise ValueError('No valid phenotypes found.')

Since there is no common chromosome number, if igc.n_phenotypes == 0: is triggered, resulting in No valid phenotypes found.
If I add chr13 to the bfile, this error will no longer occur and the results will be output, but the chr13 enhancer vs chr22 SNP results are not in output.

tensorqtl version is 1.0.9, python version is 3.6.9.

I would appreciate it if you could check it.

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

No branches or pull requests

1 participant