Ancestral paralogs and pseudoparalogs and their role in the emergence of the eukaryotic cell

DiffPaSS-high-performance differentiable pairing of protein sequences using soft scores

MOTIVATION

Identifying interacting partners from two sets of protein sequences has important applications in computational biology. Interacting partners share similarities across species due to their common evolutionary history, and feature correlations in amino acid usage due to the need to maintain complementary interaction interfaces. Thus, the problem of finding interacting pairs can be formulated as searching for a pairing of sequences that maximizes a sequence similarity or a coevolution score. Several methods have been developed to address this problem, applying different approximate optimization methods to different scores.

RESULTS

We introduce Differentiable Pairing using Soft Scores (DiffPaSS), a differentiable framework for flexible, fast, and hyperparameter-free optimization for pairing interacting biological sequences, which can be applied to a wide variety of scores. We apply it to a benchmark prokaryotic dataset, using mutual information and neighbor graph alignment scores. DiffPaSS outperforms existing algorithms for optimizing the same scores. We demonstrate the usefulness of our paired alignments for the prediction of protein complex structure. DiffPaSS does not require sequences to be aligned, and we also apply it to nonaligned sequences from T-cell receptors.

AVAILABILITY AND IMPLEMENTATION

A PyTorch implementation and installable Python package are available at https://github.com/Bitbol-Lab/DiffPaSS.

The emerging view on the origin and early evolution of eukaryotic cells

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.

Pairing interacting protein sequences using masked language modeling

Predicting which proteins interact together from amino acid sequences is an important task. We develop a method to pair interacting protein sequences which leverages the power of protein language models trained on multiple sequence alignments (MSAs), such as MSA Transformer and the EvoFormer module of AlphaFold. We formulate the problem of pairing interacting partners among the paralogs of two protein families in a differentiable way. We introduce a method called Differentiable Pairing using Alignment-based Language Models (DiffPALM) that solves it by exploiting the ability of MSA Transformer to fill in masked amino acids in multiple sequence alignments using the surrounding context. MSA Transformer encodes coevolution between functionally or structurally coupled amino acids within protein chains. It also captures inter-chain coevolution, despite being trained on single-chain data. Relying on MSA Transformer without fine-tuning, DiffPALM outperforms existing coevolution-based pairing methods on difficult benchmarks of shallow multiple sequence alignments extracted from ubiquitous prokaryotic protein datasets. It also outperforms an alternative method based on a state-of-the-art protein language model trained on single sequences. Paired alignments of interacting protein sequences are a crucial ingredient of supervised deep learning methods to predict the three-dimensional structure of protein complexes. Starting from sequences paired by DiffPALM substantially improves the structure prediction of some eukaryotic protein complexes by AlphaFold-Multimer. It also achieves competitive performance with using orthology-based pairing.

Exploration of the truncated cytosolic Hsp70 in plants - unveiling the diverse T1 lineage and the conserved T2 lineage

The 70-kDa heat shock proteins (Hsp70s) are chaperone proteins involved in protein folding processes. Truncated Hsp70 (Hsp70T) refers to the variant lacking a conserved C-terminal motif, which is crucial for co-chaperone interactions or protein retention. Despite their significance, the characteristics of Hsp70Ts in plants remain largely unexplored. In this study, we performed a comprehensive genome-wide analysis of 192 sequenced plant and green algae genomes to investigate the distribution and features of Hsp70Ts. Our findings unveil the widespread occurrence of Hsp70Ts across all four Hsp70 forms, including cytosolic, endoplasmic reticulum, mitochondrial, and chloroplast Hsp70s, with cytosolic Hsp70T being the most prevalent and abundant subtype. Cytosolic Hsp70T is characterized by two distinct lineages, referred to as T1 and T2. Among the investigated plant and green algae species, T1 genes were identified in approximately 60% of cases, showcasing a variable gene count ranging from one to several dozens. In contrast, T2 genes were prevalent across the majority of plant genomes, usually occurring in fewer than five gene copies per species. Sequence analysis highlights that the putative T1 proteins exhibit higher similarity to full-length cytosolic Hsp70s in comparison to T2 proteins. Intriguingly, the T2 lineage demonstrates a higher level of conservation within their protein sequences, whereas the T1 lineage presents a diverse range in the C-terminal and SBDα region, leading to categorization into four distinct subtypes. Furthermore, we have observed that T1-rich species characterized by the possession of 15 or more T1 genes exhibit an expansion of T1 genes into tandem gene clusters. The T1 gene clusters identified within the Laurales order display synteny with clusters found in a species of the Chloranthales order and another species within basal angiosperms, suggesting a conserved evolutionary relationship of T1 gene clusters among these plants. Additionally, T2 genes demonstrate distinct expression patterns in seeds and under heat stress, implying their potential roles in seed development and stress response.

Genome-Wide Identification, Evolution, and Expression Analysis of the WD40 Subfamily in Oryza Genus

The WD40 superfamily is widely found in eukaryotes and has essential subunits that serve as scaffolds for protein complexes. WD40 proteins play important regulatory roles in plant development and physiological processes, such as transcription regulation and signal transduction; it is also involved in anthocyanin biosynthesis. In rice, only OsTTG1 was found to be associated with anthocyanin biosynthesis, and evolutionary analysis of the WD40 gene family in multiple species is less studied. Here, a genome-wide analysis of the subfamily belonging to WD40-TTG1 was performed in nine AA genome species: Oryza sativa ssp. japonica, Oryza sativa ssp. indica, Oryza rufipogon, Oryza glaberrima, Oryza meridionalis, Oryza barthii, Oryza glumaepatula, Oryza nivara, and Oryza longistaminata. In this study, 383 WD40 genes in the Oryza genus were identified, and they were classified into four groups by phylogenetic analysis, with most members in group C and group D. They were found to be unevenly distributed across 12 chromosomes. A total of 39 collinear gene pairs were identified in the Oryza genus, and all were segmental duplications. WD40s had similar expansion patterns in the Oryza genus. Ka/Ks analyses indicated that they had undergone mainly purifying selection during evolution. Furthermore, WD40s in the Oryza genus have similar evolutionary patterns, so Oryza sativa ssp. indica was used as a model species for further analysis. The cis-acting elements analysis showed that many genes were related to jasmonic acid and light response. Among them, OsiWD40-26/37/42 contained elements of flavonoid synthesis, and OsiWD40-15 had MYB binding sites, indicating that they might be related to anthocyanin synthesis. The expression profile analysis at different stages revealed that most OsiWD40s were expressed in leaves, roots, and panicles. The expression of OsiWD40s was further analyzed by qRT-PCR in 9311 (indica) under various hormone treatments and abiotic stresses. OsiWD40-24 was found to be responsive to both phytohormones and abiotic stresses, suggesting that it might play an important role in plant stress resistance. And many OsiWD40s might be more involved in cold stress tolerance. These findings contribute to a better understanding of the evolution of the WD40 subfamily. The analyzed candidate genes can be used for the exploration of practical applications in rice, such as cultivar culture for colored rice, stress tolerance varieties, and morphological marker development.

Defining eukaryotes to dissect eukaryogenesis

