Genome-scale comparative analysis of gene fusions, gene fissions, and the fungal tree of life

Diverse origins of near-identical antifreeze proteins in unrelated fish lineages provide insights into evolutionary mechanisms of new gene birth and protein sequence convergence

Determining the origins of novel genes and the genetic mechanisms underlying the emergence of new functions is challenging yet crucial for understanding evolutionary innovations. The novel fish antifreeze proteins, exemplifying convergent evolution, represent excellent opportunities to investigate the evolutionary origins and pathways of new genes. Particularly notable is the near-identical type I antifreeze proteins (AFPI) in four phylogenetically divergent fish taxa. This study tested the hypothesis of protein sequence convergence beyond functional convergence in three unrelated AFPI-bearing fish lineages, revealing different paths by which a similar protein arose from diverse genomic resources. Comprehensive comparative analyses of de novo sequenced genome of the winter flounder and grubby sculpin, available high-quality genome of the cunner, and those of 14 other relevant species found that the near-identical AFPI originated from a distinct genetic precursor in each lineage, and independently evolved coding regions for the novel ice-binding protein while retaining sequence identity in the regulatory regions with their respective ancestor. The deduced evolutionary processes and molecular mechanisms is consistent with the Innovation-Amplification-Divergence (IAD) model applicable to AFPI formation in all three lineages, a new Duplication-Degeneration-Divergence (DDD) model we propose for the sculpin lineage, and a DDD model with gene fission for the cunner lineage. This investigation illustrates the multiple ways by which a novel functional gene with sequence convergence at the protein level could evolve across divergent species, advancing our understanding of the mechanistic intricacies in new gene formation.

A crucial role for dynamic expression of components encoding the negative arm of the circadian clock

In the Neurospora circadian system, the White Collar Complex (WCC) drives expression of the principal circadian negative arm component frequency ( frq ). FRQ interacts with FRH (FRQ-interacting helicase) and CK-1 forming a stable complex that represses its own expression by inhibiting WCC. In this study, a genetic screen identified a gene, designated as brd-8 , that encodes a conserved auxiliary subunit of the NuA4 histone acetylation complex. Loss of brd-8 reduces H4 acetylation and RNA polymerase (Pol) II occupancy at frq and other known circadian genes, and leads to a long circadian period, delayed phase, and defective overt circadian output at some temperatures. In addition to strongly associating with the NuA4 histone acetyltransferase complex, BRD-8 is also found complexed with the transcription elongation regulator BYE-1. Expression of brd-8, bye-1, histone hH2Az , and several NuA4 subunits is controlled by the circadian clock, indicating that the molecular clock both regulates the basic chromatin status and is regulated by changes in chromatin. Taken together, our data identify new auxiliary elements of the fungal NuA4 complex having homology to mammalian components, which along with conventional NuA4 subunits, are required for timely and dynamic frq expression and thereby a normal and persistent circadian rhythm.

Spermidine Synthase and Saccharopine Reductase Have Co-Expression Patterns Both in Basidiomycetes with Fusion Form and Ascomycetes with Separate Form

Gene fusion is a process through which two or more distinct genes are fused into a single chimeric gene. Unlike most harmful fusion genes in cancer cells, in this study, we first found that spermidine synthetase- (SPDS, catalyst of spermidine biosynthesis) and saccharopine reductase- (SR, catalyst of the penultimate step of lysine biosynthesis) encoding genes form a natural chimeric gene, FfSpdsSr, in Flammulina filiformis. Through the cloning of full-length ORFs in different strains and the analysis of alternative splicing in developmental stages, FfSpdsSr has only one copy and unique transcript encoding chimeric SPDS-SR in F. filiformis. By an orthologous gene search of SpdsSr in more than 80 fungi, we found that the chimeric SpdsSr exists in basidiomycetes, while the two separate Spds and Sr independently exist in ascomycetes, chytridiomycetes, and oomycetes. Further, the transcript level of FfSpdsSr was investigated in different developmental stages and under some common environmental factors and stresses by RT-qPCR. The results showed that FfSpdsSr mainly up-regulated in the elongation stage and pileus development of F. filiformis, as well as under blue light, high temperature, H2O2, and MeJA treatments. Moreover, a total of 15 sets of RNA-Seq data, including 218 samples of Neurospora crassa, were downloaded from the GEO database and used to analyze the expression correlation of NcSpds and NcSr. The results showed that the separate NcSpds and NcSr shared highly similar co-expression patterns in the samples with different strains and different nutritional and environmental condition treatments. The chimeric SpdsSr in basidiomycetes and the co-expression pattern of the Spds and Sr in N. crassa indicate the special link of spermidine and lysine in fungi, which may play an important role in the growth and development of fruiting body and in response to the multiple environmental factors and abiotic stresses.

Improving Orthologous Signal and Model Fit in Datasets Addressing the Root of the Animal Phylogeny

There is conflicting evidence as to whether Porifera (sponges) or Ctenophora (comb jellies) comprise the root of the animal phylogeny. Support for either a Porifera-sister or Ctenophore-sister tree has been extensively examined in the context of model selection, taxon sampling, and outgroup selection. The influence of dataset construction is comparatively understudied. We re-examine five animal phylogeny datasets that have supported either root hypothesis using an approach designed to enrich orthologous signal in phylogenomic datasets. We find that many component orthogroups in animal datasets fail to recover major lineages as monophyletic with the exception of Ctenophora, regardless of the supported root. Enriching these datasets to retain orthogroups recovering ≥3 major lineages reduces dataset size by up to 50% while retaining underlying phylogenetic information and taxon sampling. Site-heterogeneous phylogenomic analysis of these enriched datasets recovers both Porifera-sister and Ctenophora-sister positions, even with additional constraints on outgroup sampling. Two datasets which previously supported Ctenophora-sister support Porifera-sister upon enrichment. All enriched datasets display improved model fitness under posterior predictive analysis. While not conclusively rooting animals at either Porifera or Ctenophora, we do see an increase in signal for Porifera-sister and a decrease in signal for Ctenophore-sister when data are filtered for orthologous signal. Our results indicate that dataset size and construction as well as model fit influence animal root inference.

