Phage Origin of Mitochondrion-Localized Family A DNA Polymerases in Kinetoplastids and Diplonemids

Acknowledgments

This work was supported by the grants from the Japanese

Society for Promotion of Sciences awarded to Y.I. (numbers

18KK0203 and 19H03280). We also thank three reviewers

for their constructive suggestions and comments on the

manuscript.

Data Availability

The supplementary data are available from the Dryad Digital

Repository: https://doi.org/10.5061/dryad.9kd51c5fv.

Adriaenssens EM, et al. 2020. Taxonomy of prokaryotic viruses: 2018–

2019 update from the ICTV bacterial and archaeal viruses subcommittee. Arch Virol. 165(5):1253–1260.

Bilewitch JP, Degnan SM. 2011. A unique horizontal gene transfer event

has provided the octocoral mitochondrial genome with an active mismatch repair gene that has potential for an unusual self-contained

function. BMC Evol Biol. 11(1):228.

Burger G, Valach M. 2018. Perfection of eccentricity: mitochondrial

genomes of diplonemids. IUBMB Life. 70(12):1197–1206.

Camacho C, et al. 2009. BLASTþ: architecture and applications. BMC

Bioinformatics 10(1):421.

n T. 2009. trimAl: a tool

Capella-Guti

errez S, Silla-Martınez JM, Gabaldo

for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25(15):1972–1973.

Dob

akov

a E, Flegontov P, Skalicky T, Lukes J. 2015. Unexpectedly streamlined mitochondrial genome of the euglenozoan Euglena gracilis.

Genome Biol Evol. 7(12):3358–3367.

El-Gebali S, et al. 2019. The Pfam protein families database in 2019.

Nucleic Acids Res. 47(D1):D427–D432.

Fricova D, et al. 2010. The mitochondrial genome of the pathogenic yeast

Candida subhashii: GC-rich linear DNA with a protein covalently attached to the 5’ termini. Microbiology 156(7):2153–2163.

Fu L, Niu B, Zhu Z, Wu S, Li W. 2012. CD-HIT: accelerated for clustering the

next-generation sequencing data. Bioinformatics 28(23):3150–3152.

 n T, Huynen MA. 2007. From endosymbiont to host-controlled

Gabaldo

organelle: the hijacking of mitochondrial protein synthesis and metabolism. PLoS Comput Biol. 3(11):2209–2218.

Gray MW. 2015. Mosaic nature of the mitochondrial proteome: implications for the origin and evolution of mitochondria. Proc Natl Acad Sci U

S A. 112(33):10133–10138.

Graziewicz MA, Longley MJ, Copeland WC. 2006. DNA polymerase c in

mitochondrial DNA replication and repair. Chem Rev.

106(2):383–405.

Harada R, et al. 2020. Inventory and evolution of mitochondrion-localized

family A DNA polymerases in Euglenozoa. Pathogens 9(4):257.

Hirakawa Y, Watanabe A. 2019. Organellar DNA polymerases in complex

plastid-bearing algae. Biomolecules 9(4):140–112.

Inagaki Y, Harada R. [dataset] 2020. Data from: phage origin of

mitochondrion-localized family A DNA polymerases in kinetoplastids

and diplonemids. Dryad. doi: 10.5061/dryad.9kd51c5fv.

Genome Biol. Evol. 13(2) doi:10.1093/gbe/evab003 Advance Access publication January 12, 2021

Downloaded from https://academic.oup.com/gbe/article/13/2/evab003/6081025 by University of Tsukuba user on 16 February 2021

We retrieved 175 famA DNAP aa sequences of autographiviruses from the NCBI nr database. The details of the survey

were the same as described above. The 175 famA DNAPs

were sampled from 99 members belonging to 57 genera

and

76

unclassified

members

in

the

family

Autographiviridae. The autographivirus famA DNAPs were

found to comprise two types based on the presence/absence

of an insertion of eight aa residues (8-aa insertion; see above).

The famA DNAPs with 8-aa insertion (AGVþins famA DNAPs)

appeared to be closely related to PolIBCDþ, mitochondrionlocalized famA DNAPs in kinetoplastids (PolIB, C, D, PolIPerk1/2) and that in diplonemids (PolI-dipl). The redundancy

among the AGVþins famA DNAPs was reduced by a cluster

analysis using CD-HIT v4.7 with a threshold of 90%. Finally,

we aligned the aa sequences of 74 AGVþins famA DNAPs, 24

PolIBCDþ, and famA DNAPs of four members of Podoviridae

by MAFFT v7.455 with the L-INS-i model. Ambiguously

aligned positions were discarded manually, and gapcontaining positions were trimmed by using trimAl v1.4

with the -gt 0.9 option. The final version of the second alignment is provided as a part of the supplementary materials

(Inagaki and Harad[dataset] 2020). The final alignment containing 102 sequences with 581 unambiguously aligned aa

positions was subjected to both ML and Bayesian phylogenetic analyses. The ML and ML bootstrap analyses were performed as described above. For Bayesian analysis using

Phylobayes v4.1 (Lartillot et al. 2009), we ran four Markov

Chain Monte Carlo chains for 100,000 cycles with burn-in

of 25,000 (maxdiff ¼ 0.09472) and calculated the consensus

tree with branch lengths and BPPs from the remaining trees.

The aa substitution model was set to CAT þ GTR in

Phylobayes analysis described above.

Supplementary data are available at Genome Biology and

Evolution online.

GBE

Harada and Inagaki

