The evolution of intracellular endosymbiosis marks a major transition in the biology of the host and endosymbiont. Yet, how adaptation manifests in the genomes of the participants remains relatively understudied. We investigated this question by sequencing the genomes of Tetrahymena utriculariae, a commensal of the aquatic carnivorous bladderwort Utricularia reflexa, and its intracellular algae, Micractinium tetrahymenae. We discovered an expansion in copy number and negative selection in a TLD domain-bearing gene family in the genome of T. utriculariae, identifying it as a candidate for being an adaptive response to oxidative stress resulting from the physiology of its endosymbionts. We found that the M. tetrahymenae genome is larger than those of other Micractinium and Chlorella and contains a greater number of rapidly expanding orthogroups. These were enriched for Gene Ontology terms relevant to the regulation of intracellular signal transduction and cellular responses to stress and stimulus. Single-exon tandem repeats were overrepresented in paralogs belonging to these rapidly expanding orthogroups, which implicates long terminal repeat retrotransposons (LTRs) as potential agents of adaptation. We additionally performed a comparative transcriptomic analysis of M. tetrahymenae in a free-living state and in endosymbiosis with T. utriculariae and discovered that the genes that are differentially expressed were enriched for pathways that evidence shifts in energy generation and storage and in cellular protection strategies. Together, our results elucidate the axes along which the participants must adapt in this young endosymbiosis and highlight evolutionary responses to stress as a shared trend.

Aquatic Ecology and Evolution Department of Biology University of Konstanz Konstanz 78464 Germany

Department of Biology University of Konstanz Konstanz 78464 Germany

Department of Experimental and Functional Morphology Institute of Botany CAS Třeboň CZ 379 01 Czech Republic

Institute for Advanced Computational Science Stony Brook University Stony Brook NY 11794 USA

Zobrazit více v PubMed

Adamec  L. Oxygen concentrations inside the traps of the carnivorous plants Utricularia and Genlisea (Lentibulariaceae). Ann Bot. 2007:100(4):849–856. 10.1093/aob/mcm182. PubMed DOI PMC

Altschul  SF, Gish  W, Miller  W, Myers  EW, Lipman  DJ. Basic local alignment search tool. J Mol Biol. 1990:215(3):403–410. 10.1016/S0022-2836(05)80360-2. PubMed DOI

Arriola  MB, Velmurugan  N, Zhang  Y, Plunkett  MH, Hondzo  H, Barney  BM. Genome sequences of Chlorella sorokiniana UTEX 1602 and Micractinium conductrix SAG 241.80: implications to maltose excretion by a green alga. Plant J.  2018:93(3):566–586. 10.1111/tpj.13789. PubMed DOI

Blanc  G, Duncan  G, Agarkova  I, Borodovsky  M, Gurnon  J, Kuo  A, Lindquist  E, Lucas  S, Pangilinan  J, Polle  J, et al.  The Chlorella variabilis NC64A genome reveals adaptation to photosymbiosis, coevolution with viruses, and cryptic sex. Plant Cell. 2010:22(9):2943–2955. 10.1105/tpc.110.076406. PubMed DOI PMC

Brown  JA, Nielsen  PJ. Transfer of photosynthetically produced carbohydrate from endosymbiotic Chlorellae to Paramecium bursaria. J Protozool. 1974:21(4):569–570. 10.1111/j.1550-7408.1974.tb03702.x. PubMed DOI

Brůna  T, Lomsadze  A, Borodovsky  M. GeneMark-EP+: eukaryotic gene prediction with self-training in the space of genes and proteins. NAR Genom Bioinform. 2020:2(2):lqaa026. 10.1093/nargab/lqaa026. PubMed DOI PMC

Brunk  CF. Complete sequence of the mitochondrial genome of Tetrahymena thermophila and comparative methods for identifying highly divergent genes. Nucleic Acids Res. 2003:31(6):1673–1682. 10.1093/nar/gkg270. PubMed DOI PMC

Cantalapiedra  CP, Hernández-Plaza  A, Letunic  I, Bork  P, Huerta-Cepas  J. EggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol Biol Evol. 2021:38(12):5825–5829. 10.1093/molbev/msab293. PubMed DOI PMC

Castresana  J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000:17(4):540–552. 10.1093/oxfordjournals.molbev.a026334. PubMed DOI

Celaya  G, Fernández-Higuero  JA, Martin  I, Rivas  G, Moro  F, Muga  A. Crowding modulates the conformation, affinity, and activity of the components of the bacterial disaggregase machinery. J Mol Biol. 2016:428(11):2474–2487. 10.1016/j.jmb.2016.04.027. PubMed DOI

Chen  S. Ultrafast one-pass FASTQ data preprocessing, quality control, and deduplication using fastp. iMeta. 2023:2(2):e107. 10.1002/imt2.107. PubMed DOI PMC

Cheng  H, Concepcion  GT, Feng  X, Zhang  H, Li  H. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods. 2021:18(2):170–175. 10.1038/s41592-020-01056-5. PubMed DOI PMC