Genes that are Used Together are More Likely to be Fused Together in Evolution by Mutational Mechanisms: A Bioinformatic Test of the Used-Fused Hypothesis

Cases of parallel or recurrent gene fusions in evolution as well as in genetic disease and cancer are difficult to explain, because unlike point mutations, they can require the repetition of a similar configuration of multiple breakpoints rather than the repetition of a single point mutation. The used-together-fused-together hypothesis holds that genes that are used together repeatedly and persistently in a specific context are more likely to undergo fusion mutation in the course of evolution for mechanistic reasons. This hypothesis offers to explain gene fusion in both evolution and disease under one umbrella. Using bioinformatic data, we tested this hypothesis against alternatives, including that all gene pairs can fuse by random mutation, but among pairs thus fused, those that had interacted previously are more likely to be favored by selection. Results show that across multiple measures of gene interaction, human genes whose orthologs are fused in one or more species are more likely to interact with each other than random pairs of genes of the same genomic distance between pair members; that an overlap exists between genes that fused in the course of evolution in non-human species and genes that undergo fusion in human cancers; and that across six primate species studied, fusions predominate over fissions and exhibit substantial evolutionary parallelism. Together, these results support the used-together-fused-together hypothesis over its alternatives. Multiple implications are discussed, including the relevance of mutational mechanisms to the evolution of genome organization, to the distribution of fitness effects of mutation, to evolutionary parallelism and more.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11692-022-09579-9.

Comparative genomics and pangenomics of vancomycin-resistant and susceptible Enterococcus faecium from Irish hospitals

Introduction. Enterococcus faecium has emerged as an important nosocomial pathogen, which is increasingly difficult to treat due to the genetic acquisition of vancomycin resistance. Ireland has a recalcitrant vancomycin-resistant bloodstream infection rate compared to other developed countries.

Hypothesis/Gap statement. Vancomycin resistance rates persist amongst E. faecium isolates from Irish hospitals. The evolutionary genomics governing these trends have not been fully elucidated.

Methodology. A set of 28 vancomycin-resistant isolates was sequenced to construct a dataset alongside 61 other publicly available Irish genomes. This dataset was extensively analysed using in silico methodologies (comparative genomics, pangenomics, phylogenetics, genotypics and comparative functional analyses) to uncover distinct evolutionary, coevolutionary and clinically relevant population trends.

Results. These results suggest that a stable (in terms of genome size, GC% and number of genes), yet genetically diverse population (in terms of gene content) of E. faecium persists in Ireland with acquired resistance arising via plasmid acquisition (vanA) or, to a lesser extent, chromosomal recombination (vanB). Population analysis revealed five clusters with one cluster partitioned into four clades which transcend isolation dates. Pangenomic and recombination analyses revealed an open (whole genome and chromosomal specific) pangenome illustrating a rampant evolutionary pattern. Comparative resistomics and virulomics uncovered distinct chromosomal and mobilomal propensity for multidrug resistance, widespread chromosomal point-mutation-mediated resistance and chromosomally harboured arsenals of virulence factors. Interestingly, a potential difference in biofilm formation strategies was highlighted by coevolutionary analysis, suggesting differential biofilm genotypes between vanA and vanB isolates.

Conclusions. These results highlight the evolutionary history of Irish E. faecium isolates and may provide insight into underlying infection dynamics in a clinical setting. Due to the apparent ease of vancomycin resistance acquisition over time, susceptible E. faecium should be concurrently reduced in Irish hospitals to mitigate potential resistant infections.

A Genome Sequence Assembly of the Phototactic and Optogenetic Model Fungus Blastocladiella emersonii Reveals a Diversified Nucleotide-Cyclase Repertoire

The chytrid fungus Blastocladiella emersonii produces spores with swimming tails (zoospores); these cells can sense and swim toward light. Interest in this species stems from ongoing efforts to develop B. emersonii as a model for understanding the evolution of phototaxis and the molecular cell biology of the associated optogenetic circuits. Here, we report a highly contiguous genome assembly and gene annotation of the B. emersonii American Type Culture Collection 22665 strain. We integrate a PacBio long-read library with an Illumina paired-end genomic sequence survey leading to an assembly of 21 contigs totaling 34.27 Mb. Using these data, we assess the diversity of sensory system encoding genes. These analyses identify a rich complement of G-protein-coupled receptors, ion transporters, and nucleotide cyclases, all of which have been diversified by domain recombination and tandem duplication. In many cases, these domain combinations have led to the fusion of a protein domain to a transmembrane domain, tying a putative signaling function to the cell membrane. This pattern is consistent with the diversification of the B. emersonii sensory-signaling systems, which likely plays a varied role in the complex life cycle of this fungus.

Biogenesis of telomerase RNA from a protein-coding mRNA precursor

Telomerase is a eukaryotic ribonucleoprotein (RNP) enzyme that adds DNA repeats onto chromosome ends to maintain genomic stability and confer cellular immortality in cancer and stem cells. The telomerase RNA (TER) component is essential for telomerase catalytic activity and provides the template for telomeric DNA synthesis. The biogenesis of TERs is extremely divergent across eukaryotic kingdoms, employing distinct types of transcription machinery and processing pathways. In ciliates and plants, TERs are transcribed by RNA polymerase III (Pol III), while animal and ascomycete fungal TERs are transcribed by RNA Pol II and share biogenesis pathways with small nucleolar RNA (snoRNA) and small nuclear RNA (snRNA), respectively. Here, we report an unprecedented messenger RNA (mRNA)-derived biogenesis pathway for the 1,291 nucleotide TER from the basidiomycete fungus Ustilago maydis. The U. maydis TER (UmTER) contains a 5'-monophosphate, distinct from the 5' 2,2,7-trimethylguanosine (TMG) cap common to animal and ascomycete fungal TERs. The mature UmTER is processed from the 3'-untranslated region (3'-UTR) of a larger RNA precursor that possesses characteristics of mRNA including a 5' 7-methyl-guanosine (m7G) cap, alternative splicing of introns, and a poly(A) tail. Moreover, this mRNA transcript encodes a protein called Early meiotic induction protein 1 (Emi1) that is conserved across dikaryotic fungi. A recombinant UmTER precursor expressed from an mRNA promoter is processed correctly to yield mature UmTER, confirming an mRNA-processing pathway for producing TER. Our findings expand the plethora of TER biogenesis mechanisms and demonstrate a pathway for producing a functional long noncoding RNA from a protein-coding mRNA precursor.

