Concatenated Analysis Sheds Light on Early Metazoan Evolution and Fuels a Modern “Urmetazoon” Hypothesis

For more than a century, the origin of metazoan animals has been debated. One aspect of this debate has been centered on what the hypothetical “urmetazoon” bauplan might have been. The morphologically most simply organized metazoan animal, the placozoan Trichoplax adhaerens, resembles an intriguing model for one of several “urmetazoon” hypotheses: the placula hypothesis. Clear support for a basal position of Placozoa would aid in resolving several key issues of metazoan-specific inventions (including, for example, head–foot axis, symmetry, and coelom) and would determine a root for unraveling their evolution. Unfortunately, the phylogenetic relationships at the base of Metazoa have been controversial because of conflicting phylogenetic scenarios generated while addressing the question. Here, we analyze the sum of morphological evidence, the secondary structure of mitochondrial ribosomal genes, and molecular sequence data from mitochondrial and nuclear genes that amass over 9,400 phylogenetically informative characters from 24 to 73 taxa. Together with mitochondrial DNA genome structure and sequence analyses and Hox-like gene expression patterns, these data (1) provide evidence that Placozoa are basal relative to all other diploblast phyla and (2) spark a modernized “urmetazoon” hypothesis.     

Concatenated Analysis Sheds Light on Early Metazoan Evolution and Fuels a Modern “Urmetazoon” Hypothesis

For more than a century, the origin of metazoan animals has been debated. One aspect of this debate has been centered on what the hypothetical “urmetazoon” bauplan might have been. The morphologically most simply organized metazoan animal, the placozoan Trichoplax adhaerens, resembles an intriguing model for one of several “urmetazoon” hypotheses: the placula hypothesis. Clear support for a basal position of Placozoa would aid in resolving several key issues of metazoan-specific inventions (including, for example, head–foot axis, symmetry, and coelom) and would determine a root for unraveling their evolution. Unfortunately, the phylogenetic relationships at the base of Metazoa have been controversial because of conflicting phylogenetic scenarios generated while addressing the question. Here, we analyze the sum of morphological evidence, the secondary structure of mitochondrial ribosomal genes, and molecular sequence data from mitochondrial and nuclear genes that amass over 9,400 phylogenetically informative characters from 24 to 73 taxa. Together with mitochondrial DNA genome structure and sequence analyses and Hox-like gene expression patterns, these data (1) provide evidence that Placozoa are basal relative to all other diploblast phyla and (2) spark a modernized “urmetazoon” hypothesis.     

Modular Skeletal Evolution in Sticklebacks Is Controlled by Additive and Clustered Quantitative Trait Loci

Understanding the genetic architecture of evolutionary change remains a long-standing goal in biology. In vertebrates, skeletal evolution has contributed greatly to adaptation in body form and function in response to changing ecological variables like diet and predation. Here we use genome-wide linkage mapping in threespine stickleback fish to investigate the genetic architecture of evolved changes in many armor and trophic traits. We identify >100 quantitative trait loci (QTL) controlling the pattern of serially repeating skeletal elements, including gill rakers, teeth, branchial bones, jaws, median fin spines, and vertebrae. We use this large collection of QTL to address long-standing questions about the anatomical specificity, genetic dominance, and genomic clustering of loci controlling skeletal differences in evolving populations. We find that most QTL (76%) that influence serially repeating skeletal elements have anatomically regional effects. In addition, most QTL (71%) have at least partially additive effects, regardless of whether the QTL controls evolved loss or gain of skeletal elements. Finally, many QTL with high LOD scores cluster on chromosomes 4, 20, and 21. These results identify a modular system that can control highly specific aspects of skeletal form. Because of the general additivity and genomic clustering of major QTL, concerted changes in both protective armor and trophic traits may occur when sticklebacks inherit either marine or freshwater alleles at linked or possible “supergene” regions of the stickleback genome. Further study of these regions will help identify the molecular basis of both modular and coordinated changes in the vertebrate skeleton.     

Symbiophagy and biomineralization in the “living fossil” Astrosclera willeyana

