# Discovery of archaeal fusexins homologous to eukaryotic HAP2/GCS1 gamete fusion proteins

David Moi, Shunsuke Nishio, Xiaohui Li, Clari Valansi, Mauricio Langleib, Nicolas G. Brukman, Kateryna Flyak, Christophe Dessimoz, Daniele de Sanctis, Kathryn Tunyasuvunakool, et al.

► To cite this version:

David Moi, Shunsuke Nishio, Xiaohui Li, Clari Valansi, Mauricio Langleib, et al. . Discovery of archaeal fusexins homologous to eukaryotic HAP2/GCS1 gamete fusion proteins. Nature Communications, 2022, 13 (1), pp.3880. 10.1038/s41467-022-31564-1 . hal-05075972

HAL Id: hal-05075972

https://cnrs.hal.science/hal-05075972v1

Submitted on 20 May 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-022-31564-1

OPEN

<figure><img src="image_3.png" /></figure>

Check for updates

# Discovery of archaeal fusexins homologous to eukaryotic HAP2/GCS1 gamete fusion proteins

David Moi ID 1,2,3,14, Shunsuke Nishio ID 4,14, Xiaohui Li 5,14, Clari Valansi 5, Mauricio Langleib 6,7, Nicolas G. Brukman 5, Kateryna Flyak ID 5, Christophe Dessimoz 2,3,8,9, Daniele de Sanctis 10, Kathryn Tunyasuvunakool ID 11, John Jumper ID 11, Martin Graña ID 7✉, Héctor Romero ID 6,12✉, Pablo S. Aguilar ID 1,13✉, Luca Jovine ID 4✉ & Benjamin Podbilewicz ID 5✉

Sexual reproduction consists of genome reduction by meiosis and subsequent gamete fusion. The presence of genes homologous to eukaryotic meiotic genes in archaea and bacteria suggests that DNA repair mechanisms evolved towards meiotic recombination. However, fusogenic proteins resembling those found in gamete fusion in eukaryotes have so far not been found in prokaryotes. Here, we identify archaeal proteins that are homologs of fusexins, a superfamily of fusogens that mediate eukaryotic gamete and somatic cell fusion, as well as virus entry. The crystal structure of a trimeric archaeal fusexin (Fusexin1 or Fsx1) reveals an archetypical fusexin architecture with unique features such as a six-helix bundle and an additional globular domain. Ectopically expressed Fusexin1 can fuse mammalian cells, and this process involves the additional globular domain and a conserved fusion loop. Furthermore, archaeal fusexin genes are found within integrated mobile elements, suggesting potential roles in cell-cell fusion and gene exchange in archaea, as well as different scenarios for the evolutionary history of fusexins.

<footnote n="1">¹Instituto de Fisiología, Biología Molecular y Neurociencias (IFIBYNE-CONICET), Buenos Aires, Argentina. ²Department of Computational Biology, University of Lausanne, Lausanne, Switzerland. ³Swiss Institute of Bioinformatics, Lausanne, Switzerland. ⁴Department of Biosciences and Nutrition, Karolinska Institutet, Huddinge, Sweden. ⁵Department of Biology, Technion- Israel Institute of Technology, Haifa, Israel. ⁶Unidad de Genómica Evolutiva, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay. ⁷Unidad de Bioinformática, Institut Pasteur de Montevideo, Montevideo, Uruguay.</footnote>

<footnote n="8">⁸Department of Genetics, Evolution and Environment, Centre for Life's Origins and Evolution, University College London, London, UK. ⁹Department of Computer Science, University College London, London, UK. ¹⁰ESRF—The European Synchrotron, Grenoble, France. ¹¹DeepMind, London, UK. ¹²Centro Universitario Regional Este - CURE, Centro Interdisciplinario de Ciencia de Datos y Aprendizaje Automático - CICADA, Universidad de la República, Montevideo, Uruguay. ¹³Instituto de Investigaciones Biotecnológicas Universidad Nacional de San Martín (IIB-CONICET), San Martín, Buenos Aires, Argentina. ¹⁴These authors contributed equally: David Moi, Shunsuke Nishio, Xiaohui Li. ✉email: mgrana@pasteur.edu.uy; eletor@fcien.edu.uy; paguilar@iib.unsam.edu.ar; luca.jovine@ki.se; podbilew@technion.ac.il</footnote>
How early eukaryotes developed the capacity for gamete fusion is a central question entangled with the origin of the eukaryotic cell itself. The widespread presence of a conserved set of meiosis, gamete, and nuclear fusion proteins among extant eukaryotes indicates that meiotic sex emerged once, predating the last eukaryotic common ancestor (LECA)¹,². Two essential molecular events are required for meiotic sex: DNA recombination and plasma membrane fusion. Prokaryotic cells contain DNA repair machines that may have been precursors of the recombination machinery used in eukaryotic meiosis¹,². However, the genes encoding for proteins that are essential and sufficient to merge plasma membranes have not been identified in prokaryotes³.

In eukaryotes, different families of cellular and viral fusion proteins (fusogens) have been described³. For example, class I viral fusogens include the spike glycoproteins of Influenza, HIV, Ebola, and SARS-CoV that have similar structures that appear to have converged during evolution as a way to merge viral and eukaryotic membranes⁴⁻⁶. More recently it was shown that myoblast fusion requires two unrelated proteins to form muscles in vertebrates⁷,⁸. However, many fusogens have not been identified yet, and the molecular basis of gamete fusion in fungi and vertebrates remain unclear⁹.

The first eukaryotic cell-cell fusogen was discovered only two decades ago by genetic screens in *C. elegans*¹⁰. The fusogen EFF-1 and its paralog AFF-1 fuse one-third of all the somatic cells in the skin, excretory, reproductive, nervous and digestive systems of nematodes. The functions of these cellular fusions are to sculpt cells, tissues, and organs to restrict cellular fates for a robust development¹¹. In the absence of sequence similarity between EFF-1, AFF-1, and other known proteins, predictions of the structure of AFF-1 suggested structural similarity to class II viral fusogens¹². EFF-1 and AFF-1 can fuse cells in *C. elegans*, promote fusion between heterologous cells, and substitute a viral fusogen to mediate plasma membrane fusion only when these fusogenic proteins are expressed in both fusing membranes¹³⁻¹⁵. The crystal structure of the ectodomain of EFF-1 demonstrated structural similarity with class II viral fusogens such as the glycoproteins on the surface of Zika, dengue, and rubella viruses¹⁶. Thus, these eukaryotic and viral fusogens have remarkably similar functions and structures despite undetectable sequence similarity. However, their mechanisms of membrane fusion are different because EFF-1 and AFF-1 use a bilateral mechanism while class II viral fusogens use a unilateral mechanism¹³⁻¹⁵.

More recently, genetic screens uncovered a protein involved in gamete fusion, HAP2/GCS1, which is conserved in *Arabidopsis*, *Chlamydomonas*, *Plasmodium*, *Tetrahymena*, and *Dictyostelium*¹⁷⁻²¹. Later, structural bioinformatics, crystallographic structure elucidation and functional assays demonstrated that HAP2/GCS1 is a bona fide fusogen homologous to EFF-1, AFF-1, and class II viral fusogens²²⁻²⁴. The crystal structures of HAP2/GCS1 from *Chlamydomonas* and *Arabidopsis* showed remarkable structural conservation without sequence similarity²³,²⁵,²⁶. This superfamily of fusion proteins was named fusexins: fusion proteins essential for sexual reproduction and exoplasmic merger of plasma membranes²²⁻²⁴. Thus, the fusexin superfamily encompasses class II viral fusogens (viral fusexins) that fuse the envelope of some animal viruses with the membranes of host cells during infection⁵,⁶; EFF-1 and AFF-1 (somatic fusexins) that promote cell fusion during syncytial organ development¹⁰,¹³,¹⁵,¹⁶; and HAP2/GCS1 (sexual fusexins; hereafter referred to as HAP2) that mediate gamete fusion¹⁷⁻¹⁹.

Although it is assumed that sexual fusexins were already present in the LECA¹,²⁷, their shared ancestry with viral fusexins posed a “the virus or the egg” evolutionary dilemma²²,²⁴,²⁸. In one scenario, fusexins are proper eukaryal innovations that were

captured by some viruses and used for host invasion. Alternatively, a viral fusexin gene was transferred to an early eukaryotic cell and then repurposed for gamete fusion. Solving this evolutionary conundrum is not a trivial task because sequence-based phylogenetics cannot be applied to the whole fusexin superfamily due to lack of sequence conservation.

Here we identify a family of fusexins in genomes of Archaea and prokaryotic fractions of metagenomes from very diverse environments. We provide crystallographic and functional evidence suggesting that these proteins are cellular fusogens. Genomic analyses show that archaeal fusexins are carried by integrated mobile genetic elements. Evolutionary analyses of the whole fusexin superfamily reveal alternative working models regarding the relationships between viral, eukaryotic and archaeal fusexins and the emergence of meiotic sex during eukaryogenesis.

## Results

**Fusexin genes in Archaea.** To search for fusexins we used the crystallographic structures of *C. reinhardtii*, *A. thaliana*, and *T. cruzi* HAP2 (Cr/At/TcHAP2)²³,²⁵,²⁶ to build dedicated Hidden Markov Models (HMMs) for scanning the Uniclust30 database. We detected 24 high-confidence candidates in prokaryotes: 8 belong to isolated and cultivated archaea, and the remaining 16 to metagenome-assembled genomes (MAGs, Supplementary Table 1). We then built HMMs of the candidate ectodomains and compared them to HMMs of sexual, somatic, and viral fusexins. Figure 1a shows that the prokaryotic candidates have detectable sequence similarities with HAP2, with *E*-values below 0.001 and HHblits-derived probabilities higher than 0.95 (Supplementary Fig. 1a). We named these proteins Fusexin1 (Fx1). *fxs1* genes found in cultivated and isolated prokaryotes are restricted to the Haloarchaea class (also called Halobacteria, Euryarchaeota superphylum) whereas MAGs containing Fsx1s include all major Archaea superphyla (Supplementary Table 1). Next, we used this Fsx1 sequence set to search the Metacust database, which comprises 1.59 billion clustered proteins from over 2200 metagenomic/metatranscriptomic datasets. Performing a scan pipeline using PHMMER, PSI-BLAST, HMM-HMM comparisons and topology filtering we found 96 high-confidence *fxs1* genes. The identified *fxs1s* come from different environments (with preeminence of saline samples) and a wide temperature range (−35 to 80 °C, see Supplementary Data 1).

**Fx1s is a structural homolog of HAP2/GCS1.** To experimentally investigate the presence of fusexin-like proteins in Archaea, a selection of the candidate genes was expressed in mammalian cells (Supplementary Fig. 1b, c). High-level expression was observed for a metagenomic Fsx1 sequence from a hypersaline environment, predicted to encode a ~55 kDa ectodomain region (Fx1E) followed by three transmembrane domains (TMs) (Supplementary Data 1). Fsx1E is a monomer in solution but crystallized as a homotrimer in the presence of 2.5 M NaCl, 0.2 M CaCl₂ (Supplementary Fig. 2). These conditions precluded experimental phasing, and attempts to phase the data by molecular replacement (MR) with different kinds of homology models also failed, due to insufficient sequence identity to known fusexin structures. However, we could determine the structure of Fsx1E at 2.3 Å resolution by running MR with a combination of fragments from ab initio predictions generated by AlphaFold²⁹ (Fig. 1b, Supplementary Figs. 3, 4 and Supplementary Table 2).

The Fsx1E homotrimer has overall dimensions of 119 × 77 × 76 Å (Fig. 1b). Each protomer consists of four domains (Fig. 2a, b), the first three of which match the approximate dimensions and relative arrangement of domains I–III of known fusexins in their post-fusion conformation³⁰; accordingly, fold
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Fsx1 is a member of the fusexin protein superfamily. **a** HMM vs. HMM homology probabilities of a subset of eukaryotic, viral, and archaeal fusexin ectodomains. With exception of Fsx1, which derives from a metagenomic sequence, archaeal fusexins (red), viral fusexins (blue), EFF-1 (green) and HAP2s (black) are indicated by RefSeq/UniProt identifiers. **b** Crystal structure of the trimeric ectodomain of Fsx1. The three-fold non-crystallographic axis is indicated. Subunit A is shown as a cartoon colored by domain, with disulfides and the fusion loop (FL) colored magenta and orange, respectively; subunits B, C are in mixed cartoon/surface representation.</figcaption></figure>

and interface similarity searches identify HAP2 as the closest structural homolog of Fsx1E, followed by viral fusexins and C. elegans EFF-1 (Fig. 2c, d). Fsx1 domains I and III are relatively sequence-conserved among archaeal homologs (Supplementary Figs. 5a and 6) and closely resemble the corresponding domains of HAP2 (RMSD 2.1 Å over 218 Ca), including the invariant disulfide bond between domain III strands βC and βF²³ (C₃₃₈₉–C₄₃₂₄; Fig. 2, Supplementary Fig. 4c). On the other hand, Fsx1 domain II shares the same topology as that of HAP2 but differs significantly in its secondary structure elements and their relative orientation, as well as disulfide bonds (Fig. 2c). In particular, Fsx1 domain II is characterized by a four-helix hairpin, whose N-terminal half interacts with the same region of the other two subunits to generate a six-helix bundle around the molecule's three-fold axis (Figs. 1b and 3a–c and Supplementary Figs. 4a and 5b, c).

