Distinct functions for the paralogous RBM41 and U11/U12-65K proteins in the minor spliceosome.

Here, we identify RBM41 as a novel unique protein component of the minor spliceosome. RBM41 has no previously recognized cellular function but has been identified as a paralog of U11/U12-65K, a known unique component of the U11/U12 di-snRNP. Both proteins use their highly similar C-terminal RRMs to bind to 3'-terminal stem-loops in U12 and U6atac snRNAs with comparable affinity. Our BioID data indicate that the unique N-terminal domain of RBM41 is necessary for its association with complexes containing DHX8, an RNA helicase, which in the major spliceosome drives the release of mature mRNA from the spliceosome. Consistently, we show that RBM41 associates with excised U12-type intron lariats, is present in the U12 mono-snRNP, and is enriched in Cajal bodies, together suggesting that RBM41 functions in the post-splicing steps of the minor spliceosome assembly/disassembly cycle. This contrasts with U11/U12-65K, which uses its N-terminal region to interact with U11 snRNP during intron recognition. Finally, while RBM41 knockout cells are viable, they show alterations in U12-type 3' splice site usage. Together, our results highlight the role of the 3'-terminal stem-loop of U12 snRNA as a dynamic binding platform for the U11/U12-65K and RBM41 proteins, which function at distinct stages of the assembly/disassembly cycle.

Intracellular signaling in proto-eukaryotes evolves to alleviate regulatory conflicts of endosymbiosis.

The complex eukaryotic cell resulted from a merger between simpler prokaryotic cells, yet the role of the mitochondrial endosymbiosis with respect to other eukaryotic innovations has remained under dispute. To investigate how the regulatory challenges associated with the endosymbiotic state impacted genome and network evolution during eukaryogenesis, we study a constructive computational model where two simple cells are forced into an obligate endosymbiosis. Across multiple in silico evolutionary replicates, we observe the emergence of different mechanisms for the coordination of host and symbiont cell cycles, stabilizing the endosymbiotic relationship. In most cases, coordination is implicit, without signaling between host and symbiont. Signaling only evolves when there is leakage of regulatory products between host and symbiont. In the fittest evolutionary replicate, the host has taken full control of the symbiont cell cycle through signaling, mimicking the regulatory dominance of the nucleus over the mitochondrion that evolved during eukaryogenesis.

Uncoupled evolution of the Polycomb system and deep origin of non-canonical PRC1.

Polycomb group proteins, as part of the Polycomb repressive complexes, are essential in gene repression through chromatin compaction by canonical PRC1, mono-ubiquitylation of histone H2A by non-canonical PRC1 and tri-methylation of histone H3K27 by PRC2. Despite prevalent models emphasizing tight functional coupling between PRC1 and PRC2, it remains unclear whether this paradigm indeed reflects the evolution and functioning of these complexes. Here, we conduct a comprehensive analysis of the presence or absence of cPRC1, nPRC1 and PRC2 across the entire eukaryotic tree of life, and find that both complexes were present in the Last Eukaryotic Common Ancestor (LECA). Strikingly, ~42% of organisms contain only PRC1 or PRC2, showing that their evolution since LECA is largely uncoupled. The identification of ncPRC1-defining subunits in unicellular relatives of animals and fungi suggests ncPRC1 originated before cPRC1, and we propose a scenario for the evolution of cPRC1 from ncPRC1. Together, our results suggest that crosstalk between these complexes is a secondary development in evolution.

Obligate endosymbiosis enables genome expansion during eukaryogenesis.

The endosymbiosis of an alpha-proteobacterium that gave rise to mitochondria was one of the key events in eukaryogenesis. One striking outcome of eukaryogenesis was a much more complex cell with a large genome. Despite the existence of many alternative hypotheses for this and other patterns potentially related to endosymbiosis, a constructive evolutionary model in which these hypotheses can be studied is still lacking. Here, we present a theoretical approach in which we focus on the consequences rather than the causes of mitochondrial endosymbiosis. Using a constructive evolutionary model of cell-cycle regulation, we find that genome expansion and genome size asymmetry arise from emergent host-symbiont cell-cycle coordination. We also find that holobionts with large host and small symbiont genomes perform best on long timescales and mimic the outcome of eukaryogenesis. By designing and studying a constructive evolutionary model of obligate endosymbiosis, we uncovered some of the forces that may drive the patterns observed in nature. Our results provide a theoretical foundation for patterns related to mitochondrial endosymbiosis, such as genome size asymmetry, and reveal evolutionary outcomes that have not been considered so far, such as cell-cycle coordination without direct communication.