Representatives of all major metazoan lineages form biominerals. The molecular mechanisms that underlie this widespread and evolutionarily ancient ability are gradually being revealed for some lineages. However, until a wider range of metazoan biomineralization strategies are understood, the true diversity, and therefore the evolutionary origins of this process, will remain unknown. We have previously shown that the coralline demosponge, Astrosclera willeyana, in some way employs its endobiotic bacterial community to form its highly calcified skeleton. Here, using in situ hybridization and immunohistochemistry, we show that an ortholog of ATG8 (most likely a GABARAPL2/GATE-16 ortholog) is expressed in cells that construct the individual skeletal elements of the sponge. In TEM sections sponge cells can be observed to contain extensive populations of bacteria, and frequently possesses double-membrane structures which we interpret to be autophagosomes. In combination with our previous work, these findings support the hypothesis that the host sponge actively degrades a proportion of its bacterial community using an autophagy pathway, and uses the prokaryotic organic remains as a framework upon which calcification of the sponge skeleton is initiated.     

Contribution of sponge genes to unravel the genome of the hypothetical ancestor of Metazoa (Urmetazoa)

Recently the term Urmetazoa, as the hypothetical metazoan ancestor, was introduced to highlight the finding that all metazoan phyla including the Porifera (sponges) are derived from one common ancestor. Sponges as the evolutionarily oldest, still extant phylum, are provided with a complex network of structural and functional molecules. Analyses of sponge genomes from Demospongiae (Suberites domuncula and Geodia cydonium), Calcarea (Sycon raphanus) and Hexactinellida (Aphrocallistes vastus) have contributed also to the reconstruction of the evolutionary position of Metazoa with respect to Fungi. Furthermore, these analyses have provided evidence that the characteristic evolutionary novelties of Metazoa, such as the extracellular matrix molecules, the cell surface receptors, the nervous signal transduction molecules as well as the immune molecule existing in Porifera, share high sequence and in some aspects also functional similarities to related polypeptides found in other metazoan phyla. During the transition to Metazoa new domains occurred; as one example, the formation of the death domain from the ankyrin is outlined. In parallel, domanial proteins have been formed, such as the receptor tyrosine kinases. The metazoan essentials have been defined by analyzing and comparing the sponge sequences with the related sequences from the metazoans Homo sapiens, Caenorhabditis elegans and Drosophila melanogaster, the fungus Saccharomyces cerevisiae and the plant Arabidopsis thaliana. The data revealed that those sponge molecules grouped to cell adhesion cell recognition proteins are predominantly found in Protostomia and Deuterostomia while they are missing in Fungi and Viridiplantae. Moreover, evidence is presented allowing the conclusion that the sponge molecules are more closely related to the corresponding molecules from H. sapiens than to those of C. elegans or D. melanogaster. Especially surprising was the finding that the Demospongiae are provided with elements of adaptive immunity.     

Proteomics of Protein Secretion by Bacillus subtilis: Separating the “Secrets” of the Secretome

Secretory proteins perform a variety of important “remote-control” functions for bacterial survival in the environment. The availability of complete genome sequences has allowed us to make predictions about the composition of bacterial machinery for protein secretion as well as the extracellular complement of bacterial proteomes. Recently, the power of proteomics was successfully employed to evaluate genome-based models of these so-called secretomes. Progress in this field is well illustrated by the proteomic analysis of protein secretion by the gram-positive bacterium Bacillus subtilis, for which ~90 extracellular proteins were identified. Analysis of these proteins disclosed various “secrets of the secretome,” such as the residence of cytoplasmic and predicted cell envelope proteins in the extracellular proteome. This showed that genome-based predictions reflect only ~50% of the actual composition of the extracellular proteome of B. subtilis. Importantly, proteomics allowed the first verification of the impact of individual secretion machinery components on the total flow of proteins from the cytoplasm to the extracellular environment. In conclusion, proteomics has yielded a variety of novel leads for the analysis of protein traffic in B. subtilis and other gram-positive bacteria. Ultimately, such leads will serve to increase our understanding of virulence factor biogenesis in gram-positive pathogens, which is likely to be of high medical relevance.     

Multi-scale mineralogical characterization of the hypercalcified sponge Petrobiona massiliana (Calcarea, Calcaronea)