Divergent genomic trajectories predate the origin of animals and fungi

Animals and fungi have radically distinct morphologies, yet both evolved within the same eukaryotic supergroup: Opisthokonta^1,2. Here we reconstructed the trajectory of genetic changes that accompanied the origin of Metazoa and Fungi since the divergence of Opisthokonta with a dataset that includes four novel genomes from crucial positions in the Opisthokonta phylogeny. We show that animals arose only after the accumulation of genes functionally important for their multicellularity, a tendency that began in the pre-metazoan ancestors and later accelerated in the metazoan root. By contrast, the pre-fungal ancestors experienced net losses of most functional categories, including those gained in the path to Metazoa. On a broad-scale functional level, fungal genomes contain a higher proportion of metabolic genes and diverged less from the last common ancestor of Opisthokonta than did the gene repertoires of Metazoa. Metazoa and Fungi also show differences regarding gene gain mechanisms. Gene fusions are more prevalent in Metazoa, whereas a larger fraction of gene gains were detected as horizontal gene transfers in Fungi and protists, in agreement with the long-standing idea that transfers would be less relevant in Metazoa due to germline isolation^3-5. Together, our results indicate that animals and fungi evolved under two contrasting trajectories of genetic change that predated the origin of both groups. The gradual establishment of two clearly differentiated genomic contexts thus set the stage for the emergence of Metazoa and Fungi.

Annotation-free delineation of prokaryotic homology groups

Phylogenomic studies of prokaryotic taxa often assume conserved marker genes are homologous across their length. However, processes such as horizontal gene transfer or gene duplication and loss may disrupt this homology by recombining only parts of genes, causing gene fission or fusion. We show using simulation that it is necessary to delineate homology groups in a set of bacterial genomes without relying on gene annotations to define the boundaries of homologous regions. To solve this problem, we have developed a graph-based algorithm to partition a set of bacterial genomes into Maximal Homologous Groups of sequences (MHGs) where each MHG is a maximal set of maximum-length sequences which are homologous across the entire sequence alignment. We applied our algorithm to a dataset of 19 Enterobacteriaceae species and found that MHGs cover much greater proportions of genomes than markers and, relatedly, are less biased in terms of the functions of the genes they cover. We zoomed in on the correlation between each individual marker and their overlapping MHGs, and show that few phylogenetic splits supported by the markers are supported by the MHGs while many marker-supported splits are contradicted by the MHGs. A comparison of the species tree inferred from marker genes with the species tree inferred from MHGs suggests that the increased bias and lack of genome coverage by markers causes incorrect inferences as to the overall relationship between bacterial taxa.

Assuming genes to be the basic evolutionary unit has been commonplace in bacterial genomics. For example, when quantifying the extent of horizontal gene transfer it is common to infer gene trees and reconcile them against a species tree to account for recombination-based processes. We have developed a new method which challenges this assumption by identifying contiguous regions of true homology without regards to gene boundaries and applied it to Enterobacteriaceae, a family of bacteria containing several important human pathogens. Our results show that genes are composed of distinct homologous regions with conflicting phylogenetic histories. We further demonstrate that failing to take account of this conflict, together with the functional biases we show exist among single-copy marker genes, significantly changes the consensus evolutionary tree of Enterobacteriaceae.

A light-sensing system in the common ancestor of the fungi

Diverse light-sensing organs (i.e., eyes) have evolved across animals. Interestingly, several subcellular analogs have been found in eukaryotic microbes.¹ All of these systems have a common “recipe”: a light occluding or refractory surface juxtaposed to a membrane-layer enriched in type I rhodopsins.1, 2, 3, 4 In the fungi, several lineages have been shown to detect light using a diversity of non-homologous photo-responsive proteins.5, 6, 7 However, these systems are not associated with an eyespot-like organelle with one exception found in the zoosporic fungus Blastocladiella emersonii (Be).⁸Be possesses both elements of this recipe: an eyespot composed of lipid-filled structures (often called the side-body complex [SBC]), co-localized with a membrane enriched with a gene-fusion protein composed of a type I (microbial) rhodopsin and guanylyl cyclase enzyme domain (CyclOp-fusion protein).^8,9 Here, we identify homologous pathway components in four Chytridiomycota orders (Chytridiales, Synchytriales, Rhizophydiales, and Monoblepharidiales). To further explore the architecture of the fungal zoospore and its lipid organelles, we reviewed electron microscopy data (e.g., the works of Barr and Hartmann¹⁰ and Reichle and Fuller¹¹) and performed fluorescence-microscopy imaging of four CyclOp-carrying zoosporic fungal species, showing the presence of a variety of candidate eyespot-cytoskeletal ultrastructure systems. We then assessed the presence of canonical photoreceptors across the fungi and inferred that the last common fungal ancestor was able to sense light across a range of wavelengths using a variety of systems, including blue-green-light detection. Our data imply, independently of how the fungal tree of life is rooted, that the apparatus for a CyclOp-organelle light perception system was an ancestral feature of the fungi.

•

A wide diversity of flagellated fungi possess the CyclOp light response circuit

•

The same fungi possess the subcellular equipment to build lipid-based eyespots

•

The last common ancestor of fungi possessed the CyclOp eyespot system

•

The ancestral fungus could see a rainbow of light wavelengths

Molecular and cellular dynamics of early embryonic cell divisions in Volvox carteri

