hassansaei / VNtyper

A tool to genotype MUC1 coding VNTR in ADTKD
Other
3 stars 3 forks source link

VNtyper 2.0 - A Pipeline to genotype the MUC1-VNTR

VNtyper 2.0 is an advanced pipeline designed to genotype MUC1 coding Variable Number Tandem Repeats (VNTR) in Autosomal Dominant Tubulointerstitial Kidney Disease (ADTKD-MUC1) using Short-Read Sequencing (SRS) data. This version integrates enhanced variant calling algorithms, robust logging mechanisms, and streamlined installation processes to provide researchers with a powerful tool for VNTR analysis.


Table of Contents

  1. Features
  2. Installation
  3. Usage
  4. Pipeline Overview
  5. Dependencies
  6. Pipeline Logic Diagram
  7. Results
  8. Notes
  9. Citations
  10. Contributing
  11. License
  12. Contact

Features


Installation

VNtyper 2.0 can be installed using either pip with setup.py or via Conda environments for streamlined dependency management.

Using setup.py and pip

  1. Clone the Repository:

    mkdir vntyper
    git clone https://github.com/hassansaei/vntyper.git
    cd vntyper
    pip install .

Usage

VNtyper 2.0 offers multiple subcommands that can be used depending on your input data and requirements. Below are the main subcommands available:

1. Running the Full Pipeline

To run the entire pipeline on paired-end FASTQ files or BAM files:

vntyper pipeline \
    --config-path /path/to/config.json \
    --fastq1 /path/to/sample_R1.fastq.gz \
    --fastq2 /path/to/sample_R2.fastq.gz \
    --output-dir /path/to/output/dir \
    --threads 4

Alternatively, using a BAM file:

vntyper pipeline \
    --config-path /path/to/config.json \
    --bam /path/to/sample.bam \
    --output-dir /path/to/output/dir \
    --threads 4

2. Installing References

vntyper install-references \
    --output-dir /path/to/reference/install \
    --config-path /path/to/config.json \
    --skip-indexing  # Optional: skip BWA indexing if needed

3. Generating Reports

Generate a summary report for your VNTR genotyping analysis:

vntyper report \
    --output-dir /path/to/output/dir \
    --config-path /path/to/config.json

Process raw FASTQ files to prepare them for genotyping:

vntyper fastq \
    --fastq1 /path/to/sample_R1.fastq.gz \
    --fastq2 /path/to/sample_R2.fastq.gz \
    --output-dir /path/to/output/dir
vntyper bam \
    --alignment /path/to/sample.bam \
    --output-dir /path/to/output/dir \
    --threads 4

Pipeline Overview

VNtyper 2.0 integrates multiple steps into a streamlined pipeline. The following is an overview of the steps involved:

  1. FASTQ Quality Control: Raw FASTQ files are checked for quality.
  2. Alignment: Reads are aligned using BWA (if FASTQ files are provided).
  3. Kestrel Genotyping: Mapping-free genotyping of VNTRs.
  4. (Optional) adVNTR Genotyping: Profile-HMM based method for VNTR genotyping (requires additional setup).
  5. Summary Report Generation: A final HTML report is generated to summarize the results.

Dependencies

VNtyper 2.0 relies on several tools and Python libraries. Ensure that the following dependencies are available in your environment:

You can easily set up these dependencies via the provided Conda environment file.


Pipeline Logic Diagram

Below is a logical overview of the VNtyper pipeline:

graph TD
  A[Input: FASTQ/BAM] -->|Quality Control| B[Alignment BWA]
  B -->|Genotyping| C[Kestrel]
  C --> D[Optional: adVNTR]
  D --> E[Generate Summary Report]
  E --> F[Output: VCF, Summary HTML]

Notes

  1. This tool is for research use only.
  2. Ensure high-coverage WES data is used to genotype MUC1 VNTR accurately.
  3. For questions or issues, refer to the GitHub repository for support.

If you use VNtyper 2.0 in your research, please cite the following:

  1. Saei H, Morinière V, Heidet L, et al. VNtyper enables accurate alignment-free genotyping of MUC1 coding VNTR using short-read sequencing data. iScience. 2023.
  2. Audano PA, Ravishankar S, et al. Mapping-free variant calling using haplotype reconstruction from k-mer frequencies. Bioinformatics. 2018.
  3. Park J, Bakhtiari M, et al. Detecting tandem repeat variants in coding regions using code-adVNTR. iScience. 2022.

Contributing

We welcome contributions to VNtyper. Please refer to the CONTRIBUTING.md file for guidelines.


License

VNtyper is licensed under the BSD 3-Clause License. See the LICENSE file for more details.