SuShiE (Sum of Shared Single Effect) is a Python package for multiancestry SNP fine-mapping, estimating effect size correlations across ancestries, and computing ancestry-specific prediction weights using either individual-level or summary-level data for molecular or complex traits.

- We detest usage of our software or scientific outcome to promote racial discrimination.

SuShiE is described in

Improved multiancestry fine-mapping identifies cis-regulatory variants underlying molecular traits and disease risk.

Zeyun Lu, Xinran Wang, Matthew Carr, Artem Kim, Steven Gazal, Pejman Mohammadi, Lang Wu, James Pirruccello, Linda Kachuri, Alexander Gusev, Nicholas Mancuso.

Nature Genetics. July, 2025.

Check here for full documentation.

Installation | Example | Notes | Version History | Support | Other Software

1. Before installation, we highly recommend to create a new environment using conda so that it will not affect the software versions of the other projects. For example, use following codes:

   conda create -n env-sushie python=3.8

   We currently only support Python3.8+.

2. If you are using a Mac with an Apple M1 or newer chip, you should install cbgen package or other required pacakges from conda-forge first to ensure compatibility (see this link for previous issue). One easy workaround is to initiate your conda using miniforge. On most HPC systems, this is usually not necessary.

   conda install -c conda-forge cbgen

3. Last, users can download the latest repository and then use pip:

   git clone https://github.com/mancusolab/sushie.git
   cd sushie
   pip install .

SuShiE software is very easy to use:

For fine-mapping using individual-level data:

cd ./data/
sushie finemap --pheno EUR.pheno AFR.pheno --vcf vcf/EUR.vcf vcf/AFR.vcf --covar EUR.covar AFR.covar --output ./test_result

For fine-mapping using summary-level data:

cd ./data/
sushie finemap --summary --gwas EUR.gwas AFR.gwas --vcf vcf/EUR.vcf vcf/AFR.vcf --sample-size 489 639 --gwas-header chrom snp pos a1 a0 z --output ./test_result

It can perform:

• SuShiE: multi-ancestry fine-mapping accounting for ancestral correlation
• Single-ancestry SuSiE (Sum of Single Effect)
• Independent SuShiE: multi-ancestry SuShiE without accounting for correlation
• Meta-SuSiE: single-ancestry SuSiE followed by meta-analysis
• Mega-SuSiE: single-ancestry SuSiE on row-wise stacked data across ancestries (individual-level data only)
• cis-molQTL effect size correlation estimation
• cis-SNP heritability estimation (individual-level data only)
• Cross-validation for SuShiE prediction weights (individual-level data only)
• Convert prediction results to FUSION format, thus can be used in TWAS

See here for more details on how to use SuShiE.

If you want to use in-software SuShiE inference function, you can use following Python code as an example:

from sushie.infer import infer_sushie
# Xs is for genotype data, and it should be a list of numpy array whose length is the number of ancestry.
# ys is for phenotype data, and it should also be a list of numpy array whose length is the number of ancestry.
infer_sushie(Xs=X, ys=y)
# Or summary-level data
# lds is for LD data, and it should be a list of p by p numpy array whose length is the number of ancestry.
# zs is for GWAS data, and it should be a list of numpy array whose length is the number of ancestry/
infer_sushie_ss(lds=LD, zs=GWAS, ns=np.array([100, 100]))

You can customize this function with your own ideas!

• SuShiE currently only supports continuous phenotype fine-mapping for individual-level data.

For any questions, comments, bug reporting, and feature requests, please contact Zeyun Lu ([email protected]) and Nicholas Mancuso ([email protected]), and open a new thread in the Issue Tracker.

Feel free to use other software developed by Mancuso Lab:

• jaxQTL: a single-cell eQTL mapping tool using highly efficient count-based model (i.e., negative binomial or Poisson).
• MA-FOCUS: a Bayesian fine-mapping framework using TWAS statistics across multiple ancestries to identify the causal genes for complex traits.
• SuSiE-PCA: a scalable Bayesian variable selection technique for sparse principal component analysis
• twas_sim: a Python software to simulate TWAS statistics.
• FactorGo: a scalable variational factor analysis model that learns pleiotropic factors from GWAS summary statistics.
• HAMSTA: a Python software to estimate heritability explained by local ancestry data from admixture mapping summary statistics.
• Traceax: a Python library to perform stochastic trace estimation for linear operators.

This project has been set up using PyScaffold 4.1.1. For details and usage information on PyScaffold see https://pyscaffold.org/.