{"id":20074711,"url":"https://github.com/greenelab/pdx_exomeseq","last_synced_at":"2026-01-31T07:02:55.845Z","repository":{"id":79359775,"uuid":"66270694","full_name":"greenelab/pdx_exomeseq","owner":"greenelab","description":"Pipeline analysis for whole exome sequencing of pancreatic cancer PDX models","archived":false,"fork":false,"pushed_at":"2018-10-22T22:04:30.000Z","size":139608,"stargazers_count":22,"open_issues_count":3,"forks_count":15,"subscribers_count":5,"default_branch":"master","last_synced_at":"2025-06-22T00:06:31.084Z","etag":null,"topics":["analysis","cancer","methodology","mutation","notebook","pancreatic-cancer","pdx","pipeline","sequencing","wes"],"latest_commit_sha":null,"homepage":"","language":"HTML","has_issues":true,"has_wiki":null,"has_pages":null,"mirror_url":null,"source_name":null,"license":"mit","status":null,"scm":"git","pull_requests_enabled":true,"icon_url":"https://github.com/greenelab.png","metadata":{"files":{"readme":"README.md","changelog":null,"contributing":null,"funding":null,"license":"LICENSE.md","code_of_conduct":null,"threat_model":null,"audit":null,"citation":null,"codeowners":null,"security":null,"support":null,"governance":null}},"created_at":"2016-08-22T12:26:39.000Z","updated_at":"2025-03-10T06:17:41.000Z","dependencies_parsed_at":"2023-03-12T07:49:50.578Z","dependency_job_id":null,"html_url":"https://github.com/greenelab/pdx_exomeseq","commit_stats":null,"previous_names":[],"tags_count":0,"template":false,"template_full_name":null,"purl":"pkg:github/greenelab/pdx_exomeseq","repository_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/greenelab%2Fpdx_exomeseq","tags_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/greenelab%2Fpdx_exomeseq/tags","releases_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/greenelab%2Fpdx_exomeseq/releases","manifests_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/greenelab%2Fpdx_exomeseq/manifests","owner_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners/greenelab","download_url":"https://codeload.github.com/greenelab/pdx_exomeseq/tar.gz/refs/heads/master","sbom_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/greenelab%2Fpdx_exomeseq/sbom","scorecard":null,"host":{"name":"GitHub","url":"https://github.com","kind":"github","repositories_count":286080680,"owners_count":28932594,"icon_url":"https://github.com/github.png","version":null,"created_at":"2022-05-30T11:31:42.601Z","updated_at":"2026-01-31T04:05:25.756Z","status":"ssl_error","status_checked_at":"2026-01-31T04:02:35.005Z","response_time":128,"last_error":"SSL_read: unexpected eof while reading","robots_txt_status":"success","robots_txt_updated_at":"2025-07-24T06:49:26.215Z","robots_txt_url":"https://github.com/robots.txt","online":false,"can_crawl_api":true,"host_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub","repositories_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories","repository_names_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repository_names","owners_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners"}},"keywords":["analysis","cancer","methodology","mutation","notebook","pancreatic-cancer","pdx","pipeline","sequencing","wes"],"created_at":"2024-11-13T14:53:51.437Z","updated_at":"2026-01-31T07:02:55.831Z","avatar_url":"https://github.com/greenelab.png","language":"HTML","funding_links":[],"categories":[],"sub_categories":[],"readme":"# Whole Exome Sequencing Pipeline for JAX FNA-PDX models of Pancreatic Cancer\n\n**Gregory Way\u003csup\u003e1\u003c/sup\u003e, Casey Greene\u003csup\u003e1\u003c/sup\u003e, Yolanda Sanchez\u003csup\u003e2\u003c/sup\u003e**\n\n1. University of Pennsylvania\n2. Geisel School of Medicine at Dartmouth\n\n## Summary\n\nPatient derived xenograft (PDX) models were derived from primary and metastatic tumors from patients admitted to Dartmouth-Hitchcock Medical Center (DHMC) with pancreatic adenocarcinoma (PAAD).\nThe PDX models and tumor samples were whole exome sequenced (WES) to determine how the mutations from primary tissue and metastases propagate and evolve.\nThe following repository outlines the wes and analysis pipelines.\n\nThis is a tumor-only analysis; there were no pooled or patient-matched normal samples available.\nThe following flowchart summarizes the wes pipeline.\n\n![pdx wes flowchart](figures/pdx_wes_flowchart.png?raw=true)\n\n**Figure 1A** describes the technical replicates and data-types available across tumor and mouse passages.\n**Figure 1B** outlines our whole exome sequencing pipeline. We first apply quality control processing to raw reads, then align and remove mouse reads, and finally call and annotate variants.\n\n## WES Pipeline\n\nSee [`wes_pipeline.sh`](https://github.com/greenelab/pdx_exomeseq/blob/master/wes_pipeline.sh) for our current variant-calling pipeline for tumor-only WES.\nThis script was run step-by-step on the Dartmouth Discovery compute cluster.\n\n### WES Compute Environment\n\nAll work was performed using the Dartmouth Discovery Cluster Computer with the conda environment specified in [`environment.yml`](https://github.com/greenelab/pdx_exomeseq/blob/master/environment.yml).\n\n### Steps to Reproduce\n\nThere are 3 major steps this repository provides to get from raw sequencing reads to annotated variants.\n\n#### 1. Setup reproducible computational environment ([`setup_environment.sh`](https://github.com/greenelab/pdx_exomeseq/blob/master/setup_environment.sh), [`install.sh`](https://github.com/greenelab/pdx_exomeseq/blob/master/install.sh))\n\n```bash\n# Setup conda (version 4.5 or greater) environment\nbash setup_environment.sh\n\n# NOTE: run `conda activate pdx-exomeseq` at the beginning of each session\n\n# Install dependencies and initialize files\n# This includes downloading reference genomes and generating several index files\nbash install.sh\n```\n\n#### 2. Run data processing pipeline ([`wes_pipeline.sh`](https://github.com/greenelab/pdx_exomeseq/blob/master/wes_pipeline.sh))\n\n```bash\n# NOTE: the commands in the following script must be run sequentially\n# The script will submit several jobs per specified file that can take upwards of\n# 12 hours per sample to run _for each command_. This requires the user to specify\n# which command is being run by commenting out all others.\nbash wes_pipeline.sh\n```\n\nAlso note that the configuration file `discovery_variables.yml` includes absolute paths to each tool or resource.\nIt is sufficient to update this file only if paths to current tools change.\n\n#### 3. Visualize and summarize results ([`analysis_pipeline.sh`](https://github.com/greenelab/pdx_exomeseq/blob/master/analysis_pipeline.sh))\n\nWe use Jupyter notebooks and R scripts to visualize and summarize results.\nWe describe the analysis in the next section.\n\n## Analysis Pipeline\n\nAfter obtaining the called variants, we perform a series of analyses and visualizations.\nThese analyses use a separate conda environment which is specified in\n[`analysis_environment.yml`](https://github.com/greenelab/pdx_exomeseq/blob/master/analysis_environment.yml).\n\n### Computational Environment\n\nFollow these steps to install and begin using this conda environment:\n\n```bash\n# Using conda version 4.5 or greater\nconda env create --force --file analysis_environment.yml\nconda activate pdx-exomeseq-analysis\n```\n\n### Reproduce Results\n\nIn order to reproduce the results of the analysis pipeline perform the following steps.\n(Note that the variants are expected to be processed before running the pipeline)\n\n```bash\nbash analysis_pipeline.sh\n```\n\n### Scripts\n\nThe following notebooks perform the analysis and obtain figures and results:\n\n| Script | Output |\n| :----- | :----- |\n| [`1.read-depth-stats.ipynb`](https://github.com/greenelab/pdx_exomeseq/blob/master/1.read-depth-stats.ipynb) | Determine read depth against proportion of genome covered |\n| [`2.disambiguate-reads.ipynb`](https://github.com/greenelab/pdx_exomeseq/blob/master/2.disambiguate-reads.ipynb) | Visualizing the separation of mouse and human reads |\n| [`3.filter-variants.ipynb`](https://github.com/greenelab/pdx_exomeseq/blob/master/3.filter-variants.ipynb) | Visualize variant filtration and process filtered VCFs |\n| [`4.variant-allele-frequency.ipynb`](https://github.com/greenelab/pdx_exomeseq/blob/master/4.variant-allele-frequency.ipynb) | visualize gnomAD by SIFT scores for replicates and filtered merged files |\n|[`5.upset-plots.ipynb`](https://github.com/greenelab/pdx_exomeseq/blob/master/5.upset-plots.ipynb) | Generate UpSet plots to visualize variant overlaps across patient sets |\n|[`6.generate-oncoprint-data.ipynb`](https://github.com/greenelab/pdx_exomeseq/blob/master/5.generate-oncoprint-data.ipynb) | Wrangle variant calls to generate data for input into oncoprint visualization |\n| [`7.visualize-oncoprint.ipynb`](https://github.com/greenelab/pdx_exomeseq/blob/master/6.visualize-oncoprint.ipynb) | Visualize oncoprint diagrams and variant similarity matrices |\n","project_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fgreenelab%2Fpdx_exomeseq","html_url":"https://awesome.ecosyste.ms/projects/github.com%2Fgreenelab%2Fpdx_exomeseq","lists_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fgreenelab%2Fpdx_exomeseq/lists"}