Phylogenetic Resolution of Deep Eukaryotic and Fungal Relationships Using Highly Conserved Low-Copy Nuclear Genes

Species evolution: cryptic species and phenotypic noise with a particular focus on fungal systematics

The evolution of a species can be understood in the context of two major concepts-the cryptic species concept and the phenotypic noise concept. The former represents morphologically indistinguishable but genetically distinct evolutionary lineages, while the latter represents the phenotypic variations of an isogenic population. Although the concept of cryptic species currently represents a general topic, its effect on other aspects of biology, such as biodiversity, ecology, evolutionary biology, and taxonomy, is still unclear. In particular, cryptic species cause complications and prevent the development of a clear taxonomy. The phenotypic noise concept or phenotypic plasticity generally refers to the various expressions of phenotypes in different environments. Hence, the cryptic species concept refers to genetic variations, while the phenotypic noises concept is about non-genetic variations. Although both concepts are opposites, they each contribute significantly to the evolutionary process of an organism. Despite the extensive research studies and publications discussing those two concepts in separate accounts, a concise account that combines and compares both concepts are generally lacking. Nevertheless, these are essential to understand the evolutionary process clearly. This review addresses the available literature on this topic, intending to provide a general and overall discussion on both the cryptic species concept and the phenotypic noise concept and their effect on evolution, ecology, biodiversity, and taxonomy with a special focus on fungal systematics. hence, several fungal case studies representing the two concepts are presented, compared, and discussed for a better understanding.

Identification of hidden associations among eukaryotic genes through statistical analysis of coevolutionary transitions

Coevolution at the gene level, as reflected by correlated events of gene loss or gain, can be revealed by phylogenetic profile analysis. The optimal method and metric for comparing phylogenetic profiles, especially in eukaryotic genomes, are not yet established. Here, we describe a procedure suitable for large-scale analysis, which can reveal coevolution based on the assessment of the statistical significance of correlated presence/absence transitions between gene pairs. This metric can identify coevolution in profiles with low overall similarities and is not affected by similarities lacking coevolutionary information. We applied the procedure to a large collection of 60,912 orthologous gene groups (orthogroups) in 1,264 eukaryotic genomes extracted from OrthoDB. We found significant cotransition scores for 7,825 orthogroups associated in 2,401 coevolving modules linking known and unknown genes in protein complexes and biological pathways. To demonstrate the ability of the method to predict hidden gene associations, we validated through experiments the involvement of vertebrate malate synthase-like genes in the conversion of (S)-ureidoglycolate into glyoxylate and urea, the last step of purine catabolism. This identification explains the presence of glyoxylate cycle genes in metazoa and suggests an anaplerotic role of purine degradation in early eukaryotes.

Metagenomic Analysis Using Phylogenetic Placement-A Review of the First Decade

Phylogenetic placement refers to a family of tools and methods to analyze, visualize, and interpret the tsunami of metagenomic sequencing data generated by high-throughput sequencing. Compared to alternative (e. g., similarity-based) methods, it puts metabarcoding sequences into a phylogenetic context using a set of known reference sequences and taking evolutionary history into account. Thereby, one can increase the accuracy of metagenomic surveys and eliminate the requirement for having exact or close matches with existing sequence databases. Phylogenetic placement constitutes a valuable analysis tool per se, but also entails a plethora of downstream tools to interpret its results. A common use case is to analyze species communities obtained from metagenomic sequencing, for example via taxonomic assignment, diversity quantification, sample comparison, and identification of correlations with environmental variables. In this review, we provide an overview over the methods developed during the first 10 years. In particular, the goals of this review are 1) to motivate the usage of phylogenetic placement and illustrate some of its use cases, 2) to outline the full workflow, from raw sequences to publishable figures, including best practices, 3) to introduce the most common tools and methods and their capabilities, 4) to point out common placement pitfalls and misconceptions, 5) to showcase typical placement-based analyses, and how they can help to analyze, visualize, and interpret phylogenetic placement data.

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.

Repeated Evolution of Inactive Pseudonucleases in a Fungal Branch of the Dis3/RNase II Family of Nucleases