Cheng  Y-H, Liu  C-FJ, Yu  Y-H, Jhou  Y-T, Fujishima  M, Tsai  IJ, Leu  J-Y. Genome plasticity in Paramecium bursaria revealed by population genomics. BMC Biol. 2020:18(1):180. 10.1186/s12915-020-00912-2. PubMed DOI PMC

Cvrčková  F, De Virgilio  C, Manser  E, Pringle  JR, Nasmyth  K. Ste20-like protein kinases are required for normal localization of cell growth and for cytokinesis in budding yeast. Genes Dev. 1995:9(15):1817–1830. 10.1101/gad.9.15.1817. PubMed DOI

Dema  A, Perets  E, Schulz  MS, Deák  VA, Klussmann  E. Pharmacological targeting of AKAP-directed compartmentalized cAMP signalling. Cell Signal. 2015:27(12):2474–2487. 10.1016/j.cellsig.2015.09.008. PubMed DOI

de Silva  NDG, Murmu  J, Chabot  D, Hubbard  K, Ryser  P, Molina  I, Rowland  O. Root suberin plays important roles in reducing water loss and sodium uptake in Arabidopsis thaliana. Metabolites. 2021:11(11):735. 10.3390/metabo11110735. PubMed DOI PMC

Diaz  JM, Plummer  S. Production of extracellular reactive oxygen species by phytoplankton: past and future directions. J Plankton Res. 2018:40(6):655–666. 10.1093/plankt/fby039. PubMed DOI PMC

Dobin  A, Davis  CA, Schlesinger  F, Drenkow  J, Zaleski  C, Jha  S, Batut  P, Chaisson  M, Gingeras  TR. STAR: ultrafast universal RNA-Seq aligner. Bioinformatics. 2013:29(1):15–21. 10.1093/bioinformatics/bts635. PubMed DOI PMC

Durand  M, Kolpak  A, Farrell  T, Elliott  NA, Shao  W, Brown  M, Volkert  MR. The OXR domain defines a conserved family of eukaryotic oxidation resistance proteins. BMC Cell Biol. 2007:8(1):13. 10.1186/1471-2121-8-13. PubMed DOI PMC

Edgar  RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004:32(5):1792–1797. 10.1093/nar/gkh340. PubMed DOI PMC

Edreva  A. Generation and scavenging of reactive oxygen species in chloroplasts: a submolecular approach. Agric Ecosyst Environ. 2005:106(2-3):119–133. 10.1016/j.agee.2004.10.022. DOI

Elliott  NA, Volkert  MR. Stress induction and mitochondrial localization of Oxr1 proteins in yeast and humans. Mol Cell Biol. 2004:24(8):3180–3187. 10.1128/MCB.24.8.3180-3187.2004. PubMed DOI PMC

Emms  DM, Kelly  S. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019:20(1):238. 10.1186/s13059-019-1832-y. PubMed DOI PMC

Esteban  GF, Fenchel  T, Finlay  BJ. Mixotrophy in ciliates. Protist. 2010:161(5):621–641. 10.1016/j.protis.2010.08.002. PubMed DOI

Estrela  S, Kerr  B, Morris  JJ. Transitions in individuality through symbiosis. Curr Opin Microbiol. 2016:31:191–198. 10.1016/j.mib.2016.04.007. PubMed DOI

Fellows  I. Deducer: a data analysis GUI for R. J Stat Softw. 2012:49(8):1–15. 10.18637/jss.v049.i08. DOI

Fichman  Y, Rowland  L, Oliver  MJ, Mittler  R. ROS are evolutionary conserved cell-to-cell stress signals. Proc Natl Acad Sci U S A.  2023:120(31):e2305496120. 10.1073/pnas.2305496120. PubMed DOI PMC

Finelli  MJ, Oliver  PL. TLDc proteins: new players in the oxidative stress response and neurological disease. Mamm Genome.  2017:28(9-10):395–406. 10.1007/s00335-017-9706-7. PubMed DOI PMC

Finelli  MJ, Sanchez-Pulido  L, Liu  KX, Davies  KE, Oliver  PL. The evolutionarily conserved Tre2/Bub2/Cdc16 (TBC), lysin motif (LysM), domain catalytic (TLDc) domain is neuroprotective against oxidative stress. J Biol Chem.  2016:291(6):2751–2763. 10.1074/jbc.M115.685222. PubMed DOI PMC

Flynn  JM, Hubley  R, Goubert  C, Rosen  J, Clark  AG, Feschotte  C, Smit  AF. RepeatModeler2 for automated genomic discovery of transposable element families. Proc Natl Acad Sci U S A.  2020:117(17):9451–9457. 10.1073/pnas.1921046117. PubMed DOI PMC

Fokin  SI, Przyboś  E, Chivilev  SM, Beier  CL, Horn  M, Skotarczak  B, Wodecka  B, Fujishima  M. Morphological and molecular investigations of Paramecium schewiakoffi sp. nov. (Ciliophora, Oligohymenophorea) and current status of distribution and taxonomy of Paramecium spp. Eur J Protistol. 2004:40(3):225–243. 10.1016/j.ejop.2004.02.001. DOI