Evolution and implications of de novo genes in humans.

Genes and translated open reading frames (ORFs) that emerged de novo from previously non-coding sequences provide species with opportunities for adaptation. When aberrantly activated, some human-specific de novo genes and ORFs have disease-promoting properties-for instance, driving tumour growth. Thousands of putative de novo coding sequences have been described in humans, but we still do not know what fraction of those ORFs has readily acquired a function. Here, we discuss the challenges and controversies surrounding the detection, mechanisms of origin, annotation, validation and characterization of de novo genes and ORFs. Through manual curation of literature and databases, we provide a thorough table with most de novo genes reported for humans to date. We re-evaluate each locus by tracing the enabling mutations and list proposed disease associations, protein characteristics and supporting evidence for translation and protein detection. This work will support future explorations of de novo genes and ORFs in humans.

Increased Sampling and Intracomplex Homologies Favor Vertical Over Horizontal Inheritance of the Dam1 Complex.

Kinetochores connect chromosomes to spindle microtubules to ensure their correct segregation during cell division. Kinetochores of human and yeasts are largely homologous, their ability to track depolymerizing microtubules, however, is carried out by the nonhomologous complexes Ska1-C and Dam1-C, respectively. We previously reported the unique anti-correlating phylogenetic profiles of Dam1-C and Ska-C found among a wide variety of eukaryotes. Based on these profiles and the limited presence of Dam1-C, we speculated that horizontal gene transfer could have played a role in the evolutionary history of Dam1-C. Here, we present an expanded analysis of Dam1-C evolution, using additional genome as well as transcriptome sequences and recently published 3D structures. This analysis revealed a wider and more complete presence of Dam1-C in Cryptista, Rhizaria, Ichthyosporea, CRuMs, and Colponemidia. The fungal Dam1-C cryo-EM structure supports earlier hypothesized intracomplex homologies, which enables the reconstruction of rooted and unrooted phylogenies. The rooted tree of concatenated Dam1-C subunits is statistically consistent with the species tree of eukaryotes, suggesting that Dam1-C is ancient, and that the present-day phylogenetic distribution is best explained by multiple, independent losses and no horizontal gene transfer was involved. Furthermore, we investigated the ancient origin of Dam1-C via profile-versus-profile searches. Homology among 8 out of the 10 Dam1-C subunits suggests that the complex largely evolved from a single multimerizing subunit that diversified into a hetero-octameric core via stepwise subunit duplication and subfunctionalization of the subunits before the origin of the last eukaryotic common ancestor.

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.

The Sugar Metabolic Model of Aspergillus niger Can Only Be Reliably Transferred to Fungi of Its Phylum.

Fungi play a critical role in the global carbon cycle by degrading plant polysaccharides to small sugars and metabolizing them as carbon and energy sources. We mapped the well-established sugar metabolic network of Aspergillus niger to five taxonomically distant species (Aspergillus nidulans, Penicillium subrubescens, Trichoderma reesei, Phanerochaete chrysosporium and Dichomitus squalens) using an orthology-based approach. The diversity of sugar metabolism correlates well with the taxonomic distance of the fungi. The pathways are highly conserved between the three studied Eurotiomycetes (A. niger, A. nidulans, P. subrubescens). A higher level of diversity was observed between the T. reesei and A. niger, and even more so for the two Basidiomycetes. These results were confirmed by integrative analysis of transcriptome, proteome and metabolome, as well as growth profiles of the fungi growing on the corresponding sugars. In conclusion, the establishment of sugar pathway models in different fungi revealed the diversity of fungal sugar conversion and provided a valuable resource for the community, which would facilitate rational metabolic engineering of these fungi as microbial cell factories.

PTRN-1 (CAMSAP) and NOCA-2 (NINEIN) are required for microtubule polarity in Caenorhabditis elegans dendrites.

