Update of MmtDB: a Metazoa mitochondrial DNA variants database.

The present paper describes the improvements in MmtDB, a specialised database designed to collect Metazoa mitochondrial DNA variants. Priority in the data collection has been given to Metazoa for which a large amount of variants is available, e.g., for humans. Starting from the sequences available in the Nucleotide Sequence Databases, the redundant sequences have been removed and new sequences from other sources have been added. Value-added information is associated to each variant sequence, e.g., analysed region, experimental method, tissue and cell lines, population data, sex, age, family code and information about the variation events (nucleotide position, involved gene, restriction site gain or loss). Cross-references are introduced to the EMBL Data Library, as well as an internal cross-referencing among MmtDB entries according to tissual, heteroplasmic, familiar and aplotypical correlation. Furthermore MmtDB has a new section, AMmtDB: Aligned Metazoan mitochondrial biosequences. MmtDB can be accessed through the World Wide Web at URL http://WWW.ba.cnr.it/[symbol: see text]areamt08/MmtDBWWW.htm     

Evolutionary relationships of Metazoa within the eukaryotes based on molecular data from Porifera

Recent molecular data provide strong support for the view that all metazoan phyla, including Porifera, are of monophyletic origin. The relationship of Metazoa, including the Porifera, to Plantae, Fungi and unicellular eukaryotes has only rarely been studied by using cDNAs coding for proteins. Sequence data from rDNA suggested a relationship of Porifera to unicellular eukaryotes (choanoflagellates). However, ultrastructural studies of choanocytes did not support these findings. In the present study, we compared amino acid sequences that are found in a variety of metazoans (including sponges) with those of Plantae, Fungi and unicellular eukaryotes, to obtain an answer to this question. We used the four sequences from 70 kDa heat-shock proteins, the serine-threonine kinase domain found in protein kinases, beta-tubulin and calmodulin. The latter two sequences were deduced from cDNAs, isolated from the sponge Geodia cydonium for the phylogenetic analyses presented. These revealed that the sponge molecules were grouped into the same branch as the Metazoa, which is statistically (significantly) separated from those branches that comprise the sequences from Fungi, Plantae and unicellular eukaryotes. From our molecular data it seems evident that the unicellular eukaryotes existed at an earlier stage of evolution, and the Plantae and especially the Fungi and the Metazoa only appeared later.     

Receptor tyrosine kinase, an autapomorphic character of metazoa: identification in marine sponges

In the present review we summarize sequence data obtained from cloning of sponge receptor tyrosine kinases [RTK]. The cDNA sequences were mainly obtained from the marine sponge Geodia cydonium. RTKs (i) with immunoglobulin [Ig]-like domains in the extracellular region, (ii) of the type of insulin-like receptors, as well as (iii) RTKs with one extracellular speract domain, have been identified. The analyses revealed that the RTK genes are constructed in blocks [domains], suggesting a blockwise evolution. The phylogenetic relationships of the sequences obtained revealed that all sponge sequences fall into one branch of the evolutionary tree, while related sequences from higher Metazoa, human, mouse and rat, including also invertebrate sequences, together form a second branch. It is concluded that the RTK molecules have evolved in sponges prior to the "Cambrian Explosion" and have contributed to the rapid appearance of the higher metazoan phyla and that sponges are, as a taxon, also monophyletic. Due to the fact that protein tyrosine kinases in general and RTKs in particular have only been identified in Metazoa, they are, as a group qualified, to be considered as an autapomorphic character of all metazoan phyla.     

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.     

Variability of the mitochondrial genome in mammals at the inter-species/intra-species boundary

