the pipeline about panEDTA

[QianghuiZhu]

Hi, Oushujun!
I have some animal genomes and want to annotate them by panEDTA. But now, I am confused about the some pipeline of panEDTA.
I read some others' issues about panEDTA, and they seem that:

1. Using EDTA for each genome;
2. Using panEDTA to get the pan TE lib;
3. Using panTElib to reannotate each genome again by EDTA.

But I read the example you provided (https://github.com/HuffordLab/NAM-genomes/tree/master/te-annotation)

It seems that:

1. Using EDTA for each genome;
2. Using TElib of each genome to reannotate each genome again by EDTA;
3. Using panEDTA to get the pan TE lib;
4. Using panTElib to reannotate each genome again by EDTA.

So, my question is about,should I reannotate the genome again with their TEblib in step2?

Thanks!
Best Wishes!

[QianghuiZhu]

I used two ways of above with EDTA: pipeline1) using EDTA for genomeA; pipeline2) using the output TElib of pipeline1 as a curated library (--curatedlib), and rerun EDTA for genomeA again.
And here is the result statistics

Some types of TE with different num in these two output results. I checked some TE position and they are same.

Is it possible that reannotating genome by EDTA again is optional, or not need to reannotate again?

Best!

[oushujun]

Hello, Your pipeline 2 step 2 can be merged into step 1 if you specify —anno 1 in your step 1. This annotation is necessary to identify copy numbers of each family. Your results are pretty consistent. The variations can be considered uncertainty of TE annotation. Please check the wiki. Best, Shujun

…

On Thu, Feb 22, 2024 at 9:56 PM Hui ***@***.***> wrote: I used two ways of above with EDTA: pipeline1) using EDTA for genomeA; pipeline2) using the output TElib of pipeline1 as a curated library (--curatedlib), and rerun EDTA for genomeA again. And here is the result statistics image.png (view on web) <https://github.com/oushujun/EDTA/assets/77773766/80e110d2-b8ec-408e-8f17-f0c67101b0af> Some types of TE with different num in these two output results. I checked some TE position and they are same. Is it possible that reannotating genome by EDTA again is optional, or not need to reannotate again? Best! — Reply to this email directly, view it on GitHub <#436 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ABNX4NAQLQ44LY5FLXQL7N3YVAAMXAVCNFSM6AAAAABDV7ZO4WVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTSNRQGY3DKNRSGI> . You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>

[QianghuiZhu]

Thank you, Shujun.
I run pipeline1 and pipeline2 both with --anno 1 previously. I'll just use pipeline1 (with --anno 1) to save time.
Best!

Hello, Your pipeline 2 step 2 can be merged into step 1 if you specify —anno 1 in your step 1. This annotation is necessary to identify copy numbers of each family.
Your results are pretty consistent. The variations can be considered uncertainty of TE annotation. Please check the wiki.
Best, Shujun

[CSU-KangHu]

Hi @oushujun,
I have some confusion about the panEDTA pipeline. It seems that panEDTA runs EDTA on each genome individually and then merges the EDTA TE libraries from each genome. Considering that genomes constituting the panGenome have a significant proportion of shared conserved sequences (core genome), why not first construct the panGenome from different genomes and then use EDTA to identify and annotate TEs on the panGenome?

[oushujun]

Yes, you may do that. Be aware that the graph pangenome may present variable sequences to you without the adjacent sequences, the latter are needed for EDTA to determine the authenticity of a structurally intact TE. Best, Shujun

…

On Fri, Feb 23, 2024 at 1:14 AM Kang Hu ***@***.***> wrote: Hi @oushujun <https://github.com/oushujun>, I have some confusion about the panEDTA pipeline. It seems that panEDTA runs EDTA on each genome individually and then merges the EDTA TE libraries from each genome. Considering that genomes constituting the panGenome have a significant proportion of shared conserved sequences (core genome), why not first construct the panGenome from different genomes and then use EDTA to identify and annotate TEs on the panGenome? — Reply to this email directly, view it on GitHub <#436 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ABNX4NBJ6QDDOIPBIRSJW4TYVAXTVAVCNFSM6AAAAABDV7ZO4WVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTSNRQG44DKMJTGM> . You are receiving this because you were mentioned.Message ID: ***@***.***>

[QianghuiZhu]

Hi, Shujun.
I am sorry to bother you again.

I downloaded some genome and cds data from NCBI, and want to annotate them using EDTA.
Some homologous chr sequence in different genomes are inverted. And I plan to use the reverse complement sequence of these chr to ensure the consistency of strand. May I also need to use the reverse complement sequence of cds data in these chr? Or primary cds data is ok and strand INFO have no impacts on TE annotation?

Thanks!
Best!