The neuronal microtubule cytoskeleton is key to establish axon-dendrite polarity. Dendrites are characterized by the presence of minus-end out microtubules. However, the mechanisms that organize these microtubules with the correct orientation are still poorly understood. Using Caenorhabditis elegans as a model system for microtubule organization, we characterized the role of 2 microtubule minus-end related proteins in this process, the microtubule minus-end stabilizing protein calmodulin-regulated spectrin-associated protein (CAMSAP/PTRN-1), and the NINEIN homologue, NOCA-2 (noncentrosomal microtubule array). We found that CAMSAP and NINEIN function in parallel to mediate microtubule organization in dendrites. During dendrite outgrowth, RAB-11-positive vesicles localized to the dendrite tip to nucleate microtubules and function as a microtubule organizing center (MTOC). In the absence of either CAMSAP or NINEIN, we observed a low penetrance MTOC vesicles mislocalization to the cell body, and a nearly fully penetrant phenotype in double mutant animals. This suggests that both proteins are important for localizing the MTOC vesicles to the growing dendrite tip to organize microtubules minus-end out. Whereas NINEIN localizes to the MTOC vesicles where it is important for the recruitment of the microtubule nucleator γ-tubulin, CAMSAP localizes around the MTOC vesicles and is cotranslocated forward with the MTOC vesicles upon dendritic growth. Together, these results indicate that microtubule nucleation from the MTOC vesicles and microtubule stabilization are both important to localize the MTOC vesicles distally to organize dendritic microtubules minus-end out.

Caenorhabditis elegans MES-3 is a highly divergent ortholog of the canonical PRC2 component SUZ12.

Polycomb Repressive Complex 2 (PRC2) catalyzes the mono-, di-, and trimethylation of histone protein H3 on lysine 27 (H3K27), which is strongly associated with transcriptionally silent chromatin. The functional core of PRC2 is highly conserved in animals and consists of four subunits. One of these, SUZ12, has not been identified in the genetic model Caenorhabditis elegans, whereas C. elegans PRC2 contains the clade-specific MES-3 protein. Through unbiased sensitive sequence similarity searches complemented by high-quality structure predictions of monomers and multimers, we here demonstrate that MES-3 is a highly divergent ortholog of SUZ12. MES-3 shares protein folds and conserved residues of key domains with SUZ12 and is predicted to interact with core PRC2 members similar to SUZ12 in human PRC2. Thus, in agreement with previous genetic and biochemical studies, we provide evidence that C. elegans contains a diverged yet evolutionary conserved core PRC2, like other animals.

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.

Evolution of Complex Regulation for Cell-Cycle Control.

Many questions remain about the interplay between adaptive and neutral processes leading to genome expansion and the evolution of cellular complexity. Genome size appears to be tightly linked to the size of the regulatory repertoire of cells (van Nimwegen E. 2003. Scaling laws in the functional content of genomes. Trends Gen. 19(9):479-484). In the context of gene regulation, we here study the interplay between adaptive and nonadaptive forces on genome and regulatory network in a computational model of cell-cycle adaptation to different environments. Starting from the well-known Caulobacter crescentus network, we report on ten replicate in silico evolution experiments where cells evolve cell-cycle control by adapting to increasingly harsh spatial habitats. We find adaptive expansion of the regulatory repertoire of cells. Having a large genome is inherently costly, but also allows for improved cell-cycle behavior. Replicates traverse different evolutionary trajectories leading to distinct eco-evolutionary strategies. In four replicates, cells evolve a generalist strategy to cope with a variety of nutrient levels; in two replicates, different specialist cells evolve for specific nutrient levels; in the remaining four replicates, an intermediate strategy evolves. These diverse evolutionary outcomes reveal the role of contingency in a system under strong selective forces. This study shows that functionality of cells depends on the combination of regulatory network topology and genome organization. For example, the positions of dosage-sensitive genes are exploited to signal to the regulatory network when replication is completed, forming a de novo evolved cell cycle checkpoint. Our results underline the importance of the integration of multiple organizational levels to understand complex gene regulation and the evolution thereof.

Phylogenetic profiling in eukaryotes: The effect of species, orthologous group, and interactome selection on protein interaction prediction.