The origin of eukaryotes is among the most contentious debates in evolutionary biology, attracting multiple seemingly incompatible theories seeking to explain the sequence in which eukaryotic characteristics were acquired. Much of the controversy arises from differing views on the defining characteristics of eukaryotes. We argue that eukaryotes should be defined phylogenetically, and that doing so clarifies where competing hypotheses of eukaryogenesis agree and how we may test among aspects of disagreement. Some hypotheses make predictions about the phylogenetic origins of eukaryotic genes and are distinguishable on that basis. However, other hypotheses differ only in the order of key evolutionary steps, like mitochondrial endosymbiosis and nuclear assembly, which cannot currently be distinguished phylogenetically. Stages within eukaryogenesis may be made identifiable through the absolute dating of gene duplicates that map to eukaryotic traits, such as in genes of host or mitochondrial origin that duplicated and diverged functionally prior to emergence of the last eukaryotic common ancestor. In this way, it may finally be possible to distinguish heat from light in the debate over eukaryogenesis.

Origin and Early Diversification of the Papain Family of Cysteine Peptidases

Peptidases of the papain family play a key role in protein degradation, regulated proteolysis, and the host-pathogen arms race. Although the papain family has been the subject of many studies, knowledge about its diversity, origin, and evolution in Eukaryota, Bacteria, and Archaea is limited; thus, we aimed to address these long-standing knowledge gaps. We traced the origin and expansion of the papain family with a phylogenomic analysis, using sequence data from numerous prokaryotic and eukaryotic proteomes, transcriptomes, and genomes. We identified the full complement of the papain family in all prokaryotic and eukaryotic lineages. Analysis of the papain family provided strong evidence for its early diversification in the ancestor of eukaryotes. We found that the papain family has undergone complex and dynamic evolution through numerous gene duplications, which produced eight eukaryotic ancestral paralogous C1A lineages during eukaryogenesis. Different evolutionary forces operated on C1A peptidases, including gene duplication, horizontal gene transfer, and gene loss. This study challenges the current understanding of the origin and evolution of the papain family and provides valuable insights into their early diversification. The findings of this comprehensive study provide guidelines for future structural and functional studies of the papain family.

The symbiotic origin of the eukaryotic cell

Eukaryogenesis represented a major evolutionary transition that led to the emergence of complex cells from simpler ancestors. For several decades, the most accepted scenario involved the evolution of an independent lineage of proto-eukaryotes endowed with an endomembrane system, including a nuclear compartment, a developed cytoskeleton and phagocytosis, which engulfed the alphaproteobacterial ancestor of mitochondria. However, the recent discovery by metagenomic and cultural approaches of Asgard archaea, which harbour many genes in common with eukaryotes and are their closest relatives in phylogenomic trees, rather supports scenarios based on the symbiosis of one Asgard-like archaeon and one or more bacteria at the origin of the eukaryotic cell. Here, we review the recent discoveries that led to this conceptual shift, briefly evoking current models of eukaryogenesis and the challenges ahead to discriminate between them and to establish a detailed, plausible scenario that accounts for the evolution of eukaryotic traits from those of their prokaryotic ancestors.

Compensatory relationship between low-complexity regions and gene paralogy in the evolution of prokaryotes

The evolution of genomes in all life forms involves two distinct, dynamic types of genomic changes: gene duplication (and loss) that shape families of paralogous genes and extension (and contraction) of low-complexity regions (LCR), which occurs through dynamics of short repeats in protein-coding genes. Although the roles of each of these types of events in genome evolution have been studied, their co-evolutionary dynamics is not thoroughly understood. Here, by analyzing a wide range of genomes from diverse bacteria and archaea, we show that LCR and paralogy represent two distinct routes of evolution that are inversely correlated. The emergence of LCR is a prominent evolutionary mechanism in fast evolving, young protein families, whereas paralogy dominates the comparatively slow evolution of old protein families. The analysis of multiple prokaryotic genomes shows that the formation of LCR is likely a widespread, transient evolutionary mechanism that temporally and locally affects also ancestral functions, but apparently, fades away with time, under mutational and selective pressures, yielding to gene paralogy. We propose that compensatory relationships between short-term and longer-term evolutionary mechanisms are universal in the evolution of life.

Integrating Phylogenetics With Intron Positions Illuminates the Origin of the Complex Spliceosome

Eukaryotic genes are characterized by the presence of introns that are removed from pre-mRNA by a spliceosome. This ribonucleoprotein complex is comprised of multiple RNA molecules and over a hundred proteins, which makes it one of the most complex molecular machines that originated during the prokaryote-to-eukaryote transition. Previous works have established that these introns and the spliceosomal core originated from self-splicing introns in prokaryotes. Yet, how the spliceosomal core expanded by recruiting many additional proteins remains largely elusive. In this study, we use phylogenetic analyses to infer the evolutionary history of 145 proteins that we could trace back to the spliceosome in the last eukaryotic common ancestor. We found that an overabundance of proteins derived from ribosome-related processes was added to the prokaryote-derived core. Extensive duplications of these proteins substantially increased the complexity of the emerging spliceosome. By comparing the intron positions between spliceosomal paralogs, we infer that most spliceosomal complexity postdates the spread of introns through the proto-eukaryotic genome. The reconstruction of early spliceosomal evolution provides insight into the driving forces behind the emergence of complexes with many proteins during eukaryogenesis.

Endosymbiotic selective pressure at the origin of eukaryotic cell biology

The dichotomy that separates prokaryotic from eukaryotic cells runs deep. The transition from pro- to eukaryote evolution is poorly understood due to a lack of reliable intermediate forms and definitions regarding the nature of the first host that could no longer be considered a prokaryote, the first eukaryotic common ancestor, FECA. The last eukaryotic common ancestor, LECA, was a complex cell that united all traits characterising eukaryotic biology including a mitochondrion. The role of the endosymbiotic organelle in this radical transition towards complex life forms is, however, sometimes questioned. In particular the discovery of the asgard archaea has stimulated discussions regarding the pre-endosymbiotic complexity of FECA. Here we review differences and similarities among models that view eukaryotic traits as isolated coincidental events in asgard archaeal evolution or, on the contrary, as a result of and in response to endosymbiosis. Inspecting eukaryotic traits from the perspective of the endosymbiont uncovers that eukaryotic cell biology can be explained as having evolved as a solution to housing a semi-autonomous organelle and why the addition of another endosymbiont, the plastid, added no extra compartments. Mitochondria provided the selective pressures for the origin (and continued maintenance) of eukaryotic cell complexity. Moreover, they also provided the energetic benefit throughout eukaryogenesis for evolving thousands of gene families unique to eukaryotes. Hence, a synthesis of the current data lets us conclude that traits such as the Golgi apparatus, the nucleus, autophagosomes, and meiosis and sex evolved as a response to the selective pressures an endosymbiont imposes.

The role of mitochondrial energetics in the origin and diversification of eukaryotes

The origin of eukaryotic cell size and complexity is often thought to have required an energy excess supplied by mitochondria. Recent observations show energy demands to scale continuously with cell volume, suggesting that eukaryotes do not have higher energetic capacity. However, respiratory membrane area scales superlinearly with the cell surface area. Furthermore, the consequences of the contrasting genomic architectures between prokaryotes and eukaryotes have not been precisely quantified. Here, we investigated (1) the factors that affect the volumes at which prokaryotes become surface area-constrained, (2) the amount of energy divested to DNA due to contrasting genomic architectures and (3) the costs and benefits of respiring symbionts. Our analyses suggest that prokaryotes are not surface area-constrained at volumes of 10⁰‒10³ µm³, the genomic architecture of extant eukaryotes is only slightly advantageous at genomes sizes of 10⁶‒10⁷ base pairs and a larger host cell may have derived a greater advantage (lower cost) from harbouring ATP-producing symbionts. This suggests that eukaryotes first evolved without the need for mitochondria since these ranges hypothetically encompass the last eukaryotic common ancestor and its relatives. Our analyses also show that larger and faster-dividing prokaryotes would have a shortage of respiratory membrane area and divest more energy into DNA. Thus, we argue that although mitochondria may not have been required by the first eukaryotes, eukaryote diversification was ultimately dependent on mitochondria.