Cell division is fundamental to all organisms and the green alga used here exhibits both key animal and plant functions. Specifically, we analyzed the molecular and cellular dynamics of early embryonic divisions of the multicellular green alga Volvox carteri (Chlamydomonadales). Relevant proteins related to mitosis and cytokinesis were identified in silico, the corresponding genes were cloned, fused to yfp, and stably expressed in Volvox, and the tagged proteins were studied by live-cell imaging. We reveal rearrangements of the microtubule cytoskeleton during centrosome separation, spindle formation, establishment of the phycoplast, and generation of previously unknown structures. The centrosomes participate in initiation of spindle formation and determination of spindle orientation. Although the nuclear envelope does not break down during early mitosis, intermixing of cytoplasm and nucleoplasm results in loss of nuclear identity. Finally, we present a model for mitosis in Volvox. Our study reveals enormous dynamics, clarifies spatio-temporal relationships of subcellular structures, and provides insight into the evolution of cell division.

Hundreds of out-of-frame remodelled gene families in the E. coli pangenome

All genomes include gene families with very limited taxonomic distributions that potentially represent new genes and innovations in protein-coding sequence, raising questions on the origins of such genes. Some of these genes are hypothesized to have formed de novo, from noncoding sequences, and recent work has begun to elucidate the processes by which de novo gene formation can occur. A special case of de novo gene formation, overprinting, describes the origin of new genes from noncoding alternative reading frames of existing open reading frames (ORFs). We argue that additionally, out-of-frame gene fission/fusion events of alternative reading frames of ORFs and out-of-frame lateral gene transfers could contribute to the origin of new gene families. To demonstrate this, we developed an original pattern-search in sequence similarity networks, enhancing the use of these graphs, commonly used to detect in-frame remodeled genes. We applied this approach to gene families in 524 complete genomes of Escherichia coli. We identified 767 gene families whose evolutionary history likely included at least one out-of-frame remodeling event. These genes with out-of-frame components represent ∼2.5% of all genes in the E. coli pangenome, suggesting that alternative reading frames of existing ORFs can contribute to a significant proportion of de novo genes in bacteria.

Resolution of deep divergence of club fungi (phylum Basidiomycota)

A long-standing question about the early evolution of club fungi (phylum Basidiomycota) is the relationship between the three major groups, Pucciniomycotina, Ustilaginomycotina and Agaricomycotina. It is unresolved whether Agaricomycotina are more closely related to Ustilaginomycotina or to Pucciniomycotina. Here we reconstructed the branching order of the three subphyla through two sources of phylogenetic signals, i.e. standard phylogenomic analysis and alignment-free phylogenetic approach. Overall, beyond congruency within the frame of standard phylogenomic analysis, our results consistently and robustly supported the early divergence of Ustilaginomycotina and a closer relationship between Agaricomycotina and Pucciniomycotina.

The Early Asexual Development Regulator fluG Codes for a Putative Bifunctional Enzyme

FluG is a long recognized early regulator of asexual development in Aspergillus nidulans. fluG null mutants show profuse aerial growth and no conidial production. Initial studies reported sequence homology of FluG with a prokaryotic type I glutamine synthetase, but catalytic activity has not been demonstrated. In this study, we conducted an in-depth analysis of the FluG sequence, which revealed a single polypeptide containing a putative N-terminal amidohydrolase region linked to a putative C-terminal γ-glutamyl ligase region. Each region corresponded, separately and completely, to respective single function bacterial enzymes. Separate expression of these regions confirmed that the C-terminal region was essential for asexual development. The N-terminal region alone did not support conidial development, but contributed to increased conidial production under high nutrient availability. Point mutations directed at respective key catalytic residues in each region demonstrated that they were essential for biological function. Moreover, the substitution of the N- and C-terminal regions with homologs from Lactobacillus paracasei and Pseudomonas aeruginosa, respectively, maintained functionality, albeit with altered characteristics. Taken together, the results lead us to conclude that FluG is a bifunctional enzyme that participates in an as yet unidentified metabolic or signaling pathway involving a γ-glutamylated intermediate that contributes to developmental fate.

Genome mining and biosynthesis of a polyketide from a biofertilizer fungus that can facilitate reductive iron assimilation in plant

Fungi have the potential to produce a large repertoire of bioactive molecules, many of which can affect the growth and development of plants. Genomic survey of sequenced biofertilizer fungi showed many secondary metabolite gene clusters are anchored by iterative polyketide synthases (IPKSs), which are multidomain enzymes noted for generating diverse small molecules. Focusing on the biofertilizer Trichoderma harzianum t-22, we identified and characterized a cryptic IPKS-containing cluster that synthesizes tricholignan A, a redox-active ortho-hydroquinone. Tricholignan A is shown to reduce Fe(III) and may play a role in promoting plant growth under iron-deficient conditions. The construction of tricholignan by a pair of collaborating IPKSs was investigated using heterologous reconstitution and biochemical studies. A regioselective methylation step is shown to be a key step in formation of the ortho-hydroquinone. The responsible methyltransferase (MT) is fused with an N-terminal pseudo-acyl carrier protein (ψACP), in which the apo state of the ACP is essential for methylation of the growing polyketide chain. The ψACP is proposed to bind to the IPKS and enable the trans MT to access the growing polyketide. Our studies show that a genome-driven approach to discovering bioactive natural products from biofertilizer fungi can lead to unique compounds and biosynthetic knowledge.

Comparative genomic analysis of the 'pseudofungus' Hyphochytrium catenoides