Frickel  J, Sieber  M, Becks  L. Eco-evolutionary dynamics in a coevolving host–virus system. Ecol Lett. 2016:19(4):450–459. 10.1111/ele.12580. PubMed DOI

Fujishima  M, Kodama  Y. Mechanisms for establishing primary and secondary endosymbiosis in Paramecium. J Eukaryot Microbiol.  2022:69(5):e12901. 10.1111/jeu.12901. PubMed DOI

Gorovsky  MA, Yao  M-C, Keevert  JB, Pleger  GL. Isolation of micro- and macronuclei of Tetrahymena pyriformis. Methods Cell Biol. 1975:9:311–327. 10.1016/S0091-679X(08)60080-1. PubMed DOI

Grabherr  MG, Haas  BJ, Yassour  M, Levin  JZ, Thompson  DA, Amit  I, Adiconis  X, Fan  L, Raychowdhury  R, Zeng  Q, et al.  Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011:29(7):644–652. 10.1038/nbt.1883. PubMed DOI PMC

Gray  PC, Johnson  BD, Westenbroek  RE, Hays  LG, Yates  JR, Scheuer  T, Catterall  WA, Murphy  BJ. Primary structure and function of an A kinase anchoring protein associated with calcium channels. Neuron. 1998:20(5):1017–1026. 10.1016/S0896-6273(00)80482-1. PubMed DOI

Guerrero  MG, Vega  JM, Losada  M. The assimilatory nitrate-reducing system and its regulation. Annu Rev Plant Physiol. 1981:32(1):169–204. 10.1146/annurev.pp.32.060181.001125. DOI

Guillard  RRL, Lorenzen  CJ. Yellow-green algae with chlorophyllide C1,2. J Phycol. 1972:8:10–14. 10.1111/j.1529-8817.1972.tb03995.x. DOI

Haas  BJ. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003:31(19):5654–5666. 10.1093/nar/gkg770. PubMed DOI PMC

Haas  BJ, Salzberg  SL, Zhu  W, Pertea  M, Allen  JE, Orvis  J, White  O, Buell  CR, Wortman  JR. Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biol. 2008:9(1):R7. 10.1186/gb-2008-9-1-r7. PubMed DOI PMC

Hamilton  EP, Kapusta  A, Huvos  PE, Bidwell  SL, Zafar  N, Tang  H, Hadjithomas  M, Krishnakumar  V, Badger  JH, Caler E  V, et al.  Structure of the germline genome of Tetrahymena thermophila and relationship to the massively rearranged somatic genome. Elife. 2016:5:e19090. 10.7554/eLife.19090. PubMed DOI PMC

Han  D, Williams  E, Cadenas  E. Mitochondrial respiratory chain-dependent generation of superoxide anion and its release into the intermembrane space. Biochem J.  2001:353(2):411–416. 10.1042/bj3530411. PubMed DOI PMC

Hansel  CM, Diaz  JM. Production of extracellular reactive oxygen species by marine biota. Ann Rev Mar Sci. 2021:13(1):177–200. 10.1146/annurev-marine-041320-102550. PubMed DOI

He  M, Wang  J, Fan  X, Liu  X, Shi  W, Huang  N, Zhao  F, Miao  M. Genetic basis for the establishment of endosymbiosis in Paramecium. ISME J. 2019:13(5):1360–1369. 10.1038/s41396-018-0341-4. PubMed DOI PMC

Hennessy  F, Cheetham  ME, Dirr  HW, Blatch  GL. Analysis of the levels of conservation of the J domain among the various types of DnaJ-like proteins. Cell Stress Chaperones. 2000:5(4):347–358. 10.1379/1466-1268(2000)005<0347:AOTLOC>2.0.CO;2. PubMed DOI PMC

Hoang  DT, Chernomor  O, von Haeseler  A, Minh  BQ, Vinh  LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018:35(2):518–522. 10.1093/molbev/msx281. PubMed DOI PMC

Horas  EL, Metzger  SM, Platzer  B, Kelly  JB, Becks  L. Context-dependent costs and benefits of endosymbiotic interactions in a ciliate–algae system. Environ Microbiol. 2022:24(12):5924–5935. 10.1111/1462-2920.16112. PubMed DOI

Horner  A, Goetz  F, Tampé  R, Klussmann  E, Pohl  P. Mechanism for targeting the A-kinase anchoring protein AKAP18δ to the membrane. J Biol Chem.  2012:287(51):42495–42501. 10.1074/jbc.M112.414946. PubMed DOI PMC

Hoshina  R, Fujiwara  Y. Photobiont flexibility in Paramecium bursaria: double and triple photobiont co-habitation. Adv Microbiol. 2012:2(03):227–233. 10.4236/aim.2012.23027. DOI

Hoshina  R, Tsukii  Y, Harumoto  T, Suzaki  T. Characterization of a green Stentor with symbiotic algae growing in an extremely oligotrophic environment and storing large amounts of starch granules in its cytoplasm. Sci Rep. 2021:11(1):2865. 10.1038/s41598-021-82416-9. PubMed DOI PMC