Repeat sequences limit the effectiveness of lateral gene transfer and favored the evolution of meiotic sex in early eukaryotes

The transition from prokaryotic lateral gene transfer to eukaryotic meiotic sex is poorly understood. Phylogenetic evidence suggests that it was tightly linked to eukaryogenesis, which involved an unprecedented rise in both genome size and the density of genetic repeats. Expansion of genome size raised the severity of Muller's ratchet, while limiting the effectiveness of lateral gene transfer (LGT) at purging deleterious mutations. In principle, an increase in recombination length combined with higher rates of LGT could solve this problem. Here, we show using a computational model that this solution fails in the presence of genetic repeats prevalent in early eukaryotes. The model demonstrates that dispersed repeat sequences allow ectopic recombination, which leads to the loss of genetic information and curtails the capacity of LGT to prevent mutation accumulation. Increasing recombination length in the presence of repeat sequences exacerbates the problem. Mutational decay can only be resisted with homology along extended sequences of DNA. We conclude that the transition to homologous pairing along linear chromosomes was a key innovation in meiotic sex, which was instrumental in the expansion of eukaryotic genomes and morphological complexity.

Utilizing paralogs for phylogenetic reconstruction has the potential to increase species tree support and reduce gene tree discordance in target enrichment data

The analysis of target enrichment data in phylogenetics lacks optimization toward using paralogs for phylogenetic reconstruction. We developed a novel approach of detecting paralogs and utilizing them for phylogenetic tree inference, by retrieving both ortho- and paralogous copies and creating orthologous alignments, from which the gene trees are built. We implemented this approach in ParalogWizard and demonstrate its performance in plant groups that underwent a whole genome duplication relatively recently: the subtribe Malinae (family Rosaceae), using Angiosperms353 as well as Malinae481 probes, the genus Oritrophium (family Asteraceae), using Compositae1061 probes, and the genus Amomum (family Zingiberaceae), using Zingiberaceae1180 probes. Discriminating between orthologs and paralogs reduced gene tree discordance and increased the species tree support in the case of the Malinae, but not for Oritrophium and Amomum. This may relate to the difference in the proportion of paralogous loci between the datasets, which was highest for the Malinae. Overall, retrieving paralogs for phylogenetic reconstruction following ParalogWizard has the potential to increase the species tree support and reduce gene tree discordance in target enrichment data, particularly if the proportion of paralogous loci is high.

Production of nascent ribosome precursors within the nucleolar microenvironment of Saccharomyces cerevisiae

35S rRNA transcripts include a 5'-external transcribed spacer followed by rRNAs of the small and large ribosomal subunits. Their processing yields massive precursors that include dozens of assembly factor proteins. In Saccharomyces cerevisiae, nucleolar assembly factors form 2 coaxial layers/volumes around ribosomal DNA. Most of these factors are cyclically recruited from a latent state to an operative state, and are extensively conserved. The layers match, at least approximately, known subcompartments found in higher eukaryotic cells. ∼80% of assembly factors are essential. The number of copies of these assembly factors is comparable to the number of nascent transcripts. Moreover, they exhibit "isoelectric balance," with RNA-binding candidate "nucleator" assembly factors being notably basic. The physical properties of pre-small subunit and pre-large subunit assembly factors are similar, as are their 19 motif signatures detected by hierarchical clustering, unlike motif signatures of the 5'-external transcribed spacer rRNP. Additionally, many assembly factors lack shared motifs. Taken together with the progression of rRNP composition during subunit maturation, and the realization that the ribosomal DNA cable is initially bathed in a subunit-nonspecific assembly factor reservoir/microenvironment, we propose a "3-step subdomain assembly model": Step (1): predominantly basic assembly factors sequentially nucleate sites along nascent rRNA; Step (2): the resulting rRNPs recruit numerous less basic assembly factors along with notably basic ribosomal proteins; Step (3): rRNPs in nearby subdomains consolidate. Cleavages of rRNA then promote release of rRNPs to the nucleoplasm, likely facilitated by the persistence of assembly factors that were already associated with nucleolar precursors.

The spread of the first introns in proto-eukaryotic paralogs

Spliceosomal introns are a unique feature of eukaryotic genes. Previous studies have established that many introns were present in the protein-coding genes of the last eukaryotic common ancestor (LECA). Intron positions shared between genes that duplicated before LECA could in principle provide insight into the emergence of the first introns. In this study we use ancestral intron position reconstructions in two large sets of duplicated families to systematically identify these ancient paralogous intron positions. We found that 20-35% of introns inferred to have been present in LECA were shared between paralogs. These shared introns, which likely preceded ancient duplications, were wide spread across different functions, with the notable exception of nuclear transport. Since we observed a clear signal of pervasive intron loss prior to LECA, it is likely that substantially more introns were shared at the time of duplication than we can detect in LECA. The large extent of shared introns indicates an early origin of introns during eukaryogenesis and suggests an early origin of a nuclear structure, before most of the other complex eukaryotic features were established.

Eukarya the chimera: eukaryotes, a secondary innovation of the two domains of life?

One of the most significant events in the evolution of life is the origin of the eukaryotic cell, an increase in cellular complexity that occurred approximately 2 billion years ago. Ground-breaking research has centered around unraveling the characteristics of the Last Eukaryotic Common Ancestor (LECA) and the nuanced archaeal and bacterial contributions in eukaryogenesis, resulting in fundamental changes in our understanding of the Tree of Life. The archaeal and bacterial roles are covered by theories of endosymbiogenesis wherein an ancestral host archaeon and a bacterial endosymbiont merged to create a new complex cell type - Eukarya - and its mitochondrion. Eukarya is often regarded as a unique and distinct domain due to complex innovations not found in archaea or bacteria, despite housing a chimeric genome containing genes of both archaeal and bacterial origin. However, the discovery of complex cell machineries in recently described Asgard archaeal lineages, and the growing support for diverse bacterial gene transfers prior to and during the time of LECA, is redefining our understanding of eukaryogenesis. Indeed, the uniqueness of Eukarya, as a domain, is challenged. It is likely that many microbial syntrophies, encompassing a 'microbial village', were required to 'raise' a eukaryote during the process of eukaryogenesis.

Gene Duplications Trace Mitochondria to the Onset of Eukaryote Complexity

