This repo contains code and information related to the processing of DNA and RNA data from patient-derived xenografts (PDX).
Simply ported over from the Seven Bridges Cancer Genomics Cloud (SB CGC), PDX WES Tumor-Only (Xenome) with Variant Calling, MSI, and TMB scores.
This workflow takes reads generated from a PDX model and performs variant calling as well as calculates microsatellite instability (MSI) and tumor mutation burden (TMB).
It was developed by Jackson Labs and implemented using common workflow language (CWL).
Open Targets group has ported over this workflow code into the Seven Bridges Cavatica environment, the source if which can be found in pdx-wes-tumor-only.cwl
.
A brief outline of what the workflow does, using fastq reads generated from the PDX:
- Run QC check on reads using FastQC
- Filter out mouse contaminants using xenome tool
- Align to reference (default hg38)
- Run alignment QC using picard metrics
- Call variants with Mutect2
- Calculate microsatellite instability (MSI)
- Calculate tumor mutation burden (TMB)
- Summarize results
Detailed documentation, including specific reference inputs, generated by the authors can be found here
Currently two workflows are present as we are comparing their relative efficacy.
As above, the first was developed by Jackson labs and ported over from SB CGC PDX RNA Expression Estimation Workflow.
Open Targets group has ported over this workflow code into the Seven Bridges Cavatica environment, the source if which can be found in pdx-rna-tumor-only.cwl
.
A brief outline of what the workflow does, using fastq reads generated from the PDX:
- Run QC check on reads using FastQC
- Filter out mouse contaminants using xenome tool
- Use RSEM to align reads (via STAR one-pass) to reference (default GRCh38, gene model ENSEMBL91) and quantify expression
- Run Somalier to generate inputs for relatedness calculations
Detailed documentation, including specific reference inputs, generated by the authors can be found here
The second workflow being compared is from the Kids First Data Resource Center (KFDRC)
This workflow has been modified to add xenome filtering up front, but the rest remains the same.
It has the added benefit of also performing fusion calling.
Open Targets group has ported over this workflow code into the Seven Bridges Cavatica environment, the source if which can be found in workflows/kfdrc-pdx-RNAseq-workflow.cwl
.
A brief outline of what the workflow does, using fastq reads generated from the PDX:
- Quantify expression with kallisto
- Align using STAR two pass against reference (default GENCODE 27, a GRCh38 reference)
- Call fusions using STAR-Fusion and arriba
- Collate and filter fusions using annofuse
- Additionally quantify expression using RSEM
Detailed documentation, including specific reference inputs, generated by the authors can be found here
PDX data from the paper "Genomic Profiling of Childhood Tumor Patient-Derived Xenograft Models to Enable Rational Clinical Trial Design", DOI: https://doi.org/10.1016/j.celrep.2019.09.071 are here as well as workflows used