Genomic variations represent the molecular basis of the biodiversity of living organisms on which selection operates to generate evolution. In eukaryotes, genomic variability can be experienced in both nuclear and organellar, i.e. mitochondrial and plastid (where present), genomes, which can follow completely different evolution pathways, as revealed by comparative genomics analyses. In Metazoa, for which a substantial number of complete genome sequences are available (nuclear, but mainly mitochondrial), we are just starting to grasp the selective pressures operating on some basic features of the genome as a whole. In this brief review, we discuss the variability of the mitochondrial metazoan genome, with particular reference to mitochondrial DNA in mammals. In light of the recent assumption that a small segment of mitochondrial DNA may be used, particularly in Metazoa, as a species marker, some data on mitochondrial gene variability at the inter-species/intra-species boundary are reported. Intra-species variability has been evaluated in four mammalian species, Homo sapiens, Bos taurus, Sus scrofa and Canis familiaris, whereas the relationship between intra- and inter-species variability has been investigated in Bos taurus and Bos indicus.     

ATRIUM: Testing Untyped SNPs in Case-Control Association Studies with Related Individuals

In genome-wide association studies, only a subset of all genomic variants are typed by current, high-throughput, SNP-genotyping platforms. However, many of the untyped variants can be well predicted from typed variants, with linkage disequilibrium (LD) information among typed and untyped variants available from an external reference panel such as HapMap. Incorporation of such external information can allow one to perform tests of association between untyped variants and phenotype, thereby making more efficient use of the available genotype data. When related individuals are included in case-control samples, the dependence among their genotypes must be properly addressed for valid association testing. In the context of testing untyped variants, an additional analytical challenge is that the dependence, across related individuals, of the partial information on untyped-SNP genotypes must also be assessed and incorporated into the analysis for valid inference. We address this challenge with ATRIUM, a method for case-control association testing with untyped SNPs, based on genome screen data in samples in which some individuals are related. ATRIUM uses LD information from an external reference panel to specify a one-degree-of-freedom test of association with an untyped SNP. It properly accounts for dependence in the partial information on untyped-SNP genotypes across related individuals. We demonstrate that ATRIUM is robust in that it maintains the nominal type I error rate even when the external reference panel is not well matched to the case-control sample. We apply the method to detect association between type 2 diabetes and variants on chromosome 10 in the Framingham SHARe data.     

Plant snoRNA database