The massive basal skeleton of a few remnant living hypercalcified sponges rediscovered since the 1960s are valuable representatives of ancient calcium carbonate biomineralization mechanisms in basal Metazoa. A multi-scale mineralogical characterization of the easily accessible Mediterranean living hypercalcified sponge belonging to Calcarea, Petrobiona massiliana (Vacelet and Lévi, 1958), was conducted. Oriented observations in light and electron microscopy of mature and growing areas of the Mg-calcite basal skeleton were combined in order to describe all structural levels from the submicronic to the macroscopic scale. The smallest units produced are ca. 50–100 nm grains that are in a mushy amorphous state before their crystallization. Selected area electron diffraction (SAED) further demonstrated that submicronic grains are assembled into crystallographically coherent clusters or fibers, the latter are even laterally associated into single-crystal bundles. A model of crystallization propagation through amorphous submicronic granular units is proposed to explain the formation of coherent micron-scale structural units. Finally, XRD and EELS analyses highlighted, respectively, inter-individual variation of skeletal Mg contents and heterogeneous spatial distribution of Ca ions in skeletal fibers. All mineralogical features presented here cannot be explained by classical inorganic crystallization principles in super-saturated solutions, but rather underlined a highly biologically regulated formation of the basal skeleton. This study extending recent observations on corals, mollusk and echinoderms confirms that occurrence of submicronic granular units and a possible transient amorphous precursor phase in calcium carbonate skeletons is a common biomineralization strategy already selected by basal metazoans.     

Interplay between Chaperones and Protein Disorder Promotes the Evolution of Protein Networks

Evolution is driven by mutations, which lead to new protein functions but come at a cost to protein stability. Non-conservative substitutions are of interest in this regard because they may most profoundly affect both function and stability. Accordingly, organisms must balance the benefit of accepting advantageous substitutions with the possible cost of deleterious effects on protein folding and stability. We here examine factors that systematically promote non-conservative mutations at the proteome level. Intrinsically disordered regions in proteins play pivotal roles in protein interactions, but many questions regarding their evolution remain unanswered. Similarly, whether and how molecular chaperones, which have been shown to buffer destabilizing mutations in individual proteins, generally provide robustness during proteome evolution remains unclear. To this end, we introduce an evolutionary parameter λ that directly estimates the rate of non-conservative substitutions. Our analysis of λ in Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens sequences reveals how co- and post-translationally acting chaperones differentially promote non-conservative substitutions in their substrates, likely through buffering of their destabilizing effects. We further find that λ serves well to quantify the evolution of intrinsically disordered proteins even though the unstructured, thus generally variable regions in proteins are often flanked by very conserved sequences. Crucially, we show that both intrinsically disordered proteins and highly re-wired proteins in protein interaction networks, which have evolved new interactions and functions, exhibit a higher λ at the expense of enhanced chaperone assistance. Our findings thus highlight an intricate interplay of molecular chaperones and protein disorder in the evolvability of protein networks. Our results illuminate the role of chaperones in enabling protein evolution, and underline the importance of the cellular context and integrated approaches for understanding proteome evolution. We feel that the development of λ may be a valuable addition to the toolbox applied to understand the molecular basis of evolution.     

Distinct Types of Disorder in the Human Proteome: Functional Implications for Alternative Splicing

Intrinsically disordered regions have been associated with various cellular processes and are implicated in several human diseases, but their exact roles remain unclear. We previously defined two classes of conserved disordered regions in budding yeast, referred to as “flexible” and “constrained” conserved disorder. In flexible disorder, the property of disorder has been positionally conserved during evolution, whereas in constrained disorder, both the amino acid sequence and the property of disorder have been conserved. Here, we show that flexible and constrained disorder are widespread in the human proteome, and are particularly common in proteins with regulatory functions. Both classes of disordered sequences are highly enriched in regions of proteins that undergo tissue-specific (TS) alternative splicing (AS), but not in regions of proteins that undergo general (i.e., not tissue-regulated) AS. Flexible disorder is more highly enriched in TS alternative exons, whereas constrained disorder is more highly enriched in exons that flank TS alternative exons. These latter regions are also significantly more enriched in potential phosphosites and other short linear motifs associated with cell signaling. We further show that cancer driver mutations are significantly enriched in regions of proteins associated with TS and general AS. Collectively, our results point to distinct roles for TS alternative exons and flanking exons in the dynamic regulation of protein interaction networks in response to signaling activity, and they further suggest that alternatively spliced regions of proteins are often functionally altered by mutations responsible for cancer.     

Performance of Single and Concatenated Sets of Mitochondrial Genes at Inferring Metazoan Relationships Relative to Full Mitogenome Data