Phylogenetic profiling in eukaryotes is of continued interest to study and predict the functional relationships between proteins. This interest is likely driven by the increased number of available diverse genomes and computational methods to infer orthologies. The evaluation of phylogenetic profiles has mainly focussed on reference genome selection in prokaryotes. However, it has been proven to be challenging to obtain high prediction accuracies in eukaryotes. As part of our recent comparison of orthology inference methods for eukaryotic genomes, we observed a surprisingly high performance for predicting interacting orthologous groups. This high performance, in turn, prompted the question of what factors influence the success of phylogenetic profiling when applied to eukaryotic genomes. Here we analyse the effect of species, orthologous group and interactome selection on protein interaction prediction using phylogenetic profiles. We select species based on the diversity and quality of the genomes and compare this supervised selection with randomly generated genome subsets. We also analyse the effect on the performance of orthologous groups defined to be in the last eukaryotic common ancestor of eukaryotes to that of orthologous groups that are not. Finally, we consider the effects of reference interactome set filtering and reference interactome species. In agreement with other studies, we find an effect of genome selection based on quality, less of an effect based on genome diversity, but a more notable effect based on the amount of information contained within the genomes. Most importantly, we find it is not merely selecting the correct genomes that is important for high prediction performance. Other choices in meta parameters such as orthologous group selection, the reference species of the interaction set, and the quality of the interaction set have a much larger impact on the performance when predicting protein interactions using phylogenetic profiles. These findings shed light on the differences in reported performance amongst phylogenetic profiles approaches, and reveal on a more fundamental level for which types of protein interactions this method has most promise when applied to eukaryotes.

Chromosomal instability by mutations in the novel minor spliceosome component CENATAC.

Aneuploidy is the leading cause of miscarriage and congenital birth defects, and a hallmark of cancer. Despite this strong association with human disease, the genetic causes of aneuploidy remain largely unknown. Through exome sequencing of patients with constitutional mosaic aneuploidy, we identified biallelic truncating mutations in CENATAC (CCDC84). We show that CENATAC is a novel component of the minor (U12-dependent) spliceosome that promotes splicing of a specific, rare minor intron subtype. This subtype is characterized by AT-AN splice sites and relatively high basal levels of intron retention. CENATAC depletion or expression of disease mutants resulted in excessive retention of AT-AN minor introns in ˜ 100 genes enriched for nucleocytoplasmic transport and cell cycle regulators, and caused chromosome segregation errors. Our findings reveal selectivity in minor intron splicing and suggest a link between minor spliceosome defects and constitutional aneuploidy in humans.

Benchmarking orthology methods using phylogenetic patterns defined at the base of Eukaryotes.

Insights into the evolution of ancestral complexes and pathways are generally achieved through careful and time-intensive manual analysis often using phylogenetic profiles of the constituent proteins. This manual analysis limits the possibility of including more protein-complex components, repeating the analyses for updated genome sets or expanding the analyses to larger scales. Automated orthology inference should allow such large-scale analyses, but substantial differences between orthologous groups generated by different approaches are observed. We evaluate orthology methods for their ability to recapitulate a number of observations that have been made with regard to genome evolution in eukaryotes. Specifically, we investigate phylogenetic profile similarity (co-occurrence of complexes), the last eukaryotic common ancestor's gene content, pervasiveness of gene loss and the overlap with manually determined orthologous groups. Moreover, we compare the inferred orthologies to each other. We find that most orthology methods reconstruct a large last eukaryotic common ancestor, with substantial gene loss, and can predict interacting proteins reasonably well when applying phylogenetic co-occurrence. At the same time, derived orthologous groups show imperfect overlap with manually curated orthologous groups. There is no strong indication of which orthology method performs better than another on individual or all of these aspects. Counterintuitively, despite the orthology methods behaving similarly regarding large-scale evaluation, the obtained orthologous groups differ vastly from one another. Availability and implementation The data and code underlying this article are available in github and/or upon reasonable request to the corresponding author: https://github.com/ESDeutekom/ComparingOrthologies.

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 find 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.

Recurrent sequence evolution after independent gene duplication.

BACKGROUND: Convergent and parallel evolution provide unique insights into the mechanisms of natural selection. Some of the most striking convergent and parallel (collectively recurrent) amino acid substitutions in proteins are adaptive, but there are also many that are selectively neutral. Accordingly, genome-wide assessment has shown that recurrent sequence evolution in orthologs is chiefly explained by nearly neutral evolution. For paralogs, more frequent functional change is expected because additional copies are generally not retained if they do not acquire their own niche. Yet, it is unknown to what extent recurrent sequence differentiation is discernible after independent gene duplications in different eukaryotic taxa. RESULTS: We develop a framework that detects patterns of recurrent sequence evolution in duplicated genes. This is used to analyze the genomes of 90 diverse eukaryotes. We find a remarkable number of families with a potentially predictable functional differentiation following gene duplication. In some protein families, more than ten independent duplications show a similar sequence-level differentiation between paralogs. Based on further analysis, the sequence divergence is found to be generally asymmetric. Moreover, about 6% of the recurrent sequence evolution between paralog pairs can be attributed to recurrent differentiation of subcellular localization. Finally, we reveal the specific recurrent patterns for the gene families Hint1/Hint2, Sco1/Sco2 and vma11/vma3. CONCLUSIONS: The presented methodology provides a means to study the biochemical underpinning of functional differentiation between paralogs. For instance, two abundantly repeated substitutions are identified between independently derived Sco1 and Sco2 paralogs. Such identified substitutions allow direct experimental testing of the biological role of these residues for the repeated functional differentiation. We also uncover a diverse set of families with recurrent sequence evolution and reveal trends in the functional and evolutionary trajectories of this hitherto understudied phenomenon.