Janec̆ek  S. α-amylase family: molecular biology and evolution. Prog Biophys Mol Biol. 1997:67(1):67–97. 10.1016/S0079-6107(97)00015-1. PubMed DOI

Jaramillo-Gutierrez  G, Molina-Cruz  A, Kumar  S, Barillas-Mury  C. The Anopheles gambiae oxidation resistance 1 (OXR1) gene regulates expression of enzymes that detoxify reactive oxygen species. PLoS One. 2010:5(6):e11168. 10.1371/journal.pone.0011168. PubMed DOI PMC

Jenkins  BH, Maguire  F, Leonard  G, Eaton  JD, West  S, Housden  BE, Milner  DS, Richards  TA. Emergent RNA–RNA interactions can promote stability in a facultative phototrophic endosymbiosis. Proc Natl Acad Sci U S A.  2021:118(38):e2108874118. 10.1073/pnas.2108874118. PubMed DOI PMC

Jones  P, Binns  D, Chang  H-Y, Fraser  M, Li  W, McAnulla  C, McWilliam  H, Maslen  J, Mitchell  A, Nuka  G, et al.  InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014:30(9):1236–1240. 10.1093/bioinformatics/btu031. PubMed DOI PMC

Jumper  J, Evans  R, Pritzel  A, Green  T, Figurnov  M, Ronneberger  O, Tunyasuvunakool  K, Bates  R, Žídek  A, Potapenko  A, et al.  Highly accurate protein structure prediction with AlphaFold. Nature. 2021:596(7873):583–589. 10.1038/s41586-021-03819-2. PubMed DOI PMC

Kalyaanamoorthy  S, Minh  BQ, Wong  TKF, von Haeseler  A, Jermiin  LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017:14(6):587–589. 10.1038/nmeth.4285. PubMed DOI PMC

Kato  Y, Ueno  S, Imamura  N. Studies on the nitrogen utilization of endosymbiotic algae isolated from Japanese Paramecium bursaria. Plant Sci.  2006:170(3):481–486. 10.1016/j.plantsci.2005.09.018. DOI

Katoh  K, Standley  DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013:30(4):772–780. 10.1093/molbev/mst010. PubMed DOI PMC

Kawano  T, Irie  K, Kadono  T. Oxidative stress-mediated development of symbiosis in green paramecia. In: Seckbach  J, Grube  M, editors. Symbioses and stress. Dordrecht, Netherlands: Springer; 2010. p. 177–195.

Kawano  T, Kadono  T, Kosaka  T, Hosoya  H. Green paramecia as an evolutionary winner of oxidative symbiosis: a hypothesis and supportive data. Zeitschrift für Naturforschung C. 2004:59(7-8):538–542. 10.1515/znc-2004-7-816. PubMed DOI

Khorobrykh  S, Havurinne  V, Mattila  H, Tyystjärvi  E. Oxygen and ROS in photosynthesis. Plants. 2020:9(1):91. 10.3390/plants9010091. PubMed DOI PMC

Kim  J-S, Liu  L, Fitzsimmons  LF, Wang  Y, Crawford  MA, Mastrogiovanni  M, Trujillo  M, Till  JKA, Radi  R, Dai  S, et al.  Dksa–DnaJ redox interactions provide a signal for the activation of bacterial RNA polymerase. Proc Natl Acad Sci U S A.  2018:115(50):E11780–E11789. 10.1073/pnas.1813572115. PubMed DOI PMC

Kim  J-S, Liu  L, Vázquez-Torres  A. The DnaK/DnaJ chaperone system enables RNA polymerase-DksA complex formation in Salmonella experiencing oxidative stress. mBio. 2021:12(3):e03443–e03420. 10.1128/mBio.03443-20. PubMed DOI PMC

Kodama  Y, Suzuki  H, Dohra  H, Sugii  M, Kitazume  T, Yamaguchi  K, Shigenobu  S, Fujishima  M. Comparison of gene expression of Paramecium bursaria with and without Chlorella variabilis symbionts. BMC Genomics. 2014:15(1):183. 10.1186/1471-2164-15-183. PubMed DOI PMC

Kolmogorov  M, Yuan  J, Lin  Y, Pevzner  PA. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol. 2019:37(5):540–546. 10.1038/s41587-019-0072-8. PubMed DOI

Kondrashov  FA. Gene duplication as a mechanism of genomic adaptation to a changing environment. Proc R Soc Lond B Biol Sci.  2012:279(1749):5048–5057. 10.1098/rspb.2012.1108. PubMed DOI PMC

Kosakovsky Pond  SL, Frost  SDW. Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol. 2005:22(5):1208–1222. 10.1093/molbev/msi105. PubMed DOI

Kosakovsky Pond  SL, Poon  AFY, Velazquez  R, Weaver  S, Hepler  NL, Murrell  B, Shank  SD, Magalis  BR, Bouvier  D, Nekrutenko  A, et al.  Hyphy 2.5—a customizable platform for evolutionary hypothesis testing using phylogenies. Mol Biol Evol. 2020:37(1):295–299. 10.1093/molbev/msz197. PubMed DOI PMC

