Stephan Köstlbacher

Eukaryotic cells change their shapes, actively segregate their DNA, and contain membrane networks, facilitated by a complex cytoskeleton containing actin filaments, microtubules made from tubulin, and other components. These filaments have ancient evolutionary origins because actin- and tubulin-like proteins form prokaryotic cytoskeletons in archaea and bacteria. Bona fide eukaryotic F-actin can be traced back to crenarchaea and Asgard archaea, which are the closest known relatives of eukaryotes. A possible Asgard archaeal origin of microtubules was suggested recently with the discovery of a lokiarchaeon containing AtubAB mini microtubules that share architectural features with their eukaryotic counterparts. Using phylogenetic analyses of metagenomic data, here we report the broad occurrence of tubulins in Asgard archaea. Biochemical and structural analyses showed that one of our previously unidentified heimdallarchaeial AtubAB tubulin pairs forms four-protofilament mini microtubules that show dynamic instability and are inhibited by the tubulin drug maytansine. Our work raises the possibility that microtubule architecture and dynamics evolved in Asgard archaea prior to eukaryogenesis.

Asgard archaea played a key role in the origin of the eukaryotic cell, with extant genomes encoding relatives of diverse eukaryotic signature proteins (ESPs) involved in cellular organization. However, their often punctuated distribution and the absence of detectable homologues for many eukaryotic proteins limit our ability to reconstruct the cellular complexity of the Asgard archaeal ancestor of eukaryotes. Here we used de novo protein structure modelling and sequence similarity detection across an expanded Asgard archaeal genomic dataset to build a structural catalogue of the Asgard archaeal pangenome. We identified 908 ‘isomorphic’ ESPs—Asgard archaeal proteins with statistically enriched structural matches to eukaryotic proteins, often bridging deep sequence divergence. These isomorphic ESPs are enriched in information storage and processing roles and contain key components of the eukaryotic Vault (MVP) and Commander (COMMD) complexes, with potential roles in cellular compartmentalization and endosomal processing. These findings expand the repertoire of eukaryotic-like proteins in Asgard archaea and suggest a higher degree of eukaryote-like cellular complexity in the archaeal ancestor of eukaryotes.

Eukaryotic cells change their shapes, actively segregate their DNA and contain membrane networks, facilitated by a complex cytoskeleton containing actin filaments, microtubules made from tubulin, and other components. These filaments have ancient evolutionary origins since actin- and tubulin-like proteins form prokaryotic cytoskeletons in archaea and bacteria. Bona fide eukaryotic F-actin can be traced back to crenarchaea and Asgard archaea, which are the closest known relatives of eukaryotes. A possible Asgard archaeal origin of microtubules was suggested recently with the discovery of a lokiarchaeon containing AtubAB mini microtubules that share architectural features with their eukaryotic counterparts. Using phylogenetic analyses of metagenomic data, here we report the broad occurrence of tubulins in Asgard archaea. Biochemical and structural analyses showed that one of our newly discovered heimdallarchaeial AtubAB tubulin pairs forms four-protofilament mini-microtubules that show dynamic instability and are inhibited by the tubulin drug maytansine. Our work raises the possibility that microtubule architecture and dynamics evolved in Asgard archaea prior to eukaryogenesis.

DNA replication initiation is orchestrated in bacteria by the replication initiator DnaA. Two models for regulation of DnaA activity in Escherichia coli have been proposed: the switch between an active and inactive form, and the titration of DnaA on the chromosome. Although proposed decades ago, experimental evidence of a titration-based control mechanism is still lacking. Here, we first identified a conserved high-density region of binding motifs near the origin of replication, an advantageous trait for titration of DnaA. We then investigated the mobility of DnaA by visualising single proteins inside single cells of wild-type and deletion mutants E. coli strains, while monitoring cellular size and DNA content. Our results indicate that the chromosome of E. coli controls the free amount of DnaA in a growth rate-dependent fashion. Moreover, they address long-standing questions on the relevance of DnaA titration in stabilising DNA replication by preventing re-initiation events during slow growth.

Convergent evolution, where unrelated taxa independently evolve similar nucleotide or amino acid compositions, can introduce compositional bias into biological sequence data. Such biases distort phylogenetic inference, particularly in deep or unevenly sampled phylogenomic datasets. While composition-aware models can mitigate this issue, their computational demands often preclude their use in large-scale analyses. We present WitChi, a computationally efficient tool for identifying and removing compositionally biased alignment columns using empirical significance testing. WitChi calculates taxon-specific chi-squared (χ²) scores and compares them to null distributions derived from permutations within alignment columns that preserve the phylogenetic structure of the alignment. Sites most responsible for deviation from the expected null are iteratively pruned using one of three scoring algorithms until the bias is no longer statistically detectable. Z-scores and p-values are provided for both taxa and alignments, offering interpretable metrics of the magnitude of compositional bias. Pruning of simulated compositional heterogeneous alignments show that WitChi reliably restores correct topologies under standard, compositionally stationary models. In benchmarks, WitChi outperforms BMGE’s stationary-based trimming while scaling linearly with taxon number. Applied to the archaeal GTDB r220 dataset (5,869 taxa; 10,101 sites), WitChi completes pruning in under one hour on four CPU cores. The resulting phylogeny recovers key clades previously resolved only by in-depth analyses using complex models of sequence evolution. WitChi provides an efficient, scalable solution for detecting and removing compositional bias in phylogenomic datasets comprising thousands to tens of thousands of taxa, enabling more accurate phylogenetic inference across the tree of life.