Notably, unlike previously characterized viral and eukaryotic fusexins, Fsx1 also contains a fourth globular domain conserved among archaeal homologs (Figs. 1b, 2 and Supplementary Figs. 4d and 6). Its antiparallel β-sandwich, which includes the two C-terminal disulfides of Fsx1, resembles the carbohydrate-binding fold of dust mite allergen Der p 23 and related chitin-binding proteins³¹ (Fig. 3d); accordingly, it is also structurally similar to a high-confidence AlphaFold2 model of the C-terminal domain of acidic mammalian chitinase³². In addition to being coaxially stacked with domain III as a result of a loop/loop interaction stabilized by the C₅₄₅₇–C₆₄₇₇ disulfide, domain IV contributes to the quaternary structure of the protein by interacting with domain II of the adjacent subunit to which domain III also binds (Figs. 1b and 3c).

The Fsx1E monomer has a net charge of –67, and another feature stabilizing its homotrimeric assembly is a set of Ca²⁺ and Na⁺ ions that interacts with negatively charged residues at the interface between subunits (Fig. 3a–c and Supplementary Figs. 4a and 5b). Additional metal ions bind to sites located within individual subunits; in particular, a Ca²⁺ ion shapes the conformation of the domain II cd loop (S143-V148) so that its uncharged surface protrudes from the rest of the molecule (Fig. 3e, f and Supplementary Fig. 5d). Strikingly, the position of

this element matches that of the fusion loops (FLs) of other fusexins, including the Ca²⁺-binding fusion surface of rubella virus E1 protein³³,³⁴ (Fig. 3e). Moreover, as previously observed in the case of CrHAP2²⁶, the loops of each trimer interact with those of another trimer within the Fsx1 crystal lattice.

In summary, despite significant differences in the fold of domain II, the unprecedented presence of a domain IV and extreme electrostatic properties, the overall structural similarity between Fsx1 and viral or eukaryotic fusexins suggests that this prokaryotic molecule also functions to fuse membranes.

**Fsx1 can fuse eukaryotic cells.** To test the fusogenic activities of the candidate archaeal fusexins we studied their fusion activity upon transfection in eukaryotic cells¹⁵,¹⁶,²². Cells with either red or green nuclei are mixed with each other and fusion is measured by the formation of hybrid cells with both red and green nuclei revealing merger of their cytoplasms. For this, we co-cultured two batches of Baby Hamster Kidney (BHK) cells independently transfected with Fsx1 and co-expressing either nuclear H2B-RFP or H2B-GFP²². We then performed immunofluorescence against a V5 tag fused to the cytoplasmic tail of Fsx1 (Fig. 4a, b, and Supplementary Fig. 7a, b). We observed a five-fold increase in the mixing of the nuclear H2B-GFP and H2B-RFP compared to vector control, showing that Fsx1 is a bona fide fusogen, as efficient as the fusexin AtHAP2 (Fig. 4c). To determine whether Fsx1 expression is required in both fusing cells or, alternatively, it suffices in one of the fusing partners, we mixed BHK-Fsx1 coexpressing cytoplasmic GFP with BHK cells expressing only nuclear RFP. We found increased multinucleation of GFP⁺ cells (revealing cell–cell fusion) but very low mixing with RFP⁺ cells not expressing Fsx1. In contrast, the vesicular stomatitis virus G-glycoprotein (VSVG) fusogen-induced efficient unilateral fusion¹⁵ (Fig. 4d–f). While VSVG requires acidic pH for maximum fusion activity Fsx1-mediated multinucleation was not stimulated by low pH (Supplementary Fig. 1e). Thus, Fsx1 acts in a bilateral way, similarly to EFF-1 and AFF-1 fusexins¹⁴,¹⁵,³⁵. We then performed live imaging using spinning disk confocal microscopy and observed bilateral cell-cell fusion of BHK-Fsx1 cells (Fig. 4g, h).
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