Eukaryotic microbes have three primary mechanisms for obtaining nutrients and energy: phagotrophy, photosynthesis and osmotrophy. Traits associated with the latter two functions arose independently multiple times in the eukaryotes. The Fungi successfully coupled osmotrophy with filamentous growth, and similar traits are also manifested in the Pseudofungi (oomycetes and hyphochytriomycetes). Both the Fungi and the Pseudofungi encompass a diversity of plant and animal parasites. Genome-sequencing efforts have focused on host-associated microbes (mutualistic symbionts or parasites), providing limited comparisons with free-living relatives. Here we report the first draft genome sequence of a hyphochytriomycete ‘pseudofungus’; Hyphochytrium catenoides. Using phylogenomic approaches, we identify genes of recent viral ancestry, with related viral derived genes also present on the genomes of oomycetes, suggesting a complex history of viral coevolution and integration across the Pseudofungi. H. catenoides has a complex life cycle involving diverse filamentous structures and a flagellated zoospore with a single anterior tinselate flagellum. We use genome comparisons, drug sensitivity analysis and high-throughput culture arrays to investigate the ancestry of oomycete/pseudofungal characteristics, demonstrating that many of the genetic features associated with parasitic traits evolved specifically within the oomycete radiation. Comparative genomics also identified differences in the repertoire of genes associated with filamentous growth between the Fungi and the Pseudofungi, including differences in vesicle trafficking systems, cell-wall synthesis pathways and motor protein repertoire, demonstrating that unique cellular systems underpinned the convergent evolution of filamentous osmotrophic growth in these two eukaryotic groups.

Ubiquitin fusion proteins in algae: implications for cell biology and the spread of photosynthesis

Background

The process of gene fusion involves the formation of a single chimeric gene from multiple complete or partial gene sequences. Gene fusion is recognized as an important mechanism by which genes and their protein products can evolve new functions. The presence-absence of gene fusions can also be useful characters for inferring evolutionary relationships between organisms.

Results

Here we show that the nuclear genomes of two unrelated single-celled algae, the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans, possess an unexpected diversity of genes for ubiquitin fusion proteins, including novel arrangements in which ubiquitin occupies amino-terminal, carboxyl-terminal, and internal positions relative to its fusion partners. We explore the evolution of the ubiquitin multigene family in both genomes, and show that both algae possess a gene encoding an ubiquitin-nickel superoxide dismutase fusion protein (Ubiq-NiSOD) that is widely but patchily distributed across the eukaryotic tree of life – almost exclusively in phototrophs.

Conclusion

Our results suggest that ubiquitin fusion proteins are more common than currently appreciated; because of its small size, the ubiquitin coding region can go undetected when gene predictions are carried out in an automated fashion. The punctate distribution of the Ubiq-NiSOD fusion across the eukaryotic tree could serve as a beacon for the spread of plastids from eukaryote to eukaryote by secondary and/or tertiary endosymbiosis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5412-4) contains supplementary material, which is available to authorized users.

A bacteriophage enzyme induces bacterial metabolic perturbation that confers a novel promiscuous function

One key concept in the evolution of new functions is the ability of enzymes to perform promiscuous side-reactions that serve as a source of novelty that may become beneficial under certain conditions. Here, we identify a mechanism where a bacteriophage-encoded enzyme introduces novelty by inducing expression of a promiscuous bacterial enzyme. By screening for bacteriophage DNA that rescued an auxotrophic Escherichia coli mutant carrying a deletion of the ilvA gene, we show that bacteriophage-encoded S-adenosylmethionine (SAM) hydrolases reduce SAM levels. Through this perturbation of bacterial metabolism, expression of the promiscuous bacterial enzyme MetB is increased, which in turn complements the absence of IlvA. These results demonstrate how foreign DNA can increase the metabolic capacity of bacteria, not only by transfer of bona fide new genes, but also by bringing cryptic bacterial functions to light via perturbations of cellular physiology.

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.

What Defines the "Kingdom" Fungi?

The application of environmental DNA techniques and increased genome sequencing of microbial diversity, combined with detailed study of cellular characters, has consistently led to the reexamination of our understanding of the tree of life. This has challenged many of the definitions of taxonomic groups, especially higher taxonomic ranks such as eukaryotic kingdoms. The Fungi is an example of a kingdom which, together with the features that define it and the taxa that are grouped within it, has been in a continual state of flux. In this article we aim to summarize multiple lines of data pertinent to understanding the early evolution and definition of the Fungi. These include ongoing cellular and genomic comparisons that, we will argue, have generally undermined all attempts to identify a synapomorphic trait that defines the Fungi. This article will also summarize ongoing work focusing on taxon discovery, combined with phylogenomic analysis, which has identified novel groups that lie proximate/adjacent to the fungal clade-wherever the boundary that defines the Fungi may be. Our hope is that, by summarizing these data in the form of a discussion, we can illustrate the ongoing efforts to understand what drove the evolutionary diversification of fungi.

Dynamics of genomic innovation in the unicellular ancestry of animals

Which genomic innovations underpinned the origin of multicellular animals is still an open debate. Here, we investigate this question by reconstructing the genome architecture and gene family diversity of ancestral premetazoans, aiming to date the emergence of animal-like traits. Our comparative analysis involves genomes from animals and their closest unicellular relatives (the Holozoa), including four new genomes: three Ichthyosporea and Corallochytrium limacisporum. Here, we show that the earliest animals were shaped by dynamic changes in genome architecture before the emergence of multicellularity: an early burst of gene diversity in the ancestor of Holozoa, enriched in transcription factors and cell adhesion machinery, was followed by multiple and differently-timed episodes of synteny disruption, intron gain and genome expansions. Thus, the foundations of animal genome architecture were laid before the origin of complex multicellularity – highlighting the necessity of a unicellular perspective to understand early animal evolution.

DOI:http://dx.doi.org/10.7554/eLife.26036.001

Hundreds of millions of years ago, some single-celled organisms gained the ability to work together and form multicellular organisms. This transition was a major step in evolution and took place at separate times in several parts of the tree of life, including in animals, plants, fungi and algae.

Animals are some of the most complex organisms on Earth. Their single-celled ancestors were also quite genetically complex themselves and their genomes (the complete set of the organism’s DNA) already contained many genes that now coordinate the activity of the cells in a multicellular organism.

The genome of an animal typically has certain features: it is large, diverse and contains many segments (called introns) that are not genes. By seeing if the single-celled relatives of animals share these traits, it is possible to learn more about when specific genetic features first evolved, and whether they are linked to the origin of animals.

Now, Grau-Bové et al. have studied the genomes of several of the animal kingdom’s closest single-celled relatives using a technique called whole genome sequencing. This revealed that there was a period of rapid genetic change in the single-celled ancestors of animals during which their genes became much more diverse. Another ‘explosion’ of diversity happened after animals had evolved. Furthermore, the overall amount of the genomic content inside cells and the number of introns found in the genome rapidly increased in separate, independent events in both animals and their single-celled ancestors.