The last eukaryote common ancestor (LECA) possessed mitochondria and all key traits that make eukaryotic cells more complex than their prokaryotic ancestors, yet the timing of mitochondrial acquisition and the role of mitochondria in the origin of eukaryote complexity remain debated. Here, we report evidence from gene duplications in LECA indicating an early origin of mitochondria. Among 163,545 duplications in 24,571 gene trees spanning 150 sequenced eukaryotic genomes, we identify 713 gene duplication events that occurred in LECA. LECA's bacterial-derived genes include numerous mitochondrial functions and were duplicated significantly more often than archaeal-derived and eukaryote-specific genes. The surplus of bacterial-derived duplications in LECA most likely reflects the serial copying of genes from the mitochondrial endosymbiont to the archaeal host's chromosomes. Clustering, phylogenies and likelihood ratio tests for 22.4 million genes from 5,655 prokaryotic and 150 eukaryotic genomes reveal no evidence for lineage-specific gene acquisitions in eukaryotes, except from the plastid in the plant lineage. That finding, and the functions of bacterial genes duplicated in LECA, suggests that the bacterial genes in eukaryotes are acquisitions from the mitochondrion, followed by vertical gene evolution and differential loss across eukaryotic lineages, flanked by concomitant lateral gene transfer among prokaryotes. Overall, the data indicate that recurrent gene transfer via the copying of genes from a resident mitochondrial endosymbiont to archaeal host chromosomes preceded the onset of eukaryotic cellular complexity, favoring mitochondria-early over mitochondria-late hypotheses for eukaryote origin.

LUCA to LECA, the Lucacene: A model for the gigayear delay from the first prokaryote to eukaryogenesis

It is puzzling why life on Earth consisted of prokaryotes for up to 2.5 ± 0.5 billion years (Gy) before the appearance of the first eukaryotes. This period, from LUCA (Last Universal Common Ancestor) to LECA (Last Eucaryotic Common Ancestor), we have named the Lucacene, to suggest all prokaryotic descendants of LUCA before the appearance of LECA. Here we present a simple model based on horizontal gene transfer (HGT). It is the process of HGT from Bacteria to Archaea and its reverse that we wish to simulate and estimate its duration until eukaryogenesis. Rough quantitation of its parameters shows that the model may explain the long duration of the Lucacene.

Slaying the last unicorn: discovery of histones in the microalga Nanochlorum eucaryotum

Histones are the principal constituents of eukaryotic chromatin. The four core histones (H2A, H2B, H3 and H4) are conserved across sequenced eukaryotic genomes and therefore thought to be universal to eukaryotes. In the early 1980s, however, a series of biochemical investigations failed to find evidence for histones or nucleosomal structures in the microscopic green alga Nanochlorum eucaryotum. If true, derived histone loss in this lineage would constitute an exceptional case that might help us further understand the principles governing eukaryotic gene regulation. To substantiate these earlier reports of histone loss in N. eucaryotum, we sequenced, assembled and quantified its transcriptome. Following a systematic search for histone-fold domains in the assembled transcriptome, we detect orthologues to all four core histones. We also find histone mRNAs to be highly expressed, comparable to the situation in other eukaryotes. Finally, we obtain characteristic protection patterns when N. eucaryotum chromatin is subjected to micrococcal nuclease digestion, indicating widespread formation of nucleosomal complexes in vivo. We conclude that previous reports of missing histones in N. eucaryotum were mistaken. By all indications, Nanochlorum eucaryotum has histone-based chromatin characteristic of most eukaryotes.

Evolution of new enzymes by gene duplication and divergence

Thousands of new metabolic and regulatory enzymes have evolved by gene duplication and divergence since the dawn of life. New enzyme activities often originate from promiscuous secondary activities that have become important for fitness due to a change in the environment or a mutation. Mutations that make a promiscuous activity physiologically relevant can occur in the gene encoding the promiscuous enzyme itself, but can also occur elsewhere, resulting in increased expression of the enzyme or decreased competition between the native and novel substrates for the active site. If a newly useful activity is inefficient, gene duplication/amplification will set the stage for divergence of a new enzyme. Even a few mutations can increase the efficiency of a new activity by orders of magnitude. As efficiency increases, amplified gene arrays will shrink to provide two alleles, one encoding the original enzyme and one encoding the new enzyme. Ultimately, genomic rearrangements eliminate co-amplified genes and move newly evolved paralogs to a distant region of the genome.

Timing the origin of eukaryotic cellular complexity with ancient duplications

Eukaryogenesis is one of the most enigmatic evolutionary transitions, during which simple prokaryotic cells gave rise to complex eukaryotic cells. While evolutionary intermediates are lacking, gene duplications provide information on the order of events by which eukaryotes originated. Here we use a phylogenomics approach to reconstruct successive steps during eukaryogenesis. We found that gene duplications roughly doubled the proto-eukaryotic gene repertoire, with families inherited from the Asgard archaea-related host being duplicated most. By relatively timing events using phylogenetic distances we inferred that duplications in cytoskeletal and membrane trafficking families were among the earliest events, whereas most other families expanded predominantly after mitochondrial endosymbiosis. Altogether, we infer that the host that engulfed the proto-mitochondrion had some eukaryote-like complexity, which drastically increased upon mitochondrial acquisition. This scenario bridges the signs of complexity observed in Asgard archaeal genomes to the proposed role of mitochondria in triggering eukaryogenesis.

Tissue-Specific Requirement for the GINS Complex During Zebrafish Development

Efficient and accurate DNA replication is particularly critical in stem and progenitor cells for successful proliferation and survival. The replisome, an amalgam of protein complexes, is responsible for binding potential origins of replication, unwinding the double helix, and then synthesizing complimentary strands of DNA. According to current models, the initial steps of DNA unwinding and opening are facilitated by the CMG complex, which is composed of a GINS heterotetramer that connects Cdc45 with the mini-chromosome maintenance (Mcm) helicase. In this work, we provide evidence that in the absence of GINS function DNA replication is cell autonomously impaired, and we also show that gins1 and gins2 mutants exhibit elevated levels of apoptosis restricted to actively proliferating regions of the central nervous system (CNS). Intriguingly, our results also suggest that the rapid cell cycles during early embryonic development in zebrafish may not require the function of the canonical GINS complex as neither zygotic Gins1 nor Gins2 isoforms seem to be present during these stages.

Division of labour in a matrix, rather than phagocytosis or endosymbiosis, as a route for the origin of eukaryotic cells

Abstract

Two apparently irreconcilable models dominate research into the origin of eukaryotes. In one model, amitochondrial proto-eukaryotes emerged autogenously from the last universal common ancestor of all cells. Proto-eukaryotes subsequently acquired mitochondrial progenitors by the phagocytic capture of bacteria. In the second model, two prokaryotes, probably an archaeon and a bacterial cell, engaged in prokaryotic endosymbiosis, with the species resident within the host becoming the mitochondrial progenitor. Both models have limitations. A search was therefore undertaken for alternative routes towards the origin of eukaryotic cells. The question was addressed by considering classes of potential pathways from prokaryotic to eukaryotic cells based on considerations of cellular topology. Among the solutions identified, one, called here the “third-space model”, has not been widely explored. A version is presented in which an extracellular space (the third-space), serves as a proxy cytoplasm for mixed populations of archaea and bacteria to “merge” as a transitionary complex without obligatory endosymbiosis or phagocytosis and to form a precursor cell. Incipient nuclei and mitochondria diverge by division of labour. The third-space model can accommodate the reorganization of prokaryote-like genomes to a more eukaryote-like genome structure. Nuclei with multiple chromosomes and mitosis emerge as a natural feature of the model. The model is compatible with the loss of archaeal lipid biochemistry while retaining archaeal genes and provides a route for the development of membranous organelles such as the Golgi apparatus and endoplasmic reticulum. Advantages, limitations and variations of the “third-space” models are discussed.

Reviewers

This article was reviewed by Damien Devos, Buzz Baum and Michael Gray.