The RNase II family of 3'-5' exoribonucleases is present in all domains of life, and eukaryotic family members Dis3 and Dis3L2 play essential roles in RNA degradation. Ascomycete yeasts contain both Dis3 and inactive RNase II-like "pseudonucleases." The latter function as RNA-binding proteins that affect cell growth, cytokinesis, and fungal pathogenicity. However, the evolutionary origins of these pseudonucleases are unknown: What sequence of events led to their novel function, and when did these events occur? Here, we show how RNase II pseudonuclease homologs, including Saccharomyces cerevisiae Ssd1, are descended from active Dis3L2 enzymes. During fungal evolution, active site mutations in Dis3L2 homologs have arisen at least four times, in some cases following gene duplication. In contrast, N-terminal cold-shock domains and regulatory features are conserved across diverse dikarya and mucoromycota, suggesting that the nonnuclease function requires these regions. In the basidiomycete pathogenic yeast Cryptococcus neoformans, the single Ssd1/Dis3L2 homolog is required for cytokinesis from polyploid "titan" growth stages. This phenotype of C. neoformans Ssd1/Dis3L2 deletion is consistent with those of inactive fungal pseudonucleases, yet the protein retains an active site sequence signature. We propose that a nuclease-independent function for Dis3L2 arose in an ancestral hyphae-forming fungus. This second function has been conserved across hundreds of millions of years, whereas the RNase activity was lost repeatedly in independent lineages.

Tracking gene expression, metabolic profiles, and biochemical analysis in the halotolerant basidiomycetous yeast Rhodotorula mucilaginosa EXF-1630 during benzo[a]pyrene and phenanthrene biodegradation under hypersaline conditions

Polyaromatic phenanthrene (Phe) and benzo[a]pyrene (BaP) are highly toxic, mutagenic, and carcinogenic contaminants widely dispersed in nature, including saline environments. Polyextremotolerant Rhodotorula mucilaginosa EXF-1630, isolated from Arctic sea ice, was grown on a huge concentration range -10 to 500 ppm- of Phe and BaP as sole carbon sources at hypersaline conditions (1 M NaCl). Selected polycyclic aromatic hydrocarbons (PAHs) supported growth as well as glucose, even at high PAH concentrations. Initially, up to 40% of Phe and BaP were adsorbed, followed by biodegradation, resulting in 80% removal in 10 days. While extracellular laccase, peroxidase, and un-specific peroxygenase activities were not detected, NADPH-cytochrome c reductase activity peaked at 4 days. The successful removal of PAHs and the absence of toxic metabolites were confirmed by toxicological tests on moss Physcomitrium patens, bacterium Aliivibrio fischeri, human erythrocytes, and pulmonary epithelial cells (A549). Metabolic profiles were determined at the midpoint of the biodegradation exponential phase, with added Phe and BaP (100 ppm) and 1 M NaCl. Different hydroxylated products were found in the culture medium, while the conjugative metabolite 1-phenanthryl-β-D-glucopyranose was detected in the medium and in the cells. Transcriptome analysis resulted in 870 upregulated and 2,288 downregulated transcripts on PAHs, in comparison to glucose. Genomic mining of 61 available yeast genomes showed a widespread distribution of 31 xenobiotic degradation pathways in different yeast lineages. Two distributions with similar metabolic capacities included black yeasts and mainly members of the Sporidiobolaceae family (including EXF-1630), respectively. This is the first work describing a metabolic profile and transcriptomic analysis of PAH degradation by yeast.

SUMOylation in fungi: A potential target for intervention

SUMOylation is a post-translational, reversible modification process which occurs in eukaryotes. Small Ubiquitin like MOdifier or (SUMO) proteins are a family of small proteins that are covalently attached to and detached from other proteins to modify the target protein function. In pathogenic fungi, SUMO has been identified and preliminary studies indicate its importance either for survival and/or for virulence. In this review we provide an overview of the current state of knowledge of SUMOylation in fungi and the effects on pathogenesis. Subsequently we identify the orthologs of the SUMOylation pathway components across fungi. We also show the level of conservation of the proteins involved and identify the similarities/differences in the orthologs across fungi and the human and plant hosts to identify potential targets of intervention.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

Transcriptomic analysis of polyaromatic hydrocarbon degradation by the halophilic fungus Aspergillus sydowii at hypersaline conditions

