0
$\begingroup$

The attached figure is generated based on Illumina reads from multiple individuals compared to genome assembly. Looks like there are a lot of kmers are reads only (grey colored). Also, a blue peak (2x kmers), does it represents high heterozygosity? GenomeScope gives 1.5% of heterozygosity result. But it seems the blue peak normally should be near the 70x on the x-axisenter image description here. Thanks.

$\endgroup$

1 Answer 1

1
$\begingroup$

Assuming a diploid genome, usually, we aim to build a monoploid reference, representing a single haplotype.

Practically speaking, we need to include only one of the two k-mers that are representing heterozygous sites, however, we would desire all the homozygous k-mers to be included exactly once.

Looking at your plot, I would guess 1n coverage is somewhere around 18x, and 2n about 36x. If that's right, your plot is suggesting that your assembly contains lots of uncollapsed haplotypes, which is quite understandable given very high heterozygosity you estimated with genomescope. This is also supported by many homozygous k-mers are there twice (the blue bump). These are probably regions in between heterozygous sites, or flanking them.

Anyway, I would recommend something like purge haplotigs to collapse haplotypes into single monoploid reference and redo the plot. The blue bump should disappear and the black and red distributions should have about the same height around 18x, while nearly all the 36x k-mers should be red (once).

Also, if the interpretation of the plot still would not be clear, you can check this part of a lecture by Lucía Campos-Dominguez on the topic (it's a different software called KAT she talks about, but it is using the same principle as Mercury).

$\endgroup$
2
  • $\begingroup$ Hi Kamil, thank you for your informative reply. You are exactly right about the 1n = 18, although not sure how you get that so accurate ;-). I do have another Merqury kmer plot with a very low blue peak near the 4n location rather than the 2n location, while the red peak at the same 2n. And this happens for another species genome assembly. Is that normal? Because that makes the blue peak at 2n look weired. $\endgroup$ Dec 14, 2021 at 23:26
  • $\begingroup$ I will post another question. Thank you! $\endgroup$ Dec 15, 2021 at 4:26

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.