Evolution of P2A and P5A ATPases: ancient gene duplications and the red algal connection to green plants revisited

In a search for slowly evolving nuclear genes that may cast light on the deep evolution of plants, we carried out phylogenetic analyses of two well-characterized subfamilies of P-type pumps (P2A and P5A ATPases) from representative branches of the eukaryotic tree of life. Both P-type ATPase genes were duplicated very early in eukaryotic evolution and before the divergence of the present eukaryotic supergroups. Synapomorphies identified in the sequences provide evidence that green plants and red algae are more distantly related than are green plants and eukaryotic supergroups in which secondary or tertiary plastids are common, such as several groups belonging to the clade that includes Stramenopiles, Alveolata, Rhizaria, Cryptophyta and Haptophyta (SAR). We propose that red algae branched off soon after the first photosynthesizing eukaryote had acquired a primary plastid, while in another lineage that led to SAR, the primary plastid was lost but, in some cases, regained as a secondary or tertiary plastid.

Paralogization and New Protein Architectures in Planctomycetes Bacteria with Complex Cell Structures

Bacteria of the phylum Planctomycetes have a unique cell plan with an elaborate intracellular membrane system, thereby resembling eukaryotic cells. The origin and evolution of these remarkable features is debated. To study the evolutionary genomics of bacteria with complex cell architectures, we have resequenced the 9.2-Mb genome of the model organism Gemmata obscuriglobus and sequenced the 10-Mb genome of G. massiliana Soil9, the 7.9-Mb genome of CJuql4, and the 6.7-Mb genome of Tuwongella immobilis, all of which belong to the family Gemmataceae. A gene flux analysis of the Planctomycetes revealed a massive emergence of novel protein families at multiple nodes within the Gemmataceae. The expanded protein families have unique multidomain architectures composed of domains that are characteristic of prokaryotes, such as the sigma factor domain of extracytoplasmic sigma factors, and domains that have proliferated in eukaryotes, such as the WD40, leucine-rich repeat, tetratricopeptide repeat and Ser/Thr kinase domains. Proteins with identifiable domains in the Gemmataceae have longer lengths and linkers than proteins in most other bacteria, and the analyses suggest that these traits were ancestrally present in the Planctomycetales. A broad comparison of protein length distribution profiles revealed an overlap between the longest proteins in prokaryotes and the shortest proteins in eukaryotes. We conclude that the many similarities between proteins in the Planctomycetales and the eukaryotes are due to convergent evolution and that there is no strict boundary between prokaryotes and eukaryotes with regard to features such as gene paralogy, protein length, and protein domain composition patterns.

Characterization of Aminoacyl-tRNA Synthetases in Chromerids

Aminoacyl-tRNA synthetases (AaRSs) are enzymes that catalyze the ligation of tRNAs to amino acids. There are AaRSs specific for each amino acid in the cell. Each cellular compartment in which translation takes place (the cytosol, mitochondria, and plastids in most cases), needs the full set of AaRSs; however, individual AaRSs can function in multiple compartments due to dual (or even multiple) targeting of nuclear-encoded proteins to various destinations in the cell. We searched the genomes of the chromerids, Chromera velia and Vitrella brassicaformis, for AaRS genes: 48 genes encoding AaRSs were identified in C. velia, while only 39 AaRS genes were found in V. brassicaformis. In the latter alga, ArgRS and GluRS were each encoded by a single gene occurring in a single copy; only PheRS was found in three genes, while the remaining AaRSs were encoded by two genes. In contrast, there were nine cases for which C. velia contained three genes of a given AaRS (45% of the AaRSs), all of them representing duplicated genes, except AsnRS and PheRS, which are more likely pseudoparalogs (acquired via horizontal or endosymbiotic gene transfer). Targeting predictions indicated that AaRSs are not (or not exclusively), in most cases, used in the cellular compartment from which their gene originates. The molecular phylogenies of the AaRSs are variable between the specific types, and similar between the two investigated chromerids. While genes with eukaryotic origin are more frequently retained, there is no clear pattern of orthologous pairs between C. velia and V. brassicaformis.

Mosaic origin of the eukaryotic kinetochore

The emergence of eukaryotes from ancient prokaryotic lineages embodied a remarkable increase in cellular complexity. While prokaryotes operate simple systems to connect DNA to the segregation machinery during cell division, eukaryotes use a highly complex protein assembly known as the kinetochore. Although conceptually similar, prokaryotic segregation systems and the eukaryotic kinetochore are not homologous. Here we investigate the origins of the kinetochore before the last eukaryotic common ancestor (LECA) using phylogenetic trees, sensitive profile-versus-profile homology detection, and structural comparisons of its protein components. We show that LECA's kinetochore proteins share deep evolutionary histories with proteins involved in a few prokaryotic systems and a multitude of eukaryotic processes, including ubiquitination, transcription, and flagellar and vesicular transport systems. We find that gene duplications played a major role in shaping the kinetochore; more than half of LECA's kinetochore proteins have other kinetochore proteins as closest homologs. Some of these have no detectable homology to any other eukaryotic protein, suggesting that they arose as kinetochore-specific folds before LECA. We propose that the primordial kinetochore evolved from proteins involved in various (pre)eukaryotic systems as well as evolutionarily novel folds, after which a subset duplicated to give rise to the complex kinetochore of LECA.

Grammar of protein domain architectures

Genomes appear similar to natural language texts, and protein domains can be treated as analogs of words. To investigate the linguistic properties of genomes further, we calculated the complexity of the “protein languages” in all major branches of life and identified a nearly universal value of information gain associated with the transition from a random domain arrangement to the current protein domain architecture. An exploration of the evolutionary relationship of the protein languages identified the domain combinations that discriminate between the major branches of cellular life. We conclude that there exists a “quasi-universal grammar” of protein domains and that the nearly constant information gain we identified corresponds to the minimal complexity required to maintain a functional cell.

From an abstract, informational perspective, protein domains appear analogous to words in natural languages in which the rules of word association are dictated by linguistic rules, or grammar. Such rules exist for protein domains as well, because only a small fraction of all possible domain combinations is viable in evolution. We employ a popular linguistic technique, n-gram analysis, to probe the “proteome grammar”—that is, the rules of association of domains that generate various domain architectures of proteins. Comparison of the complexity measures of “protein languages” in major branches of life shows that the relative entropy difference (information gain) between the observed domain architectures and random domain combinations is highly conserved in evolution and is close to being a universal constant, at ∼1.2 bits. Substantial deviations from this constant are observed in only two major groups of organisms: a subset of Archaea that appears to be cells simplified to the limit, and animals that display extreme complexity. We also identify the n-grams that represent signatures of the major branches of cellular life. The results of this analysis bolster the analogy between genomes and natural language and show that a “quasi-universal grammar” underlies the evolution of domain architectures in all divisions of cellular life. The nearly universal value of information gain by the domain architectures could reflect the minimum complexity of signal processing that is required to maintain a functioning cell.

Archaeal DNA replication and repair: new genetic, biophysical and molecular tools for discovering and characterizing enzymes, pathways and mechanisms