Future research is needed to investigate whether other multicellular life forms – such as plants, fungi and algae – originated in the same way as animal life. Understanding how the genetic material of animals evolved also helps us to understand the genetic structures that affect our health. For example, genes that coordinate the behavior of cells (and so are important for multicellular organisms) also play a role in cancer, where cells break free of this regulation to divide uncontrollably.

DOI:http://dx.doi.org/10.7554/eLife.26036.002

Comparative Analysis of Type IV Pilin in Desulfuromonadales

During anaerobic respiration, the bacteria Geobacter sulfurreducens can transfer electrons to extracellular electron accepters through its pilus. G. sulfurreducens pili have been reported to have metallic-like conductivity that is similar to doped organic semiconductors. To study the characteristics and origin of conductive pilin proteins found in the pilus structure, their genetic, structural, and phylogenetic properties were analyzed. The genetic relationships, and conserved structures and sequences that were obtained were used to predict the evolution of the pilins. Homologous genes that encode conductive pilin were found using PilFind and Cluster. Sequence characteristics and protein tertiary structures were analyzed with MAFFT and QUARK, respectively. The origin of conductive pilins was explored by building a phylogenetic tree. Truncation is a characteristic of conductive pilin. The structures of truncated pilins and their accompanying proteins were found to be similar to the N-terminal and C-terminal ends of full-length pilins respectively. The emergence of the truncated pilins can probably be ascribed to the evolutionary pressure of their extracellular electron transporting function. Genes encoding truncated pilins and proteins similar to the C-terminal of full-length pilins, which contain a group of consecutive anti-parallel beta-sheets, are adjacent in bacterial genomes. According to the genetic, structure, and phylogenetic analyses performed in this study, we inferred that the truncated pilins and their accompanying proteins probably evolved from full-length pilins by gene fission through duplication, degeneration, and separation. These findings provide new insights about the molecular mechanisms involved in long-range electron transport along the conductive pili of Geobacter species.

Proteomic profiling of the mitochondrial ribosome identifies Atp25 as a composite mitochondrial precursor protein

Whereas the structure and function of cytosolic ribosomes are well characterized, we only have a limited understanding of the mitochondrial translation apparatus. Using SILAC-based proteomic profiling, we identified 13 proteins that cofractionated with the mitochondrial ribosome, most of which play a role in translation or ribosomal biogenesis. One of these proteins is a homologue of the bacterial ribosome-silencing factor (Rsf). This protein is generated from the composite precursor protein Atp25 upon internal cleavage by the matrix processing peptidase MPP, and in this respect, it differs from all other characterized mitochondrial proteins of baker's yeast. We observed that cytosolic expression of Rsf, but not of noncleaved Atp25 protein, is toxic. Our results suggest that eukaryotic cells face the challenge of avoiding negative interference from the biogenesis of their two distinct translation machineries.

Protein networks identify novel symbiogenetic genes resulting from plastid endosymbiosis

The integration of foreign genetic information is central to the evolution of eukaryotes, as has been demonstrated for the origin of the Calvin cycle and of the heme and carotenoid biosynthesis pathways in algae and plants. For photosynthetic lineages, this coordination involved three genomes of divergent phylogenetic origins (the nucleus, plastid, and mitochondrion). Major hurdles overcome by the ancestor of these lineages were harnessing the oxygen-evolving organelle, optimizing the use of light, and stabilizing the partnership between the plastid endosymbiont and host through retargeting of proteins to the nascent organelle. Here we used protein similarity networks that can disentangle reticulate gene histories to explore how these significant challenges were met. We discovered a previously hidden component of algal and plant nuclear genomes that originated from the plastid endosymbiont: symbiogenetic genes (S genes). These composite proteins, exclusive to photosynthetic eukaryotes, encode a cyanobacterium-derived domain fused to one of cyanobacterial or another prokaryotic origin and have emerged multiple, independent times during evolution. Transcriptome data demonstrate the existence and expression of S genes across a wide swath of algae and plants, and functional data indicate their involvement in tolerance to oxidative stress, phototropism, and adaptation to nitrogen limitation. Our research demonstrates the "recycling" of genetic information by photosynthetic eukaryotes to generate novel composite genes, many of which function in plastid maintenance.

Fungal traits that drive ecosystem dynamics on land

Fungi contribute extensively to a wide range of ecosystem processes, including decomposition of organic carbon, deposition of recalcitrant carbon, and transformations of nitrogen and phosphorus. In this review, we discuss the current knowledge about physiological and morphological traits of fungi that directly influence these processes, and we describe the functional genes that encode these traits. In addition, we synthesize information from 157 whole fungal genomes in order to determine relationships among selected functional genes within fungal taxa. Ecosystem-related traits varied most at relatively coarse taxonomic levels. For example, we found that the maximum amount of variance for traits associated with carbon mineralization, nitrogen and phosphorus cycling, and stress tolerance could be explained at the levels of order to phylum. Moreover, suites of traits tended to co-occur within taxa. Specifically, the genetic capacities for traits that improve stress tolerance-β-glucan synthesis, trehalose production, and cold-induced RNA helicases-were positively related to one another, and they were more evident in yeasts. Traits that regulate the decomposition of complex organic matter-lignin peroxidases, cellobiohydrolases, and crystalline cellulases-were also positively related, but they were more strongly associated with free-living filamentous fungi. Altogether, these relationships provide evidence for two functional groups: stress tolerators, which may contribute to soil carbon accumulation via the production of recalcitrant compounds; and decomposers, which may reduce soil carbon stocks. It is possible that ecosystem functions, such as soil carbon storage, may be mediated by shifts in the fungal community between stress tolerators and decomposers in response to environmental changes, such as drought and warming.

Bacterial proteins pinpoint a single eukaryotic root