Nguyen LT, Schmidt HA, Von Haeseler A, Minh BQ. 2015. IQ-TREE: a fast

and effective stochastic algorithm for estimating maximum-likelihood

phylogenies. Mol Biol Evol. 32(1):268–274.

Nishimura Y, et al. 2019. Horizontally-acquired genetic elements in the

mitochondrial genome of a centrohelid Marophrys sp. Sci Rep.

9(1):4850.

Rodrıguez-Ezpeleta N, et al. 2007. Toward resolving the eukaryotic tree:

the phylogenetic positions of Jakobids and Cercozoans. Curr Biol.

17(16):1420–1425.

~ oz-Go

mez SA, Kamikawa R. 2017. The origin and diversiRoger AJ, Mun

fication of mitochondria. Curr Biol. 27(21):R1177–R1192.

Roy J, Faktorov

a D, Lukes J, Burger G. 2007. Unusual mitochondrial genome structures throughout the Euglenozoa. Protist 158(3):385–396.

Sayers EW, et al. 2020. Database resources of the national center for

biotechnology information. Nucleic Acids Res. 48(D1):D9–D16.

Shutt TE, Gray MW. 2006. Bacteriophage origins of mitochondrial replication and transcription proteins. Trends Genet. 22(2):90–95.

Spencer DF, Gray MW. 2011. Ribosomal RNA genes in Euglena gracilis

mitochondrial DNA: fragmented genes in a seemingly fragmented

genome. Mol Genet Genomics. 285(1):19–31.

Swart EC, et al. 2012. The Oxytricha trifallax mitochondrial genome.

Genome Biol Evol. 4(2):136–154.

Wang HC, Minh BQ, Susko E, Roger AJ. 2018. Modeling site heterogeneity

with posterior mean site frequency profiles accelerates accurate phylogenomic estimation. Syst Biol. 67(2):216–235.

Wang Z, Wu M. 2014. Phylogenomic reconstruction indicates mitochondrial ancestor was an energy parasite. PLoS One 9(10):e110685.

Yabuki A, et al. 2011. Tsukubamonas globosa n. gen., n. sp., a novel

excavate flagellate possibly holding a key for the early evolution in

“Discoba”. J Eukaryot Microbiol. 58(4):319–331.

Yabuki A, Gyaltshen Y, Heiss AA, Fujikura K, Kim E. 2018. Ophirina

amphinema n. gen., n. sp., a new deeply branching discobid with

phylogenetic affinity to jakobids. Sci Rep. 8(1):1–14.

Yabuki A, Tanifuji G, Kusaka C, Takishita K, Fujikura K. 2016. Hyper-eccentric structural genes in the mitochondrial genome of the algal parasite Hemistasia phaeocysticola. Genome Biol Evol. 8(9):2870–2878.

Associate editor: Martin Embley

Genome Biol. Evol. 13(2) doi:10.1093/gbe/evab003 Advance Access publication January 12, 2021

Downloaded from https://academic.oup.com/gbe/article/13/2/evab003/6081025 by University of Tsukuba user on 16 February 2021

Janouskovec J, et al. 2015. Factors mediating plastid dependency and the

origins of parasitism in apicomplexans and their close relatives. Proc

Natl Acad Sci U S A. 112(33):10200–10207.

Jones P, et al. 2014. InterProScan 5: genome-scale protein function classification. Bioinformatics 30(9):1236–1240.

Jung GH, Leavitt MC, Hsieh JC, Ito J. 1987. Bacteriophage PRD1 DNA

polymerase: evolution of DNA polymerases. Proc Natl Acad Sci U S

A. 84(23):8287–8291.

Kalyaanamoorthy S, Minh BQ, Wong TKF, Von Haeseler A, Jermiin LS.

2017. ModelFinder: fast model selection for accurate phylogenetic

estimates. Nat Methods. 14(6):587–589.

Kamikawa R, et al. 2014. Gene content evolution in discobid mitochondria

deduced from the phylogenetic position and complete mitochondrial

genome of Tsukubamonas globosa. Genome Biol Evol. 6(2):306–315.

Katoh K, Standley DM. 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol

Evol. 30(4):772–780.

Klingbeil MM, Motyka SA, Englund PT. 2002. Multiple mitochondrial DNA

polymerases in Trypanosoma brucei. Mol Cell. 10(1):175–186.

Krasich R, Copeland WC. 2017. DNA polymerases in the mitochondria: a

critical review of the evidence. Physiol Behav. 22(1):692–709.

Lartillot N, Lepage T, Blanquart S. 2009. PhyloBayes 3: a Bayesian software

package for phylogenetic reconstruction and molecular dating.

Bioinformatics 25(17):2286–2288.

Lavigne R, Seto D, Mahadevan P, Ackermann HW, Kropinski AM. 2008.

Unifying classical and molecular taxonomic classification: analysis of

the Podoviridae using BLASTP-based tools. Res Microbiol.

159(5):406–414.

Li W, Godzik A. 2006. Cd-hit: a fast program for clustering and comparing

large sets of protein or nucleotide sequences. Bioinformatics

22(13):1658–1659.

Lukes J, et al. 2002. Kinetoplast DNA network: evolution of an improbable

structure. Eukaryot Cell. 1(4):495–502.

Moriyama T, Sato N. 2014. Enzymes involved in organellar DNA replication

in photosynthetic eukaryotes. Front Plant Sci. 5:480.

Moriyama T, Terasawa K, Sato N. 2011. Conservation of POPs, the plant

organellar DNA polymerases, in eukaryotes. Protist 162(1):177–187.

...