**Fig. 2** Domain architecture of Fsx1 and topological comparison with HAP2. **a** Schematic diagram of the domains of Fsx1. SP, signal peptide; TM, predicted transmembrane helices. Note that Fsx1 is predicted to contain three C-terminal TMs, with a cytoplasmic loop between TM1 and TM2 that lacks Cys residues; on the other hand, HAP2 homologs are characterized by having a single TM, followed by a cytoplasmic tail that often contains Cys implicated in fusion⁶⁴. **b** Crystal structure of the ectodomain of Fsx1 and predicted topology of the full-length protein relative to the plasma membrane. Domains I-IV are shown in red, yellow, blue, and green, respectively; disulfide bonds are indicated and colored magenta. **c** Topology diagrams of the ectodomains of Fsx1, C. reinhardtii HAP2 (PDB 6E18 [https://www.rcsb.org/structure/6E18])⁶⁵, A. thaliana HAP2 (PDB 5OW3 [https://www.rcsb.org/structure/5OW3])²⁵, and T. cruzi HAP2 (PDB 5OW4 [https://www.rcsb.org/structure/5OW4])²⁵. Domains and disulfide bonds are colored as in panel **b**. Note how, although domain II of Fsx1 has the same topology as the corresponding domain of HAP2, it contains only one of its conserved disulfide bonds (C₁₁₂₅–C₂₁₅₅, corresponding to C147–C210/disulfide bond 3 of CrHAP2²³). **d** Side-by-side comparison of Fsx1E and known class II fusogens. Elements are colored as in panel **b**; the stem region and the linker between domains I and III are shown in pink and cyan, respectively.
<figure><img src="image_9.png" /><figcaption>**Fig. 3** Distinct structural features of Fsx1. **a** Fsx1E surface colored by electrostatic potential from red (−10 kT/e) to blue (+10 kT/e) through white (0 kT/e). Orange arrows indicate the FLs. **b** Location of the Ca²⁺ and Na⁺ ions (depicted as yellow and purple spheres, respectively; see also Supplementary Figs. 4 and 5) stabilizing the Fsx1E trimer. The molecular surface of a protomer is shaded green and the outline of the other two subunits is colored cyan and wheat. **c** Details of interactions at the level of the six-helical bundle made by domain II of the Fsx1 subunits (right subpanels) and the domain IV/domain II interface (bottom left subpanel). Selected side chains are colored by the type of inter-chain contacts in which they are involved (gray: hydrophobic interaction; blue bell: hydrogen bonding; yellow: Ca²⁺ coordination; fuchsa: Na⁺ coordination), with dashed lines indicating hydrogen bonds. Note that the helical bundle of Fsx1 is not a leucine-zipper coiled-coil structure, such as those found in class I/III viral fusion proteins or in the SNARE four-helix bundles, and see also Supplementary Fig. 5b. **d** Superposition of Fsx1 domain IV (green) and Der p 23 (PDB 4ZCE [https://www.rcsb.org/structure/4ZCE])³¹, blue) (Dali Z-score 3.6, RMSD 2.2 Å). **e** Comparison of the Fsx1 region that includes the FL and the corresponding parts of CrHAP2 and Rubella virus E1 protein (PDB 4B3V [https://www.rcsb.org/structure/4B3V])³³. Residues coordinating the Ca²⁺ ion that stabilizes the Fsx1 FL are underlined, and compared to the E1 protein Ca²⁺-binding region in the boxed panels on the far right. **f** The Fsx1 FL adopts a highly ordered conformation stabilized by a Ca²⁺ ion. Presence and identity of the latter, indicated by a yellow arrow, are supported by two other maps shown in addition to the 2mFo-DFc map (blue mesh, contoured at 1.0 σ): a difference map calculated upon omitting all metal ions from the model (thick green(+)/red(−) mesh, 6.0 σ) and a phased anomalous difference map calculated from a 2.9 Å-resolution dataset collected at 7.1 keV (thick magenta mesh, 3.2 σ).</figcaption></figure>

**Structure-function analysis of Fsx1.** To compare archaeal Fsx1 activity with fusexins from eukaryotes and viruses, we introduced mutations into three different structural domains of Fsx1 and tested surface expression and fusogenic activity in mammalian cells.

First, to test whether the putative FL of Fsx1 (143-SVTSPV-148) is involved in fusion, we replaced it with a linker of 4G between Y142A and Y149A (Figs. 3e, 5a, and Supplementary Fig. 5d; ΔFL → AG₄A). This replacement does not affect surface expression yet reduces cell-cell fusion to levels similar to those of the negative control (Fig. 5b-f).

Second, we asked whether domain IV, which is only present in archaeal fusexins, has a function in the fusion process. For this, we replaced the entire domain with the stem region of EFF-1 (Figs. 3d and 5a; ΔDIV → EFF-1 stem). While this mutant Fsx1 reaches the cell surface, suggesting that it folds normally, it

shows a significantly reduced activity compared to wild-type Fsx1 (Fig. 5b-f).

Third, to test whether the three TMs of Fsx1 are required for fusion, we replaced them with the TM and cytoplasmic domains of EFF-1 (Fig. 5a; ΔTMs → EFF-1 TM) or a glycosylphosphatidylinositol (GPI) anchor signal (Fig. 5a; ΔTMs → GPI). We found that both Fsx1 mutants remained active (Fig. 5b), indicating that the Fsx1 TMs are not essential for fusion. Finally, we also replaced the TM and cytoplasmic domains of AtHAP2 with a signal for GPI and found that the protein also maintained its fusogenic activity (Figs. 4c, and 5a, b). Thus, contrary to some viral fusogens in which the GPI-anchored glycoproteins fail to drive complete fusion³⁶⁻³⁸, lipid-anchored Fsx1 or eukaryotic HAP2 promote syncytia formation when expressed on the surface of BHK cells.
<figure><img src="image_10.png" /><figcaption>**Fig. 4 Fsx1 mediates bilateral cell-cell fusion.** **a–c** Cell-cell fusion was measured by content-mixing, indicated by the appearance of multinucleated cells containing green nuclei (H2B-GFP) and magenta nuclei (H2B-RFP). Immunofluorescence against the V5 tag was also performed (gray). **a** Images of mixed cells. DAPI, blue. Scale bars, 20 µm. See also Supplementary Fig. 7a. **b** Scheme of experimental design. **c** Quantification of content-mixing. The mixing indices presented as means ± SEM of four independent experiments. Comparisons by one-way ANOVA followed by Bonferroni's test. ns = non-significant, \*\*\*p < 0.001, \*\*\*\*p < 0.0001. Source data are provided as a Source Data file. **d–f** Unilateral fusion was evaluated by mixing control cells expressing nuclear H2B-RFP (magenta) with cells expressing GFP with a nuclear export signal (GFPnes, green cytoplasm) only or together with Fsx1, EFF-1 or VSV G. **d** Images of cells transfected with empty::GFPnes vector, Fsx1::GFPnes, EFF-1::GFPnes or VSV G::GFPnes. Fsx1 and EFF-1 show multinucleated GFPnes positive cells (arrowheads). VSV G multinucleated cells are found with GFP only (arrowheads) or with both markers (arrows). Scale bars, 20 µm. See also Supplementary Fig. 7e. **e** Scheme of experimental design. **f** Quantification of content-mixing experiments in which only the GFP population of cells express Fsx1, EFF-1, VSV G, or none of them (vector). Bar chart showing means ± SEM of three independent experiments. Comparisons by one-way ANOVA followed by Dunett's test against the vector. ns = non-significant, \*\*p < 0.01, \*\*\*p < 0.001, \*\*\*\*p < 0.0001. Source data are provided as a Source Data file. **g** Spinning disk microscopy time-lapse images indicating the merging of two cells expressing myr-mCherry and Fsx1. Time in hours:minutes. The red channel (mCherry, white) and the DIC are shown. Refer to Supplementary Movie 1. **h** For the last point a Z-projection showing the myr-mCherry fluorescence (white) and the nuclei Hoechst (blue; Supplementary Movie 2). Scale bar, 20 µm.</figcaption></figure>

**Fsx1s are ancient fusogens associated with integrated mobile elements.** The *fsx1* genes here identified are present in a wide physicochemical landscape (Fig. 6a). We observed that the branching pattern of Fsx1 sequences from complete genomes is incompatible with their species tree (Fig. 6a, b). This and the sparse pattern of Fsx1 presence in Archaea led us to perform genomic comparisons of related species with and without the *fsx1* gene. These revealed >50 kbp DNA insertions in the genomes of species with *fsx1* genes (Fig. 7a). To investigate them, we performed k-mer spectrum analysis on *fsx1*-containing Pure Culture

Genomes (PCGs) and found divergent regions containing the *fsx1* ORF (Fig. 7b). Gene content analyses of *fsx1*-containing regions show that they share a portion of their genes (Supplementary Fig. 9) and display conserved synteny (Fig. 8), suggesting common ancestry. These regions are enriched in ORFs homologous to proteins involved in DNA mobilization and integration (Fig. 8 and Supplementary Table 3). Thus, our results indicate that *fsx1* genes are contained in integrated mobile elements (IMEs) that can be mobilized by a conjugative-like, cell fusion-dependent mechanism.
<figure><img src="image_11.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_21.png" /></figure>

To describe Fsx1's tempo and mode of evolution we first compared archaeal and sexual fusexins, which share enough sequence conservation to apply standard phylogenetic analyses, not possible for somatic and viral fusexins, as these methods are not able to cope with such amount of divergence (Fig. 1a

and Supplementary Fig. 1a). We built maximum likelihood (ML) phylogenies for a set of Fsx1 sequences derived from isolated species and metagenomes, and a subset of HAP2s which capture the full phylogenetic diversity present in eukaryotic lineages (Supplementary Fig. 10a). A major finding
**Fig. 5** *Structure-function analysis of Fsx1.* **a** Schematic diagram of wild-type Fsx1, four mutants and AtHAP2ΔTM → GPI. SP signal peptide, FL fusion loop. For colors and abbreviations see legend of Fig. 2. **b** Quantification of content-mixing (cell-cell fusion) in populations of cells expressing vectors (n = 7), Fsx1 (wt) (n = 7), its mutants (ΔFL → AG₄A (n = 6), ΔDIV → EFF-1 stem (n = 4), ΔTM5 → EFF-1 TM (n = 4), ΔTM5 → GPI (n = 3), or AtHAP2ΔTM → GPI (n = 3). Bar chart showing means ± SEM. Comparisons by one-way ANOVA followed by Bonferroni's test against the vector (black) and against Fsx1 (red). ns = non-significant, \*\*\*p < 0.001, \*\*\*\*p < 0.0001. Source data are provided as a Source Data file. **c** Representative merged images from the experiments in (b): magenta (RFP); green (GFP) and blue (DAPI). Fused cells with RFP and GFP (arrows). Scale bars, 20 µm. See also Supplementary Fig. 7f. **d** Immunoblot of EFF-1-V5, control (untransfected cells) and Fsx1-V5 expressing cells. "Surface" indicates surface biotinylation followed by affinity purification using neutravidin agarose beads; "Total" indicates the expression in whole cell extracts. Actin is used as a loading control. The amount of initial cells for Fsx1 is 4 times higher than EFF-1. n = 3. **e** Surface biotinylation as explained in panel d for cells expressing Fsx1-V5 (WT), ΔFL → AG₄A-V5, ΔDIV → EFF-1 stem-V5 or ΔTM5s → EFF-1 TM-V5. n = 3. **f** Immunofluorescence images on non-permeabilized cells expressing Fsx1-FLAG (WT), AFF-1-FLAG (negative control, cytotail), Fsx1-ΔFL → AG₄A-FLAG, Fsx1-ΔDIV → EFF-1 stem-FLAG, AFF-1-FLAG (permeabilized), Fsx1-ΔTM5s → EFF-1 TM-FLAG, Fsx1-ΔTM5s → GPI and AtHAP2-ΔTM → GPI. The FLAG tag was inserted before the first TM or the GPI signal of each construct except for C. elegans AFF-1 in which the FLAG is at C-terminal after the cytoplasmic tail. Transfected BHK cells were incubated with anti-FLAG antibody on ice before fixation. Non-permeabilized staining of FLAG antibody showed the surface expression of Fsx1 and the mutants. C. elegans AFF-1 tagged with a cytoplasmic FLAG is a negative control for non-permeabilized staining. Permeabilized staining of CeAFF-1-FLAG shows the localization on plasma membrane and internal compartments (see also Supplementary Fig. 7g). Scale bars, 10 µm.*

comes from these phylogenies: eukaryal and archaeal fusexins cluster into strongly supported clades suggesting they diverged before LECA.

To place *fsx1* in the broader fusexin superfamily context, we performed structural phylogenetic analysis comparing crystal structures from viral, somatic, and eukaryotic gamete fusogens (Supplementary Fig. 10b). This structure-based tree supports a viral origin of somatic fusexins (EFF-1)¹⁶ and is also compatible with archaeal fusexins appearing before the radiation of eukaryotes.

Discussion

All *fsx1* genes found in cultured and isolated genomes are restricted to the Haloarchaea clade. Although 83% of *fsx1* genes were found in saline environments, they are not restricted to a particular niche, neither geographically nor environmentally and, by extension, potentially not restricted to halophilic archaea (Supplementary Data 1 and Fig. 6a). Only 16 out of 88 *fsx1* metagenomic genes have taxonomy assigned. The 16 MAG-containing fusexins are distributed in other archaeal clades (including Asgard) and also in bacteria but we take those taxonomic assignments with caution. *fsx1*-containing MAGs are highly fragmented and were assembled using methods that include sequence composition (k-mer) criteria, allocating scaffolds with similar k-mer spectra into MAGs. As *fsx1* genes from PCGs are located within distinctive k-mer regions and the metagenomic contigs containing the *fsx1* gene are compositionally homogeneous, and plausibly are also IMEs, it is quite possible that MAG-fusexins are misassigned. The association of haloarchaeal *fsx1s* with IMEs (Fig. 8), genomic comparisons of close species (Fig. 7a) and the incompatibility between their sequence phylogeny and cognate species tree (Fig. 6), indicate lateral mobility within the Haloarchaea class. This evidence suggests that Fsx1s mediate a cell fusion-dependent genetic exchange in archaea. This hypothesis is consistent with the genetic structure and lifestyle of halophilic archaea which are notorious for being polyploid³⁹ and undergoing massive genetic exchanges that overcome species and genera barriers⁴⁰,⁴¹. Moreover, compelling evidence of archaeal cell fusion comes from studies showing bilateral DNA exchange that correlates with cytoplasmic bridges made up of fused lipid bilayers connecting haloarchaeal cells⁴²⁻⁴⁴. Thus, it is plausible that Haloarchaea evolved HGT mechanisms based on conjugative-like DNA mobilization and cell-cell fusion⁴⁵.

The “virus or the egg” dilemma²⁸ posits that fusexins may have been either viral innovations (class II fusogens), later acquired by eukaryotes, or vice versa. Archaeal fusexins expand this dilemma:

gamete fusogens may have prokaryotic origins. Both structure- and sequence-based trees (Supplementary Fig. 10) do not solve but provide insights to address this conundrum, in which we distinguish three main hypotheses indicating alternative origins: Virus, Eukarya, and Archaea (Fig. 9). For all three scenarios we assume that sexual fusexins (HAP2) were present in the LECA¹,²⁷.

A Virus-first scenario has circumstantial evidence favoring it. Exaptation of viral genes is documented for all three domains of life⁴⁶. To be consistent with the observed basal divergence between archaeal and eukaryal (HAP2) fusexins (Supplementary Fig. 10a) the Virus-first hypothesis must include two HGT events before the eukaryal radiation (Fig. 9a). Thus, an archaeon could have exapted a fusexin from an enveloped archaeal virus and then transferred it to an early eukaryote. The reciprocal is also formally possible: an eukaryal viral fusexin was exapted by a pre-LECA eukaryotic cell and then transferred to an archaeon, before the eukaryal radiation. Alternatively, the ancestral fusexin-containing virus was able to infect both archaea and eukaryal cells (Fig. 9a). These putative events are at odds with the distribution of fusexins in extant viruses. All currently known viral fusexins belong to RNA viruses that are confined to a few multicellular hosts: vertebrates, arthropods, and flowering plants⁴⁷,⁴⁸. This distribution favors a scenario where viral fusexins, like many other eukaryotic viral proteins, have eukaryotic cellular origins⁴⁹. Unlike currently known fusexin-containing viruses, all isolated archaeal viruses to date have DNA genomes. The recently elucidated structure of VP5, a haloarchaeal virus envelope protein that mediates cell invasion shows a fold that differs from all previously known viral fusion proteins, including fusexins⁵⁰. Although there is no evidence for the presence of fusexins in archaeal viruses, upcoming metagenomic and structural analyses may provide support to the Virus-first hypothesis.

The widespread presence of sexual fusexins in Eukarya indicates evolutionary success, in line with the Eukarya-first hypothesis (Fig. 9b). However, introduction of an eukaryal fusexin into the Archaeal domain is less supported by currently available evidence as interdomain gene transfers from eukaryotes to archaea are hardly documented and thought to be scarce⁴⁵,⁵¹.

The presence of *fsx1* genes in Haloarchaea IMEs is consistent with gene transfer in the opposite direction. Eukaryogenesis, and by extension the emergence of sex, is marked by massive horizontal gene transfer events to the archaeal ancestor of eukaryotes⁵², disregarding if it belonged to the Asgard superphylum or to a sister group of Archaea⁵³. During the First
a

<figure><img src="image_22.png" /></figure>

b

<figure><img src="image_23.png" /><figcaption>**Fig. 6 Environmental distribution of Fsx1s.** a Archaeal fusexins unrooted phylogeny, environmental details, and trimeric models⁶⁶ based on the Fsx1 X-ray structure (PDB 7P4L [https://www.rcsb.org/structure/7P4L], enlarged). Tree tip IDs coming from metagenomic data coded by sampling site, sample id and sequence (see Supplementary Data 1). Fsx1s from cultivated genomes are shown as Natriabales and Haloferacales orders in blue and red, respectively. Surfaces colored and oriented as in Fig. 3a. Found in diverse environments and often differing in global surface properties, modeled trimers share hydrophobic tips. b Panel with phylogenomic tree for *fsx1*-containing cultured archaeal species, showing incongruences between Fsx1s in Haloferacales.</figcaption></figure>

Eukaryotic Common Ancestor (FECA) to LECA transition, in addition to the alphaproteobacterial endosymbiont-related inherited genes, the pre-LECA received hundreds of archaeal genes from other lineages, including Euryarchaeota⁵⁴. Thus, an Archaea-first hypothesis (Fig. 9c) considers horizontal transfer of archaeal fusexins into the ancestor of eukaryotes. A weakness of this hypothesis is the sparse distribution of fusexins in archaeal genomes. This distribution is not a product of a recent

HGT from Eukarya, as phylogenetic analysis indicates basal divergence between eukaryotic and archaeal fusexins (Supplementary Fig. 10a). Lateral mobility of extant *fsx1* genes within Haloarchaea, their relative confinement to few archaeal lineages and their basal divergence from sexual fusexins suggests they are molecular relics, and that cell fusion-based HGT might have declined during archaeal evolution in favor of conjugation, transduction and transformation.
a

<figure><img src="image_24.png" /></figure>

b

<figure><img src="image_25.png" /><figcaption>**Fig. 7** Genomic features of Fsx1s. **a** Whole genome comparison of species with and without *fsx1*. Each blue dot represents a segment of 500 bp with more than 80% identity between the species harboring *fsx1* (e.g., *Haloplanus natans* DSM 17083) and the species with no *fsx1* (e.g., *Haloplanus* sp. CBA1112, see Supplementary Table 4). Species with *fsx1* are in the x-axis, the base of the green rectangles represent the detected IME carrying the *fsx1* gene, locus of *fsx1* is in red vertical line and pointed with a red arrowhead. **b** K-mer spectra deviation of *fsx1*-containing IMEs. K-mer spectrum deviation from centroid is shown for each of the Pure Culture Genomes (PCGs) where *fsx1* was detected. Blue region shows the standard deviation. Locus of *fsx1* is in red vertical line and pointed with a red arrowhead, the mobile element containing *fsx1* is in green. Dashed vertical white lines indicate the end of a contig in the genome assembly. *fsx1* is consistently found within regions that deviate from the core genome's spectrum, indicating they belong to a mobile element. Kbp, kilo base pairs.</figcaption></figure>

This third scenario, an Archaeal origin of fusexins, poses new challenges to both sex origins and eukaryogenesis models. Discovery of the Asgard superphylum⁵⁵ and the recent cultivation of one of its members⁵⁶ support eukaryogenesis scenarios where populations of bacteria and archaea lived in syntrophy, transferring metabolites and genes⁵⁷. Acquisition of a *fsx1* gene during the FECA to LECA transition could have enabled pre-LECA cells to undergo genome expansion, explore syncytial forms⁵⁸ and evolve into mononucleated cells fully equipped for meiosis and gamete fusion⁵⁹. Our findings raise the possibility that gamete fusion is the product of over two billion years of evolution of this ancient archaeal cell fusion machine.

The archaeal proteins herein identified place fusexins in yet another domain of life, with different membrane chemistries and along a broad niche landscape, from cold hypersaline lakes to hot springs and hydrothermal vents (Fig. 6a). Our structural and functional analyses show that Fsx1 has both conserved and divergent properties when compared to eukaryotic and viral fusexins (Figs. 3 and 4). Like its viral counterparts, Fsx1 has an uncharged loop required for fusion. However, unlike previously known fusexins, Fsx1 harbors an additional domain (IV) involved in

fusogenic activity that may bind sugars (Figs. 2d and 3d). Considering that cell surface glycosylation was found to be important for fusion-based mating of halophilic archaea⁶⁰, this domain may actively promote fusion by interacting with carbohydrates attached to lipids or proteins such as S-layer glycoproteins⁴². Unlike HAP2s, Fsx1 homologs have 1–4 TMs and a variable Cys number (5–30, see Supplementary Data 1). Like eukaryotic fusexins, Fsx1 mediates BHK cell fusion in a bilateral fashion (Fig. 4f). However, in contrast to viral fusogens³⁶⁻³⁸,⁶¹,⁶², the fusion activity of Fsx1 is maintained following substitution of its three TMs with a single TM or a GPI anchor. The retention of fusion activity when the transmembrane domains are replaced by a GPI anchor supports the model in which Fsx1 mediates homotypic fusion with fusogens required on both target membranes. These findings also suggest that interaction between the TMs during fusion is not essential for trimerization and expansion of the fusion pores. Since GPI-anchored AtHAP2 is also fusogenic, other fusexins may also drive complete cell fusion without a specific involvement of TMs. Future studies will address the function of the six-helix bundle formed by Fsx1 domain II, which is unprecedented among fusexins and raises an unexpected structural connection with class I viral fusogens⁵,⁶.
<figure><img src="image_26.png" /><figcaption>**Fig. 8 Fsx1s are embedded in integrated mobile elements (IMEs).** Synteny plots for IMEs from PCGs and metagenomic data. Annotated regions plus inferred ORFs belonging to homologous clusters identified by the workflow depicted in Supplementary Fig. 8. Homology relationships are represented by gray links. *fsx1* genes are marked in red and selected ORFs homologous to IME signature genes are labeled and color-coded. XerC/XerD recombinases (green); HerA helicase (dark blue); VirB4, Type IV secretion system (T4SS) pathway (cyan); TraG/TraD/VirD4 family enzyme, ATPase, T4SS (see Supplementary Table 3 and Supplementary Data 2, 3 for details). The 11 segments analyzed correspond to the cluster marked in Supplementary Fig. 9.</figcaption></figure>

What are the limitations of this study? First, the description of the cellular and molecular functions of the *fsx1* family presented in this work is reliant upon expression within a heterologous system designed to probe fusion activity. Despite membrane chemistry differing from the type that would be associated with haloarchaeal hosts of the IME, Fsx1 was able to promote membrane fusion reliably. Second, evolutionary analysis and comparative genomics confirm *fsx1*'s link to IMEs. These sequences appear to be spread across a wide variety of niches all across the globe despite a relatively sparse distribution in sequenced archaeal genomes. Third, in this study we focused on describing the Fsx1 family evolutionarily, structurally and functionally to place it within its context in the Fusexin superfamily, but future studies will be needed to elucidate which biological processes *fsx1* is involved in as well as its relationship to the rest of the archaeal mobilome and virome. Additionally, structural features, such as the presence of a fourth domain, were modified in order to observe their effect on fusion activity within our experimental system, but it remains to be seen what importance they have in their native context. Future work will focus on studying *fsx1* in archaeal experimental systems as well as leverage metagenomic sampling and assembly techniques to exhaustively detect possible *fsx1* homologs in environments where it may be facilitating horizontal gene transfer.

## Methods

**Initial fusexin search using structurally guided MSAs.** HMMs were prepared using structurally guided multiple sequence alignments (MSAs) of known eukaryotic HAP2 sequences (ectodomains only). Structural MSAs were derived using I-TASSER⁶⁸-generated models of HAP2 homologs for *Erythranthe guttata* (A0A022QRC8), *Phytomonas* sp. isolate EM1 (W6KUI1), *Plasmodium falciparum* (A0A1C3KGX6), *Chlorella variabilis* (E1Z455) and the HAP2 crystal structures for

*Chlamydomonas reinhardtii* (PDB 6E18 [https://www.rcsb.org/structure/6E18]⁶⁵, 6DBS [https://www.rcsb.org/structure/6DBS]²⁶) and *Arabidopsis thaliana* (PDB 5OW3 [https://www.rcsb.org/structure/5OW3]²⁵).

Searches for fusexin homologs using structurally guided MSAs were performed for 3 iterations on the Uniclust database⁶⁹ using default HHblits⁷⁰ parameters.

**HMM-based distance matrices.** A taxonomically representative list of known viral and eukaryotic fusexin homologs, covering major lineages, was manually curated. A MSA was built for each homolog by using the sequence as a query on the Uniclust database with HHblits for three iterations. This set of MSAs was compiled into an HH-suite database and each MSA was used as a query against this database to establish a profile-based distance matrix using the probability of homology (Fig. 1a and Supplementary Fig. 1).

**Metaclust database search pipeline.** We searched the Metaclust⁷¹ dataset (nr50) using an HMM made of Fsx1 sequences found in PCGs and MAGs (Supplementary Data 1; see also codes, notebooks and datasets available at Zenodo⁶³). Fsx1 sequences were aligned using ClustalO⁷² with default settings for 3 iterations and the resulting MSA was used as a query with HMMER hmmsearch⁷³ against the Metaclust50 dataset⁷¹. All returned sequences with an E-value < 0.0001 with a match length greater than 100 residues were selected for further analysis. PSI-BLAST⁷⁴ was also used on the Metaclust (nr90) with Fsx1 sequences found in PCGs and MAGs with default parameters for 3 iterations. All returned sequences with an E-value < 0.0001 and an alignment length greater than 100 were added to the pool of candidates. Manual curation was performed using membrane protein topology predictor TOPCONS⁷⁵ and distant homology searches using HHblits⁷⁶ against PDB70.

**DNA constructs.** Ten archaeal genes were synthesized (GenScript) and cloned into pGene/V5-His vectors (Supplementary Table 5). Details of nucleotides used for synthesis and protein sequences are described in Supplementary Data 4.

For structural studies, a synthetic gene fragment encoding the extracellular region of a metagenomic Fsx1 ORF (IMG genome 3300000868, scaffold JGI12330J12834_1000008, ORF 8; Supplementary Data 1) (GenScript) was subcloned by PCR in frame with the 5' chicken Crypa signal peptide- and 3' 8xHis-tag-encoding sequences of pLJ6, a mammalian expression vector derived from
<figure><img src="image_27.png" /><figcaption>**Fig. 9** Possible evolutionary trajectories. **a** Virus-first. Originating in the virosphere (1), ancestral fusexins were transfered to an ancient eukariotic organism (2) before the eukaryal radiation. Also before the eukaryal radiation, either viral (3) or early eukaryal (not numbered) fusexins were horizontally transferred to Archaea where they became integrated mobile elements. **b** Eukarya-first. Originating before the eukaryal radiation (1), fusexins were horizontally transferred to Archaea (2) where they became integrated in mobile elements. Eukaryal fusexins were captured by eukaryal viruses (3) leading to extant viral (class II) fusexins. **c** Archaea-first. Originating in Euryarchaeota (1), fusexins were horizontally transfered (2) and became fixed during eukaryogenesis and the emergence of meiotic sex. Viral capture of fusexin genes from early eukaryotic cells and further evolution within the virosphere led to extant viral (class II) fusexins (3). Common to all models, viral fusexin genes were captured by different eukaryotic lineages (4), leading to phlebovirus-like integrated fusogens and EFF-1 somatic fusogens⁴⁸,⁶⁷. Solid lines represent evolutionary trajectories of Archaea (red), viruses (blue), Eukarya (black), and eukaryogenesis (red to black gradient). Dashed arrows represent HGT events. Question marks denote uncertainty regarding the presence of *fsx1*-related genes in the respective lineages. LACA Last Archaeal Common Ancestor.</figcaption></figure>

pHLsec3⁷⁷. The protein construct that yielded the final high-resolution dataset included residues D25-S535 and contained a T369C substitution, introduced by PCR mutagenesis with the aim of facilitating heavy atom derivatization for experimental phasing. Oligonucleotides were from Sigma-Aldrich or IDT and all constructs were verified by DNA sequencing (Eurofins Genomics or Macrogen).

To generate pCI::GFPnes plasmid (see list of plasmids in Supplementary Table 6), an oligo DNA encoding for the nuclear export signal (LQKKLEELELD) was cloned downstream the region encoding EGFP of the pCAGIG plasmid using the enzyme BsrGI. Then, the GFPnes coding sequence was amplified, cut with BmgBI and BglII and used to replace the H2B-GFP coding sequence of the pCI::H2B-GFP plasmid (see list of primers in Supplementary Table 7). Fsx1-V5, AtHAP2-V5²², EFF-1-V5, V5V-G¹⁵ and other archaeal fusexins (NaFsx1, HQ22Fsx1, HnFsx1) were subcloned into corresponding pCI::H2B-RFP/H2B-GFP/GFPnes vectors separately. For mutagenesis of Fsx1, (i) Fsx1-ΔFL-AG₄A: The mutation of Y142A, Y149A and four glycines inserted between them were achieved using PCR with overlapping primers. (ii) Fsx1-ΔDIV-EFF-1 stem: The stem region of EFF-1 (E510-D561) was amplified from pGene::EFF-1-V5 and fused to the upstream and downstream regions of Fsx1-DIV with overlapping primers. (iii) Fsx1ΔTM → EFF-1 TM: The TM and cytoplasmic tail of EFF-1 (I562-I658) were amplified from pGene::EFF-1-V5 and fused to the ectodomain of Fsx1 to replace its TMs. (iv) Fsx1ΔTM → GPI: The Fsx1 TMs were replaced with the carboxy-terminal 37 amino acids of decay accelerating factor (DAF) which confer the signal for GPI anchor⁷⁸. Similarly, the TM and cytoplasmic tail of AtHAP2 were replaced with the GPI signal from DAF to get AtHAP2ΔTM → GPI. All mutants were ligated into pCI::H2B-RFP and pCI::GFPnes vectors for mixing assay. Additional details are found in Supplementary Tables 6 and 7.

**Protein expression and purification.** HEK293T cells (ATCC CRL-3216)⁷⁹ were transiently transfected using 25 kDa branched polyethyleneimine and cultured in DMEM media (Invitrogen) supplemented with 2% (v/v) fetal bovine serum (Biological Industries). 90–96 h after transfection, the conditioned media from HEK293T cells was harvested, 0.2 µm-filtered (Pall) and adjusted to 20 mM Na-HEPES pH 7.8, 2.5 M NaCl, 5 mM imidazole. 10 ml Ni Sepharose excel beads (GE Healthcare) pre-equilibrated with immobilized metal affinity chromatography (IMAC) buffer (20 mM Na-HEPES pH 7.8, 2.5 M NaCl, 10 mM imidazole) were added to 1 l adjusted conditioned media and incubated overnight at 4 °C. After washing the beads with 100 column volumes IMAC buffer, captured Fsx1E was

batch-eluted with 30 mL 20 mM Na-HEPES pH 7.8, 2.5 M NaCl, 500 mM imidazole, and concentrated with 30 kDa-cutoff centrifugal filtration devices (Amicon). The material was then further purified by SEC at 4 °C, using an ÄKTAflpc chromatography system (GE Healthcare) equipped with a Superdex 200 Increase 10/300 GL column (GE Healthcare) pre-equilibrated with 20 mM Na-HEPES pH 7.8, 2.5 M NaCl. Peak fractions were pooled and concentrated to 5 mg mL⁻¹ (Supplementary Fig. 2a, b).

**Size-exclusion chromatography-multiangle light scattering (SEC-MALS).**

Purified Fsx1E samples (120–150 µg) were measured using an Etvan LC high-performance liquid chromatography system with a UV-900 detector (Amersham Pharmacia Biotech; λ = 280 nm), coupled with a miniDawn Treos MALS detector (Wyatt Technology; λ = 658 nm) and an Optilab T-rEX dRI detector (Wyatt Technology; λ = 660 nm). Separation was performed at 20 °C using a Superdex 200 Increase 10/300 GL column (GE Healthcare) with a flow rate of 0.5 mL min⁻¹ and mobile phases consisting of 20 mM Tris-HCl pH 8.5, 150 mM NaCl (normal salt condition) or 20 mM Tris-HCl pH 8.5, 2.0 M NaCl and 0.2 M CaCl₂ (high salt condition) (Supplementary Fig. 2c). Data processing and weight-averaged molecular mass calculations were performed using ASTRA (Wyatt Technology). BSA (150 µg) was used as a control.

**Small-angle X-ray scattering (SAXS).** SAXS experiments were performed at beamline BM29 of the European Synchrotron Radiation Facility (ESRF)⁸⁰, using Fsx1E (4.5 mg mL⁻¹) in 20 mM Na-HEPES pH 7.8, 150 mM NaCl. Sample delivery and measurements were performed using a 1 mm-thick quartz capillary, which is part of the BM29 BioSAXS automated sample changer unit⁸¹. Data were collected at 1 Å wavelength in 10 frames of 1 s at 20 °C, using an estimated beam size of 1 mm × 100 µm; buffer blank measurements were carried out under the same conditions, both before and after sample measurement. Data were averaged and subtracted using PRIMUS⁸² from the ATSAS package⁸³, which was also used to calculate the pair-distance distribution function, as well as the radius of gyration and the Porod volume. Theoretical scattering curves for monomeric and trimeric Fsx1E were calculated and compared with the experimental data using CRY SOL⁸⁴. Ab initio envelope reconstruction was performed with DAMMIF⁸⁵, resulting in 20 models that were superimposed and averaged with DAMAVER⁸⁶. Chain A of the refined Fsx1E model was either rigidly fitted with UCSF ChimeraX⁸⁷ into the
envelope generated by averaging all 20 independent ab initio SAXS models (Supplementary Fig. 2d, top envelope), or flexibly fit with Namdatora⁸⁸ into the average envelopes generated from the two most abundant clusters of SAXS models (accounting for 4 and 5 of the 20 SAXS models, respectively; Supplementary Fig. 2d, middle and bottom envelopes).

**Crystallization and X-ray diffraction data collection.** Two similar initial hits obtained from extensive screening using a mosquito crystallization robot (TTP Labtech) were manually optimized by setting up vapor diffusion experiments at 20 °C in 24-well plates. To grow diffraction-quality crystals, 1 µl purified Fsx1E was mixed with 1 µl 23% (w/v) PEG 4000, 0.1 M Tris-HCl pH 8.5, 0.2 M CaCl₂, and equilibrated against 1 mL of the same solution. Rhomboidal plates of Fsx1E grew in 1–3 months from protein precipitate that appeared after overnight equilibration of the crystallization drops (Supplementary Fig. 2e). For data collection, specimens were freed from the precipitate by micromanipulation with MicroMounts (MiTeGen) and flash frozen in liquid nitrogen. More than a hundred crystals were screened at beamlines ID23-1 of the ESRF⁸⁹ and 104 of Diamond Light Source, yielding datasets of highly variable quality. The final X-ray diffraction dataset at 2.3 Å resolution was collected at ESRF ID23-1.

**Data reduction and non-crystallographic symmetry analysis.** Datasets were processed in space group C2 with XDS⁹⁰ (Supplementary Table 2). By revealing a strong non-origin peak at chi = 120 (Supplementary Fig. 2f), self rotation functions calculated with MOLREP⁹¹ or POLARRFN⁹² clearly indicated the presence of three-fold non-crystallographic symmetry (NCS) within the asymmetric unit of the centered monoclinic crystals. Combined with Matthews coefficient calculations⁹³,⁹⁴, this strongly suggested that Fsx1E crystallized as a homotrimer.

**Structure determination by molecular replacement with AlphaFold2 models.**

Multiple attempts to experimentally determine the structure of Fsx1E using a variety of heavy atoms failed, probably because the high-salt mother liquor composition hindered heavy atom binding. Because molecular replacement (MR) with HAP2-derived homology models also failed, we phased the data by taking advantage of the recent significant advances in protein 3D structure prediction using machine learning⁹⁵ to phase the data by MR⁹⁶ (Supplementary Fig. 3). To do so, we used AlphaFold²⁹⁷ (with default monomer prediction parameters) to generate five independent models of Fsx1 ectodomain residues D25-S535, with per-residue pseudo-B factors corresponding to 100-(per-residue confidence (pLDDT²⁹⁷)). These models had relative root-mean-square deviations (RMSD) of 1.4–3.3, or 0.7–1.9 Å after excluding 26 C-terminal residues predicted with low-confidence. Initial attempts to solve the structure by MR with Phaser⁹⁷, using an ensemble including these models (further truncated to Q453, the predicted C-terminal end of domain III), yielded 4 solutions (with top Log Likelihood Gain (LLG) 188, final translation function Z score (TFZ) 9.6) that were retrospectively correct in terms of domain I/II placement, but completely wrong in the positioning of domain III. Because of the latter, automatic refinement of these solutions did not progress beyond Rfree ~ 0.53. On the other hand, a parallel consecutive search for three copies of a domain I/II ensemble (D25-A335; RMSD 0.3–0.9 Å) followed by three copies of domain III (P350-Q453; RMSD 0.1–0.3 Å), using a model RMSD variance of 1 Å, yielded a clear single solution (LLG 876, TFZ 23.1) that could be automatically refined to initial R 0.45, Rfree 0.46.

Remarkably, although a single copy of domain 3 corresponds to only 7% of the total scattering mass in the asymmetric unit of the Fsx1E crystal, the very high accuracy of its AlphaFold2 model (reflected by a posteriori-calculated global RMSD and Distance Test Total Score (GDT_TS) of 0.7 Å and 97.6, respectively) allowed Phaser to also find a correct MR solution using just this part of the structure. Specifically, a consecutive search for three copies of the domain resulted in a trimeric model with LLG 275 and TFZ 15.1, which could be refined to starting R 0.51, Rfree 0.51.

Also worth mentioning is the observation that the same domain I/II + domain III MR strategy used to phase the 2.3 Å resolution data could also be successfully applied to an initial dataset at much lower resolution (3.5 Å, with outer shell mean I/σI 0.6 and CC₁/₂ 0.31); in this case, the Phaser LLG and TFZ values for the solution were 361 and 13.5, respectively, and initial automatic refinement of the corresponding model yielded R 0.44, Rfree 0.48.

**Model building, refinement, and validation.** The initial model of Fsx1E was first automatically rebuilt using PHENIX AutoBuild⁹⁸ (1083 residues; R 0.34, Rfree 0.38) and then significantly improved with the machine learning-based sequence-docking method of ARP/wARP⁹⁹, as implemented in CCP4⁹² (1390 residues; REFMAC¹⁰⁰ R 0.23). The resulting set of coordinates was subsequently subjected to alternating cycles of manual rebuilding with Coot¹⁰¹/SOLIDE¹⁰² and refinement with phenix refine¹⁰³, using torsion-based NCS restraints and three Translation-Libration-Screw-rotation groups per chain. Putative identities of the metal ions were assigned based on electron density level; difference Fourier maps generated using alternative atom types; correspondence with peaks in phased anomalous difference maps, calculated with PHENIX¹⁰⁴ or ANODE¹⁰⁵ from data collected at low energy; and coordination properties¹⁰⁶. Protein geometry was validated using MolProbity¹⁰⁷ (Supplementary Table 2).

**Sequence-structure analysis.** Sequence alignments were rendered with ESPript¹⁰⁸ and manually annotated. Transmembrane helices were predicted using TMHMM¹⁰⁹. GDT_TS scores were calculated using LGA¹¹⁰ and structural similarities were assessed with Dali¹¹¹ and PDBeFold¹¹². Secondary structure was assigned using DSSP¹¹³. Subunit interfaces were analyzed using PDBsum¹¹⁴, PIC¹¹⁵ and PDBePISA¹¹⁶. Molecular charge was calculated using the YASARA2 force field of YASARA Structure¹¹⁷ and electrostatic surface potential calculations were performed with PDB2PQR¹¹⁸ and APBS¹¹⁹, via the APBS Tools plugin of PyMOL (Schrödinger, LLC). Mapping of amino acid conservation onto the 3D structure of Fsx1E was carried out by analyzing a sequence alignment of archaeal homologs with ConSurf¹²⁰. Structural figures were generated with PyMOL.

**Structural modeling of trimeric Fsx1.** Models of homotrimeric Fsx1 were generated using a local copy of AlphaFold-Multimer¹²¹, installed using the open source code and instructions available at https://github.com/deepmind/alphafold.

**Cells and reagents.** Baby Hamster Kidney (BHK-21) cells (kindly obtained from Judith White, University of Virginia) were maintained in DMEM supplemented with 10% FBS (Biological Industries), 100 U/ml penicillin, 100 µg/ml streptomycin (Biological Industries), 2 mM l-glutamine (Biological Industries), 1 mM sodium pyruvate (Gibco), and 30 mM HEPES buffer, pH 7.3, at 37 °C with 5% CO₂. Transfections were performed using Fugene HD (Promega) or jetPRIME (PolyPlus) according to the manufacturer's instructions.

**Immunofluorescence.** BHK cells were grown on 24-well tissue-culture plates with glass coverslips. Permeabilized cells were fixed with 4% paraformaldehyde (EM grade, Bar Naor, Israel) in PBS, followed by incubation in 40 mM NH₄Cl to block free aldehydes, permeabilized in 0.1% Triton X-100 in PBS and blocked in 1% FBS in PBS. After fixation, the coverslips were incubated 1 h with mouse anti-V5 antibody (Invitrogen, 1:500) and 1 h with the secondary antibody which was donkey anti-mouse coupled to Alexa Fluor 488 (Invitrogen, 1:500). Alternatively, for immunofluorescence without permeabilization, cells were blocked on ice in PBS with 1% FBS for 20 min, and then stained with Monoclonal ANTI-FLAG M2 antibody (Sigma, 1:1000) on ice for 1 h. After anti-FLAG staining, cells were washed and fixed with 4% PFA in PBS. Cells were blocked again and stained with the secondary antibody (donkey anti-mouse coupled to Alexa Fluor 488; Invitrogen) diluted 1:500 in PBS for 1 h. In all cases, nuclei were stained with 1 µg/ml DAPI. Images were captured using a Nikon Eclipse E800 with a 60X/1.40 Plan Apochromat objective and an optical zoom lens (Nikon) using a Hamamatsu ORCA-ER camera controlled by Micro-Manager software¹²² (Fig. 5f).

**Western blots.** 24 h post-transfection, cells were treated with lysis buffer (50 mM Tris-HCl pH 8.0, 100 mM NaCl, 5 mM EDTA, 1% Triton X-100 supplemented with chymostatin, leupeptin, antipain and pepstatin) on ice for 10 min. After centrifugation at 21,000 × g for 10 min at 4 °C, supernatants of lysates were mixed with reducing sample buffer (+DTT) and incubated 5 min at 95 °C. Samples were loaded on a 10% SDS-PAGE gel and transferred to PVDF membrane. After blocking, membranes were incubated with primary antibody anti-V5 mouse monoclonal antibody (1:5000; Invitrogen) or anti-actin (1:2000; MP Biomedicals) at 4 °C overnight and HRP-conjugated goat anti-mouse secondary antibody 1 h at room temperature. Membranes were imaged by the ECL detection system using FUSION-PULSE.6 (VILBER).

**Content mixing assays with immunofluorescence.** BHK-21 cells at 70% confluence were transfected (using JetPrime; Polyplus at a ratio of 1:2 DNA-transfection reagent) with 1 µg pCI::Fsx1-V5::H2B-eGFP, pCI::Fsx1-V5::H2B-RFP, pCI::AtHAP2-V5::H2B-eGFP, pCI::AtHAP2-V5::H2B-RFP, respectively. Control cells were co-transfected with pCI::H2B-eGFP and pRFPnes or pCI::H2B-RFP and pRFPnes. 4 h after transfection, the cells were washed 4 times with DMEM with 10% serum (Invitrogen), 4 times with PBS and detached using Trypsin (Biological Industries). The transfected cells were collected in Eppendorf tubes, resuspended in DMEM with 10% serum, and counted. Equal amounts of H2B-RFP and H2B-eGFP cells were mixed and seeded on glass-bottom plates (12-well black, glass-bottom #1.5H; Cellvis) and incubated at 37 °C and 5% CO₂. 18 h after mixing, 20 µM 5-fluoro-2'-deoxyuridine (FdUrd) was added to the plates to arrest the cell cycle and 24 h later, the cells were fixed with 4% PFA in PBS and processed for immunofluorescence. To assay mixed cells and detect the transfected proteins (Fsx1-V5 or AtHAP2-V5), we stained cells with anti-V5 mAb (Life Science). The secondary antibody was Alexa Fluor 488 goat anti-mouse, with 1 µg/ml DAPI²². Micrographs were obtained using wide-field illumination using an ELYRA system S1 microscope (Plan-Apochromat ×20 NA 0.8; Zeiss) and recorded with a iXon+ EMCCD camera (Andor). The GFP + RFP mixing index was calculated as the number of Red and Green nuclei in mixed cells out of the total number of nuclei of fluorescent cells in contact (Fig. 4b).

**Cell fusion assay by content mixing with nuclear and cytoplasmic markers.** For the unilateral setup, BHK-21 cells were transfected (as explained above) with 1 µg pCI::GFPnes; pCI::Fsx1-V5::GFPnes; 0.25 µg pCI::EFF-1-V5::GFPnes; 1 µg
pCI::VSV-G::GFPnese in respective 35 mm plates. The cells were incubated, washed, and mixed with pCI::H2B-RFP (empty vector) transfected cells (Fig. 4e). For evaluating the mutants, BHK-21 cells were transfected with 1 µg pCI::Fsx1-V5::GFPnese or pCI::Fsx1-V5::H2B-RFP or the plasmids encoding for each mutant: ΔFL → AG4A, ΔDIV → EFF-1 stem, ΔTMs → EFF-1 TM, Fsx1ΔTMs → GPI, or AtHAP2ΔTM → GPI (Fig. 5a). Empty pCI::GFPnese or pCI::H2B-RFP were used as negative controls. 4 h after transfection, the cells were washed, counted, mixed, and incubated as previously described. In all cases, 18 h after mixing, 20 µM FdUrd was added to the plates, and 24 h later, the cells were fixed with 4% paraformaldehyde diluted in PBS. Nuclei were stained with 1 µg/ml DAPI. Images were obtained using wide-field illumination with an ELYRA system S1 microscope as described above.

The GFP + RFP mixing index was calculated as the number of nuclei in mixed cells, green cytoplasm (GFPnese) with red (H2B-RFP) and blue (DAPI) nuclei out of the total number of nuclei in fluorescent cells in contact (Figs. 4f, 5b). The multinucleation indexes were defined as the ratio between the number of nuclei in multinucleated cells ($N_m$) and the total number of nuclei in multinucleated cells and expressing cells that were in contact ($N_c$) but did not fuse, using the following equation: $% multinucleation = N_m/(N_c + N_m) \times 100$. The percentage of multinucleation was calculated for GFPnese cells with RFP and DAPI nuclei. For the unilateral assay, multinucleation was determined as the ratio between the number of nuclei in multinucleated green cells and the total number of nuclei in green multinucleated cells and GFPnese expressing cells that were in contact but did not fuse (Fig. 4f).

**Live imaging of fusing cells.** BHK cells were plated on 15 mm glass bottom plates (Wuxi NEST Biotechnology Co., Ltd.) and transfected with 1 µg pCI::Fsx1-V5::H2B-FGP together with 0.5 µg myristoylated-mCherry (myr-palm-mCherry; kindly provided by Valentin Dunsing and Salvatore Chiantia¹²³), 18 h after transfection, the cells were incubated with 2 µg/ml Hoechst dye for 10 min at 37 °C and washed once with fresh medium. Time-lapse microscopy to identify fusing cells was performed using a spinning disc confocal microscope (CSU-X; Yokogawa Electric Corporation) with an Eclipse Ti and a Plan-Apochromat ×20 (NA, 0.75; Nikon) objective. Images in differential interference contrast and red channels were recorded every 4 min in different positions of the plate using high gain and minimum laser exposure. Time lapse images were captured with an iXon 3 EMCCD camera (Andor Technology). After 5 h, confocal z-series, including detection of the DAPI channel, were obtained to confirm the formation of multinucleated cells. Image analyses were performed in MetaMorph (Molecular Devices) and ImageJ¹²⁴ (National Institutes of Health).

**Surface biotinylation.** Proteins localizing on the surface were detected as previously described²². Briefly, BHK cells were transfected with 1 µg pCAGGS, pCAGGS::EFF-1-V5, pCAGGS::Fsx1-V5, pCAGGS::ΔFL → AG4A-V5, pCAGGS::ΔDIV → EFF-1 stem-V5 or pCAGGS::fsx1ΔTMs → EFF-1 TM-V5. 24 h later, cells were washed twice with ice-cold PBS²⁺ (with Ca²⁺ and Mg²⁺) and incubated with 0.5 mg/ml EZ-Link Sulfo NHS-Biotin (Thermo Fisher Scientific) for 30 min on ice. The cells were washed four times with ice-cold PBS²⁺, once with DMEM with 10% FBS (to quench residual biotin), followed by two more washes with PBS²⁺. To each plate 300 µl of lysis buffer supplemented with 10 mM iodoacetamide were added and the cells detached using a scrapper. The insoluble debris was separated by centrifugation (10 min at 21,000 × g), and the lysate was mixed with NeutrAvidin Agarose Resin (Thermo Fisher Scientific) and 0.3% SDS. After an incubation of 12 h at 4 °C the resin was separated by centrifugation (2 min at 21,000 × g), washed three times with lysis buffer and then mixed with SDS-PAGE loading solution with freshly added 5% β-mercaptoethanol and incubated 5 min at 100 °C. After pelleting by centrifugation, the samples were separated by SDS-PAGE gel and analyzed by Western blotting as described above using anti-V5 mouse monoclonal antibody. Loading was controlled using anti-actin C4 monoclonal (1:2000; MP Biomedicals).

**Integrated mobile element (IME) identification by k-mer spectra analysis and comparative genomics.** Comparison between close species with presence (fsx1+) or absence (fsx1−) of archaeal fuxins to detect insertion sites was done performing sequence similarity searches in complete genomes from the closest relatives available in the PATRIC database¹²⁵ (Fig. 7a and Supplementary Table 4). Coordinates of fsx1-containing IMEs present in PCGs are annotated in Supplementary Table 4.

Among different methodologies that rely on DNA composition to identify horizontally transferred genomic regions¹²⁶, k-mer spectrum analysis is a standard tool for this purpose¹²⁷,¹²⁸. Normalized k-mer spectra for DNA sequences of arbitrary length were generated by counting occurrences of all k-mers and normalizing by the total amount of words counted. k-mer sizes from 3 to 8 bp were tested with no effect on results. A length of 4 bp was selected. To detect possible horizontally transferred regions, an average spectrum for each genome was calculated. A spectrum was calculated for a sliding window of 1 kb using 500 bp steps and subtracted from the genomic average at each window position (Fig. 7b). The absolute value of the difference between the genomic average and window spectra is represented over the entire genome. Gaussian mixture models using two distributions were fitted¹²⁹ to the k-mer content of all windows, to classify these as

belonging to either the core genome or transferred elements. This deviation in k-mer spectra has been explored in the context of the archaeal mobilome and contains information on the ecological niche and evolutionary history of DNA sequences¹³⁰.

**IME gene content and homology analyses.** We followed the pipeline depicted in Supplementary Fig. 8. Briefly, PCGs' IMEs were determined by a combination of k-mer spectra and genomic alignments (see Supplementary Table 4). We initially inspected fsx1-containing scaffolds and kept only sequences that were 20 kb or longer for downstream analyses. We generated an enriched annotation for each IME. Then, we obtained an initial set of groups of homologous sequences, and each of these groups was enriched by means of HMM searches. Subsequently, the enriched homology groups showing similarity between them, as judged by HMM-HMM comparisons, were collapsed into unique groups.

In detail, first, we re-annotated the identified mobile elements combining the corresponding segment of the PATRIC¹²⁵ GFP annotation file with in-house ORF predictions (minimum ORF length of 30 nucleotides, option by default). ORF inference was done by means of getorf of the EMBOSS package¹³¹, specifying genetic code by Table 11 (Bacteria and Archaea) and other parameters running by default. The similarity of inferred ORFs and annotated features in these mobile elements (i.e., features in their GFP annotation file) was established by means of BLASTP reciprocal searches⁷⁴. We kept all the predicted ORFs and homologs that were at least annotated in one genome, in this way we tried to recover misannotated conserved ORFs.

Initial sets of homologs were generated with get_homologs¹³². Sequence identity and query coverage thresholds were set to 35% and 70%, respectively. In-paralogues were not allowed within these groups (option '-e'), and remaining parameters were run by default.

HMM profiles were constructed for each homolog group. To this aim, homologous sequences were retrieved for members of each group from the UniRef50 database¹³³ with jackhmm from the HMMER package⁷³ running with one iteration ('-N 1' parameter). MSAs were then generated for each group and its relevant hits with MAFFT¹³⁴ running under '-auto' parameter, and HMMs were created with HMMER hmmbuild. Homolog groups were enriched by means of HMM searches with HMMER hmmsearch, using each HMM as a query against a database comprising all predicted ORFs described above. Hits showing an e-value < 1e-10 and covering at least 50% of the HMM were added to the groups.

Enriched homology groups showing homology were collapsed. For this purpose, HMM-vs.-HMM comparisons were performed with HHalign from the HHsuite¹³⁵. A graph was created with the networkx Python library (https://networkx.org), each node being an enriched group of homologs. An edge was established between nodes if their HMM-HMM alignment was significant (i.e., e-value < 1e-10, HMM coverage of longest HMM > = 50%). Groups of interconnected nodes were established with the 'connected_components' routine, creating a collapsed homology group in each case.

Finally, we assessed the gene content similarity between mobile elements using a Jaccard Index based on the homology groups defined above. Usual Jaccard index of two sets is defined as (# of the intersection)/(# of the union). In this case:

$$J(\text{MEA}, \text{MEB}) = \frac{\text{Nhomology groups shared between ME A:ME B}}{\text{N homol. groups MEA} + \text{N homol. groups MEB} - \text{Nhomol. groups shared between ME A:ME B}}$$

We performed a hierarchical clustering of the MEs based on a distance matrix obtained from the pairwise Jaccard Indexes (distance(A,B) = $1-J_{A,B}$). This was done in Python with seaborn¹³⁶, employing the clustermap function. A subset of 11 mobile elements (red cluster in Supplementary Fig. 9), which included ME from PCGs and JGI12330/12834-1000008 (Supplementary Data 1, 3 and 4), was selected for synteny conservation analysis. Plots depicting synteny in gene content between homolog groups were generated employing the MCscan tool¹³⁷.

HMMER and Pfam¹³⁸ were used on default parameters to assign domains and their associated arCOG¹³⁹,¹⁴⁰ identifiers to ORFs (Supplementary Data 2).

These analyses, including collapsed clusters, can be found in Zenodo⁶³.

**Sequence and structure phylogenies.** Maximum-likelihood phylogenetic trees were generated with sequences aligned with MAFFT (L-INS-i option) as input for IQ-TREE¹⁴¹ and selecting the best evolutionary model with ModelFinder¹⁴². Homology trimeric models of archaeal homologs of Fsx1E were built with MODELLER using our crystal structure as template.

Protein folds preserve deeper evolutionary signals than sequences¹⁴³⁻¹⁴⁵. Fsx1 models and crystal structures of Fsx1E and eukaryotic and viral fuxins were all-lll compared with FATCAT¹⁴⁶ to establish their structural distances between them. The following experimental crystal structures from other works were used: Flavivirus E: West Nile virus (2I69)¹⁴⁷; Dengue virus serotype 1 (4GSX)¹⁴⁸; Alphavirus E1: Semliki Forest virus (1RER)¹⁴⁹; Chikungunya virus (3N43)¹⁵⁰; C. elegans EFF-1 (4OJC)¹⁶; Bunyavirus Gc Rift Valley fever virus (6EGU)¹⁵¹; eukaryotic HAP2/GCS1 from A. thaliana (5OW3)²⁵ and C. reinhardtii (6E18)⁶⁵. The PDB files produced by flexible alignment with FATCAT were compared with TMalign¹⁵² to build a TMscore¹⁵³ distance matrix (distance = 1-TMscore). This distance matrix was the basis to compute a minimum evolution tree with FastME¹⁵⁴ on default parameters (Supplementary Fig. 10b).
**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

Crystallographic structure factors and atomic coordinates have been deposited in the Protein Data Bank under accession code 7P4L.

Source Data for this paper's structure-function analyses (Figs. 4 and 5, and Supplementary Fig. 1) are provided in the Source Data file; sequences of synthesized *fsx1* genes are in Supplementary Data 4. Source data are provided with this paper.

## Code availability

All relevant codes, notebooks and datasets necessary for: HHblits and HMMER searches and comparisons (Fig. 1a, Supplementary Fig. 1a); genomic comparisons and K-mer spectra analyses (Fig. 7); bioinformatic workflow for IMEs, IMEs clustering, content and synteny analyses (Supplementary Figs. 8 and 9, Fig. 8 and Supplementary Data 2 and 3); and structure- and sequence-based comparisons (Supplementary Fig. 10 and Supplementary Data 5), are provided on GitHub (https://github.com/DessimozLab/Archaea-Fusexins) and Zenodo (https://doi.org/10.5281/zenodo.6677729)⁶³.

Received: 27 May 2022; Accepted: 22 June 2022;
Published online: 06 July 2022

## References

1. Speijer, D., Lukeš, J. & Eliáš, M. Sex is a ubiquitous, ancient, and inherent attribute of eukaryotic life. *Proc. Natl Acad. Sci. USA* **112**, 8827–8834 (2015).

2. Ramesh, M. A., Malik, S.-B. & Logsdon, J. M. Jr. A phylogenomic inventory of meiotic genes; evidence for sex in Giardia and an early eukaryotic origin of meiosis. *Curr. Biol.* **15**, 185–191 (2005).

3. Podbilewicz, B. Virus and cell fusion mechanisms. *Annu. Rev. Cell Dev. Biol.* **30**, 111–139 (2014).

4. White, J. M., Delos, S. E., Brecher, M. & Schornberg, K. Structures and mechanisms of viral membrane fusion proteins: multiple variations on a common theme. *Crit. Rev. Biochem. Mol. Biol.* **43**, 189–219 (2008).

5. Kielian, M. Mechanisms of virus membrane fusion proteins. *Annu. Rev. Virol.* **1**, 171–189 (2014).

6. Harrison, S. C. Viral membrane fusion. *Virology* **479–480**, 498–507 (2015).

7. Millay, D. P. et al. Myomaker is a membrane activator of myoblast fusion and muscle formation. *Nature* **499**, 301–305 (2013).

8. Leikina, E. et al. Myomaker and myomerger work independently to control distinct steps of membrane remodeling during myoblast fusion. *Dev. Cell* **46**, 767–780.e7 (2018).

9. Brukman, N. G., Li, X. & Podbilewicz, B. Fusexins, HAP2/GCS1 and evolution of gamete fusion. *Front. Cell Dev. Biol.* **9**, 824024 (2021).

10. Mohler, W. A. et al. The type I membrane protein EFF-1 is essential for developmental cell fusion. *Dev. Cell* **2**, 355–362 (2002).

11. Iosilevskii, Y. & Podbilewicz, B. Programmed cell fusion in development and homeostasis. *Curr. Top. Dev. Biol.* **144**, 215–244 (2021).

12. Avinoam, O. Conserved eukaryotic fusogens can fuse viral envelopes to cells. Ph.D. thesis, Technion-Israel Institute of Technology (2012).

13. Podbilewicz, B. et al. The *C. elegans* developmental fusogen EFF-1 mediates homotypic fusion in heterologous cells and in vivo. *Dev. Cell* **11**, 471–481 (2006).

14. Sapir, A. et al. AFF-1, a FOS-1-regulated fusogen, mediates fusion of the anchor cell in *C. elegans*. *Dev. Cell* **12**, 683–698 (2002).

15. Avinoam, O. et al. Conserved eukaryotic fusogens can fuse viral envelopes to cells. *Science* **332**, 589–592 (2011).

16. Pérez-Vargas, J. et al. Structural basis of eukaryotic cell–cell fusion. *Cell* **157**, 407–419 (2014).

17. Johnson, M. A. et al. Arabidopsis hapless mutations define essential gametophytic functions. *Genetics* **168**, 971–982 (2004).

18. Mori, T., Kuroiwa, H., Higashiyama, T. & Kuroiwa, T. GENERATIVE CELL SPECIFIC 1 is essential for angiosperm fertilization. *Nat. Cell Biol.* **8**, 64–71 (2006).

19. Liu, Y. et al. The conserved plant sterility gene HAP2 functions after attachment of fusogenic membranes in *Chlamydomonas* and *Plasmodium* gametes. *Genes Dev.* **22**, 1051–1068 (2008).

20. Cole, E. S. et al. Function of the male-gamete-specific fusion protein HAP2 in a seven-sexed ciliate. *Curr. Biol.* **24**, 2168–2173 (2014).

21. Okamoto, M. et al. Two HAP2-GCS1 homologs responsible for gamete interactions in the cellular slime mold with multiple mating types: Implication for common mechanisms of sexual reproduction shared by plants and protozoa and for male-female differentiation. *Dev. Biol.* **415**, 6–13 (2016).

22. Valansi, C. et al. Arabidopsis HAP2/GCS1 is a gamete fusion protein homologous to somatic and viral fusogens. *J. Cell Biol.* **216**, 571–581 (2017).

23. Fédry, J. et al. The ancient gamete fusogen HAP2 is a Eukaryotic Class II fusion protein. *Cell* **168**, 904–915.e10 (2017).

24. Pinello, J. F. et al. Structure–function studies link Class II viral fusogens with the ancestral gamete fusion protein HAP2. *Curr. Biol.* **27**, 651–660 (2017).

25. Fédry, J. et al. Evolutionary diversification of the HAP2 membrane insertion motifs to drive gamete fusion across eukaryotes. *PLoS Biol.* **16**, e2006357 (2018).

26. Feng, J. et al. Fusion surface structure, function, and dynamics of gamete fusogen HAP2. *Elife* **7**, e39772 (2018).

27. Wong, J. L. & Johnson, M. A. Is HAP2-GCS1 an ancestral gamete fusogen? *Trends Cell Biol.* **20**, 134–141 (2010).

28. Doms, R. W. What came first—the virus or the egg? *Cell* **168**, 755–757 (2017).

29. Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. *Nature* **596**, 583–589 (2021).

30. Modis, Y., Ogata, S., Clements, D. & Harrison, S. C. Structure of the dengue virus envelope protein after membrane fusion. *Nature* **427**, 313–319 (2004).

31. Mueller, G. A. et al. Serological, genomic and structural analyses of the major mite allergen Der p 23. *Clin. Exp. Allergy* **46**, 365–376 (2016).

32. Tunyasuvunakool, K. et al. Highly accurate protein structure prediction for the human proteome. *Nature* **596**, 590–596 (2021).

33. DuBois, R. M. et al. Functional and evolutionary insight from the crystal structure of rubella virus protein E1. *Nature* **493**, 552–556 (2013).

34. Dubé, M., Etienne, L., Fels, M. & Kielian, M. Calcium-dependent Rubella virus fusion occurs in early endosomes. *J. Virol.* **90**, 6303–6313 (2016).

35. Gattegno, T. et al. Genetic control of fusion pore expansion in the epidermis of *Caenorhabditis elegans*. *Mol. Biol. Cell* **18**, 1153–1166 (2007).

36. Kemble, G. W., Danieli, T. & White, J. M. Lipid-anchored influenza hemagglutinin promotes hemifusion, not complete fusion. *Cell* **76**, 383–391 (1994).

37. Markosyan, R. M., Cohen, F. S. & Melikyan, G. B. The lipid-anchored ectodomain of influenza virus hemagglutinin (GPI-HA) is capable of inducing nonenlarging fusion pores. *Mol. Biol. Cell* **11**, 1143–1152 (2000).

38. Jones, N. A. & Geraghty, R. J. Fusion activity of lipid-anchored envelope glycoproteins of herpes simplex virus type 1. *Virology* **324**, 213–228 (2004).

39. Ludt, K. & Soppa, J. Polyploidy in halophilic archaea: regulation, evolutionary advantages, and gene conversion. *Biochem. Soc. Trans.* **47**, 933–944 (2019).

40. Turgeman-Grott, I. et al. Pervasive acquisition of CRISPR memory driven by inter-species mating of archaea can limit gene transfer and influence speciation. *Nat. Microbiol.* **4**, 177–186 (2019).

41. DeMaere, M. Z. et al. High level of intergenera gene exchange shapes the evolution of haloarchaea in an isolated Antarctic lake. *Proc. Natl Acad. Sci. USA* **110**, 16939–16944 (2013).

42. Sivabalasarma, S. et al. Analysis of cell–cell bridges in *Haloferax volcanii* using electron cryo-tomography reveal a continuous cytoplasm and S-layer. *Front. Microbiol.* **11**, 612239 (2020).

43. Rosenshine, I., Tchelet, R. & Mevarech, M. The mechanism of DNA transfer in the mating system of an archaebacterium. *Science* **245**, 1387–1389 (1989).

44. Naor, A., Lapierre, P., Mevarech, M., Papke, R. T. & Gophna, U. Low species barriers in halophilic archaea and the formation of recombinant hybrids. *Curr. Biol.* **22**, 1444–1448 (2012).

45. Wagner, A. et al. Mechanisms of gene flow in archaea. *Nat. Rev. Microbiol.* **15**, 492–501 (2017).

46. Koonin, E. V. & Krupovic, M. The depths of virus exaptation. *Curr. Opin. Virol.* **31**, 1–8 (2018).

47. Modis, Y. Relating structure to evolution in class II viral membrane fusion proteins. *Curr. Opin. Virol.* **5**, 34–41 (2014).

48. Merchant, M. et al. A bioactive phlebovirus-like envelope protein in a hornworm endogenous virus. *Sci. Adv.* **8**, eabj6894 (2022).

49. Krupovic, M. & Koonin, E. V. Multiple origins of viral capsid proteins from cellular ancestors. *Proc. Natl Acad. Sci. USA* **114**, E2401–E2410 (2017).

50. El Omari, K. et al. The structure of a prokaryotic viral envelope protein expands the landscape of membrane fusion proteins. *Nat. Commun.* **10**, 846 (2019).

51. Brochier-Armanet, C. et al. Complete-fosmid and fosmid-end sequences reveal frequent horizontal gene transfers in marine uncultured planktonic archaea. *ISME J.* **5**, 1291–1302 (2011).

52. Gabaldón, T. Origin and early evolution of the Eukaryotic cell. *Annu. Rev. Microbiol.* **75**, 631–647 (2021).

53. Liu, Y. et al. Expanded diversity of Asgard archaea and their relationships with eukaryotes. *Nature* **593**, 553–557 (2021).

54. Vosseberg, J. et al. Timing the origin of eukaryotic cellular complexity with ancient duplications. *Nat. Ecol. Evol.* **5**, 92–100 (2021).

55. Zaremba-Niedzwiedzka, K. et al. Asgard archaea illuminate the origin of eukaryotic cellular complexity. *Nature* **541**, 353–358 (2017).
56. Imachi, H. et al. Isolation of an archaeon at the prokaryote-eukaryote interface. *Nature* **577**, 519–525 (2020).

57. O'Malley, M. A., Leger, M. M., Wideman, J. G. & Ruiz-Trillo, I. Concepts of the last eukaryotic common ancestor. *Nat. Ecol. Evol.* **3**, 338–344 (2019).

58. Skejo, J. et al. Evidence for a syncytial origin of Eukaryotes from ancestral state reconstruction. *Genome Biol. Evol.* **13**, evab096 (2021).

59. Koonin, E. V. Origin of eukaryotes from within archaea, archaeal eukaryome and bursts of gene gain: eukaryogenesis just made easier? *Philos. Trans. R. Soc. Lond. B Biol. Sci.* **370**, 20140333 (2015).

60. Shalev, Y., Turgeman-Grott, I., Tamir, A., Eichler, J. & Gophna, U. Cell surface glycosylation is required for efficient fitting of *Haloferax volcanii*. *Front. Microbiol.* **8**, 1253 (2017).

61. Odell, D., Wanas, E., Yan, J. & Ghosh, H. P. Influence of membrane anchoring and cytoplasmic domains on the fusogenic activity of vesicular stomatitis virus glycoprotein G. *J. Virol.* **71**, 7996–8000 (1997).

62. Tong, S. & Compans, R. W. Oligomerization, secretion, and biological function of an anchor-free parainfluenza virus type 2 (PI2) fusion protein. *Virology* **270**, 368–376 (2000).

63. Moi, D., Langleib, M., Graña, M. & Romero, H. Discovery of archaeal fusexins homologous to eukaryotic HAP2/GCS1 gamete fusion proteins. github.com/ DessimoLab/Archaeal-Fusexins https://doi.org/10.5281/zenodo.6677729 (2022).

64. Liu, Y., Pei, J., Grishin, N. & Snell, W. J. The cytoplasmic domain of the gamete membrane fusion protein HAP2 targets the protein to the fusion site in *Chlamydomonas* and regulates the fusion reaction. *Development* **142**, 962–971 (2015).

65. Baquero, E., Fedry, J., Legrand, P., Krey, T. & Rey, F. A. Species-specific functional regions of the green alga gamete fusion protein HAP2 revealed by structural studies. *Structure* **27**, 113–124.e4 (2019).

66. Sali, A. & Blundell, T. L. Comparative protein modelling by satisfaction of spatial restraints. *J. Mol. Biol.* **234**, 779–815 (1993).

67. Malik, H. S., Henikoff, S. & Eickbush, T. H. Poised for contagion: evolutionary origins of the infectious abilities of invertebrate retroviruses. *Genome Res.* **10**, 1307–1318 (2000).

68. Yang, J. et al. The I-TASSER Suite: protein structure and function prediction. *Nat. Methods* **12**, 7–8 (2015).

69. Mirdita, M. et al. Uniclust databases of clustered and deeply annotated protein sequences and alignments. *Nucleic Acids Res.* **45**, D170–D176 (2017).

70. Remmert, M., Biegert, A., Hauser, A. & Söding, J. HHblits: lightning-fast iterative protein sequence searching by HMM-HMM alignment. *Nat. Methods* **9**, 173–175 (2011).

71. Steinegger, M. & Söding, J. Clustering huge protein sequence sets in linear time. *Nat. Commun.* **9**, 2542 (2018).

72. Sievers, F. & Higgins, D. G. Clustal Omega, accurate alignment of very large numbers of sequences. *Methods Mol. Biol.* **1079**, 105–116 (2014).

73. Eddy, S. R. Accelerated profile HMM searches. *PLoS Comput. Biol.* **7**, e1002195 (2011).

74. Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. *Nucleic Acids Res.* **25**, 3389–3402 (1997).

75. Tsirigos, K. D., Peters, C., Shu, N., Kall, L. & Elofsson, A. The TOPCONS web server for consensus prediction of membrane protein topology and signal peptides. *Nucleic Acids Res.* **43**, W401–W407 (2015).

76. Söding, J. Protein homology detection by HMM-HMM comparison. *Bioinformatics* **21**, 951–960 (2005).

77. Raj, I. et al. Structural basis of egg coat-sperm recognition at fertilization. *Cell* **169**, 1315–1326.e17 (2017).

78. Kemble, G. W., Henis, Y. I. & White, J. M. GPI- and transmembrane-anchored influenza hemagglutinin differ in structure and receptor binding activity. *J. Cell Biol.* **122**, 1253–1265 (1993).

79. DuBridge, R. B. et al. Analysis of mutation in human cells by using an Epstein-Barr virus shuttle system. *Mol. Cell. Biol.* **7**, 379–387 (1987).

80. Pernot, P. et al. Upgraded ESRF BM29 beamline for SAXS on macromolecules in solution. *J. Synchrotron Radiat.* **20**, 660–664 (2013).

81. Round, A. et al. BioSAXS Sample Changer: a robotic sample changer for rapid and reliable high-throughput X-ray solution scattering experiments. *Acta Crystallogr. D Biol. Crystallogr.* **71**, 67–75 (2015).

82. Konarev, P. V., Volkov, V. V., Sokolova, A. V., Koch, M. H. J. & Svergun, D. I. PRIMUS: a Windows PC-based system for small-angle scattering data analysis. *J. Appl. Crystallogr.* **36**, 1277–1282 (2003).

83. Manalastas-Cantos, K. et al. ATSAS 3.0: expanded functionality and new tools for small-angle scattering data analysis. *J. Appl. Crystallogr.* **54**, 343–355 (2021).

84. Svergun, D., Barberato, C. & Koch, M. H. J. CRY SOL—a program to evaluate X-ray solution scattering of biological macromolecules from atomic coordinates. *J. Appl. Crystallogr.* **28**, 768–773 (1995).

85. Franke, D. & Svergun, D. I. DAMMIF, a program for rapid ab-initio shape determination in small-angle scattering. *J. Appl. Crystallogr.* **42**, 342–346 (2009).

86. Volkov, V. V. & Svergun, D. I., IUCr. Uniqueness of ab initio shape determination in small-angle scattering. *J. Appl. Crystallogr.* **36**, 860–864 (2003).

87. Pettersen, E. F. et al. UCSF ChimeraX: structure visualization for researchers, educators, and developers. *Protein Sci.* **30**, 70–82 (2021).

88. Kidmose, R. T. et al. Namdinator—automatic molecular dynamics flexible fitting of structural models into cryo-EM and crystallography experimental maps. *IUCrJ* **6**, 526–531 (2019).

89. Nurizzo, D. et al. The ID23-1 structural biology beamline at the ESRF. *J. Synchrotron Radiat.* **13**, 227–238 (2006).

90. Kabsch, W. XDS. *Acta Crystallogr. D. Biol. Crystallogr.* **66**, 125–132 (2010).

91. Vagin, A. & Teplyakov, A. MOLREP: an automated program for molecular replacement. *Mol. Appl. Crystallogr.* **30**, 1022–1025 (1997).

92. Winn, M. D. et al. Overview of the CCP4 suite and current developments. *Acta Crystallogr. D. Biol. Crystallogr.* **67**, 235–242 (2011).

93. Matthews, B. W. Solvent content of protein crystals. *J. Mol. Biol.* **33**, 491–497 (1968).

94. Kantardjieff, K. A. & Rupp, B. Matthews coefficient probabilities: improved estimates for unit cell contents of proteins, DNA, and protein-nucleic acid complex crystals. *Protein Sci.* **12**, 1865–1871 (2003).

95. AlQuraishi, M. Machine learning in protein structure prediction. *Curr. Opin. Chem. Biol.* **65**, 1–8 (2021).

96. Millán, C. et al. Assessing the utility of CASP14 models for molecular replacement. *Proteins* **89**, 1752–1769 (2021).

97. McCoy, A. J. et al. Phaser crystallographic software. *J. Appl. Crystallogr.* **40**, 658–674 (2007).

98. Terwilliger, T. C. et al. Iterative model building, structure refinement and density modification with the PHENIX AutoBuild wizard. *Acta Crystallogr. D Biol. Crystallogr.* **64**, 61–69 (2008).

99. Chojnowski, G., Pereira, J. & Lamzin, V. S. Sequence assignment for low-resolution modelling of protein crystal structures. *Acta Crystallogr. D Struct. Biol.* **75**, 753–763 (2019).

100. Murshudov, G. N. et al. REFMAC5 for the refinement of macromolecular crystal structures. *Acta Crystallogr. D Biol. Crystallogr.* **67**, 355–367 (2011).

101. Casañal, A., Lohkamp, B. & Emsley, P. Current developments in Coot for macromolecular model building of Electron Cryo-microscopy and Crystallographic Data. *Protein Sci.* **29**, 1069–1078 (2020).

102. Croll, T. I. *ISOLDE*: a physically realistic environment for model building into low-resolution electron-density maps. *Acta Crystallogr. D Struct. Biol.* **74**, 519–530 (2018).

103. Afonine, P. V. et al. Towards automated crystallographic structure refinement with phenix.refine. *Acta Crystallogr. D Biol. Crystallogr.* **68**, 352–367 (2012).

104. Liebschner, D. et al. Macromolecular structure determination using X-rays, neutrons and electrons: recent developments in *Phenix*. *Acta Crystallogr. D Struct. Biol.* **75**, 861–877 (2019).

105. Thorn, A. D. & Sheldrick, G. M. ANODE: anomalous and heavy-atom density calculation. *J. Appl. Crystallogr.* **44**, 1285–1287 (2011).

106. Zheng, H. et al. CheckMyMetal: a macromolecular metal-binding validation tool. *Acta Crystallogr. D Biol. Crystallogr.* **73**, 223–233 (2017).

107. Williams, C. J. et al. MolProbity: more and better reference data for improved all-atom structure validation. *Protein Sci.* **27**, 293–315 (2018).

108. Gouet, P., Robert, X. & Courcelle, E. ESPrift/ENDscript: extracting and rendering sequence and 3D information from atomic structures of proteins. *Nucleic Acids Res.* **31**, 3320–3323 (2003).

109. Sonnhammer, E. L., von Heijne, G. & Krogh, A. A hidden Markov model for predicting transmembrane helices in protein sequences. *Proc. Int. Conf. Intell. Syst. Mol. Biol.* **6**, 175–182 (1998).

110. Zemla, A. LGA: a method for finding 3D similarities in protein structures. *Nucleic Acids Res.* **31**, 3370–3374 (2003).

111. Holm, L. Using Dali for protein structure comparison. *Methods Mol. Biol.* **2112**, 29–42 (2020).

112. Krissinel, E. & Henrick, K. Secondary-structure matching (SSM), a new tool for fast protein structure alignment in three dimensions. *Acta Crystallogr. D Biol. Crystallogr.* **60**, 2256–2268 (2004).

113. Kabsch, W. & Sander, C. Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. *Biopolymers* **22**, 2577–2637 (1983).

114. de Beer, T. A. P., Berka, K., Thornton, J. M. & Laskowski, R. A. PDBsum additions. *Nucleic Acids Res.* **42**, D292–D296 (2014).

115. Tina, K. G., Bhadra, R. & Srinivasan, N. PIC: Protein Interactions Calculator. *Nucleic Acids Res.* **35**, W473–W476 (2007).

116. Krissinel, E. & Henrick, K. Inference of macromolecular assemblies from crystalline state. *J. Mol. Biol.* **372**, 774–797 (2007).

117. Krieger, E. et al. Improving physical realism, stereochemistry, and side-chain accuracy in homology modeling: Four approaches that performed well in CASP8. *Proteins* **77**, 114–122 (2009).

118. Dolinsky, T. J. et al. PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations. *Nucleic Acids Res.* **35**, W522–W525 (2007).
119. Baker, N. A., Sept, D., Joseph, S., Holst, M. J. & McCammon, J. A. Electrostatics of nanostructures: application to microtubules and the ribosome. *Proc. Natl Acad. Sci. USA* **98**, 10037–10041 (2001).

120. Ashkenazy, H. et al. ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules. *Nucleic Acids Res.* **44**, W344–W350 (2016).

121. Evans, R. et al. Protein complex prediction with AlphaFold-Multimer. Preprint at bioRxiv https://doi.org/10.1101/2021.10.04.463034 (2022).

122. Edelstein, A. D. et al. Advanced methods of microscope control using µManager software. *J. Biol. Methods* **1**, e10 (2014).

123. Dunsing, V. et al. Optimal fluorescent protein tags for quantifying protein oligomerization in living cells. *Sci. Rep.* **8**, 1–12 (2018).

124. Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. *Nat. Methods* **9**, 671–675 (2012).

125. Davis, J. J. et al. The PATRIC Bioinformatics Resource Center: expanding data and analysis capabilities. *Nucleic Acids Res.* **48**, D606–D612 (2020).

126. Lu, B. & Leong, H. W. Computational methods for predicting genomic islands in microbial genomes. *Comput. Struct. Biotechnol. J.* **14**, 200–206 (2016).

127. Sievers, A. et al. K-mer content, correlation, and position analysis of genome DNA sequences for the identification of function and evolutionary features. *Genes* **8**, 122 (2017).

128. Zhou, F., Olman, V. & Xu, Y. Barcodes for genomes and applications. *BMC Bioinform.* **9**, 546 (2008).

129. Dempster, A. P., Laird, N. M. & Rubin, D. B. Maximum likelihood from incomplete data via the EM algorithm. *J. R. Stat. Soc. Ser. B Stat. Methodol.* **39**, 1–38 (1977).

130. Bize, A. et al. Exploring short k-mer profiles in cells and mobile elements from Archaea highlights the major influence of both the ecological niche and evolutionary history. *BMC Genom.* **22**, 186 (2021).

131. Rice, P., Longden, I. & Bleasby, A. EMBOSS: the European Molecular Biology Open Software Suite. *Trends Genet.* **16**, 276–277 (2000).

132. Contreras-Moreira, B. & Vinuesa, P. GET_HOMOLOGUES, a versatile software package for scalable and robust microbial pangenome analysis. *Appl. Environ. Microbiol.* **79**, 7696–7701 (2013).

133. Suzek, B. E. et al. UniRef clusters: a comprehensive and scalable alternative for improving sequence similarity searches. *Bioinformatics* **31**, 926–932 (2015).

134. Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. *Mol. Biol. Evol.* **30**, 772–780 (2013).

135. Steinegger, M. et al. HH-suite3 for fast remote homology detection and deep protein annotation. *BMC Bioinform.* **20**, 473 (2019).

136. Waskom, M. seaborn: statistical data visualization. *J. Open Source Softw.* **6**, 3021 (2021).

137. Tang, H. et al. Synteny and collinearity in plant genomes. *Science* **320**, 486–488 (2008).

138. Bateman, A. et al. The Pfam protein families database. *Nucleic Acids Res.* **32**, D138–D141 (2004).

139. Makarova, K. S., Sorokin, A. V., Novichkov, P. S., Wolf, Y. I. & Koonin, E. V. Clusters of orthologous genes for 41 archaeal genomes and implications for evolutionary genomics of archaea. *Biol. Direct* **2**, 1–20 (2007).

140. Makarova, K. S., Wolf, Y. I. & Koonin, E. V. Archaeal Clusters of Orthologous Genes (arCOGs): an update and application for analysis of shared features between thermococcales, methanococcales, and methanobacteriales. *Life* **5**, 818–840 (2015).

141. Nguyen, L.-T., Schmidt, H. A., von Haeseler, A. & Minh, B. Q. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. *Mol. Biol. Evol.* **32**, 268–274 (2015).

142. Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K. F., von Haeseler, A. & Jermihn, L. S. ModelFinder: fast model selection for accurate phylogenetic estimates. *Nat. Methods* **14**, 587–589 (2017).

143. Chothia, C. & Lesk, A. M. The relation between the divergence of sequence and structure in proteins. *EMBO J.* **5**, 823–826 (1986).

144. Holm, L. & Sander, C. Mapping the protein universe. *Science* **273**, 595–603 (1996).

145. Lundin, D., Poole, A. M., Sjöberg, B.-M. & Högbom, M. Use of structural phylogenetic networks for classification of the ferritin-like superfamily. *J. Biol. Chem.* **287**, 20565–20575 (2012).

146. Ye, Y. & Godzik, A. FATCAT: a web server for flexible structure comparison and structure similarity searching. *Nucleic Acids Res.* **32**, W582–W585 (2004).

147. Kanai, R. et al. Crystal structure of West Nile Virus envelope glycoprotein reveals viral surface epitopes. *J. Virol.* **80**, 11000–11008 (2006).

148. Klein, D. E., Choi, J. L. & Harrison, S. C. Structure of a Dengue virus envelope protein late-stage fusion intermediate. *J. Virol.* **87**, 2287–2293 (2013).

149. Gibbons, D. L. et al. Conformational change and protein-protein interactions of the fusion protein of Semliki Forest virus. *Nature* **427**, 320–325 (2004).

150. Voss, J. E. et al. Glycoprotein organization of Chikungunya virus particles revealed by X-ray crystallography. *Nature* **468**, 709–712 (2010).

151. Guardado-Calvo, P. et al. A glycerophospholipid-specific pocket in the RVFV class II fusion protein drives target membrane insertion. *Science* **358**, 663–667 (2017).

152. Zhang, Y. & Skolnick, J. TM-align: a protein structure alignment algorithm based on the TM-score. *Nucleic Acids Res.* **33**, 2302–2309 (2005).

153. Xu, J. & Zhang, Y. How significant is a protein structure similarity with TM-score = 0.5? *Bioinformatics* **26**, 889–895 (2010).

154. Lefort, V., Desper, R. & Gascuel, O. FastME 2.0: a comprehensive, accurate, and fast distance-based phylogeny inference program. *Mol. Biol. Evol.* **32**, 2798–2800 (2015).

## Acknowledgements

We thank Sonja-Verena Albers, Dan Cassel, Alejandro Colman-Lerner, Uri Gophna, Yael Iosilevskii, Shahar Lavid, Peter Walter, as well as members of our laboratories for discussion and comments on the manuscript; Jose Flores for advice on searches of fusexins in genomes; Olivier Gascuel for discussions on phylogenetics; Yoav Henis for the GPI-BHA plasmid; Kira Makarova for sharing insights on mobile elements. Part of the computations for this work was performed at the Vital-IT Center for high-performance computing of the SIB Swiss Institute of Bioinformatics. This work was supported by Beca de Doctorado Consejo Nacional de Investigaciones Científicas y Técnicas, Argentina (D.M.); Comisión Sectorial de Investigación Científica grant CSIC I + D-2020-682, Uruguay (H.R., M.G.); European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No. 844807 (N.G.B.); FOCEM, Fondo para la Convergencia Estructural del Mercosur grant COF 03/11 (M.G.); Fondo para la Investigación Científica y Tecnológica grant, Argentina (PICT-2017-0854 (P.S.A.); Israel Science Foundation grants 257/17, 2462/18, 2327/19, and 178/20 (B.P.); Knut and Alice Wallenberg Foundation grant 2018.0042, Sweden (L.J.); Novartis Foundation for Medical-Biological Research grant 17B111, Switzerland (C.D.); Programa de Desarrollo de las Ciencias Básicas, Uruguay (H.R., M.G., and M.L.); Swedish Research Council grants 2016-03999 and 2020-04936 (L.J.); Swiss Leading House for the Latin American Region (C.D. and P.S.A.); Swiss National Science Foundation grant 183723 (C.D.).

## Author contributions

B.P. conceived the experiments; performed some imaging work; designed, supervised and analyzed cell fusion experiments. C.D. helped devise analysis strategies for k-mer and phylogenetic analysis. C.V. designed, performed and analyzed cell fusion assays. Děs collected X-ray data, took part in structural analysis and validated metal substructure. D.M. carried out deep homology detection of Fsx1; designed gene content analysis strategies and phylogenetic analysis, collected sequence data; performed k-mer, functional, structural and phylogenomics analyses, built homology models; coded analysis routine pipelines. H.R. supervised the bioinformatics part of the work; designed and performed sequence phylogenetic analyses and phylogenomic surveys. K.F. made constructs of Fsx1 mutants. K.T. and J.J. generated AlphaFold2 models of monomeric Fsx1E. L.J. supervised the biochemical and structural part of the work; collected X-ray data; solved the Fsx1 structure, refined and analyzed it; generated AlphaFold-Multimer models of trimeric Fsx1. M.G. designed bioinformatic strategies; supervised bioinformatic aspects of the work; analyzed sequence and structural data. M.L. performed I.M.E. synteny analyses, phylogenetic and phylogenomic surveys. N.G.B. carried out live imaging and surface biotinylation experiments; assisted with the preparation and design of plasmids; analyzed data. P.S.A. supervised the bioinformatics part of the work and analyzed data. S.N. expressed, purified and crystallized Fsx1; performed SEC-MALS experiments; analyzed SAXS data; collected X-ray data and took part in structure determination, model building and structure analysis. X.L. designed and constructed plasmids; carried out immunofluorescence and western blots for archaeal fusexins in mammalian cells; performed imaging work. D.M., S.N., X.L., N.G.B., M.G., H.R., P.S.A., L.J., and B.P. made figures and tables. D.M., S.N., X.L., M.G., H.R., P.S.A., L.J., and B.P. wrote the manuscript, which was reviewed by all authors.

## Competing interests

J.J. has filed provisional patent applications relating to machine learning for predicting protein structures. The other authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-022-31564-1.

**Correspondence** and requests for materials should be addressed to Martin Graña, Héctor Romero, Pablo S. Aguilar, Luca Jovine or Benjamin Podbilewicz.

**Peer review information** *Nature Communications* thanks the anonymous reviewers for their contribution to the peer review of this work.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_28.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022, corrected publication 2022