The large phylogenetic distance separating eukaryotic genes and their archaeal orthologs has prevented identification of the position of the eukaryotic root in phylogenomic studies. Recently, an innovative approach has been proposed to circumvent this issue: the use as phylogenetic markers of proteins that have been transferred from bacterial donor sources to eukaryotes, after their emergence from Archaea. Using this approach, two recent independent studies have built phylogenomic datasets based on bacterial sequences, leading to different predictions of the eukaryotic root. Taking advantage of additional genome sequences from the jakobid Andalucia godoyi and the two known malawimonad species (Malawimonas jakobiformis and Malawimonas californiana), we reanalyzed these two phylogenomic datasets. We show that both datasets pinpoint the same phylogenetic position of the eukaryotic root that is between "Unikonta" and "Bikonta," with malawimonad and collodictyonid lineages on the Unikonta side of the root. Our results firmly indicate that (i) the supergroup Excavata is not monophyletic and (ii) the last common ancestor of eukaryotes was a biflagellate organism. Based on our results, we propose to rename the two major eukaryotic groups Unikonta and Bikonta as Opimoda and Diphoda, respectively.

Evolution and classification of myosins, a paneukaryotic whole-genome approach

Myosins are key components of the eukaryotic cytoskeleton, providing motility for a broad diversity of cargoes. Therefore, understanding the origin and evolutionary history of myosin classes is crucial to address the evolution of eukaryote cell biology. Here, we revise the classification of myosins using an updated taxon sampling that includes newly or recently sequenced genomes and transcriptomes from key taxa. We performed a survey of eukaryotic genomes and phylogenetic analyses of the myosin gene family, reconstructing the myosin toolkit at different key nodes in the eukaryotic tree of life. We also identified the phylogenetic distribution of myosin diversity in terms of number of genes, associated protein domains and number of classes in each taxa. Our analyses show that new classes (i.e., paralogs) and domain architectures were continuously generated throughout eukaryote evolution, with a significant expansion of myosin abundance and domain architectural diversity at the stem of Holozoa, predating the origin of animal multicellularity. Indeed, single-celled holozoans have the most complex myosin complement among eukaryotes, with paralogs of most myosins previously considered animal specific. We recover a dynamic evolutionary history, with several lineage-specific expansions (e.g., the myosin III-like gene family diversification in choanoflagellates), convergence in protein domain architectures (e.g., fungal and animal chitin synthase myosins), and important secondary losses. Overall, our evolutionary scheme demonstrates that the ancestral eukaryote likely had a complex myosin repertoire that included six genes with different protein domain architectures. Finally, we provide an integrative and robust classification, useful for future genomic and functional studies on this crucial eukaryotic gene family.

Complex patterns of gene fission in the eukaryotic folate biosynthesis pathway

Shared derived genomic characters can be useful for polarizing phylogenetic relationships, for example, gene fusions have been used to identify deep-branching relationships in the eukaryotes. Here, we report the evolutionary analysis of a three-gene fusion of folB, folK, and folP, which encode enzymes that catalyze consecutive steps in de novo folate biosynthesis. The folK-folP fusion was found across the eukaryotes and a sparse collection of prokaryotes. This suggests an ancient derivation with a number of gene losses in the eukaryotes potentially as a consequence of adaptation to heterotrophic lifestyles. In contrast, the folB-folK-folP gene is specific to a mosaic collection of Amorphea taxa (a group encompassing: Amoebozoa, Apusomonadida, Breviatea, and Opisthokonta). Next, we investigated the stability of this character. We identified numerous gene losses and a total of nine gene fission events, either by break up of an open reading frame (four events identified) or loss of a component domain (five events identified). This indicates that this three gene fusion is highly labile. These data are consistent with a growing body of data indicating gene fission events occur at high relative rates. Accounting for these sources of homoplasy, our data suggest that the folB-folK-folP gene fusion was present in the last common ancestor of Amoebozoa and Opisthokonta but absent in the Metazoa including the human genome. Comparative genomic data of these genes provides an important resource for designing therapeutic strategies targeting the de novo folate biosynthesis pathway of a variety of eukaryotic pathogens such as Acanthamoeba castellanii.

Genome sequence and virulence variation-related transcriptome profiles of Curvularia lunata, an important maize pathogenic fungus

BACKGROUND

Curvularia lunata is an important maize foliar fungal pathogen that distributes widely in maize growing area in China. Genome sequencing of the pathogen will provide important information for globally understanding its virulence mechanism.

RESULTS

We report the genome sequences of a highly virulent C. lunata strain. Phylogenomic analysis indicates that C. lunata was evolved from Bipolaris maydis (Cochliobolus heterostrophus). The highly virulent strain has a high potential to evolve into other pathogenic stains based on analyses on transposases and repeat-induced point mutations. C. lunata has a smaller proportion of secreted proteins as well as B. maydis than entomopathogenic fungi. C. lunata and B. maydis have a similar proportion of protein-encoding genes highly homologous to experimentally proven pathogenic genes from pathogen-host interaction database. However, relative to B. maydis, C. lunata possesses not only many expanded protein families including MFS transporters, G-protein coupled receptors, protein kinases and proteases for transport, signal transduction or degradation, but also many contracted families including cytochrome P450, lipases, glycoside hydrolases and polyketide synthases for detoxification, hydrolysis or secondary metabolites biosynthesis, which are expected to be crucial for the fungal survival in varied stress environments. Comparative transcriptome analysis between a lowly virulent C. lunata strain and its virulence-increased variant induced by resistant host selection reveals that the virulence increase of the pathogen is related to pathways of toxin and melanin biosynthesis in stress environments, and that the two pathways probably have some overlaps.

CONCLUSIONS

The data will facilitate a full revelation of pathogenic mechanism and a better understanding of virulence differentiation of C. lunata.

A rhodopsin-guanylyl cyclase gene fusion functions in visual perception in a fungus