DNA replication and repair are essential biological processes needed for the survival of all organisms. Although these processes are fundamentally conserved in the three domains, archaea, bacteria and eukarya, the proteins and complexes involved differ. The genetic and biophysical tools developed for archaea in the last several years have accelerated the study of DNA replication and repair in this domain. In this review, the current knowledge of DNA replication and repair processes in archaea will be summarized, with emphasis on the contribution of genetics and other recently developed biophysical and molecular tools, including capillary gel electrophoresis, next-generation sequencing and single-molecule approaches. How these new tools will continue to drive archaeal DNA replication and repair research will also be discussed.

Extreme Deviations from Expected Evolutionary Rates in Archaeal Protein Families

Origin of new biological functions is a complex phenomenon ranging from single-nucleotide substitutions to the gain of new genes via horizontal gene transfer or duplication. Neofunctionalization and subfunctionalization of proteins is often attributed to the emergence of paralogs that are subject to relaxed purifying selection or positive selection and thus evolve at accelerated rates. Such phenomena potentially could be detected as anomalies in the phylogenies of the respective gene families. We developed a computational pipeline to search for such anomalies in 1,834 orthologous clusters of archaeal genes, focusing on lineage-specific subfamilies that significantly deviate from the expected rate of evolution. Multiple potential cases of neofunctionalization and subfunctionalization were identified, including some ancient, house-keeping gene families, such as ribosomal protein S10, general transcription factor TFIIB and chaperone Hsp20. As expected, many cases of apparent acceleration of evolution are associated with lineage-specific gene duplication. On other occasions, long branches in phylogenetic trees correspond to horizontal gene transfer across long evolutionary distances. Significant deceleration of evolution is less common than acceleration, and the underlying causes are not well understood; functional shifts accompanied by increased constraints could be involved. Many gene families appear to be “highly evolvable,” that is, include both long and short branches. Even in the absence of precise functional predictions, this approach allows one to select targets for experimentation in search of new biology.

The effects of repeated whole genome duplication events on the evolution of cytokinin signaling pathway

BACKGROUND

It is thought that after whole-genome duplications (WGDs), a large fraction of the duplicated gene copies is lost over time while few duplicates are retained. Which factors promote survival or death of a duplicate remains unclear and the underlying mechanisms are poorly understood. According to the model of gene dosage balance, genes encoding interacting proteins are predicted to be preferentially co-retained after WGDs. Among these are genes encoding proteins involved in complexes or in signal transduction.

RESULTS

We have investigated the way that repeated WGDs during land plant evolution have affected cytokinin signaling to study patterns of gene duplicability and co-retention in this important signal transduction pathway. Through the integration of phylogenetic analyses with comparisons of genome collinearity, we have found that signal input mediated by cytokinin receptors proved to be highly conserved over long evolutionary time-scales, with receptors showing predominantly gene loss after repeated WGDs. However, the downstream elements, e,g. response regulators, were mainly retained after WGDs and thereby formed gene families in most plant lineages.

CONCLUSIONS

Gene dosage balance between the interacting components indicated by co-retention after WGDs seems to play a minor role in the evolution of cytokinin signaling pathway. Overall, core genes of cytokinin signaling show a highly heterogeneous pattern of gene retention after WGD, reflecting complex relationships between the various factors that shape the long-term fate of a duplicated gene.

Comparative in Silico Analysis of Ferric Reduction Oxidase (FRO) Genes Expression Patterns in Response to Abiotic Stresses, Metal and Hormone Applications

The ferric reduction oxidase (FRO) gene family is involved in various biological processes widely found in plants and may play an essential role in metal homeostasis, tolerance and intricate signaling networks in response to a number of abiotic stresses. Our study describes the identification, characterization and evolutionary relationships of FRO genes families. Here, total 50 FRO genes in Plantae and 15 ‘FRO like’ genes in non-Plantae were retrieved from 16 different species. The entire FRO genes have been divided into seven clades according to close similarity in biological and functional behavior. Three conserved domains were common in FRO genes while in two FROs sub genome have an extra NADPH-Ox domain, separating the function of plant FROs. OsFRO1 and OsFRO7 genes were expressed constitutively in rice plant. Real-time RT-PCR analysis demonstrated that the expression of OsFRO1 was high in flag leaf, and OsFRO7 gene expression was maximum in leaf blade and flag leaf. Both genes showed vigorous expressions level in response to different abiotic and hormones treatments. Moreover, the expression of both genes was also substantial under heavy metal stresses. OsFRO1 gene expression was triggered following 6 h under Zn, Pb, Co and Ni treatments, whereas OsFRO7 gene expression under Fe, Pb and Ni after 12 h, Zn and Cr after 6 h, and Mn and Co after 3 h treatments. These findings suggest the possible involvement of both the genes under abiotic and metal stress and the regulation of phytohormones. Therefore, our current work may provide the foundation for further functional characterization of rice FRO genes family.

Eukaryote specific folds: Part of the whole

The origin of eukaryotes is one of the central transitions in the history of life; without eukaryotes there would be no complex multicellular life. The most accepted scenarios suggest the endosymbiosis of a mitochondrial ancestor with a complex archaeon, even though the details regarding the host and the triggering factors are still being discussed. Accordingly, phylogenetic analyses have demonstrated archaeal affiliations with key informational systems, while metabolic genes are often related to bacteria, mostly to the mitochondrial ancestor. Despite of this, there exists a large number of protein families and folds found only in eukaryotes. In this study, we have analyzed structural superfamilies and folds that probably appeared during eukaryogenesis. These folds typically represent relatively small binding domains of larger multidomain proteins. They are commonly involved in biological processes that are particularly complex in eukaryotes, such as signaling, trafficking/cytoskeleton, ubiquitination, transcription and RNA processing, but according to recent studies, these processes also have prokaryotic roots. Thus the folds originating from an eukaryotic stem seem to represent accessory parts that have contributed in the expansion of several prokaryotic processes to a new level of complexity. This might have taken place as a co-evolutionary process where increasing complexity and fold innovations have supported each other.

Formation of chimeric genes with essential functions at the origin of eukaryotes

BACKGROUND

Eukaryotes evolved from the symbiotic association of at least two prokaryotic partners, and a good deal is known about the timings, mechanisms, and dynamics of these evolutionary steps. Recently, it was shown that a new class of nuclear genes, symbiogenetic genes (S-genes), was formed concomitant with endosymbiosis and the subsequent evolution of eukaryotic photosynthetic lineages. Understanding their origins and contributions to eukaryogenesis would provide insights into the ways in which cellular complexity has evolved.

RESULTS

Here, we show that chimeric nuclear genes (S-genes), built from prokaryotic domains, are critical for explaining the leap forward in cellular complexity achieved during eukaryogenesis. A total of 282 S-gene families contributed solutions to many of the challenges faced by early eukaryotes, including enhancing the informational machinery, processing spliceosomal introns, tackling genotoxicity within the cell, and ensuring functional protein interactions in a larger, more compartmentalized cell. For hundreds of S-genes, we confirmed the origins of their components (bacterial, archaeal, or generally prokaryotic) by maximum likelihood phylogenies. Remarkably, Bacteria contributed nine-fold more S-genes than Archaea, including a two-fold greater contribution to informational functions. Therefore, there is an additional, large bacterial contribution to the evolution of eukaryotes, implying that fundamental eukaryotic properties do not strictly follow the traditional informational/operational divide for archaeal/bacterial contributions to eukaryogenesis.

CONCLUSION

This study demonstrates the extent and process through which prokaryotic fragments from bacterial and archaeal genes inherited during eukaryogenesis underly the creation of novel chimeric genes with important functions.