Krasileva  KV. The role of transposable elements and DNA damage repair mechanisms in gene duplications and gene fusions in plant genomes. Curr Opin Plant Biol. 2019:48:18–25. 10.1016/j.pbi.2019.01.004. PubMed DOI

Kusuhara  H, Sugiyama  Y. ATP-binding cassette, subfamily G (ABCG family). Pflugers Arch. 2007:453(5):735–744. 10.1007/s00424-006-0134-x. PubMed DOI

Laloi  C, Apel  K, Danon  A. Reactive oxygen signalling: the latest news. Curr Opin Plant Biol. 2004:7(3):323–328. 10.1016/j.pbi.2004.03.005. PubMed DOI

Lane  N, Martin  W. The energetics of genome complexity. Nature. 2010:467(7318):929–934. 10.1038/nature09486. PubMed DOI

Lanzoni  O, Fokin  SI, Lebedeva  N, Migunova  A, Petroni  G, Potekhin  A. Rare freshwater ciliate Paramecium chlorelligerum Kahl, 1935 and its macronuclear symbiotic bacterium “Candidatus Holospora parva.”. PLoS One. 2016:11(12):e0167928. 10.1371/journal.pone.0167928. PubMed DOI PMC

Lesser  MP. Oxidative stress in marine environments: biochemistry and physiological ecology. Annu Rev Physiol. 2006:68(1):253–278. 10.1146/annurev.physiol.68.040104.110001. PubMed DOI

Li  H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018:34(18):3094–3100. 10.1093/bioinformatics/bty191. PubMed DOI PMC

Li  H. Protein-to-genome alignment with miniprot. Bioinformatics. 2023:39(1):btad014. 10.1093/bioinformatics/btad014. PubMed DOI PMC

Liao  Y, Smyth  GK, Shi  W. FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014:30(7):923–930. 10.1093/bioinformatics/btt656. PubMed DOI

Love  MI, Huber  W, Anders  S. Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2. Genome Biol. 2014:15(12):550. 10.1186/s13059-014-0550-8. PubMed DOI PMC

Lowe  CD, Minter  EJ, Cameron  DD, Brockhurst  MA. Shining a light on exploitative host control in a photosynthetic endosymbiosis. Curr Biol.  2016:26(2):207–211. 10.1016/j.cub.2015.11.052. PubMed DOI

Manni  M, Berkeley  MR, Seppey  M, Simão  FA, Zdobnov  EM. BUSCO update: novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. Mol Biol Evol. 2021:38(10):4647–4654. 10.1093/molbev/msab199. PubMed DOI PMC

Marçais  G, Delcher  AL, Phillippy  AM, Coston  R, Salzberg  SL, Zimin  A. MUMmer4: a fast and versatile genome alignment system. PLoS Comput Biol. 2018:14(1):e1005944. 10.1371/journal.pcbi.1005944. PubMed DOI PMC

Martin  C, Smith  AM. Starch biosynthesis. Plant Cell. 1995:7(7):971–985. 10.1105/tpc.7.7.971. PubMed DOI PMC

Meier  R, Wiessner  W. Infection of algae-free Paramecium bursaria with symbiotic Chlorella sp. Isolated from green paramecia. Eur J Protistol. 1988:24(1):69–74. 10.1016/S0932-4739(88)80011-7. PubMed DOI

Mendes  FK, Vanderpool  D, Fulton  B, Hahn  MW. CAFE 5 models variation in evolutionary rates among gene families. Bioinformatics. 2021:36(22-23):5516–5518. 10.1093/bioinformatics/btaa1022. PubMed DOI

Mikheenko  A, Prjibelski  A, Saveliev  V, Antipov  D, Gurevich  A. Versatile genome assembly evaluation with QUAST-LG. Bioinformatics. 2018:34(13):i142–i150. 10.1093/bioinformatics/bty266. PubMed DOI PMC

Miller  G, Schlauch  K, Tam  R, Cortes  D, Torres  MA, Shulaev  V, Dangl  JL, Mittler  R. The plant NADPH oxidase RBOHD mediates rapid systemic signaling in response to diverse stimuli. Sci Signal. 2009:2(84):ra45–ra45. 10.1126/scisignal.2000448. PubMed DOI

Minh  BQ, Schmidt  HA, Chernomor  O, Schrempf  D, Woodhams  MD, von Haeseler  A, Lanfear  R. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020:37(5):1530–1534. 10.1093/molbev/msaa015. PubMed DOI PMC

Mirdita  M, Schütze  K, Moriwaki  Y, Heo  L, Ovchinnikov  S, Steinegger  M. ColabFold: making protein folding accessible to all. Nat Methods. 2022:19(6):679–682. 10.1038/s41592-022-01488-1. PubMed DOI PMC

