I would like to ask a methodological question regarding a specific scenario for somatic variant calling concerning whole exome sequencing data:
briefly, based on a biological protocol regarding the isolation of circulating tumor cells (CTCs) from patients with lung cancer, in vivo xenograft experiments were performed- CTCs isolated from specific patients were isolated and injected to mice, and approximately after 6 months, xenografts were removed, cut in small pieces, frozen and gDNA was extracted, which was subjected to WES. At the same time, peripheral blood cell samples (normal) are available from the same patients.
My main question is if there is a suggested pipeline that can be utilized to perform somatic variant calling on these PDX WES data ? using also as reference matched samples the normal PBMC samples from the patients ? I'm aware that one challenge is how to discriminate between mouse and human reads, as also any other infiltration due to stromal infiltration from mouse cells...
Thank you in advance,