# How to transfer gff annotations in genome with extensive duplications?

Microbial genomes can contain extensive duplications. Often we'd like to transfer annotations from an annotated species to one that is newly sequenced.

Existing tools (e.g. RATT, LiftOver, Kraken) either make specific assumptions about how closely related the species are or fail to transfer when multiple matches are found in the new genome, especially if the sequences are highly similar.

Specifically, I have a synthetic biology application where genes can duplicate extensively. They are identical in sequence but duplicated many times and be relocated (i.e., not just adjacent to each other). None of the above mentioned tools are able to transfer coordinates of annotations to genomes with multiple copies of features.

Are there any pre-existing tools or software that transfer annotations in this scenario? Ideas for ways to do this robustly?

• Can you provide more detail? What's your input (e.g. raw reads, assemblies, open-reading frames)? what kind of sequencing? what is the exact output you want? ' I'm not sure I understand your point about the assumptions. Also, LiftOver and Kraken are completely different tools with different uses May 31 '17 at 11:19
• Would this even be theoretically possible? How can you assume the annotations are transferable if there are extensive duplications? It's probably better to look for homologs instead. May 31 '17 at 11:22
• @terdon do you mean orthologs? homologs = orthologs (non-duplicated) + paralogs (duplicated) May 31 '17 at 11:25
• @Chris_Rands No, I meant homologs. Precisely because we can't know whether they are ortho- or para- (I have a nice post on the difference between the two here, by the way), so all you can do first is find homologs and then try and figure out whether they are similar enough to carry any annotations over. May 31 '17 at 11:54
• @terdon I see well resolving orthologs/paralogs is not easy of course, but it can be done, depending on the exact data (I don't know what the OP's data looks like), for example some of my colleagues maintain orthodb orthodb.org May 31 '17 at 12:03

There is one very simplistic way I use which might work for what you are doing, it is similar to what terdon proposed.

Take a de-novo microbial genome annotation tool (I have my own, but you could use/modify prokka). Tools like these often first predict gene boundaries (with other tools like prodigal or glimmer) and then try to assign a function to found genes. This function assignment is often done with BLAST and other tools ... and that is where you can go in and modify to do what you need.

I use a "knowledge" protein database of genes I want to have very strictly annotated as a first line of annotation (e.g. in your case: the annotated genomes). For that I loop through very strict identity/similarity parameters which get gradually relaxed.

E.g.: Loop 0: only transfer annotations at 100% DNA identity, same length. Loop 1: only transfer annotations at 100% similarity, same length. Loop 2: only transfer annotations at 99% similarity, length +/- 1%. ... Loop n: only transfer annotations at 100-(n-1)% similarity, length +/- (n-1)%.

In each loop, obviously only annotate what has not been annotated in previous loops.

After that, use "normal" annotation pipeline of the tool to annotate the rest.

• Doesn't that require the target genome's genes to have been found first? Or can your tool also do de-novo gene prediction? (sounds like a very useful tool, by the way, kudos!) May 31 '17 at 16:33
• Prokaryotik gene finding/prediction is a more or less solved problem, existing tools work reasonably well. See prodigal.ornl.gov and prodigal.ornl.gov (just to name two).
– BaCh
May 31 '17 at 17:27
• Yes, I know, I was just surprised you didn't mention that in your answer. If I understand correctly, the first step would be for the OP to find the list of putative genes in their newly sequenced genome, right? May 31 '17 at 17:31
• Correct. Prokka (vicbioinformatics.com/software.prokka.shtml) uses a whole battery of third-party tools (including prodigal) to annotate a genome de-novo, I started by modifying prokka before I wrote my own one, which uses some ideas from the prokka pipeline.
– BaCh
May 31 '17 at 17:37

I think you will have to first identify the regions homologous to the ones defined in your GFF and then transfer the annotations. Of course, the assumption there is that the homolog will also have the same annotation which is often not true. However, I don't see how you can do it in any other way since you cannot use genomic coordinates (and you would still be making the same assumption even if you could, anyway) when the genomes are so different.

For a very simplistic approach (which might be enough if, as you say, your sequences are almost identical), you can do something like:

2. Use a tool like genewise or exonerate to map these into the target genome. Both tools can return gff-formatted output and both can find multiple hits in the target genome. For what you want, I would suggest using a very high threshold of sequence similarity and query coverage (where the target sequence found covers all or most of the query sequence used).