Regarding the Use of Incomplete Plant Mitochondrial Genome Assembly Outputs for Subsequent Analysis

[twdawd]

Thank you very much for making this software. I have a question: When I assembled the plant mitochondrial genome, I got a file named embplant_mt.K105.scaffolds.graph1.1.path_sequence.fasta. Can I use the generated file(fasfa file) in the next step of analysis such as gene annotation and phylogeny?
Looking forward to your reply

[JianjunJin]

The assembly quality of incomplete organelle genomes (scaffolds) cannot be guaranteed in the current version of GetOrganelle (1.7.7.0 or earlier). However, enhancements were made in version 1.8.0 (although not officially released yet, see https://github.com/Kinggerm/GetOrganelle/tree/intermediate_graph) for chloroplast and animal mitochondrial genomes.

The primary concern of directly using the plant mitochondrial assembly lies in the potential inclusion of chloroplast loci, due to the similarity between pt and mt and the false hits of blast. This becomes severe when the assembly is incomplete because depth and graph information can barely be utilized to help identify contigs. Therefore, it is advisable to run get_organelle_from_assembly.py using the current output, implementing stricter contig filtering measures, which can be quick and easy. The following command is an example proposed for this purpose. Still, you are encouraged to carefully assess the possibility of chloroplast inclusion in mitochondrial assemblies even after this.

 get_organelle_from_assembly.py -F embplant_mt --slim-options "--perc-identity 99.0 -e 1e-250" -g your_old_output/filtered_K105.assembly_graph.fastg.extend-embplant_mt-embplant_pt.fastg -o further_purified_result

[twdawd]

I understand !!! Thank you !!!