Polycyclic aromatic hydrocarbons (PAHs) are among the most persistent xenobiotic compounds, with high toxicity effects. Mycoremediation with halophilic Aspergillus sydowii was used for their removal from a hypersaline medium (1 M NaCl). A. sydowii metabolized PAHs as sole carbon sources, resulting in the removal of up to 90% for both PAHs [benzo [a] pyrene (BaP) and phenanthrene (Phe)] after 10 days. Elimination of Phe and BaP was almost exclusively due to biotransformation and not adsorption by dead mycelium and did not correlate with the activity of lignin modifying enzymes (LME). Transcriptomes of A. sydowii grown on PAHs, or on glucose as control, both at hypersaline conditions, revealed 170 upregulated and 76 downregulated genes. Upregulated genes were related to starvation, cell wall remodelling, degradation and metabolism of xenobiotics, DNA/RNA metabolism, energy generation, signalling and general stress responses. Changes of LME expression levels were not detected, while the chloroperoxidase gene, possibly related to detoxification processes in fungi, was strongly upregulated. We propose that two parallel metabolic pathways (mitochondrial and cytosolic) are involved in degradation and detoxification of PAHs in A. sydowii resulting in intracellular oxidation of PAHs. To the best of our knowledge, this is the most comprehensive transcriptomic analysis on fungal degradation of PAHs.

MetaEuk-sensitive, high-throughput gene discovery, and annotation for large-scale eukaryotic metagenomics

BACKGROUND

Metagenomics is revolutionizing the study of microorganisms and their involvement in biological, biomedical, and geochemical processes, allowing us to investigate by direct sequencing a tremendous diversity of organisms without the need for prior cultivation. Unicellular eukaryotes play essential roles in most microbial communities as chief predators, decomposers, phototrophs, bacterial hosts, symbionts, and parasites to plants and animals. Investigating their roles is therefore of great interest to ecology, biotechnology, human health, and evolution. However, the generally lower sequencing coverage, their more complex gene and genome architectures, and a lack of eukaryote-specific experimental and computational procedures have kept them on the sidelines of metagenomics.

RESULTS

MetaEuk is a toolkit for high-throughput, reference-based discovery, and annotation of protein-coding genes in eukaryotic metagenomic contigs. It performs fast searches with 6-frame-translated fragments covering all possible exons and optimally combines matches into multi-exon proteins. We used a benchmark of seven diverse, annotated genomes to show that MetaEuk is highly sensitive even under conditions of low sequence similarity to the reference database. To demonstrate MetaEuk's power to discover novel eukaryotic proteins in large-scale metagenomic data, we assembled contigs from 912 samples of the Tara Oceans project. MetaEuk predicted >12,000,000 protein-coding genes in 8 days on ten 16-core servers. Most of the discovered proteins are highly diverged from known proteins and originate from very sparsely sampled eukaryotic supergroups.

CONCLUSION

