Bioinformatics Stack Exchange is a question and answer site for researchers, developers, students, teachers, and end users interested in bioinformatics. It only takes a minute to sign up.
Sign up to join this community
I'm evaluating possibilites for somatic tumor variant calling without paired-normal samples. I'm aware of the consequences without a normal sample.
All the popular variant calling tools such as Strelka, VarScan etc require a normal sample.
What are the some of the possibilites I have? I'm looking for good sensitivity in exchange for higher false positives.
There is a recent paper that attempts to do this:
ISOWN: accurate somatic mutation identification in the absence of normal tissue controls.
In this work, we describe the development, implementation, and validation of ISOWN, an accurate algorithm for predicting somatic mutations in cancer tissues in the absence of matching normal tissues. ISOWN is available as Open Source under Apache License 2.0 from https://github.com/ikalatskaya/ISOWN
I've been using the LoFreq* caller for exactly this. It is designed to find variants with very low frequency, so is well suited for this type of analysis.
LoFreq* (i.e. LoFreq version 2) is a fast and sensitive variant-caller for inferring SNVs and indels from next-generation sequencing data. It makes full use of base-call qualities and other sources of errors inherent in sequencing (e.g. mapping or base/indel alignment uncertainty), which are usually ignored by other methods or only used for filtering.
LoFreq* is very sensitive; most notably, it is able to predict variants below the average base-call quality (i.e. sequencing error rate). Each variant call is assigned a p-value which allows for rigorous false positive control. Even though it uses no approximations or heuristics, it is very efficient due to several runtime optimizations and also provides a (pseudo-)parallel implementation.
(In my experience it is not fast at all, but that's what they claim).
I have used it many times to call somatic variants using the germline caller (and not its somatic one which requires a tumor/normal pair). Because the tool is so sensitive, its germline caller actually deals quite well with somatic samples.
You can find instructions for its use here.
All the popular variant calling tools such as Strelka, VarScan etc require a normal sample.
Strelka and VarScan require a normal sample in somatic mode, but they both have a germline mode for unpaired analysis.
One of the most popular somatic variant callers is MuTect and that does not require a normal.