Domestication of self-splicing introns during eukaryogenesis: the rise of the complex spliceosomal machinery

ᅟ

The spliceosome is a eukaryote-specific complex that is essential for the removal of introns from pre-mRNA. It consists of five small nuclear RNAs (snRNAs) and over a hundred proteins, making it one of the most complex molecular machineries. Most of this complexity has emerged during eukaryogenesis, a period that is characterised by a drastic increase in cellular and genomic complexity. Although not fully resolved, recent findings have started to shed some light on how and why the spliceosome originated.

In this paper we review how the spliceosome has evolved and discuss its origin and subsequent evolution in light of different general hypotheses on the evolution of complexity. Comparative analyses have established that the catalytic core of this ribonucleoprotein (RNP) complex, as well as the spliceosomal introns, evolved from self-splicing group II introns. Most snRNAs evolved from intron fragments and the essential Prp8 protein originated from the protein that is encoded by group II introns. Proteins that functioned in other RNA processes were added to this core and extensive duplications of these proteins substantially increased the complexity of the spliceosome prior to the eukaryotic diversification. The splicing machinery became even more complex in animals and plants, yet was simplified in eukaryotes with streamlined genomes. Apparently, the spliceosome did not evolve its complexity gradually, but in rapid bursts, followed by stagnation or even simplification. We argue that although both adaptive and neutral evolution have been involved in the evolution of the spliceosome, especially the latter was responsible for the emergence of an enormously complex eukaryotic splicing machinery from simple self-splicing sequences.

Reviewers

This article was reviewed by W. Ford Doolittle, Eugene V. Koonin and Vivek Anantharaman.

Conservation of Nonsense-Mediated mRNA Decay Complex Components Throughout Eukaryotic Evolution

Nonsense-mediated mRNA decay (NMD) is an essential eukaryotic process regulating transcript quality and abundance, and is involved in diverse processes including brain development and plant defenses. Although some of the NMD machinery is conserved between kingdoms, little is known about its evolution. Phosphorylation of the core NMD component UPF1 is critical for NMD and is regulated in mammals by the SURF complex (UPF1, SMG1 kinase, SMG8, SMG9 and eukaryotic release factors). However, since SMG1 is reportedly missing from the genomes of fungi and the plant Arabidopsis thaliana, it remains unclear how UPF1 is activated outside the metazoa. We used comparative genomics to determine the conservation of the NMD pathway across eukaryotic evolution. We show that SURF components are present in all major eukaryotic lineages, including fungi, suggesting that in addition to UPF1 and SMG1, SMG8 and SMG9 also existed in the last eukaryotic common ancestor, 1.8 billion years ago. However, despite the ancient origins of the SURF complex, we also found that SURF factors have been independently lost across the Eukarya, pointing to genetic buffering within the essential NMD pathway. We infer an ancient role for SURF in regulating UPF1, and the intriguing possibility of undiscovered NMD regulatory pathways.

Sugar Lego: gene composition of bacterial carbohydrate metabolism genomic loci

Background

Bacterial carbohydrate metabolism is extremely diverse, since carbohydrates serve as a major energy source and are involved in a variety of cellular processes. Bacterial genes belonging to same metabolic pathway are often co-localized in the chromosome, but it is not a strict rule. Gene co-localization in linked to co-evolution and co-regulation. This study focuses on a large-scale analysis of bacterial genomic loci related to the carbohydrate metabolism.

Results

We demonstrate that only 53% of 148,000 studied genes from over six hundred bacterial genomes are co-localized in bacterial genomes with other carbohydrate metabolism genes, which points to a significant role of singleton genes. Co-localized genes form cassettes, ranging in size from two to fifteen genes. Two major factors influencing the cassette-forming tendency are gene function and bacterial phylogeny. We have obtained a comprehensive picture of co-localization preferences of genes for nineteen major carbohydrate metabolism functional classes, over two hundred gene orthologous clusters, and thirty bacterial classes, and characterized the cassette variety in size and content among different species, highlighting a significant role of short cassettes. The preference towards co-localization of carbohydrate metabolism genes varies between 40 and 76% for bacterial taxa. Analysis of frequently co-localized genes yielded forty-five significant pairwise links between genes belonging to different functional classes. The number of such links per class range from zero to eight, demonstrating varying preferences of respective genes towards a specific chromosomal neighborhood. Genes from eleven functional classes tend to co-localize with genes from the same class, indicating an important role of clustering of genes with similar functions. At that, in most cases such co-localization does not originate from local duplication events.

Conclusions

Overall, we describe a complex web formed by evolutionary relationships of bacterial carbohydrate metabolism genes, manifested as co-localization patterns.

Reviewers

This article was reviewed by Daria V. Dibrova (A.N. Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow, Russia), nominated by Armen Mulkidjanian (University of Osnabrück, Germany), Igor Rogozin (NCBI, NLM, NIH, USA) and Yuri Wolf (NCBI, NLM, NIH, USA).

Electronic supplementary material

The online version of this article (10.1186/s13062-017-0200-7) contains supplementary material, which is available to authorized users.

The complete structure of the small-subunit processome

The small-subunit processome represents the earliest stable precursor of the eukaryotic small ribosomal subunit. Here we present the cryo-EM structure of the Saccharomyces cerevisiae small-subunit processome at an overall resolution of 3.8 Å, which provides an essentially complete near-atomic model of this assembly. In this nucleolar superstructure, 51 ribosome-assembly factors and two RNAs encapsulate the 18S rRNA precursor and 15 ribosomal proteins in a state that precedes pre-rRNA cleavage at site A1. Extended flexible proteins are employed to connect distant sites in this particle. Molecular mimicry and steric hindrance, as well as protein- and RNA-mediated RNA remodeling, are used in a concerted fashion to prevent the premature formation of the central pseudoknot and its surrounding elements within the small ribosomal subunit.

Evolution Analysis of the Aux/IAA Gene Family in Plants Shows Dual Origins and Variable Nuclear Localization Signals

The plant hormone auxin plays pivotal roles in many aspects of plant growth and development. The auxin/indole-3-acetic acid (Aux/IAA) gene family encodes short-lived nuclear proteins acting on auxin perception and signaling, but the evolutionary history of this gene family remains to be elucidated. In this study, the Aux/IAA gene family in 17 plant species covering all major lineages of plants is identified and analyzed by using multiple bioinformatics methods. A total of 434 Aux/IAA genes was found among these plant species, and the gene copy number ranges from three (Physcomitrella patens) to 63 (Glycine max). The phylogenetic analysis shows that the canonical Aux/IAA proteins can be generally divided into five major clades, and the origin of Aux/IAA proteins could be traced back to the common ancestor of land plants and green algae. Many truncated Aux/IAA proteins were found, and some of these truncated Aux/IAA proteins may be generated from the C-terminal truncation of auxin response factor (ARF) proteins. Our results indicate that tandem and segmental duplications play dominant roles for the expansion of the Aux/IAA gene family mainly under purifying selection. The putative nuclear localization signals (NLSs) in Aux/IAA proteins are conservative, and two kinds of new primordial bipartite NLSs in P. patens and Selaginella moellendorffii were discovered. Our findings not only give insights into the origin and expansion of the Aux/IAA gene family, but also provide a basis for understanding their functions during the course of evolution.

The RecJ2 protein in the thermophilic archaeon Thermoplasma acidophilum is a 3'-5' exonuclease that associates with a DNA replication complex