Moffatt  NSC, Bruinsma  E, Uhl  C, Obermann  WMJ, Toft  D. Role of the cochaperone tpr2 in Hsp90 chaperoning. Biochemistry. 2008:47(31):8203–8213. 10.1021/bi800770g. PubMed DOI

Moradian  MM, Beglaryan  D, Skozylas  JM, Kerikorian  V. Complete mitochondrial genome sequence of three Tetrahymena species reveals mutation hot spots and accelerated nonsynonymous substitutions in Ymf Genes. PLoS One. 2007:2(7):e650. 10.1371/journal.pone.0000650. PubMed DOI PMC

Morris  JL, Puttick  MN, Clark  JW, Edwards  D, Kenrick  P, Pressel  S, Wellman  CH, Yang  Z, Schneider  H, Donoghue  PCJ. The timescale of early land plant evolution. Proc Natl Acad Sci U S A.  2018:115(10):E2274–E2283. 10.1073/pnas.1719588115. PubMed DOI PMC

Moy  GW, Mendoza  LM, Schulz  JR, Swanson  WJ, Glabe  CG, Vacquier  VD. The sea urchin sperm receptor for egg jelly is a modular protein with extensive homology to the human polycystic kidney disease protein, PKD1. J Cell Biol. 1996:133(4):809–817. 10.1083/jcb.133.4.809. PubMed DOI PMC

Muller  F. The nature and mechanism of superoxide production by the electron transport chain: its relevance to aging. J Am Aging Assoc. 2000:23(4):227–253. 10.1007/s11357-000-0022-9. PubMed DOI PMC

Nakjang  S, Williams  TA, Heinz  E, Watson  AK, Foster  PG, Sendra  KM, Heaps  SE, Hirt  RP, Martin Embley  T. Reduction and expansion in microsporidian genome evolution: new insights from comparative genomics. Genome Biol Evol. 2013:5(12):2285–2303. 10.1093/gbe/evt184. PubMed DOI PMC

Oliver  PL, Finelli  MJ, Edwards  B, Bitoun  E, Butts  DL, Becker  EBE, Cheeseman  MT, Davies  B, Davies  KE. Oxr1 is essential for protection against oxidative stress-induced neurodegeneration. PLoS Genet. 2011:7(10):e1002338. 10.1371/journal.pgen.1002338. PubMed DOI PMC

Ou  S, Su  W, Liao  Y, Chougule  K, Agda  JRA, Hellinga  AJ, Lugo  CSB, Elliott  TA, Ware  D, Peterson  T, et al.  Benchmarking transposable element annotation methods for creation of a streamlined, comprehensive pipeline. Genome Biol. 2019:20(1):275. 10.1186/s13059-019-1905-y. PubMed DOI PMC

Pan  G, Xu  J, Li  T, Xia  Q, Liu  S-L, Zhang  G, Li  S, Li  C, Liu  H, Yang  L, et al.  Comparative genomics of parasitic silkworm microsporidia reveal an association between genome expansion and host adaptation. BMC Genomics. 2013:14(1):186. 10.1186/1471-2164-14-186. PubMed DOI PMC

Pardy  RL, Spargo  B, Crowe  JH. Release of trehalose by symbiotic algae. Symbiosis. 1989:7:149–158.

Pettersen  EF, Goddard  TD, Huang  CC, Meng  EC, Couch  GS, Croll  TI, Morris  JH, Ferrin  TE. UCSF ChimeraX: structure visualization for researchers, educators, and developers. Protein Sci.  2021:30(1):70–82. 10.1002/pro.3943. PubMed DOI PMC

Pitsch  G, Adamec  L, Dirren  S, Nitsche  F, Šimek  K, Sirová  D, Posch  T. The green Tetrahymena utriculariae n. sp. (Ciliophora, Oligohymenophorea) with its endosymbiotic algae (Micractinium sp.), living in traps of a carnivorous aquatic plant. J Eukaryot Microbiol.  2017:64(3):322–335. 10.1111/jeu.12369. PubMed DOI

Prescott  DM. The DNA of ciliated protozoa. Microbiol Rev. 1994:58(2):233–267. 10.1128/mr.58.2.233-267.1994. PubMed DOI PMC

Proost  S, Fostier  J, De Witte  D, Dhoedt  B, Demeester  P, Van de Peer  Y, Vandepoele  K. . i-ADHoRe 3.0—fast and sensitive detection of genomic homology in extremely large data sets. Nucleic Acids Res. 2012:40(2):e11–e11. 10.1093/nar/gkr955. PubMed DOI PMC

Pröschold  T, Darienko  T, Silva  PC, Reisser  W, Krienitz  L. The systematics of Zoochlorella revisited employing an integrative approach. Environ Microbiol. 2011:13(2):350–364. 10.1111/j.1462-2920.2010.02333.x. PubMed DOI

Pröschold  T, Rieser  D, Darienko  T, Nachbaur  L, Kammerlander  B, Qian  K, Pitsch  G, Bruni  EP, Qu  Z, Forster  D, et al.  An integrative approach sheds new light onto the systematics and ecology of the widespread ciliate genus Coleps (Ciliophora, Prostomatea). Sci Rep. 2021:11(1):5916. 10.1038/s41598-021-84265-y. PubMed DOI PMC