The open-source (GPLv3) MetaEuk software (https://github.com/soedinglab/metaeuk) enables large-scale eukaryotic metagenomics through reference-based, sensitive taxonomic and functional annotation. Video abstract.

Fungal evolution: diversity, taxonomy and phylogeny of the Fungi

The fungal kingdom comprises a hyperdiverse clade of heterotrophic eukaryotes characterized by the presence of a chitinous cell wall, the loss of phagotrophic capabilities and cell organizations that range from completely unicellular monopolar organisms to highly complex syncitial filaments that may form macroscopic structures. Fungi emerged as a 'Third Kingdom', embracing organisms that were outside the classical dichotomy of animals versus vegetals. The taxonomy of this group has a turbulent history that is only now starting to be settled with the advent of genomics and phylogenomics. We here review the current status of the phylogeny and taxonomy of fungi, providing an overview of the main defined groups. Based on current knowledge, nine phylum-level clades can be defined: Opisthosporidia, Chytridiomycota, Neocallimastigomycota, Blastocladiomycota, Zoopagomycota, Mucoromycota, Glomeromycota, Basidiomycota and Ascomycota. For each group, we discuss their main traits and their diversity, focusing on the evolutionary relationships among the main fungal clades. We also explore the diversity and phylogeny of several groups of uncertain affinities and the main phylogenetic and taxonomical controversies and hypotheses in the field.

Inferring Ancient Relationships with Genomic Data: A Commentary on Current Practices

Contemporary phylogeneticists enjoy an embarrassment of riches, not only in the volumes of data now available, but also in the diversity of bioinformatic tools for handling these data. Here, I discuss a subset of these tools I consider well-suited to the task of inferring ancient relationships with coding sequence data in particular, encompassing data generation, orthology assignment, alignment and gene tree inference, supermatrix construction, and analysis under the best-fitting models applicable to large-scale datasets. Throughout, I compare and critique methods, considering both their theoretical principles and the details of their implementation, and offering practical tips on usage where appropriate. I also entertain different motivations for analyzing what are almost always originally DNA sequence data as codons, amino acids, and higher-order recodings. Although presented in a linear order, I see value in using the diversity of tools available to us to assess the sensitivity of clades of biological interest to different gene and taxon sets and analytical modes, which can be an indication of the presence of systematic error, of which a few forms remain poorly controlled by even the best available inference methods.

Evolutionary, structural and functional analysis of the caleosin/peroxygenase gene family in the Fungi

BACKGROUND

Caleosin/peroxygenases, CLO/PXG, (designated PF05042 in Pfam) are a group of genes/proteins with anomalous distributions in eukaryotic taxa. We have previously characterised CLO/PXGs in the Viridiplantae. The aim of this study was to investigate the evolution and functions of the CLO/PXGs in the Fungi and other non-plant clades and to elucidate the overall origin of this gene family.

RESULTS

CLO/PXG-like genes are distributed across the full range of fungal groups from the basal clades, Cryptomycota and Microsporidia, to the largest and most complex Dikarya species. However, the genes were only present in 243 out of 844 analysed fungal genomes. CLO/PXG-like genes have been retained in many pathogenic or parasitic fungi that have undergone considerable genomic and structural simplification, indicating that they have important functions in these species. Structural and functional analyses demonstrate that CLO/PXGs are multifunctional proteins closely related to similar proteins found in all major taxa of the Chlorophyte Division of the Viridiplantae. Transcriptome and physiological data show that fungal CLO/PXG-like genes have complex patterns of developmental and tissue-specific expression and are upregulated in response to a range of biotic and abiotic stresses as well as participating in key metabolic and developmental processes such as lipid metabolism, signalling, reproduction and pathogenesis. Biochemical data also reveal that the Aspergillus flavus CLO/PXG has specific functions in sporulation and aflatoxin production as well as playing roles in lipid droplet function.

CONCLUSIONS

In contrast to plants, CLO/PXGs only occur in about 30% of sequenced fungal genomes but are present in all major taxa. Fungal CLO/PXGs have similar but not identical roles to those in plants, including stress-related oxylipin signalling, lipid metabolism, reproduction and pathogenesis. While the presence of CLO/PXG orthologs in all plant genomes sequenced to date would suggest that they have core housekeeping functions in plants, the selective loss of CLO/PXGs in many fungal genomes suggests more restricted functions in fungi as accessory genes useful in particular environments or niches. We suggest an ancient origin of CLO/PXG-like genes in the 'last eukaryotic common ancestor' (LECA) and their subsequent loss in ancestors of the Metazoa, after the latter had diverged from the ancestral fungal lineage.

Optimal Rates for Phylogenetic Inference and Experimental Design in the Era of Genome-Scale Data Sets

With the rise of genome-scale data sets, there has been a call for increased data scrutiny and careful selection of loci that are appropriate to use in an attempt to resolve a phylogenetic problem. Such loci should maximize phylogenetic information content while minimizing the risk of homoplasy. Theory posits the existence of characters that evolve at an optimum rate, and efforts to determine optimal rates of inference have been a cornerstone of phylogenetic experimental design for over two decades. However, both theoretical and empirical investigations of optimal rates have varied dramatically in their conclusions: spanning no relationship to a tight relationship between the rate of change and phylogenetic utility. Herein, we synthesize these apparently contradictory views, demonstrating both empirical and theoretical conditions under which each is correct. We find that optimal rates of characters-not genes-are generally robust to most experimental design decisions. Moreover, consideration of site rate heterogeneity within a given locus is critical to accurate predictions of utility. Factors such as taxon sampling or the targeted number of characters providing support for a topology are additionally critical to the predictions of phylogenetic utility based on the rate of character change. Further, optimality of rates and predictions of phylogenetic utility are not equivalent, demonstrating the need for further development of comprehensive theory of phylogenetic experimental design. [Divergence time; GC bias; homoplasy; incongruence; information content; internode length; optimal rates; phylogenetic informativeness; phylogenetic theory; phylogenetic utility; phylogenomics; signal and noise; subtending branch length; state space; taxon and character sampling.].

Comparative genomics of nuclear envelope proteins

BACKGROUND

The nuclear envelope (NE) that encapsulates the nuclear genome is a double lipid bilayer with several integral and peripherally associated proteins. It is a characteristic feature of the eukaryotes and acts as a hub for a number of important nuclear events including transcription, repair, and regulated gene expression. The proteins associated with the nuclear envelope mediate the NE functions and maintain its structural integrity, which is crucial for survival. In spite of the importance of this structure, knowledge of the protein composition of the nuclear envelope and their function, are limited to very few organisms belonging to Opisthokonta and Archaeplastida supergroups. The NE composition is largely unknown in organisms outside these two supergroups.

RESULTS

In this study, we have taken a comparative sequence analysis approach to identify the NE proteome that is present across all five eukaryotic supergroups. We identified 22 proteins involved in various nuclear functions to be part of the core NE proteome. The presence of these proteins across eukaryotes, suggests that they are traceable to the Last Eukaryotic Common Ancestor (LECA). Additionally, we also identified the NE proteins that have evolved in a lineage specific manner and those that have been preserved only in a subset of organisms.

CONCLUSIONS

Our study identifies the conserved features of the nuclear envelope across eukaryotes and provides insights into the potential composition and the functionalities that were constituents of the LECA NE.

Multigene phylogeny and cell evolution of chromist infrakingdom Rhizaria: contrasting cell organisation of sister phyla Cercozoa and Retaria

Infrakingdom Rhizaria is one of four major subgroups with distinct cell body plans that comprise eukaryotic kingdom Chromista. Unlike other chromists, Rhizaria are mostly heterotrophic flagellates, amoebae or amoeboflagellates, commonly with reticulose (net-like) or filose (thread-like) feeding pseudopodia; uniquely for eukaryotes, cilia have proximal ciliary transition-zone hub-lattices. They comprise predominantly flagellate phylum Cercozoa and reticulopodial phylum Retaria, whose exact phylogenetic relationship has been uncertain. Given even less clear relationships amongst cercozoan classes, we sequenced partial transcriptomes of seven Cercozoa representing five classes and endomyxan retarian Filoreta marina to establish 187-gene multiprotein phylogenies. Ectoreta (retarian infraphyla Foraminifera, Radiozoa) branch within classical Cercozoa as sister to reticulose Endomyxa. This supports recent transfer of subphylum Endomyxa from Cercozoa to Retaria alongside subphylum Ectoreta which embraces classical retarians where capsules or tests subdivide cells into organelle-containing endoplasm and anastomosing pseudopodial net-like ectoplasm. Cercozoa are more homogeneously filose, often with filose pseudopodia and/or posterior ciliary gliding motility: zooflagellate Helkesimastix and amoeboid Guttulinopsis form a strongly supported clade, order Helkesida. Cercomonads are polyphyletic (Cercomonadida sister to glissomonads; Paracercomonadida deeper). Thecofilosea are a clade, whereas Imbricatea may not be; Sarcomonadea may be paraphyletic. Helkesea and Metromonadea are successively deeper outgroups within cercozoan subphylum Monadofilosa; subphylum Reticulofilosa (paraphyletic on site-heterogeneous trees) branches earliest, Granofilosea before Chlorarachnea. Our multiprotein trees confirm that Rhizaria are sisters of infrakingdom Halvaria (Alveolata, Heterokonta) within chromist subkingdom Harosa (= SAR); they further support holophyly of chromist subkingdom Hacrobia, and are consistent with holophyly of Chromista as sister of kingdom Plantae. Site-heterogeneous rDNA trees group Kraken with environmental DNA clade 'eSarcomonad', not Paracercomonadida. Ectoretan fossil dates evidence ultrarapid episodic stem sequence evolution. We discuss early rhizarian cell evolution and multigene tree coevolutionary patterns, gene-paralogue evidence for chromist monophyly, and integrate this with fossil evidence for the age of Rhizaria and eukaryote cells, and revise rhizarian classification.

Phylogenetic taxon definitions for Fungi, Dikarya, Ascomycota and Basidiomycota

Phylogenetic taxon definitions (PTDs) are explicit, phylogeny-based statements that specify clades. PTDs are central to the system of rank-free classification that is governed by the PhyloCode, but they can also be used to clarify the meanings of ranked names. We present PTDs for four major groups: Fungi, Dikarya, Ascomycota, and Basidiomycota.

Origin and evolution of fungal HECT ubiquitin ligases

Ubiquitin ligases (E3s) are basic components of the eukaryotic ubiquitination system. In this work, the emergence and diversification of fungal HECT ubiquitin ligases is described. Phylogenetic and structural data indicate that six HECT subfamilies (RSP5, TOM1, UFD4, HUL4, HUL4A and HUL5) existed in the common ancestor of all fungi. These six subfamilies have evolved very conservatively, with only occasional losses and duplications in particular fungal lineages. However, an early, drastic reduction in the number of HECT genes occurred in microsporidians, in parallel to the reduction of their genomes. A significant correlation between the total number of genes and the number of HECT-encoding genes present in fungi has been observed. However, transitions from unicellularity to multicellularity or vice versa apparently had no effect on the evolution of this family. Likely orthologs or co-orthologs of all fungal HECT genes have been detected in animals. Four genes are deduced to be present in the common ancestor of fungi, animals and plants. Protein-protein interactions detected in both the yeast Saccharomyces cerevisiae and humans suggest that some ancient functions of HECT proteins have been conserved since the animals/fungi split.

Kingdom Chromista and its eight phyla: a new synthesis emphasising periplastid protein targeting, cytoskeletal and periplastid evolution, and ancient divergences

In 1981 I established kingdom Chromista, distinguished from Plantae because of its more complex chloroplast-associated membrane topology and rigid tubular multipartite ciliary hairs. Plantae originated by converting a cyanobacterium to chloroplasts with Toc/Tic translocons; most evolved cell walls early, thereby losing phagotrophy. Chromists originated by enslaving a phagocytosed red alga, surrounding plastids by two extra membranes, placing them within the endomembrane system, necessitating novel protein import machineries. Early chromists retained phagotrophy, remaining naked and repeatedly reverted to heterotrophy by losing chloroplasts. Therefore, Chromista include secondary phagoheterotrophs (notably ciliates, many dinoflagellates, Opalozoa, Rhizaria, heliozoans) or walled osmotrophs (Pseudofungi, Labyrinthulea), formerly considered protozoa or fungi respectively, plus endoparasites (e.g. Sporozoa) and all chromophyte algae (other dinoflagellates, chromeroids, ochrophytes, haptophytes, cryptophytes). I discuss their origin, evolutionary diversification, and reasons for making chromists one kingdom despite highly divergent cytoskeletons and trophic modes, including improved explanations for periplastid/chloroplast protein targeting, derlin evolution, and ciliary/cytoskeletal diversification. I conjecture that transit-peptide-receptor-mediated 'endocytosis' from periplastid membranes generates periplastid vesicles that fuse with the arguably derlin-translocon-containing periplastid reticulum (putative red algal trans-Golgi network homologue; present in all chromophytes except dinoflagellates). I explain chromist origin from ancestral corticates and neokaryotes, reappraising tertiary symbiogenesis; a chromist cytoskeletal synapomorphy, a bypassing microtubule band dextral to both centrioles, favoured multiple axopodial origins. I revise chromist higher classification by transferring rhizarian subphylum Endomyxa from Cercozoa to Retaria; establishing retarian subphylum Ectoreta for Foraminifera plus Radiozoa, apicomonad subclasses, new dinozoan classes Myzodinea (grouping Colpovora gen. n., Psammosa), Endodinea, Sulcodinea, and subclass Karlodinia; and ranking heterokont Gyrista as phylum not superphylum.

Maximizing Power in Phylogenetics and Phylogenomics: A Perspective Illuminated by Fungal Big Data

Since its original inception over 150 years ago by Darwin, we have made tremendous progress toward the reconstruction of the Tree of Life. In particular, the transition from analyzing datasets comprised of small numbers of loci to those comprised of hundreds of loci, if not entire genomes, has aided in resolving some of the most vexing of evolutionary problems while giving us a new perspective on biodiversity. Correspondingly, phylogenetic trees have taken a central role in fields that span ecology, conservation, and medicine. However, the rise of big data has also presented phylogenomicists with a new set of challenges to experimental design, quantitative analyses, and computation. The sequencing of a number of very first genomes presented significant challenges to phylogenetic inference, leading fungal phylogenomicists to begin addressing pitfalls and postulating solutions to the issues that arise from genome-scale analyses relevant to any lineage across the Tree of Life. Here we highlight insights from fungal phylogenomics for topics including systematics and species delimitation, ecological and phenotypic diversification, and biogeography while providing an overview of progress made on the reconstruction of the fungal Tree of Life. Finally, we provide a review of considerations to phylogenomic experimental design for robust tree inference. We hope that this special issue of Advances in Genetics not only excites the continued progress of fungal evolutionary biology but also motivates the interdisciplinary development of new theory and methods designed to maximize the power of genomic scale data in phylogenetic analyses.

The Conservation of the Germline Multipotency Program, from Sponges to Vertebrates: A Stepping Stone to Understanding the Somatic and Germline Origins

The germline definition in metazoans was first based on few bilaterian models. As a result, gene function interpretations were often based on phenotypes observed in those models and led to the definition of a set of genes, considered as specific of the germline, named the “germline core”. However, some of these genes were shown to also be involved in somatic stem cells, thus leading to the notion of germline multipotency program (GMP). Because Porifera and Ctenophora are currently the best candidates to be the sister-group to all other animals, the comparative analysis of gene contents and functions between these phyla, Cnidaria and Bilateria is expected to provide clues on early animal evolution and on the links between somatic and germ lineages. Our present bioinformatic analyses at the metazoan scale show that a set of 18 GMP genes was already present in the last common ancestor of metazoans and indicate more precisely the evolution of some of them in the animal lineage. The expression patterns and levels of 11 of these genes in the homoscleromorph sponge Oscarella lobularis show that they are expressed throughout their life cycle, in pluri/multipotent progenitors, during gametogenesis, embryogenesis and during wound healing. This new study in a nonbilaterian species reinforces the hypothesis of an ancestral multipotency program.

Cytochrome P450 diversity in the tree of life

Sequencing in all areas of the tree of life has produced >300,000 cytochrome P450 (CYP) sequences that have been mined and collected. Nomenclature has been assigned to >41,000 CYP sequences and the majority of the remainder has been sorted by BLAST searches into clans, families and subfamilies in preparation for naming. The P450 sequence space is being systematically explored and filled in. Well-studied groups like vertebrates are covered in greater depth while new insights are being added into uncharted territories like horseshoe crab (Limulus polyphemus), tardigrades (Hypsibius dujardini), velvet worm (Euperipatoides_rowelli), and basal land plants like hornworts, liverworts and mosses. CYPs from the fungi, one of the most diverse groups, are being explored and organized as nearly 800 fungal species are now sequenced. The CYP clan structure in fungi is emerging with 805 CYP families sorting into 32 CYP clans. >3000 bacterial sequences are named, mostly from terrestrial or freshwater sources. Of 18,379 bacterial sequences downloaded from the CYPED database, all are >43% identical to named CYPs. Therefore, they fit in the 602 named P450 prokaryotic families. Diversity in this group is becoming saturated, however 25% of 3305 seawater bacterial P450s did not match known P450 families, indicating marine bacterial CYPs are not as well sampled as land/freshwater based bacterial CYPs. Future sequencing plans of the Genome 10K project, i5k and GIGA (Global Invertebrate Genomics Alliance) are expected to produce more than one million cytochrome P450 sequences by 2020. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.