RecJ/cell division cycle 45 (Cdc45) proteins are widely conserved in the three domains of life, i.e. in bacteria, Eukarya, and Archaea. Bacterial RecJ is a 5'-3' exonuclease and functions in DNA repair pathways by using its 5'-3' exonuclease activity. Eukaryotic Cdc45 has no identified enzymatic activity but participates in the CMG complex, so named because it is composed of Cdc45, minichromosome maintenance protein complex (MCM) proteins 2-7, and GINS complex proteins (Sld5, Psf11-3). Eukaryotic Cdc45 and bacterial/archaeal RecJ share similar amino acid sequences and are considered functional counterparts. In Archaea, a RecJ homolog in Thermococcus kodakarensis was shown to associate with GINS and accelerate its nuclease activity and was, therefore, designated GAN (GINS-associated nuclease); however, to date, no archaeal RecJ·MCM·GINS complex has been isolated. The thermophilic archaeon Thermoplasma acidophilum has two RecJ-like proteins, designated TaRecJ1 and TaRecJ2. TaRecJ1 exhibited DNA-specific 5'-3' exonuclease activity, whereas TaRecJ2 had 3'-5' exonuclease activity and preferred RNA over DNA. TaRecJ2, but not TaRecJ1, formed a stable complex with TaGINS in a 2:1 molar ratio. Furthermore, the TaRecJ2·TaGINS complex stimulated activity of TaMCM (T. acidophilum MCM) helicase in vitro, and the TaRecJ2·TaMCM·TaGINS complex was also observed in vivo However, TaRecJ2 did not interact with TaMCM directly and was not required for the helicase activation in vitro These findings suggest that the function of archaeal RecJ in DNA replication evolved divergently from Cdc45 despite conservation of the CMG-like complex formation between Archaea and Eukarya.

Mechanisms for initiating cellular DNA replication

Cellular DNA replication factories depend on ring-shaped hexameric helicases to aid DNA synthesis by processively unzipping the parental DNA helix. Replicative helicases are loaded onto DNA by dedicated initiator, loader, and accessory proteins during the initiation of DNA replication in a tightly regulated, multistep process. We discuss here the molecular choreography of DNA replication initiation across the three domains of life, highlighting similarities and differences in the strategies used to deposit replicative helicases onto DNA and to melt the DNA helix in preparation for replisome assembly. Although initiators and loaders are phylogenetically related, the mechanisms they use for accomplishing similar tasks have diverged considerably and in an unpredictable manner.

Splendor and misery of adaptation, or the importance of neutral null for understanding evolution

The study of any biological features, including genomic sequences, typically revolves around the question: what is this for? However, population genetic theory, combined with the data of comparative genomics, clearly indicates that such a "pan-adaptationist" approach is a fallacy. The proper question is: how has this sequence evolved? And the proper null hypothesis posits that it is a result of neutral evolution: that is, it survives by sheer chance provided that it is not deleterious enough to be efficiently purged by purifying selection. To claim adaptation, the neutral null has to be falsified. The adaptationist fallacy can be costly, inducing biologists to relentlessly seek function where there is none.

A natural barrier to lateral gene transfer from prokaryotes to eukaryotes revealed from genomes: the 70 % rule

BACKGROUND

The literature harbors many claims for lateral gene transfer (LGT) from prokaryotes to eukaryotes. Such claims are typically founded in analyses of genome sequences. It is undisputed that many genes entered the eukaryotic lineage via the origin of mitochondria and the origin of plastids. Claims for lineage-specific LGT to eukaryotes outside the context of organelle origins and claims of continuous LGT to eukaryotic lineages are more problematic. If eukaryotes acquire genes from prokaryotes continuously during evolution, then sequenced eukaryote genomes should harbor evidence for recent LGT, like prokaryotic genomes do.

RESULTS

Here we devise an approach to investigate 30,358 eukaryotic sequences in the context of 1,035,375 prokaryotic homologs among 2585 phylogenetic trees containing homologs from prokaryotes and eukaryotes. Prokaryote genomes reflect a continuous process of gene acquisition and inheritance, with abundant recent acquisitions showing 80-100 % amino acid sequence identity to their phylogenetic sister-group homologs from other phyla. By contrast, eukaryote genomes show no evidence for either continuous or recent gene acquisitions from prokaryotes. We find that, in general, genes in eukaryotic genomes that share ≥70 % amino acid identity to prokaryotic homologs are genome-specific; that is, they are not found outside individual genome assemblies.

CONCLUSIONS

Our analyses indicate that eukaryotes do not acquire genes through continual LGT like prokaryotes do. We propose a 70 % rule: Coding sequences in eukaryotic genomes that share more than 70 % amino acid sequence identity to prokaryotic homologs are most likely assembly or annotation artifacts. The findings further uncover that the role of differential loss in eukaryote genome evolution has been vastly underestimated.

Evolutionary consequences of polyploidy in prokaryotes and the origin of mitosis and meiosis

Background

The origin of eukaryote-specific traits such as mitosis and sexual reproduction remains disputable. There is growing evidence that both mitosis and eukaryotic sex (i.e., the alternation of syngamy and meiosis) may have already existed in the basal eukaryotes. The mating system of the halophilic archaeon Haloferax volcanii probably represents an intermediate stage between typical prokaryotic and eukaryotic sex. H. volcanii is highly polyploid, as well as many other Archaea. Here, we use computer simulation to explore genetic and evolutionary outcomes of polyploidy in amitotic prokaryotes and its possible role in the origin of mitosis, meiosis and eukaryotic sex.

Results

Modeling suggests that polyploidy can confer strong short-term evolutionary advantage to amitotic prokaryotes. However, it also promotes the accumulation of recessive deleterious mutations and the risk of extinction in the long term, especially in highly mutagenic environment. There are several possible strategies that amitotic polyploids can use in order to reduce the genetic costs of polyploidy while retaining its benefits. Interestingly, most of these strategies resemble different components or aspects of eukaryotic sex. They include asexual ploidy cycles, equalization of genome copies by gene conversion, high-frequency lateral gene transfer between relatives, chromosome exchange coupled with homologous recombination, and the evolution of more accurate chromosome distribution during cell division (mitosis). Acquisition of mitosis by an amitotic polyploid results in chromosome diversification and specialization. Ultimately, it transforms a polyploid cell into a functionally monoploid one with multiple unique, highly redundant chromosomes. Specialization of chromosomes makes the previously evolved modes of promiscuous chromosome shuffling deleterious. This can result in selective pressure to develop accurate mechanisms of homolog pairing, and, ultimately, meiosis.

Conclusion

Emergence of mitosis and the first evolutionary steps towards eukaryotic sex could have taken place in the ancestral polyploid, amitotic proto-eukaryotes, as they were struggling to survive in the highly mutagenic environment of the Early Proterozoic shallow water microbial communities, through the succession of the following stages: (1) acquisition of high-frequency between-individual genetic exchange coupled with homologous recombination; (2) acquisition of mitosis, followed by rapid chromosome diversification and specialization; (3) evolution of homolog synapsis and meiosis. Additional evidence compatible with this scenario includes mass acquisition of new families of paralogous genes by the basal eukaryotes, and recently discovered correlation between polyploidy and the presence of histones in Archaea.

Reviewer

This article was reviewed by Eugene Koonin, Uri Gophna and Armen Mulkidjanian. For the full reviews, please go to the Reviewers' comments section.