The Plant snoRNA database (http://www.scri.sari.ac.uk/plant_snoRNA/) provides information on small nucleolar RNAs from Arabidopsis and eighteen other plant species. Information includes sequences, expression data, methylation and pseudouridylation target modification sites, initial gene organization (polycistronic, single gene and intronic) and the number of gene variants. The Arabidopsis information is divided into box C/D and box H/ACA snoRNAs, and within each of these groups, by target sites in rRNA, snRNA or unknown. Alignments of orthologous genes and gene variants from different plant species are available for many snoRNA genes. Plant snoRNA genes have been given a standard nomenclature, designed wherever possible, to provide a consistent identity with yeast and human orthologues.     

Conflicting phylogenetic signals at the base of the metazoan tree

A phylogenetic framework is essential for under-standing the origin and evolution of metazoan development. Despite a number of recent molecular studies and a rich fossil record of sponges and cnidarians, the evolutionary relationships of the early branching metazoan groups to each other and to a putative outgroup, the choanoflagellates, remain uncertain. This situation may be the result of the limited amount of phylogenetic information found in single genes and the small number of relevant taxa surveyed. To alleviate the effect of these analytical factors in the phylogenetic recons-truction of early branching metazoan lineages, we cloned multiple protein-coding genes from two choanoflagellates and diverse sponges, cnidarians, and a ctenophore. Comparisons of sequences for alpha-tubulin, beta-tubulin, elongation factor 2, HSP90, and HSP70 robustly support the hypothesis that choanoflagellates are closely affiliated with animals. However, analyses of single and concatenated amino acid sequences fail to resolve the relationships either between early branching metazoan groups or between Metazoa and choano-flagellates. We demonstrate that variable rates of evolution among lineages, sensitivity of the analyses to taxon selection, and conflicts in the phylogenetic signal contained in different amino acid sequences obscure the phylogenetic associations among the early branching Metazoa. These factors raise concerns about the ability to resolve the phylogenetic history of animals with molecular sequences. A consensus view of animal evolution may require investigations of genome-scale characters.     

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.     

The closest unicellular relatives of animals

Molecular phylogenies support a common ancestry between animals (Metazoa) and Fungi, but the evolutionary descent of the Metazoa from single-celled eukaryotes (protists) and the nature and taxonomic affiliation of these ancestral protists remain elusive. We addressed this question by sequencing complete mitochondrial genomes from taxonomically diverse protists to generate a large body of molecular data for phylogenetic analyses. Trees inferred from multiple concatenated mitochondrial protein sequences demonstrate that animals are specifically affiliated with two morphologically dissimilar unicellular protist taxa: Monosiga brevicollis (Choanoflagellata), a flagellate, and Amoebidium parasiticum (Ichthyosporea), a fungus-like organism. Statistical evaluation of competing evolutionary hypotheses confirms beyond a doubt that Choanoflagellata and multicellular animals share a close sister group relationship, originally proposed more than a century ago on morphological grounds. For the first time, our trees convincingly resolve the currently controversial phylogenetic position of the Ichthyosporea, which the trees place basal to Choanoflagellata and Metazoa but after the divergence of Fungi. Considering these results, we propose the new taxonomic group Holozoa, comprising Ichthyosporea, Choanoflagellata, and Metazoa. Our findings provide insight into the nature of the animal ancestor and have broad implications for our understanding of the evolutionary transition from unicellular protists to multicellular animals.     

PhyreRisk: A Dynamic Web Application to Bridge Genomics,Proteomics and 3D Structural Data to Guide Interpretation of Human Genetic Variants

PhyreRisk is an open-access, publicly accessible web application for interactively bridging genomic, proteomic and structural data facilitating the mapping of human variants onto protein structures. A major advance over other tools for sequence-structure variant mapping is that PhyreRisk provides information on 20,214 human canonical proteins and an additional 22,271 alternative protein sequences (isoforms). Specifically, PhyreRisk provides structural coverage (partial or complete) for 70% (14,035 of 20,214 canonical proteins) of the human proteome, by storing 18,874 experimental structures and 84,818 pre-built models of canonical proteins and their isoforms generated using our in house Phyre2. PhyreRisk reports 55,732 experimentally, multi-validated protein interactions from IntAct and 24,260 experimental structures of protein complexes.Another major feature of PhyreRisk is that, rather than presenting a limited set of precomputed variant-structure mapping of known genetic variants, it allows the user to explore novel variants using, as input, genomic coordinates formats (Ensembl, VCF, reference SNP ID and HGVS notations) and Human Build GRCh37 and GRCh38. PhyreRisk also supports mapping variants using amino acid coordinates and searching for genes or proteins of interest.PhyreRisk is designed to empower researchers to translate genetic data into protein structural information, thereby providing a more comprehensive appreciation of the functional impact of variants. PhyreRisk is freely available at http://phyrerisk.bc.ic.ac.uk     

Structure of the choanocyte kinetid in the sponge <Emphasis Type="Italic">Haliclona</Emphasis> sp. (Demospongiae: Haplosclerida) and its implication for taxonomy and phylogeny of Demospongiae

Sponges are a group of primitive animals that are of particular interest in terms of evolutionary theory. Molecular and phylogenic studies have revealed the common origin of sponges and Metazoa and the relationship of Metazoa to Choanoflagellata. Comparative studies of the morphology have allowed reconstruction of the structure of ancestors common to Metazoa and Choanoflagellata, as well as for Porifera and other Metazoa. For this purpose, a comparative analysis of phylogenically informative structures, such as the kinetid, is needed. This article presents data on the structure of the choanocyte kinetid in Haliclona sp., which was investigated for the first time. The reconstruction of the kinetid based on ultrathin serial sections shows that it lacks an accessory centriole; it has fibrillar roots despite the basal position of the nucleus; two-to-three triplets of the kinetosome are longer than the others; bundles of lateral microtubules arise from fibrillar roots and reach the cytoplasmic membrane in the region of intercellular contact. The studies of the kinetid in the few representatives of Demospongiae investigated showed the compatibility of the data on the ultrastructure and molecular phylogeny.     

Phylogenetic analysis of the Hsp70 sequences reveals the monophyly of Metazoa and specific phylogenetic relationships between animals and fungi

To understand the early evolution of the Metazoa, it is necessary to determine the correct phylogenetic status of diploblastic animals. Despite cladistic studies of morphological characters and recent molecular phylogenetic studies, it remains uncertain whether diploblasts are monophyletic or paraphyletic, and how the phyla of diploblasts are phylogenetically related. The heat shock protein 70 (Hsp70) sequences, because of their ubiquity and high degree of conservation, could provide a useful model for phylogenetic analysis. We have sequenced almost the entire nucleic acid sequence of cytoplasmic Hsp70 from eight diploblastic species. Our data support the monophyly of diploblastic animals. However, the phylogenetic relationships of the diploblast groups were not significantly resolved. Our phylogenetic trees also support the monophyly of Metazoa with high bootstrap values, indicating that animals form an extremely robust clade.      

Annotating functional effects of non-coding variants in neuropsychiatric cell types by deep transfer learning

Genomewide association studies (GWAS) have identified a large number of loci associated with neuropsychiatric traits, however, understanding the molecular mechanisms underlying these loci remains difficult. To help prioritize causal variants and interpret their functions, computational methods have been developed to predict regulatory effects of non-coding variants. An emerging approach to variant annotation is deep learning models that predict regulatory functions from DNA sequences alone. While such models have been trained on large publicly available dataset such as ENCODE, neuropsychiatric trait-related cell types are under-represented in these datasets, thus there is an urgent need of better tools and resources to annotate variant functions in such cellular contexts. To fill this gap, we collected a large collection of neurodevelopment-related cell/tissue types, and trained deep Convolutional Neural Networks (ResNet) using such data. Furthermore, our model, called MetaChrom, borrows information from public epigenomic consortium to improve the accuracy via transfer learning. We show that MetaChrom is substantially better in predicting experimentally determined chromatin accessibility variants than popular variant annotation tools such as CADD and delta-SVM. By combining GWAS data with MetaChrom predictions, we prioritized 31 SNPs for Schizophrenia, suggesting potential risk genes and the biological contexts where they act. In summary, MetaChrom provides functional annotations of any DNA variants in the neuro-development context and the general method of MetaChrom can also be extended to other disease-related cell or tissue types.     

Review: How was metazoan threshold crossed? The hypothetical Urmetazoa.

The origin of Metazoa remained--until recently--the most enigmatic of all phylogenetic problems. Sponges [Porifera] as "living fossils", positioned at the base of multicellular animals, have been used to answer basic questions in metazoan evolution by molecular biological techniques. During the last few years, cDNAs/genes coding for informative proteins have been isolated and characterized from sponges, especially from the marine demosponges Suberites domuncula and Geodia cydonium. The analyses of their deduced amino acid sequences allowed a molecular biological resolution of the monophyly of Metazoa. Molecules of the extracellular matrix/basal lamina, with the integrin receptor, fibronectin and galectin as prominent examples, cell-surface receptors (tyrosine kinase receptors), elements of nerve system/sensory cells (metabotropic glutamate receptor), homologs/modules of an immune system [immunoglobulin-like molecules, SRCR- and SCR-repeats, cytokines, (2-5)A synthetase], as well as morphogens (myotrophin) classify the Porifera as true Metazoa. As "living fossils", provided with simple, primordial molecules allowing cell-cell and cell-matrix adhesion, as well as processes of signal transduction as known in a more complex manner from higher Metazoa, sponges also show peculiarities. Tissues of sponges are rich in telomerase activity, suggesting a high plasticity in the determination of cell lineages. It is concluded that molecular biological studies with sponges as models will not only help to understand the evolution to the Metazoa but also the complex, hierarchical regulatory network of cells in higher Metazoa [reviewed in Progress in Molecular Subcellular Biology, vols. 19, 21 (1998) Springer Verlag]. The hypothetical ancestral animal, the Urmetazoa, from which the metazoan lineages diverged (more than 600 MYA), may have had the following characteristics: cell adhesion molecules with intracellular signal transduction pathways, morphogens/growth factors forming gradients, a functional immune system, and a primordial nerve cell/receptor system.