Puttick  MN. MCMCtreeR: functions to prepare MCMCtree analyses and visualize posterior ages on trees. Bioinformatics. 2019:35(24):5321–5322. 10.1093/bioinformatics/btz554. PubMed DOI

Raitt  DC. Yeast Cdc42 GTPase and Ste20 PAK-like kinase regulate Sho1-dependent activation of the Hog1 MAPK pathway. EMBO J. 2000:19(17):4623–4631. 10.1093/emboj/19.17.4623. PubMed DOI PMC

Rambaut  A, Drummond  AJ, Xie  D, Baele  G, Suchard  MA. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst Biol. 2018:67(5):901–904. 10.1093/sysbio/syy032. PubMed DOI PMC

Roach  MJ, Schmidt  SA, Borneman  AR. Purge Haplotigs: allelic contig reassignment for third-gen diploid genome assemblies. BMC Bioinformatics. 2018:19(1):460. 10.1186/s12859-018-2485-7. PubMed DOI PMC

Ronquist  F, Teslenko  M, van der Mark  P, Ayres  DL, Darling  A, Höhna  S, Larget  B, Liu  L, Suchard  MA, Huelsenbeck  JP. Mrbayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012:61(3):539–542. 10.1093/sysbio/sys029. PubMed DOI PMC

Šimek  K, Pitsch  G, Salcher  MM, Sirová  D, Shabarova  T, Adamec  L, Posch  T. Ecological traits of the algae-bearing Tetrahymena utriculariae (Ciliophora) from traps of the aquatic carnivorous plant Utricularia reflexa. J Eukaryot Microbiol.  2017:64(3):336–348. 10.1111/jeu.12368. PubMed DOI

Siozios  S, Jimenez  PN, Azagi  T, Sprong  H, Frost  CL, Parratt  SR, Taylor  G, Brettell  L, Liew  KC, Croft  L, et al.  Genome dynamics across the evolutionary transition to endosymbiosis. Curr Biol  2024:34(24):5659–5670.e7. 10.1016/j.cub.2024.10.044. PubMed DOI

Sirová  D, Adamec  L, Vrba  J. Enzymatic activities in traps of four aquatic species of the carnivorous genus Utricularia. New Phytol.  2003:159(3):669–675. 10.1046/j.1469-8137.2003.00834.x. PubMed DOI

Sirová  D, Borovec  J, Picek  T, Adamec  L, Nedbalová  L, Vrba  J. Ecological implications of organic carbon dynamics in the traps of aquatic carnivorous Utricularia plants. Funct Plant Biol.  2011:38(7):583–593. 10.1071/FP11023. PubMed DOI

Sirová  D, Kreidlová  V, Adamec  L, Vrba  J. The ability of Tetrahymena utriculariae (Ciliophora, Oligohymenophorea) to colonize traps of different species of aquatic carnivorous Utricularia. J Eukaryot Microbiol.  2020:67(5):608–611. 10.1111/jeu.12812. PubMed DOI

Spanaki  C, Plaitakis  A. The role of glutamate dehydrogenase in mammalian ammonia metabolism. Neurotox Res. 2012:21(1):117–127. 10.1007/s12640-011-9285-4. PubMed DOI

Spanner  C, Darienko  T, Filker  S, Sonntag  B, Pröschold  T. Morphological diversity and molecular phylogeny of five Paramecium bursaria (Alveolata, Ciliophora, Oligohymenophorea) syngens and the identification of their green algal endosymbionts. Sci Rep. 2022:12(1):18089. 10.1038/s41598-022-22284-z. PubMed DOI PMC

Stamatakis  A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014:30(9):1312–1313. 10.1093/bioinformatics/btu033. PubMed DOI PMC

Stanke  M, Diekhans  M, Baertsch  R, Haussler  D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008:24(5):637–644. 10.1093/bioinformatics/btn013. PubMed DOI

Strassert  JFH, Irisarri  I, Williams  TA, Burki  F. A molecular timescale for eukaryote evolution with implications for the origin of red algal-derived plastids. Nat Commun. 2021:12(1):1879. 10.1038/s41467-021-22044-z. PubMed DOI PMC

Susin  MF, Baldini  RL, Gueiros-Filho  F, Gomes  SL. GroES/GroEL and DnaK/DnaJ have distinct roles in stress responses and during cell cycle progression in Caulobacter crescentus. J Bacteriol. 2006:188(23):8044–8053. 10.1128/JB.00824-06. PubMed DOI PMC

Szabo  A, Korszun  R, Hartl  FU, Flanagan  J. A zinc finger-like domain of the molecular chaperone DnaJ is involved in binding to denatured protein substrates. EMBO J. 1996:15(2):408–417. 10.1002/j.1460-2075.1996.tb00371.x. PubMed DOI PMC

