DOI: 10.5281/zenodo.7807739
title: ATP synthase evolution on a cross-braced dated tree of life

description: <strong>Abstract</strong> The timing of early cellular evolution from the divergence of Archaea and Bacteria to the origin of eukaryotes remains poorly constrained. The ATP synthase complex is thought to have originated prior to the Last Universal Common Ancestor (LUCA) and analyses of ATP synthase genes, together with ribosomes, have played a key role in inferring and rooting the tree of life. Here we reconstruct the evolutionary history of ATP synthases using an expanded sampling of Archaea, Bacteria, and eukaryotes. We developed a phylogenetic cross-bracing approach making use of endosymbioses and ancient gene duplications of the major ATP synthase subunits to infer a highly resolved, dated species tree and establish an absolute timeline for ATP synthase evolution. Our analyses show that the divergence of the ATP synthase into F- and A/V-type lineages, was a very early event in cellular evolution dating back to more than 4Ga potentially predating the diversification of Archaea and Bacteria. Our cross-braced, dated tree of life also provides insight into more recent evolutionary transitions including eukaryogenesis, showing that the eukaryotic nuclear and mitochondrial lineages diverged from their closest archaeal (2.67-2.19Ga) and bacterial (2.58-2.12Ga) relatives at roughly the same time, with the nuclear stem being moderately longer. <strong>Repository Contents</strong> <strong>1_100Eukaryote_genomes.tar.gz</strong>: includes all protein sequence files for the 100 Eukaryotes sampled in this study. <strong>2_Phylogenies.tar.gz</strong>: includes all files used for phylogenetic analyses. Folders are organized as follows: <strong>1_ATPsynthase_gene_trees</strong>: this folder contains all sequence, alignment, and tree files for the ATP synthase gene trees. Files are organized as follows and are associated with the corresponding parts of the manuscript: Figure 3, Figure 5B, Supplementary Figures 5-10, Supplementary Figures 18-19 Folder '1_sequences' includes all unaligned fasta sequence files for each ATP synthase gene tree (see Methods) Folder ‘2_alignments’ includes all alignments generated using MAFFT L-INS-i (subdirectory: 1_untrimmed) and trimmed with BMGE (subdirectory: 2_trimmed) Folder '3_treefiles' includes all IQ-TREE2 output files for all ATP synthase gene phylogenies. Any files with suffix *taxa.treefile contain the full taxonomic string for each accession. Folder '4_pdfs' includes PDF files for each ATP synthase gene tree <strong>2_Eukaryotic_subsets</strong>: this folder contains all sequence, alignment, and tree files for ATP synthase Eukaryotic subset gene trees. Files are organized as follows and are associated with the corresponding parts of the manuscript: Supplementary Figure 11 Folder ‘1_sequences’ includes all unaligned fasta sequence files for the eukaryotic subsets. Folder ‘2_alignments’ includes all alignments generated using MAFFT L-INS-i (subdirectory: 1_untrimmed) and trimmed with BMGE (subdirectory: 2_trimmed). Folder ‘3_treefiles’ includes all Bayesian trees inferred for eukaryotic subsets. Folder '4_pdfs' includes PDF files for each eukaryotic subset tree <strong>3_21eLife_concatenated_species_tree</strong>: this folder contains all sequence, alignment, and tree files for the single gene tree and concatenated phylogeny analyses (inferred using 21 single-copy marker genes, see Methods). Files are organized as follows and are associated with the following parts of the manuscript: Figure 1, Supplementary Figure 20 Folder ‘1_inspection_start’ corresponds to the initial manual inspection of the single gene trees and includes the following subdirectories: Folder ‘1_sequences’ includes all protein sequence fasta files corresponding to the 27 original single-copy marker genes Folder ‘2_alignments’ includes all alignment files generated using MAFFT L-INS-i (subdirectory: 1_untrimmed) and trimmed with BMGE (subdirectory: 2_untrimmed) Folder ‘3_treefiles’ includes all IQ-TREE2 output files for all phylogenies (27 single-copy marker genes) Folder '4_pdfs' includes PDF files for each single gene tree Folder ‘2_inspection_final’ corresponds to the final manual inspection of the single gene trees and includes the following subdirectories: Folder ‘1_sequences’ includes all protein sequence fasta files corresponding to the final 21 single-copy marker genes Folder ‘2_alignments’ includes all alignment files generated using MAFFT L-INS-i (subdirectory: 1_untrimmed) and trimmed with BMGE (subdirectory: 2_untrimmed) Folder ‘3_treefiles’ includes all IQ-TREE2 output files for all phylogenies (21 single-copy marker genes) Folder '4_pdfs' includes PDF files for each single gene tree Folder ‘3_concatenated_phylogeny’ contains concatenated alignment generated from the final 21 single-copy marker gene alignments Folder ‘1_alignment’ includes the concatenated alignment generated from the 21 trimmed alignments from the final inspection Folder ‘2_treefiles’ includes all IQ-TREE2 output files for trees inferred using the two different models (subdirectories: LG+C20+R+F and LG+C60+R+F) Folder '4_Eukaryote_only_phylogeny' contains sequence, alignment, and tree files for 21 single-copy marker genes used to infer a Eukaryote-only phylogeny. Folder is organized as follows and files correspond to Supplementary Figure 3: Folder ‘1_sequences’ includes all protein sequence fasta files corresponding to the 21 single-copy marker genes with only Eukaryotes Folder ‘2_alignments’ includes all alignment files generated using MAFFT L-INS-i (subdirectory: 1_untrimmed) and trimmed with BMGE (subdirectory: 2_untrimmed) Folder ‘3_concatenated_phylogeny’ includes concatenated alignment generated from 21 single-copy markers with only Eukaryotes (subdirectory: 1_alignment) and all IQ-TREE2 output files for the concatenated phylogeny (subdirectory: 2_treefiles) Folder '4_pdfs' includes PDF files for the concatenated Eukaryote tree <strong>4_Ribosomal_species_tree</strong>: this folder contains all sequence, alignment, and tree files for the single gene tree and concatenated phylogeny analyses (inferred using 12 ribosomal marker genes, see Methods). Files are organized as follows and are associated with the corresponding parts of the manuscript: Figure 5A, Figure 5C, Supplementary Figures 12-16, Supplementary Figure 21 Folder ‘1_sequences’ includes all protein sequence fasta files for the original 15 ribosomal proteins. Sequence sets include the best-hit Archaea and Bacteria, and nuclear, mitochondrial, and plastid eukaryotic homologs Folder ‘2_alignments’ includes all alignment files generated using MAFFT L-INS-i (subdirectory: 1_untrimmed) and trimmed with TRIMAL (gappy-out) (subdirectory: 2_trimmed) Folder ‘3_treefiles’ includes all original FastTree tree files, tree files with highlighted sequences to remove (*blue-to-rem = eukaryotic nuclear homolog only; *colored-to-rem = eukaryotic nuclear, mitochondrial, and plastid homologs). PDFs of each marker gene tree are also included that depict highlighting of sequences to keep and/or remove. Folder ‘4_concatenated_phylogeny’ contains concatenated alignment generated from the final 12 ribosomal marker genes Folder ‘1_alignment’ includes the concatenated alignment generated with 12 ribosomal marker proteins in MAFFT L-INS-i and trimmed with TRIMAL (gappy-out) Folder ‘2_phylogeny’ includes all IQ-TREE2 output files for the species tree inferred using the LG+C60+R+F model <strong>5_Dating_analysis</strong>: includes all Mcmcdate output files for the dating analyses (species tree and ATP synthase gene tree, see Methods). Folder '1_Edited1_dating' includes all dated tree files and monitor files for braced and unbraced analyses of the Edited1 species tree topology. Data corresponds to Supplementary Figure 12, Supplementary Figure 14-15 Folder '2_Edited2_dating' includes all dated tree files and monitor files for braced and unbraced analyses of the Edited2 (focal) species tree topology. Data corresponds to Figure 5A, Figure 5C, Supplementary Figure 13, Supplementary Figure 16. Folder '3_ATP_synthase_dating' includes all dated tree files and monitor files for braced and unbraced analyses of the ATP synthase gene tree. Data corresponds to Figure 5B, Supplementary Figures 18-19). <strong>3_Scripts.tar.gz</strong>: includes all workflows and scripts used for phylogenetic analyses. <strong>1_workflows</strong>: includes bash workflows for phylogenetic analyses (details on software versions are included in each workflow summary): Workflow_ATPsynthase_gene_trees.sh: generation of the ATP synthase phylogenies Workflow_21eLife_marker_phylogeny.sh: inferring the 21 marker-gene species tree Workflow_Ribosomal_species_tree.sh: inferring the 12 ribosomal marker-gene species tree Workflow_Database_annotations.sh: workflow for gene annotation for 800 sampled Archaea, Bacteria, and Eukaryota <strong>2_R_scripts</strong>: includes R scripts used to plot tree figures (Supplementary Figures 5-10) and for the presence-absence analyses (Figure 1, Figure 2, Supplementary Figure 2). Input mapping files are included Folder ‘1_Presence_absence’ includes the treefile used for ordering the plots (‘1_tree’), the taxonomic and COG mapping files (‘2_input_files’), and script 'Presence-absence.R' Folder ‘2_plotting_trees’ includes the rectangular and radial trees generated for each ATP synthase trees (see Supplementary Figures 5-10). Trees were generated from the treefiles for the ATP synthase gene trees (see above), and script 'Plotting_trees.Rmd' 'Marker_gene_counts.R' script used to count marker genes per genome (see Methods)
description: This work was supported by the Gordon and Betty Moore Foundation through grant GBMF9741 to TAW, AS, and GJSz. Furthermore, A.S. has received funding from the Swedish Research Council (VR starting grant 2016-03559), the NWO-I foundation of the Netherlands Organisation for Scientific Research (WISE fellowship) and the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 947317, ASymbEL). Further, this work was supported by a Royal Society University Research Fellowship to TAW. ERRM, DP, PCJD and TAW were supported by the John Templeton Foundation (62220). The opinions expressed in this publication are those of the author(s) and do not necessarily reflect the views of the John Templeton Foundation. GJSz, DS and LLSz received funding from the European Union's Horizon 2020 research and innovation programme (grant agreement No. 714774, GENECLOCKS).