Skip to content

RUN_GISTIC:ASSEMBLE_GISTIC_OUTPUT - KeyError: 'Unique Name' #33

Closed

Description

@pageale

Description of the bug

Hello, good afternoon !

I encountered the following error while trying to run GISTIC along with ichorCNA. The error message refers to a missing column named 'Unique Name', but I am unsure in which file this column is supposed to be present.

This issue appeared after I made the latest pull to version v1.0.6. Before this update, GISTIC was running without any problems. I also had the same issue with qpdf, as described in #32.

I am attaching the .nextflow.log, command.sh, and .command.run files for reference.

Thank you in advance for your help, and congratulations on your latest publication of SAMURAI! :D

executor >  local (2)
[71/96fd69] process > DINCALCILAB_SAMURAI:SAMURAI:SAMTOOLS_INDEX (16_11418_R)                   [100%] 6 of 6, cached: 6 ✔
[b3/ec30fc] process > DINCALCILAB_SAMURAI:SAMURAI:BAM_QC_PICARD:PICARD_COLLECTMULTIPLEMETRIC... [100%] 6 of 6, cached: 6 ✔
[-        ] process > DINCALCILAB_SAMURAI:SAMURAI:BAM_QC_PICARD:PICARD_COLLECTHSMETRICS         -
[e7/5de44c] process > DINCALCILAB_SAMURAI:SAMURAI:BAM_QC_PICARD:PICARD_COLLECTWGSMETRICS (VH... [100%] 6 of 6, cached: 6 ✔
[67/736e2a] process > DINCALCILAB_SAMURAI:SAMURAI:LIQUID_BIOPSY:HMMCOPY_READCOUNTER_ICHORCNA... [100%] 6 of 6, cached: 6 ✔
[5f/f979ef] process > DINCALCILAB_SAMURAI:SAMURAI:LIQUID_BIOPSY:RUN_ICHORCNA (M16_1190_R)       [100%] 6 of 6, cached: 6 ✔
[c4/ff4d01] process > DINCALCILAB_SAMURAI:SAMURAI:LIQUID_BIOPSY:AGGREGATE_ICHORCNA_TABLE (As... [100%] 1 of 1, cached: 1 ✔
[93/45a501] process > DINCALCILAB_SAMURAI:SAMURAI:LIQUID_BIOPSY:CORRECT_LOGR_ICHORCNA (Corre... [100%] 1 of 1, cached: 1 ✔
[62/bb08b1] process > DINCALCILAB_SAMURAI:SAMURAI:LIQUID_BIOPSY:CONCATENATE_BIN_PLOTS (Conca... [100%] 1 of 1, cached: 1 ✔
[d5/6721e1] process > DINCALCILAB_SAMURAI:SAMURAI:RUN_GISTIC:GISTIC2 (Run Gistic2)              [100%] 1 of 1, cached: 1 ✔
[52/316646] process > DINCALCILAB_SAMURAI:SAMURAI:RUN_GISTIC:ASSEMBLE_GISTIC_OUTPUT (Rearran... [  0%] 0 of 1
[fc/fcc57e] process > DINCALCILAB_SAMURAI:SAMURAI:RUN_GISTIC:MAFTOOLS (Maftools Plot)           [  0%] 0 of 1
[-        ] process > DINCALCILAB_SAMURAI:SAMURAI:MULTIQC 
-[dincalcilab/samurai] Pipeline completed with errors-
ERROR ~ Error executing process > 'DINCALCILAB_SAMURAI:SAMURAI:RUN_GISTIC:ASSEMBLE_GISTIC_OUTPUT (Rearrange GISTIC2 output)'

Caused by:
  Process `DINCALCILAB_SAMURAI:SAMURAI:RUN_GISTIC:ASSEMBLE_GISTIC_OUTPUT (Rearrange GISTIC2 output)` terminated with an error exit status (1)

Command executed:

  gistic lesion_table -t all --genome hg38 gistic_results gistic_summary_mqc
  mv gistic_summary_mqc.focal_lesions.txt gistic_summary_mqc.txt
  
  # Only do the rename if the file was produced, i.e. when broad analysis runs
  if [ -f gistic_summary_mqc.broad_lesions.txt ]; then
      mv gistic_summary_mqc.broad_lesions.txt gistic_broad_lesions_mqc.txt
  fi
  
  gistic gene_table --genome hg38 gistic_results gistic
  gistic copy_number_state_table --genome hg38 gistic_results gistic_cn_states.txt
  gistic log2ratio_table --genome hg38 gistic_results gistic_log2R.txt
  
  cat <<-END_VERSIONS > versions.yml
  "DINCALCILAB_SAMURAI:SAMURAI:RUN_GISTIC:ASSEMBLE_GISTIC_OUTPUT":
      gistic-cli: $(gistic --version | cut -d ' ' -f2)
  END_VERSIONS

Command used and terminal output

nextflow run dincalcilab/samurai --input short_ver_test_new_version_samurai.csv  -profile docker --outdir SAMURAI_SHORT_VER_RUN_FROM_BAM --genome hg38 --fasta ~/Reference/bwa_GRCh38/Homo_sapiens.GRCh38.dna.primary_assembly.fa --fai ~/Reference/bwa_GRCh38/Homo_sapiens.GRCh38.dna.primary_assembly.fa.fai --dict ~/Reference/bwa_GRCh38/Homo_sapiens.GRCh38.dna.primary_assembly.fa.dict --analysis_type liquid_biopsy --caller ichorcna --binsize 1000 --normal_panel ~/Projects/ichorcna_u/ichorCNA_U/src_ichorCNA/inst/extdata/HD_ULP_PoN_1Mb_median_normAutosome_mapScoreFiltered_median.rds --ichorcna_readcounter_chrs "1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,X,Y" --ichorcna_chrs_to_use "c(1:22)" --ichorcna_chrs_to_train "c(1:22)" --ichorcna_chrs_to_normalize "c(1:22)" --ichorcna_estimate_sc TRUE --ichorcna_gc_wig ~/Projects/ichorcna_u/ichorCNA_U/src_ichorCNA/inst/extdata/gc_hg38_1000kb.wig --ichorcna_map_wig ~/Projects/ichorcna_u/ichorCNA_U/src_ichorCNA/inst/extdata/map_hg38_1000kb.wig --ichorcna_centromere_file ~/Projects/ichorcna_u/ichorCNA_U/src_ichorCNA/inst/extdata/GRCh38.GCA_000001405.2_centromere_acen.txt   --ichorcna_normal_states "c(0.5,0.6,0.7,0.8,0.9)" --ichorcna_plotylim "c(-2, 2)" --max_memory 28.GB --compute_signatures TRUE --run_gistic TRUE --gistic_broad_analysis TRUE -resume

Relevant files

nextflow.log

command_files.zip

System information

N E X T F L O W ~ version 23.10.1
Container: Docker
OS: Linux
Version: dincalcilab/samurai v1.0.6-gb82e24f

Activity

  1. lbeltrame commented on Feb 4, 2025

    @lbeltrame
    Contributor

    Can you check the header output files from GISTIC directly (in particular all_lesions_conf_xx.txt where xx is your confidence level)? Unique Name should be a column in the GISTIC output: it may mean that no regions where identified.

  2. SaraPotente commented on Jul 2, 2026

    @SaraPotente
    Collaborator

    No response. Closed.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions