2
$\begingroup$

I have what at first sight appears to be a high-quality MAG (~10 pieces, high completion%) that I built from a hybrid assembly (Illumina + Nanopore data) from a cyanobacterium.

Workflow: Quality control (BBDuk) > Hybrid assembly (Unicycler) > co-binning multiple different samples (vamb) > Quality control (CheckM, anvi'o)

*The Nanopore data is pretty old, and therefore I'd thought Unicycler might be a safer bet rather than long-read-first assemblers, due to the higher error rate

The same strain also has 1-2 very fragmented genomes (100s of pieces) available in NCBI at scaffold level. The reads themselves are not available, so I can't just add them to the assembly.

Is there a way to use the available NCBI genome(s) to try to check/polish mine? (I've checked out Improve a reference genome with sequencing data but it's already a few years old so I imagine different tools might be applicable)

Thank you for your time!

$\endgroup$
7
  • 2
    $\begingroup$ In bacteria de novo assemblies are the norm because they shift their genomes about a lot (which @Laura will know). A reference genome would be easy - even Spades would do this and that is likely a long way from the optimal approach. First question, which genome? Secondly, this is an impressive data set - long read, short read why not simply use advanced gap closing algorithms and the latest de novo assemblers? $\endgroup$
    – M__
    Jan 9, 2023 at 16:49
  • $\begingroup$ Can you please explain how you built this genome in the first place? That will help to avoid answers that duplicate things you've already tried. $\endgroup$
    – gringer
    Jan 9, 2023 at 19:27
  • 2
    $\begingroup$ @M__ Thank you! Edited the question with more details, including workflow. I've used de novo assembly (never used reference-guided assembly before). But I can't help but wonder whether I can close some of the gaps still remaining by considering the NCBI genomes as well. Perhaps not - I imagine the gaps are characterized by repetitive sequences that would break any assembly...but maybe? One of the NCBI genomes is the same strain from the same culture collection that I used, so I thought there's a chance that rearrangements would not have happened. I'm not sure if this is a good idea. $\endgroup$
    – Laura
    Jan 10, 2023 at 12:31
  • 1
    $\begingroup$ @Laura I will put an answer together this week. Mine is old technology in this area for this technique. However, I personally don't think any reference genome approach will be an improvement, I think it will result in alternative gaps. I'll highlight an alternative modern de novo approach as well. $\endgroup$
    – M__
    Jan 10, 2023 at 13:03
  • 1
    $\begingroup$ What is your input data? Is it isolate or metagenome? On one hand, you mentioned MAG (metagenome-assembled genome) and vamb. On the other hand, you are running Unicycler – I thought Unicycler wouldn't work well with metagenome samples? Also, how "old" is your Nanopore data? Flye/Canu can handle >10% error rate. If you have FAST5 files for the Nanopore data, redo base calling with the latest basecaller. This might dramatically improve the base accuracy. Your assembly seems to have much higher quality than the NCBI one. I would ignore the NCBI assembly. $\endgroup$
    – user172818
    Jan 11, 2023 at 18:58

2 Answers 2

2
$\begingroup$

If your goal is to compare to the NCBI version I'd suggest using QUAST to compare the two assemblies, supplying the NCBI assembly as the reference genome (-r option I think). This will not guarantee anything, but you would expect the two assemblies to be fairly syntenic and to have good overlap of sequences. The tool includes visualization, but you could also use an approach like circoletto or dot plots for something a little more intuitive/simpler.

I'd expect your genome to be quite a bit higher quality even as a metagenome-derived assembly due to nanopore contiguity, I'm not sure that using the NCBI reference to improve yours is a great idea. I'd instead suggest following @M__'s suggestion to scaffold with your nanopore reads (you will have to map them onto your full metagenome first, though approach may differ based on what the scaffolder wants). I'm not up on the current methods here but here is a recent tool. If you have 10 pieces I'd expect to be able to close at least some of that using nanopore (even old).

If you do choose to use the reference to improve, the tool that people seem to be using now is RagTag.

$\endgroup$
2
  • 1
    $\begingroup$ I've tried a couple of things that were suggested: in the end, LongStitch, which you've linked to, proved the best option, and it looks like it got the main genome to 3-4 pieces with some potential plasmids on the side (will have to check). Thank you! $\endgroup$
    – Laura
    Jan 19, 2023 at 13:51
  • $\begingroup$ @Laura great to know, glad that this seems to have been helpful! $\endgroup$ Jan 19, 2023 at 16:57
2
$\begingroup$

Summary What I am saying in this post is there approximates to 3 "generations" of in silico gap-closers/polishes:

  1. GapFiller - amongst other "early" generations (2012);
  2. medaka/racon/Nanopolish/Pilon (2014 - 2018 inclusively)
  3. Nextpolish/PEPPER/Apollo/Homopolish/NeuralPolish (2020 onwards).

The trendy approach is combining 2. and 3. I can't guarantee it will work for cyanobacteria gap-closing, but the shear volume of algorithms dedicated to this purpose, infers they must be delivering results.

Background Whats the issue with a reference genome? The key thing about many bacteria is their ability to "dump" genes and then regain them via plasmid associated recombination. This can occur extensively within a single species. This makes reference associated assembly essentially impossible in many bacteria - particularly against the official NCBI reference genome. Hence, the official NCBI reference gene, is not like a human genome reference - which is a gold standard. With due to respect to those who have dedicated time to presenting the official bacterial reference genomes are not very useful.

This BioStars post here demonstrates that even Mycobacteria tuberculosis would not assemble via a reference genome and TB in context is one of the most conserved genomes in bacteria, the SNP behaviour is famously limited.

If attempting a reference genome it is worth considering the following:

  1. Any amount of genetic divergence can result genetic island phenomena and non-contiguous genomes that will make assembly difficult, therefore a complete genome as close to the target genome as possible is a good reference.
  • For example, simply Blasting sections of you genome to look for a consistent reference that is close to yours. A custom database of complete genomes for your species my be useful, but just sending bits of your target onto the NCBI website will be useful. The idea is to get lots of bits of the genome matching the closest strain.
  1. QUAST was going to be my suggestion too.
  2. Mapping via e.g. Mummer4 against the reference genome as the "scaffold".

Key question:

The same strain also has 1-2 very fragmented genomes (100s of pieces) available in NCBI at scaffold level. The reads themselves are not available, so I can't just add them to the assembly.

I'm not sure what "not available" means.

Key suggestion: Polishing The other approach are recent de novo assemblers and in silico "polishing". What I'm recommending is to first try "polishing". This doesn't exclude subsequent reference genome assemblies - not at all - but it's a way to improve the de novo assemblies and if there are still issues then consider a reference genome.

Firstly I would read someone who is really up to date, e.g.

Comparative evaluation of Nanopore polishing tools for microbial genome assembly and polishing strategies for downstream analysis

They are medaka fans and are using this as part of a combination approach. medaka was part of the established stable of "polishers"

  • GapFiller - I don't think this is useful.
  • racon is specifically for long reads
  • Nanopolish here

The authors prefer medaka then goes on to describe

Nextpolish, PEPPER, Apollo, Homopolish, and NeuralPolish

Much more recent from 2020 ... thats a lot of polishing tools. Note, PEPPER and Apollo are often used for algorithms names - make sure its 2020 onwards, with PEPPER being 2021. They suggest medaka x PEPPER.

I really think this combination polishing approach with an established (medaka) and one the very latest polishing algorithms, such as PEPPER (2021) is really worth trying. If it doesn't work for the specific purpose it will certainly improve the contig quality. However, given the shear volume of work dedicated to this area then there has to be good reason.

By the same token that

Money doesnt give happiness, but sure gives a better standard of misery.

Then ...

Polishing might not give gap closure, but sure gives a better standard of quality to the same old bunch of miserable contigs


de novo Assembly Maybe Miniasm? I'd look at "polishing" the de novo assembly you have.

RagTag as a SPAdes alternative is interesting @MaximilianPress

$\endgroup$
2
  • 2
    $\begingroup$ I'm not convinced long-read polishing is as useful for hybrid data, but I've given medaka and pilon (short-read polishing) a try. It didn't work in this particular case (at best moved some indels around, at worst it made more), but that's more because I've only tested things on my absolutely worst sample. Either way, it didn't work in this case, but now I know about these tools, so I wanted to say thank you! $\endgroup$
    – Laura
    Jan 19, 2023 at 13:57
  • 1
    $\begingroup$ @Laura one day the algorithms will catch up with the problem, they always do, but clearly not today. So, I'm really sorry it's not worked. Keep in mind the 2020+ tools. From a distance its seems an exciting future, where all de novo assemblies would be resolved via combination gap-filling and polishing. I was too optimistic ... for now. That means complex mapping - and a lot of work, via QUAST, Mummer and Daniel Standage's recommendations (his suggests looked cool). $\endgroup$
    – M__
    Jan 19, 2023 at 15:16

Your Answer

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

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