Questions Regarding Heterozygous Variants, Somatic Mutations, and Phasing in ClairS Usage

[sloth-eat-pudding]

Questions Regarding Heterozygous Variants, Somatic Mutations, and Phasing in ClairS Usage

[zhengzhenxian]

Hi, @sloth-eat-pudding,

Thank you for your interest in ClairS.

1. Our design is to select heterozygous SNPs from both normal and tumor samples for phasing and haplotagging. These signals are more likely to represent true germline heterozygous SNPs rather than somatic mutations.

2. We should only use heterozygous germline SNPs variants for phasing. In some cases, some somatic variants may have similar patterns to germline variants(might be due to the quality or high normal AF), which were identified as germline by Clair3. We are also actively working on excluding the rest of them from the phasing process to avoid confusing the model.

3. It would be highly beneficial to have a longer phaseset and an improved haplotagging ratio. We believe that having more phased alignments will significantly enhance performance. Any hints on parameter settings to have an improved phasing performance? Thanks!

4. Currently, we categorize the haplotypes into germline H1 and H2 only. However, it would be beneficial to include somatic haplotypes (H3) as well. However, linking distant somatic variants to obtain somatic haplotypes would be challenging. After analyzing the data, we have observed that the distance between two somatic variants can range from 10k to 100k, which presents a challenge in acquiring somatic haplotypes even in ONT reads.

Look forward to having a new LongPhase version for somatic variant calling!

Zhenxian

[sloth-eat-pudding]

I apologize for not being clear earlier.

Q1. I searched for confident heterozygous germline variants identified by ClairS in the normal.bam file, but found that they actually contain homozygous variants.(as shown in Figure 2).

Q3. This is our development goal.Additionally, I noticed that Clair3's Makefile uses LongPhase v1.3. Our current version is v1.6, and I suggest upgrading to this version for improvements in accuracy and processing time.

[zhengzhenxian]

For Q1, seems there are no homozygous variants in normal BAM in Figure 2, are you referring to a homozygous reference(that is the same allele as the reference base)? But thanks for reporting this, we will check the details then.

For Q3, thanks for the suggestion, sure, we will update LongPhase to v1.6 in our next release.

[sloth-eat-pudding]

Q1. Your understanding is correct. Thank you for your confirmation.