Mitochondrial (mt) genes are some of the most popular and widely-utilized genetic loci in phylogenetic studies of metazoan taxa. However, their linked nature has raised questions on whether using the entire mitogenome for phylogenetics is overkill (at best) or pseudoreplication (at worst). Moreover, no studies have addressed the comparative phylogenetic utility of mitochondrial genes across individual lineages within the entire Metazoa. To comment on the phylogenetic utility of individual mt genes as well as concatenated subsets of genes, we analyzed mitogenomic data from 1865 metazoan taxa in 372 separate lineages spanning genera to subphyla. Specifically, phylogenies inferred from these datasets were statistically compared to ones generated from all 13 mt protein-coding (PC) genes (i.e., the “supergene” set) to determine which single genes performed “best” at, and the minimum number of genes required to, recover the “supergene” topology. Surprisingly, the popular marker COX1 performed poorest, while ND5, ND4, and ND2 were most likely to reproduce the “supergene” topology. Averaged across all lineages, the longest ∼2 mt PC genes were sufficient to recreate the “supergene” topology, although this average increased to ∼5 genes for datasets with 40 or more taxa. Furthermore, concatenation of the three “best” performing mt PC genes outperformed that of the three longest mt PC genes (i.e, ND5, COX1, and ND4). Taken together, while not all mt PC genes are equally interchangeable in phylogenetic studies of the metazoans, some subset can serve as a proxy for the 13 mt PC genes. However, the exact number and identity of these genes is specific to the lineage in question and cannot be applied indiscriminately across the Metazoa.     

The Pre-Endosymbiont Hypothesis: A New Perspective on the Origin and Evolution of Mitochondria

Mitochondrial DNA (mtDNA) is unquestionably the remnant of an α-proteobacterial genome, yet only ∼10%–20% of mitochondrial proteins are demonstrably α-proteobacterial in origin (the “α-proteobacterial component,” or APC). The evolutionary ancestry of the non-α-proteobacterial component (NPC) is obscure and not adequately accounted for in current models of mitochondrial origin. I propose that in the host cell that accommodated an α-proteobacterial endosymbiont, much of the NPC was already present, in the form of a membrane-bound metabolic organelle (the premitochondrion) that compartmentalized many of the non-energy-generating functions of the contemporary mitochondrion. I suggest that this organelle also possessed a protein import system and various ion and small-molecule transporters. In such a scenario, an α-proteobacterial endosymbiont could have been converted relatively directly and rapidly into an energy-generating organelle that incorporated the extant metabolic functions of the premitochondrion. This model (the “pre-endosymbiont hypothesis”) effectively represents a synthesis of previous, contending mitochondrial origin hypotheses, with the bulk of the mitochondrial proteome (much of the NPC) having an endogenous origin and the minority component (the APC) having a xenogenous origin.Considering the central role played in all eukaryotic cells by mitochondria or mitochondrion-related organelles (MROs, such as hydrogenosomes and mitosomes) (Hjort et al. 2010; Shiflett and Johnson 2010; Müller et al. 2012), the question of the origin and subsequent evolution of the mitochondrion has long captivated and challenged biologists. In a recent article in this series (Gray 2012), I discussed in detail several aspects of mitochondrial evolution, focusing particularly on how well the accumulating molecular data can be accommodated in current models of mitochondrial origin. In this context, the origin and evolution of the mitochondrial proteome, as opposed to the origin and evolution of the mitochondrial genome, were examined from the perspective of comparative mitochondrial proteomics. Somewhat disconcertingly, as more data have become available, we find ourselves considerably less certain about key aspects of how mitochondria originated than we were (or thought we were) several decades ago.Here, I summarize key points discussed in more detail in the previous article before presenting a novel perspective on how the mitochondrion might have originated. The new model proposed here, which represents a synthesis of both endogenous (“origin from within”) and xenogenous (“origin from outside”) modes, is advanced in an attempt to account for the inability of a purely endosymbiotic model, whose strongest support has come from studies of the mitochondrial genome, to adequately accommodate data on the mitochondrial proteome.     

Intrinsically Disordered Regions May Lower the Hydration Free Energy in Proteins: A Case Study of Nudix Hydrolase in the Bacterium Deinococcus radiodurans