Evolutionary Dynamics of the Spindle Assembly Checkpoint in Eukaryotes.

The tremendous diversity in eukaryotic life forms can ultimately be traced back to evolutionary modifications at the level of molecular networks. Deep understanding of these modifications will not only explain cellular diversity, but will also uncover different ways to execute similar processes and expose the evolutionary 'rules' that shape the molecular networks. Here, we review the evolutionary dynamics of the spindle assembly checkpoint (SAC), a signaling network that guards fidelity of chromosome segregation. We illustrate how the interpretation of divergent SAC systems in eukaryotic species is facilitated by combining detailed molecular knowledge of the SAC and extensive comparative genome analyses. Ultimately, expanding this to other core cellular systems and experimentally interrogating such systems in organisms from all major lineages may start outlining the routes to and eventual manifestation of the cellular diversity of eukaryotic life.

Measuring the impact of gene prediction on gene loss estimates in Eukaryotes by quantifying falsely inferred absences.

In recent years it became clear that in eukaryotic genome evolution gene loss is prevalent over gene gain. However, the absence of genes in an annotated genome is not always equivalent to the loss of genes. Due to sequencing issues, or incorrect gene prediction, genes can be falsely inferred as absent. This implies that loss estimates are overestimated and, more generally, that falsely inferred absences impact genomic comparative studies. However, reliable estimates of how prevalent this issue is are lacking. Here we quantified the impact of gene prediction on gene loss estimates in eukaryotes by analysing 209 phylogenetically diverse eukaryotic organisms and comparing their predicted proteomes to that of their respective six-frame translated genomes. We observe that 4.61% of domains per species were falsely inferred to be absent for Pfam domains predicted to have been present in the last eukaryotic common ancestor. Between phylogenetically different categories this estimate varies substantially: for clade-specific loss (ancestral loss) we found 1.30% and for species-specific loss 16.88% to be falsely inferred as absent. For BUSCO 1-to-1 orthologous families, 18.30% were falsely inferred to be absent. Finally, we showed that falsely inferred absences indeed impact loss estimates, with the number of losses decreasing by 11.78%. Our work strengthens the increasing number of studies showing that gene loss is an important factor in eukaryotic genome evolution. However, while we demonstrate that on average inferring gene absences from predicted proteomes is reliable, caution is warranted when inferring species-specific absences.

The molecular basis of monopolin recruitment to the kinetochore.

The monopolin complex is a multifunctional molecular crosslinker, which in S. pombe binds and organises mitotic kinetochores to prevent aberrant kinetochore-microtubule interactions. In the budding yeast S. cerevisiae, whose kinetochores bind a single microtubule, the monopolin complex crosslinks and mono-orients sister kinetochores in meiosis I, enabling the biorientation and segregation of homologs. Here, we show that both the monopolin complex subunit Csm1 and its binding site on the kinetochore protein Dsn1 are broadly distributed throughout eukaryotes, suggesting a conserved role in kinetochore organisation and function. We find that budding yeast Csm1 binds two conserved motifs in Dsn1, one (termed Box 1) representing the ancestral, widely conserved monopolin binding motif and a second (termed Box 2-3) with a likely role in enforcing specificity of sister kinetochore crosslinking. We find that Box 1 and Box 2-3 bind the same conserved hydrophobic cavity on Csm1, suggesting competition or handoff between these motifs. Using structure-based mutants, we also find that both Box 1 and Box 2-3 are critical for monopolin function in meiosis. We identify two conserved serine residues in Box 2-3 that are phosphorylated in meiosis and whose mutation to aspartate stabilises Csm1-Dsn1 binding, suggesting that regulated phosphorylation of these residues may play a role in sister kinetochore crosslinking specificity. Overall, our results reveal the monopolin complex as a broadly conserved kinetochore organiser in eukaryotes, which budding yeast have co-opted to mediate sister kinetochore crosslinking through the addition of a second, regulatable monopolin binding interface.