GISTA

April 23, 2025 ยท View on GitHub

A light-weight tool for finding Group-, Individual-sample Specific TADs (GISTA) GISTA is composed of three steps: 1). Partitioning the whole genome into a series of TADs arrays across all samples; 2). Constructing a feature vector for each TADs array; 3). Annotating the type of TAD change to each TADs arrays and statistical inference.

WorkFlow

Installation

Step1: install mamba from https://mamba.readthedocs.io/en/latest/installation/mamba-installation.html
Step2: mamba create -n gista python pandas numpy seaborn scipy tqdm logomaker click openpyxl xlsxwriter natsort
Step3: mamba activate gista
Step4: cd <path to the GISTA>
Step5: pip install --editable .

Quick Start

Currently GISTA only accept TopDom output formats (please check the exmaple data). If you need help, please reach out ot me :)

Multi-Samples Mode

# 'RT,PT;PT-RT,NT' mean process RT.vs.PT and PT+RT.vs.NT
# The comparison separate by ';' and treat and control separate by ',',The second sample is control
# see SampleList_Multi.xlsx exmaple in data
GISTA multi -sf SampleList_Multi.xlsx -c 'RT,PT;PT-RT,NT' --binsize 40000

Two-Samples Mode

# 'RT,PT;PT-RT,NT' mean process RT.vs.PT and PT+RT.vs.NT
# The comparison separate by ';' and treat and control separate by ',',The second sample is control
# see SampleList_Multi.xlsx exmaple in data
GISTA two -sf SampleList_Two.xlsx -c 'RT,PT;PT-RT,NT' --binsize 40000

Manual

Please use 'GISTA multi --help' or 'GISTA two --help' to check detailed instruction

Welcome to use GISTA :)
Usage: GISTA [OPTIONS] COMMAND [ARGS]...

Options:
  --version  Show the version and exit.
  --help     Show this message and exit.

Commands:
  multi  Multi-samples mode for GISTA
  two    Two-samples mode for GISTA

Multi

Usage: GISTA multi [OPTIONS]

  Multi-samples mode for GISTA

Options:
  -sf, --samplesfile PATH     sample metadata sheet, contain TopDom (like)
                              files, .csv or .xlsx  [required]
  -c, --comparison TEXT       Comparison string in the format 'RT,PT;PT-
                              RT,NT', The comparison separate by ';' and treat
                              and control separate by ',',The second sample is
                              control  [required]
  -bs, --binsize INTEGER      resolution/binsize of the TADs
  -gc, --groupcut FLOAT       Group level high variation cutoff
  -ic, --individualcut FLOAT  Individual level high variation cutoff
  -od, --outdir PATH          Output folder
  --help                      Show this message and exit.

Two

Usage: GISTA two [OPTIONS]

  Two-samples mode for GISTA

Options:
  -sf, --samplesfile PATH     sample metadata sheet, contain TopDom (like)
                              file, .csv or .xlsx  [required]
  -c, --comparison TEXT       Comparison string in the format 'RT,PT;PT-
                              RT,NT', The comparison separate by ';' and treat
                              and control separate by ',',The second sample is
                              control  [required]
  -bs, --binsize INTEGER      resolution/binsize of the TADs
  -gc, --groupcut FLOAT       Group level high variation cutoff
  -ic, --individualcut FLOAT  Individual level high variation cutoff
  -pr, --pseudorep INTEGER    The number of Pseudo-replication
  -od, --outdir PATH          Output folder
  --help                      Show this message and exit.

Citation

If you find this tool useful, please consider cite. Thank you!

Choppavarapu, L., Fang, K., Liu, T., Ohihoin, A. G., & Jin, V. X. (2025). Hi-C profiling in tissues reveals 3D chromatin-regulated breast tumor heterogeneity informing a looping-mediated therapeutic avenue. Cell reports, 44(4).