The proteome of the radiation- and desiccation-resistant bacterium D. radiodurans features a group of proteins that contain significant intrinsically disordered regions that are not present in non-extremophile homologues. Interestingly, this group includes a number of housekeeping and repair proteins such as DNA polymerase III, nudix hydrolase and rotamase. Here, we focus on a member of the nudix hydrolase family from D. radiodurans possessing low-complexity N- and C-terminal tails, which exhibit sequence signatures of intrinsic disorder and have unknown function. The enzyme catalyzes the hydrolysis of oxidatively damaged and mutagenic nucleotides, and it is thought to play an important role in D. radiodurans during the recovery phase after exposure to ionizing radiation or desiccation. We use molecular dynamics simulations to study the dynamics of the protein, and study its hydration free energy using the GB/SA formalism. We show that the presence of disordered tails significantly decreases the hydration free energy of the whole protein. We hypothesize that the tails increase the chances of the protein to be located in the remaining water patches in the desiccated cell, where it is protected from the desiccation effects and can function normally. We extrapolate this to other intrinsically disordered regions in proteins, and propose a novel function for them: intrinsically disordered regions increase the “surface-properties” of the folded domains they are attached to, making them on the whole more hydrophilic and potentially influencing, in this way, their localization and cellular activity.     

The Oxytricha trifallax Macronuclear Genome: A Complex Eukaryotic Genome with 16,000 Tiny Chromosomes

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor “silent” germline micronuclear genome by a process of “unscrambling” and fragmentation. The tiny macronuclear “nanochromosomes” typically encode single, protein-coding genes (a small portion, 10%, encode 2–8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.     

Meiofauna and the origins of the Metazoa*

Possible alternative habitats and life-styles of the original metazoan are considered. It is argued from the dominance of the benthic habitat in present-day groups that the original metazoan habitat was benthic rather than planktonic. Similarly, plesiomorphic metazoan taxa tend to be holobenthic rather than pelago-benthic. It is therefore probable that the early Metazoa were holobenthic. The concept of plesiohabitats and apohabitats in the evolution of taxa is presented. This leads to the proposition that the early metazoans were interstitial bionts of fine sand. Finally, the controversy concerning the aerobic or anaerobic origin of the Metazoa is considered. It is shown that competition theory predicts that plesiomorphic taxa are likely to remain in plesiohabitats. Diagrams showing the possible evolution of major taxa in relation to available habitats are presented. It is concluded that the earliest Metazoa could have evolved in anaerobic marine sand and that the early Plathelminthomorpha and Aschelminthes did so.     

An evolutionary fast‐track to biocalcification

The ability to construct mineralized shells, spicules, spines and skeletons is thought to be a key factor that fuelled the expansion of multicellular animal life during the early Cambrian. The genes and molecular mechanisms that control the process of biomineralization in disparate phyla are gradually being revealed, and it is broadly recognized that an insoluble matrix of proteins, carbohydrates and other organic molecules are required for the initiation, regulation and inhibition of crystal growth. Here, we show that Astrosclera willeyana, a living representative of the now largely extinct stromatoporid sponges (a polyphyletic grade of poriferan bauplan), has apparently bypassed the requirement to evolve many of these mineral‐regulating matrix proteins by using the degraded remains of bacteria to seed CaCO₃ crystal growth. Because stromatoporid sponges formed extensive reefs during the Paelozoic and Mesozoic eras (fulfilling the role that stony corals play in modern coral reefs), and fossil evidence suggests that the same process of bacterial skeleton formation occurred in these stromatoporid ancestors, we infer that some ancient reef ecosystems might have been founded on this microbial–metazoan relationship.     

A Comprehensive Phylogenetic Analysis of the Serpin Superfamily

Serine protease inhibitors (serpins) are found in all kingdoms of life and play essential roles in multiple physiological processes. Owing to the diversity of the superfamily, phylogenetic analysis is challenging and prokaryotic serpins have been speculated to have been acquired from Metazoa through horizontal gene transfer due to their unexpectedly high homology. Here, we have leveraged a structural alignment of diverse serpins to generate a comprehensive 6,000-sequence phylogeny that encompasses serpins from all kingdoms of life. We show that in addition to a central “hub” of highly conserved serpins, there has been extensive diversification of the superfamily into many novel functional clades. Our analysis indicates that the hub proteins are ancient and are similar because of convergent evolution, rather than the alternative hypothesis of horizontal gene transfer. This work clarifies longstanding questions in the evolution of serpins and provides new directions for research in the field of serpin biology.