Sensing light is the fundamental property of visual systems, with vision in animals being based almost exclusively on opsin photopigments [1]. Rhodopsin also acts as a photoreceptor linked to phototaxis in green algae [2, 3] and has been implicated by chemical means as a light sensor in the flagellated swimming zoospores of the fungus Allomyces reticulatus [4]; however, the signaling mechanism in these fungi remains unknown. Here we use a combination of genome sequencing and molecular inhibition experiments with light-sensing phenotype studies to examine the signaling pathway involved in visual perception in the closely related fungus Blastocladiella emersonii. Our data show that in these fungi, light perception is accomplished by the function of a novel gene fusion (BeGC1) of a type I (microbial) rhodopsin domain and guanylyl cyclase catalytic domain. Photobleaching of rhodopsin function prevents accumulation of cGMP levels and phototaxis of fungal zoospores exposed to green light, whereas inhibition of guanylyl cyclase activity negatively affects fungal phototaxis. Immunofluorescence microscopy localizes the BeGC1 protein to the external surface of the zoospore eyespot positioned close to the base of the swimming flagellum [4, 5], demonstrating this is a photoreceptive organelle composed of lipid droplets. Taken together, these data indicate that Blastocladiomycota fungi have a cGMP signaling pathway involved in phototaxis similar to the vertebrate vision-signaling cascade but composed of protein domain components arranged as a novel gene fusion architecture and of distant evolutionary ancestry to type II rhodopsins of animals.

An alternative root for the eukaryote tree of life

The root of the eukaryote tree of life defines some of the most fundamental relationships among species. It is also critical for defining the last eukaryote common ancestor (LECA), the shared heritage of all extant species. The unikont-bikont root has been the reigning paradigm for eukaryotes for more than 10 years but is becoming increasingly controversial. We developed a carefully vetted data set, consisting of 37 nuclear-encoded proteins of close bacterial ancestry (euBacs) and their closest bacterial relatives, augmented by deep sequencing of the Acrasis kona (Heterolobosea, Discoba) transcriptome. Phylogenetic analysis of these data produces a highly robust, fully resolved global phylogeny of eukaryotes. The tree sorts all examined eukaryotes into three megagroups and identifies the Discoba, and potentially its parent taxon Excavata, as the sister group to the bulk of known eukaryote diversity, the proposed Neozoa (Amorphea + Stramenopila+Alveolata+Rhizaria+Plantae [SARP]). All major alternative hypotheses are rejected with as little as ∼50% of the data, and this resolution is unaffected by the presence of fast-evolving alignment positions or distant outgroup sequences. This "neozoan-excavate" root revises hypotheses of early eukaryote evolution and highlights the importance of the poorly studied Discoba for understanding the evolution of eukaryotic diversity and basic cellular processes.

The Schizosaccharomyces pombe fusion gene hal3 encodes three distinct activities

Saccharomyces cerevisiae Hal3 and Vhs3 are moonlighting proteins, forming an atypical heterotrimeric decarboxylase (PPCDC) required for CoA biosynthesis, and regulating cation homeostasis by inhibition of the Ppz1 phosphatase. The Schizosaccharomyces pombe ORF SPAC15E1.04 (renamed as Sp hal3) encodes a protein whose amino-terminal half is similar to Sc Hal3 whereas its carboxyl-terminal half is related to thymidylate synthase (TS). We show that Sp Hal3 and/or its N-terminal domain retain the ability to bind to and modestly inhibit in vitro S. cerevisiae Ppz1 as well as its S. pombe homolog Pzh1, and also exhibit PPCDC activity in vitro and provide PPCDC function in vivo, indicating that Sp Hal3 is a monogenic PPCDC in fission yeast. Whereas the Sp Hal3 N-terminal domain partially mimics Sc Hal3 functions, the entire protein and its carboxyl-terminal domain rescue the S. cerevisiae cdc21 mutant, thus proving TS function. Additionally, we show that the 70 kDa Sp Hal3 protein is not proteolytically processed under diverse forms of stress and that, as predicted, Sp hal3 is an essential gene. Therefore, Sp hal3 represents a fusion event that joined three different functional activities in the same gene. The possible advantage derived from this surprising combination of essential proteins is discussed.

Toward genome-enabled mycology

Genome-enabled mycology is a rapidly expanding field that is characterized by the pervasive use of genome-scale data and associated computational tools in all aspects of fungal biology. Genome-enabled mycology is integrative and often requires teams of researchers with diverse skills in organismal mycology, bioinformatics and molecular biology. This issue of Mycologia presents the first complete fungal genomes in the history of the journal, reflecting the ongoing transformation of mycology into a genome-enabled science. Here, we consider the prospects for genome-enabled mycology and the technical and social challenges that will need to be overcome to grow the database of complete fungal genomes and enable all fungal biologists to make use of the new data.

Extensive intra-kingdom horizontal gene transfer converging on a fungal fructose transporter gene

Comparative genomics revealed in the last decade a scenario of rampant horizontal gene transfer (HGT) among prokaryotes, but for fungi a clearly dominant pattern of vertical inheritance still stands, punctuated however by an increasing number of exceptions. In the present work, we studied the phylogenetic distribution and pattern of inheritance of a fungal gene encoding a fructose transporter (FSY1) with unique substrate selectivity. 109 FSY1 homologues were identified in two sub-phyla of the Ascomycota, in a survey that included 241 available fungal genomes. At least 10 independent inter-species instances of horizontal gene transfer (HGT) involving FSY1 were identified, supported by strong phylogenetic evidence and synteny analyses. The acquisition of FSY1 through HGT was sometimes suggestive of xenolog gene displacement, but several cases of pseudoparalogy were also uncovered. Moreover, evidence was found for successive HGT events, possibly including those responsible for transmission of the gene among yeast lineages. These occurrences do not seem to be driven by functional diversification of the Fsy1 proteins because Fsy1 homologues from widely distant lineages, including at least one acquired by HGT, appear to have similar biochemical properties. In summary, retracing the evolutionary path of the FSY1 gene brought to light an unparalleled number of independent HGT events involving a single fungal gene. We propose that the turbulent evolutionary history of the gene may be linked to the unique biochemical properties of the encoded transporter, whose predictable effect on fitness may be highly variable. In general, our results support the most recent views suggesting that inter-species HGT may have contributed much more substantially to shape fungal genomes than heretofore assumed.