somatic_SNVs_pipeline/workflow_SEQC2.wdl

import "./tasks/BQSR.wdl" as Sentieon_BQSR
import "./tasks/TNseq.wdl" as TNseq
import "./tasks/TNscope.wdl" as TNscope
import "./tasks/manta.wdl" as manta_calling
import "./tasks/strelka.wdl" as strelka_calling
import "./tasks/bcftools.wdl" as bcftools
import "./tasks/bcftools_concat.wdl" as bcftools_concat
import "./tasks/annovar.wdl" as annovar
import "./tasks/Haplotyper.wdl" as Haplotyper
import "./tasks/MSIsensor.wdl" as MSIsensor
import "./tasks/tmb.wdl" as tmb


workflow {{ project_name }} {
	#sample info
  File? normal_deduped_bam
  File? normal_deduped_bam_bai
	File? tumor_deduped_bam
	File? tumor_deduped_bam_bai
  String sample_id
	String Seq_platform
	#reference fasta
	String ref_fasta
	File ref_fasta_dir
	File annot_gff
  File dbsnp_dir
  File dbmills_dir
  String dbsnp
  String db_mills
	File germline_resource
	File germline_resource_tbi
	File annovar_database
  File? regions
  File baseline
  File hg38_CDS

	#docker
  String docker_fastp
	String docker_sentieon
	String SENTIEON_LICENSE
	String docker_qualimap
	String docker_manta
	String docker_strelka
	String docker_bcftools
	String docker_annovar
  String docker_tmb
  String docker_MSIsensor

	#config
  String cluster_config
  String disk_size
	Boolean? fastqc
	Boolean bamqc
	Boolean Annovar
  Boolean AF_filter

#SEQC2应该是_N_1和T_1
if (normal_deduped_bam != ""){

		call Sentieon_BQSR.Sentieon_BQSR as Sentieon_BQSR_normal{
			input:
				ref_dir=ref_fasta_dir,
				dbsnp_dir=dbsnp_dir,
				dbmills_dir = dbmills_dir,
				sample_id=sample_id+'_N',
				ref_fasta=ref_fasta,
				dbsnp=dbsnp,
				db_mills=db_mills,
				deduped_bam=normal_deduped_bam,
				deduped_bam_index=normal_deduped_bam_bai,
				docker=docker_sentieon,
				cluster_config=cluster_config,
				disk_size=disk_size,
				SENTIEON_LICENSE=SENTIEON_LICENSE
		}

	}

if (tumor_deduped_bam != ""){

		call Sentieon_BQSR.Sentieon_BQSR as Sentieon_BQSR_tumor{
			input:
				ref_dir=ref_fasta_dir,
				dbsnp_dir=dbsnp_dir,
				dbmills_dir = dbmills_dir,
				sample_id=sample_id+'_T',
				ref_fasta=ref_fasta,
				dbsnp=dbsnp,
				db_mills=db_mills,
				deduped_bam=tumor_deduped_bam,
				deduped_bam_index=tumor_deduped_bam_bai,
				docker=docker_sentieon,
				cluster_config=cluster_config,
				disk_size=disk_size,
				SENTIEON_LICENSE=SENTIEON_LICENSE
		}

	}

call MSIsensor.MSIsensor as MSIsensor {
    input:
    sample=sample_id,
    fasta=ref_fasta,
    ref_dir=ref_fasta_dir,
    normal_bam=Sentieon_BQSR_normal.recaled_bam,
    normal_bam_index=Sentieon_BQSR_normal.recaled_bam_index,
    tumor_bam=Sentieon_BQSR_tumor.recaled_bam,
    tumor_bam_index=Sentieon_BQSR_tumor.recaled_bam_index,
    baseline=baseline,
    docker=docker_MSIsensor,
    cluster_config="OnDemand bcs.ps.g.large img-ubuntu-vpc",
    disk_size=disk_size
  }


call Haplotyper.Haplotyper as Haplotyper {
      input:
        SENTIEON_LICENSE=SENTIEON_LICENSE,
        sample=sample_id + '_N',
        fasta=ref_fasta,
        ref_dir=ref_fasta_dir,
        recaled_bam=Sentieon_BQSR_normal.recaled_bam,
        recaled_bam_index=Sentieon_BQSR_normal.recaled_bam_index,
        dbsnp=dbsnp,
        dbsnp_dir=dbsnp_dir,
        regions=regions,
        docker=docker_sentieon,
        disk_size=disk_size,
        cluster_config=cluster_config
    }

call bcftools.bcftools as Haplotyper_bcftools {
      input:
        vcf=Haplotyper.vcf,
        fasta=ref_fasta,
        ref_dir=ref_fasta_dir,
        docker=docker_bcftools,
        cluster_config="OnDemand bcs.ps.g.large img-ubuntu-vpc",
        disk_size=disk_size
    }

call TNseq.sentieon_TNseq as sentieon_TNseq{
	input:
		sample_id = sample_id,
		tumor_bam = Sentieon_BQSR_tumor.recaled_bam,
		tumor_bam_bai = Sentieon_BQSR_tumor.recaled_bam_index,
		normal_bam = Sentieon_BQSR_normal.recaled_bam,
		normal_bam_bai = Sentieon_BQSR_normal.recaled_bam_index,
		tumor_name = sample_id+'_T_1',
		normal_name = sample_id+'_N_1',

		ref_dir = ref_fasta_dir,
		ref_fasta = ref_fasta,
		germline_resource = germline_resource,
		germline_resource_tbi = germline_resource_tbi,

		docker = docker_sentieon,
		cluster_config = cluster_config,
		disk_size = disk_size,
		SENTIEON_LICENSE = SENTIEON_LICENSE
}



call bcftools.bcftools as TNseq_bcftools {
	input:
	vcf=sentieon_TNseq.vcf,
	fasta=ref_fasta,
	ref_dir=ref_fasta_dir,
	docker=docker_bcftools,
	cluster_config="OnDemand bcs.ps.g.large img-ubuntu-vpc",
	disk_size=disk_size
}

call TNscope.sentieon_TNscope as sentieon_TNscope{
	input:
		sample_id = sample_id,
		tumor_bam = tumor_deduped_bam,
		tumor_bam_bai = tumor_deduped_bam_bai,
		normal_bam = normal_deduped_bam,
		normal_bam_bai = normal_deduped_bam_bai,
		tumor_name = sample_id+'_T_1',
		normal_name = sample_id+'_N_1',
		tumor_recall_data = Sentieon_BQSR_tumor.recal_table,
		normal_recall_data = Sentieon_BQSR_normal.recal_table,

		ref_dir = ref_fasta_dir,
		ref_fasta = ref_fasta,
		dbsnp_dir = dbsnp_dir,
		dbsnp = dbsnp,
		# excute env
		docker = docker_sentieon,
		cluster_config = cluster_config,
		disk_size = disk_size,
		SENTIEON_LICENSE = SENTIEON_LICENSE
}

call bcftools.bcftools as TNscope_bcftools {
	input:
	vcf=sentieon_TNscope.vcf,
	fasta=ref_fasta,
	ref_dir=ref_fasta_dir,
	docker=docker_bcftools,
	cluster_config="OnDemand bcs.ps.g.large img-ubuntu-vpc",
	disk_size=disk_size
}


call strelka_calling.strelka_calling as strelka_calling{
	input:
		tumor_bam = Sentieon_BQSR_tumor.recaled_bam,
		tumor_bam_bai = Sentieon_BQSR_tumor.recaled_bam_index,
		normal_bam = Sentieon_BQSR_normal.recaled_bam,
		normal_bam_bai = Sentieon_BQSR_normal.recaled_bam_index,
		ref_fasta = ref_fasta,
		ref_dir = ref_fasta_dir,
		sample_id = sample_id,
		
		docker=docker_strelka,
		cluster_config="OnDemand bcs.ps.g.6xlarge img-ubuntu-vpc",
		disk_size=disk_size
}

call bcftools_concat.bcftools as bcftools_concat{
        input:
            ref_dir=ref_fasta_dir,
            fasta=ref_fasta,
            vcf_indels=strelka_calling.indel_vcf,
            vcf_snvs=strelka_calling.snv_vcf,
            sample_id=sample_id,
            docker=docker_bcftools,
            cluster_config="OnDemand bcs.ps.g.large img-ubuntu-vpc",
            disk_size=disk_size

    }

if (Annovar){
	call annovar.ANNOVAR as TNseq_ANNOVAR {
		input:
			vcf=TNseq_bcftools.norm_vcf,
			annovar_database=annovar_database,
			docker=docker_annovar,
			cluster_config="OnDemand bcs.ps.g.6xlarge img-ubuntu-vpc",
			disk_size=disk_size
  }

    call tmb.TMB as TMB{
        input:
            sample=sample_id,
            hg38_CDS=hg38_CDS,
            regions=regions,
            snpindel_txt=TNseq_ANNOVAR.multianno_txt,
            docker=docker_tmb,
            cluster_config="OnDemand bcs.ps.g.large img-ubuntu-vpc",
            disk_size=disk_size,
            AF_filter=AF_filter
    }


}

}