Skip to main content

A collection of tools for genotype quality control and analysis

Project description

GenoTools

The preprint is out: https://www.biorxiv.org/content/10.1101/2024.03.26.586362v1.full.pdf

DOI

Documentation

You can find the full documentation with the following links:

Getting Started

GenoTools is a suite of automated genotype data processing steps written in Python. The core pipeline was built for Quality Control and Ancestry estimation of data in the Global Parkinson's Genetics Program (GP2)

To download the most current version from pip:

pip install the_real_genotools

Alternatively, if you'd like to download from github:

git clone https://github.com/dvitale199/GenoTools.git
cd GenoTools
pip install .

you can pull the most current references by running:

genotools-download

By default, the reference panel will be downloaded to ~/.genotools/ref. but can be download to a location of choice with --destination.

To download specific references/models, you can run the download with the following options:

genotools-download --ref 1kg_30x_hgdp_ashk_ref_panel --model nba_v1 --destination /path/to/download_directory/

Currently, 1kg_30x_hgdp_ashk_ref_panel is the only available reference panel. Available models are nba_v1 for the NeuroBooster array and neurochip_v1 for the NeuroChip Array and both are in GRCh38. If using a different array, we would suggest training a new model by running the standard command below. Please ensure the reference panel and your genotypes are in the same build. If you're using our reference panel, your genotypes must be in GRCh38.

Modify the paths in the following command to run the standard GP2 pipeline:

genotools \
  --pfile /path/to/genotypes/for/qc \
  --out /path/to/qc/output \
  --ancestry \
  --ref_panel /path/to/reference/panel \
  --ref_labels /path/to/reference/ancestry/labels \
  --all_sample \
  --all_variant

This will find common snps between your genotype data and the reference panel, run PCA, UMAP-transform PCs, and train a new XGBoost classifier specific to your data/ref panel.

if you'd like to run the pipeline using an existing model, you can do that like so (take note of the --model option):

genotools \
  --pfile /path/to/genotypes/for/qc \
  --out /path/to/qc/output \
  --ancestry \
  --ref_panel /path/to/reference/panel \
  --ref_labels /path/to/reference/ancestry/labels \
  --all_sample \
  --all_variant
  --model /path/to/nba_v1/model

if you'd like to run the pipeline using the default nba_v1 model in a Docker container, you can do that like so:

genotools \
  --pfile /path/to/genotypes/for/qc \
  --out /path/to/qc/output \
  --ancestry \
  --ref_panel /path/to/reference/panel \
  --ref_labels /path/to/reference/ancestry/labels \
  --container \
  --all_sample \
  --all_variant

Note: add the --singularity flag to run containerized ancestry predictions on HPC

genotools accept --pfile, --bfile, or --vcf. Any bfile or vcf will be converted to a pfile before running any steps.

Please consult the docs links listed at the top of the README for the full argument guide, function guide, Default pipeline overview, and guide for navigating the output JSON.

Acknowledgements

GenoTools was developed as the core genotype and wgs processing pipeline for the Global Parkinson's Genetics Program (GP2) at the Center for Alzheimer's and Related Dementias (CARD) at the National Institutes of Health.

This tool relies on PLINK, a whole genome association analysis toolset, for various genetic data processing functionalities. We gratefully acknowledge the developers of PLINK for their foundational contributions to the field of genetics. More about PLINK can be found at their website.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

the_real_genotools-1.3.0.tar.gz (1.1 MB view details)

Uploaded Source

Built Distribution

the_real_genotools-1.3.0-py3-none-any.whl (1.1 MB view details)

Uploaded Python 3

File details

Details for the file the_real_genotools-1.3.0.tar.gz.

File metadata

  • Download URL: the_real_genotools-1.3.0.tar.gz
  • Upload date:
  • Size: 1.1 MB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/5.1.1 CPython/3.11.3

File hashes

Hashes for the_real_genotools-1.3.0.tar.gz
Algorithm Hash digest
SHA256 5c7678d43ba14f7c3439fdf5cb5261340fb253670650d436ae4e6bec8f4e0154
MD5 fbc8cd8f4922826e1b8835c78f50daf5
BLAKE2b-256 ff45e12a43b1a2b14299bc296c3e22507e7236c13cbb6cb966760cf273a5e3ca

See more details on using hashes here.

File details

Details for the file the_real_genotools-1.3.0-py3-none-any.whl.

File metadata

File hashes

Hashes for the_real_genotools-1.3.0-py3-none-any.whl
Algorithm Hash digest
SHA256 ab9d0a402663ed8c60a6c84c50ae2757bbb797e53473d2136da255c75e812748
MD5 bc8d4ce64ac8abdeab2a3a80d9f0509c
BLAKE2b-256 1711331137a4417b11ae4060136daff03fe93f8abc809933ab2db932ad85f540

See more details on using hashes here.

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page