The origin of the eukaryotic cell, with its compartmentalized nature and generally large size compared with bacterial and archaeal cells, represents a cornerstone event in the evolution of complex life on Earth. In a process referred to as eukaryogenesis, the eukaryotic cell is believed to have evolved between approximately 1.8 and 2.7 billion years ago from its archaeal ancestors, with a symbiosis with a bacterial (proto-mitochondrial) partner being a key event. In the tree of life, the branch separating the first from the last common ancestor of all eukaryotes is long and lacks evolutionary intermediates. As a result, the timing and driving forces of the emergence of complex eukaryotic features remain poorly understood. During the past decade, environmental and comparative genomic studies have revealed vital details about the identity and nature of the host cell and the proto-mitochondrial endosymbiont, enabling a critical reappraisal of hypotheses underlying the symbiotic origin of the eukaryotic cell. Here we outline our current understanding of the key players and events underlying the emergence of cellular complexity during the prokaryote-to-eukaryote transition and discuss potential avenues of future research that might provide new insights into the enigmatic origin of the eukaryotic cell.

Argonaute proteins are the central effectors of RNA-guided RNA silencing pathways in eukaryotes, playing crucial roles in gene repression and defense against viruses and transposons. Eukaryotic Argonautes are subdivided into two clades: AGOs generally facilitate miRNA- or siRNA-mediated silencing, while PIWIs generally facilitate piRNA-mediated silencing. It is currently unclear when and how Argonaute-based RNA silencing mechanisms arose and diverged during the emergence and early evolution of eukaryotes. Here, we show that in Asgard archaea, the closest prokaryotic relatives of eukaryotes, an evolutionary expansion of Argonaute proteins took place. In particular, a deep-branching PIWI protein (HrAgo1) encoded by the genome of the Lokiarchaeon ‘Candidatus Harpocratesius repetitus’ shares a common origin with eukaryotic PIWI proteins. Contrasting known prokaryotic Argonautes that use single-stranded DNA as guides and/or targets, HrAgo1 mediates RNA-guided RNA cleavage, and facilitates gene silencing when expressed in human cells and supplied with miRNA precursors. A cryo-EM structure of HrAgo1, combined with quantitative single-molecule experiments, reveals that the protein displays structural features and target-binding modes that are a mix of those of eukaryotic AGO and PIWI proteins. Thus, this deep-branching archaeal PIWI may have retained an ancestral molecular architecture that preceded the functional and mechanistic divergence of eukaryotic AGOs and PIWIs.

Methane is a potent greenhouse gas contributing to global warming. Microorganisms largely drive the biogeochemical cycling of methane, yet little is known about viral contributions to methane metabolism (MM). We analyzed 982 publicly available metagenomes from host-associated and environmental habitats containing microbial MM genes, expanding the known MM auxiliary metabolic genes (AMGs) from three to 24, including seven genes exclusive to MM pathways. These AMGs are recovered on 911 viral contigs predicted to infect 14 prokaryotic phyla including Halobacteriota, Methanobacteriota, and Thermoproteota. Of those 24, most were encoded by viruses from rumen (16/24), with substantially fewer by viruses from environmental habitats (0–7/24). To search for additional MM AMGs from an environmental habitat, we generate metagenomes from methane-rich sediments in Vrana Lake, Croatia. Therein, we find diverse viral communities, with most viruses predicted to infect methanogens and methanotrophs and some encoding 13 AMGs that can modulate host metabolisms. However, none of these AMGs directly participate in MM pathways. Together these findings suggest that the extent to which viruses use AMGs to modulate host metabolic processes (e.g., MM) varies depending on the ecological properties of the habitat in which they dwell and is not always predictable by habitat biogeochemical properties.

The generation of cyclic oligoadenylates and subsequent allosteric activation of proteins that carry sensory domains is a distinctive feature of type III CRISPR-Cas systems. In this work, we characterize a set of associated genes of a type III-B system from Haliangium ochraceum that contains two caspase-like proteases, SAVED-CHAT and PCaspase (prokaryotic caspase), co-opted from a cyclic oligonucleotide–based antiphage signaling system (CBASS). Cyclic tri–adenosine monophosphate (AMP)–induced oligomerization of SAVED-CHAT activates proteolytic activity of the CHAT domains, which specifically cleave and activate PCaspase. Subsequently, activated PCaspase cleaves a multitude of proteins, which results in a strong interference phenotype in vivo in Escherichia coli. Taken together, our findings reveal how a CRISPR-Cas–based detection of a target RNA triggers a cascade of caspase-associated proteolytic activities.

Eukaryotic Argonaute proteins achieve gene repression and defense against viruses and transposons by RNA-guided RNA silencing. By contrast, known prokaryotic Argonautes adopt single-stranded DNA as guides and/or targets, leaving the evolutionary origin of RNA-guided RNA silencing elusive. Here, we show an evolutionary expansion of Asgard archaeal Argonautes (asAgos), including the discovery of HrAgo1 from the Lokiarchaeon ‘ Candidatus Harpocratesius repetitus’ that shares a common origin with eukaryotic PIWI proteins. HrAgo1 exhibits RNA-guided RNA cleavage in vitro and RNA silencing in human cells. The cryo-EM structure of HrAgo1 combined with quantitative single-molecule experiments reveals that HrAgo1 possesses hybrid structural features and target binding modes bridging those of the eukaryotic AGO and PIWI clades. Finally, genomic evidence suggests that eukaryotic Dicer-like processing of double-stranded RNA likely emerged as a mechanism of generating guide RNA for asAgos prior to eukaryogenesis. Our study provides new insights into the evolutionary origin and plasticity of Argonaute-based RNA silencing.