Takahashi  M, Sasaki  Y, Ida  S, Morikawa  H. Nitrite reductase gene enrichment improves assimilation of NO2 in Arabidopsis. Plant Physiol. 2001:126(2):731–741. 10.1104/pp.126.2.731. PubMed DOI PMC

Tibbs  VC, Gray  PC, Catterall  WA, Murphy  BJ. AKAP15 Anchors cAMP-dependent protein kinase to brain sodium channels. J Biol Chem.  1998:273(40):25783–25788. 10.1074/jbc.273.40.25783. PubMed DOI

Trotter  KW, Fraser  IDC, Scott  GK, Stutts  MJ, Scott  JD, Milgram  SL. Alternative splicing regulates the subcellular localization of a-kinase anchoring protein 18 Isoforms. J Cell Biol. 1999:147(7):1481–1492. 10.1083/jcb.147.7.1481. PubMed DOI PMC

Tsai  M-F, Wang  C-C, Chang  G-C, Chen  C-Y, Chen  H-Y, Cheng  C-L, Yang  Y-P, Wu  C-Y, Shih  F-Y, Liu  C-C, et al.  A new tumor suppressor DnaJ-like heat shock protein, HLJ1, and survival of patients with non-small-cell lung carcinoma. J Natl Cancer Inst.  2006:98(12):825–838. 10.1093/jnci/djj229. PubMed DOI

Viéitez  C, Martínez-Cebrián  G, Solé  C, Böttcher  R, Potel  CM, Savitski  MM, Onnebo  S, Fabregat  M, Shilatifard  A, Posas  F, et al.  A genetic analysis reveals novel histone residues required for transcriptional reprogramming upon stress. Nucleic Acids Res. 2020:48(7):3455–3475. 10.1093/nar/gkaa081. PubMed DOI PMC

Volkert  MR, Elliott  NA, Housman  DE. Functional genomics reveals a family of eukaryotic oxidation protection genes. Proc Natl Acad Sci U S A.  2000:97(26):14530–14535. 10.1073/pnas.260495897. PubMed DOI PMC

Winter  J, Jakob  U. Beyond transcription—new mechanisms for the regulation of molecular chaperones. Crit Rev Biochem Mol Biol. 2004:39(5-6):297–317. 10.1080/10409230490900658. PubMed DOI

Wood  DE, Lu  J, Langmead  B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019:20(1):257. 10.1186/s13059-019-1891-0. PubMed DOI PMC

Wu  T, Hu  E, Xu  S, Chen  M, Guo  P, Dai  Z, Feng  T, Zhou  L, Tang  W, Zhan  L, et al.  ClusterProfiler 4.0: a universal enrichment tool for interpreting omics data. The Innovation. 2021:2(3):100141. 10.1016/j.xinn.2021.100141. PubMed DOI PMC

Xian  L, Zhang  Y, Cao  Y, Wan  T, Gong  Y, Dai  C, Ochieng  WA, Nasimiyu  AT, Li  W, Liu  F. Glutamate dehydrogenase plays an important role in ammonium detoxification by submerged macrophytes. Sci Total Environ.  2020:722:137859. 10.1016/j.scitotenv.2020.137859. PubMed DOI

Xu  Y, Zhu  Y, Wang  Y, Chang  Y-F, Zhang  Y, Jiang  X, Zhuang  X, Zhu  Y, Zhang  J, Zeng  L, et al.  Whole genome sequencing revealed host adaptation-focused genomic plasticity of pathogenic Leptospira. Sci Rep. 2016:6(1):20020. 10.1038/srep20020. PubMed DOI PMC

Yadav  V, Molina  I, Ranathunge  K, Castillo  IQ, Rothstein  SJ, Reed  JW. ABCG transporters are required for suberin and pollen wall extracellular barriers in Arabidopsis. Plant Cell. 2014:26(9):3569–3588. 10.1105/tpc.114.129049. PubMed DOI PMC

Yang  Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007:24(8):1586–1591. 10.1093/molbev/msm088. PubMed DOI

Yao  M-C, Zheng  K, Yao  C-H. A conserved nucleotide sequence at the sites of developmentally regulated chromosomal breakage in Tetrahymena. Cell. 1987:48(5):779–788. 10.1016/0092-8674(87)90075-4. PubMed DOI

Yokobayashi  K, Misaki  A, Harada  T. Purification and properties of Pseudomonas isoamylase. Biochim Biophys Acta. 1970:212(3):458–469. 10.1016/0005-2744(70)90252-4. PubMed DOI

Yoshida  K, Komae  K. A rice family 9 glycoside hydrolase isozyme with broad substrate specificity for hemicelluloses in type II cell walls. Plant Cell Physiol. 2006:47(11):1541–1554. 10.1093/pcp/pcl020. PubMed DOI

Zhang  C, Mirarab  S. Weighting by gene tree uncertainty improves accuracy of quartet-based species trees. Mol Biol Evol. 2022:39(12):msac215. 10.1093/molbev/msac215. PubMed DOI PMC

Zhu  A, Ibrahim  JG, Love  MI. Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics. 2019:35(12):2084–2092. 10.1093/bioinformatics/bty895. PubMed DOI PMC