Reference / Project information
About Metazome
A comparative genomics resource for twenty four metazoan proteomes clustered at nine ancestral nodes, with public annotations from RefSeq, UniProt, Ensembl, and JGI.
Overview
Metazome is a joint project of the U.S. Department of Energy's Joint Genome Institute and the University of California to facilitate comparative genomic studies among metazoans. Clusters of orthologous and paralogous genes that represent the modern descendants of ancestral gene sets are constructed at key phylogenetic nodes. These clusters allow easy access to clade specific orthology and paralogy relationships as well as clade specific genes and gene expansions.
As of version 2.0.4, Metazome provides access to twenty four sequenced and annotated metazoan genomes, clustered at nine evolutionarily significant nodes. Where possible, each gene has been annotated with PFAM, KOG, KEGG, and PANTHER assignments, and publicly available annotations from RefSeq, UniProt, Ensembl, and JGI are hyperlinked and searchable.
Included organisms
Genes from the following organisms are clustered in release 2.0.4 of Metazome. Each organism name links to its genome record, which lists the assembly source, clade path, and related nodes.
| Organism | Common name | Source | Node |
|---|---|---|---|
| Aedes aegypti | Yellow fever mosquito | AaegL 1 | Arthropod |
| Anopheles gambiae | Mosquito | AgamP3 | Arthropod |
| Bombyx mori | Silkworm | As distributed with Metazome 2.0.4 | Arthropod |
| Branchiostoma floridae | Amphioxus | Assembly released 11 April 2006 | Chordate |
| Caenorhabditis briggsae | Briggsae worm | WormBase release WS164 | Protostome |
| Caenorhabditis elegans | Worm | WormBase release WS164 | Protostome |
| Canis familiaris | Dog | NCBI gene models build 2 version 1, assembly released by Broad, 10 May 2005 | Tetrapod |
| Ciona intestinalis | Sea squirt intestinalis | JGI 2 | Chordate |
| Ciona savignyi | Sea squirt savignyi | CSAV 2.0 | Chordate |
| Danio rerio | Zebrafish | Zv6 | Teleost |
| Drosophila melanogaster | Fruit fly | BDGP 4 | Arthropod |
| Gallus gallus | Chicken | Assembly WASHUC 1, March 2004, gene build December 2005 | Tetrapod |
| Gasterosteus aculeatus | Stickleback | Broad S1 | Teleost |
| Homo sapiens | Human | NCBI 36 build | Tetrapod |
| Lottia gigantea | Owl limpet | Early p02S assembly | Protostome |
| Monodelphis domestica | Opossum | MonDom 4 | Tetrapod |
| Mus musculus | Mouse | NCBI m36 build | Tetrapod |
| Nematostella vectensis | Sea anemone | Assembly and annotation distributed with Metazome 2.0.4 | Eumetazoa |
| Oryzias latipes | Medaka | MEDAKA 1 | Teleost |
| Rattus norvegicus | Rat | RGSC 3.4 | Tetrapod |
| Strongylocentrotus purpuratus | Sea urchin | NCBI gene build 2 version 1 | Deuterostome |
| Takifugu rubripes | Fugu pufferfish | Assembly and annotation distributed with Metazome 2.0.4 | Teleost |
| Tribolium castaneum | Red flour beetle | NCBI gene build 1 version 1 | Arthropod |
| Xenopus tropicalis | Frog | Assembly and annotation distributed with Metazome 2.0.4 | Tetrapod |
Nodes
Clustering is used to group extant genes into sets representing the ancestral genes that existed just prior to various significant evolutionary events, called nodes. Extant genes have been clustered at nodes representing the following speciation events.
| Node | Gene families | Singletons | Largest family | Definition |
|---|---|---|---|---|
| Eumetazoa | 31,902 | 203,761 | 496, Notch 1 related | Gene families representing the most recent common ancestor of bilateria and cnidaria. |
| Bilateria | 30,367 | 192,975 | 489, Helix loop helix DNA binding domain containing protein | Gene families representing the most recent common ancestor of deuterostomes and protostomes. |
| Protostome | 16,314 | 75,876 | 141, UDP glucuronosyltransferase related | Gene families representing the most recent common ancestor of lophotrochozoans and ecdysozoans. |
| Deuterostome | 21,031 | 126,047 | 230, Major Facilitator Superfamily domain containing protein | Gene families representing the most recent common ancestor of hemichordates and echinoderms. |
| Arthropod | 8,663 | 32,389 | 49, Odorant binding protein AgamOBP3 | Gene families representing the most recent common ancestor of coleoptera and diptera. |
| Chordate | 17,527 | 91,781 | 189, RAB5A | Gene families representing the most recent common ancestor of cephalochordates and olfactores. |
| Vertebrate | 18,382 | 87,487 | 121, G protein signaling regulator | Gene families representing the most recent common ancestor of teleosts and tetrapods. |
| Teleost | 16,213 | 33,235 | 26, Glicacolin related and C1q domain containing protein | Gene families representing the most recent common ancestor of euteleostei and otocephala. |
| Tetrapod | 15,654 | 62,031 | 30, ETS domain containing predicted transcription factor | Gene families representing the most recent common ancestor of amniota and amphibia. |
Clustering methodology
Clustering is performed hierarchically, from the crown nodes to the root. At each node, ingroup paralog clustering is followed by outgroup ortholog clustering. Ingroup and outgroup relationships are defined as follows. All descendant nodes, internal as well as crown, that connect through the same deepest branch to a given node are considered ingroup nodes with respect to that node. Descendant nodes that connect through the opposite deepest branch are considered outgroup nodes.
For example, human, mouse, rat, dog, opossum, chicken, frog, and the teleosts are all ingroup with respect to the chordate node. With respect to the vertebrate node, teleosts are an outgroup for tetrapods, and tetrapods are an outgroup for teleosts. Gene families are created at a given node by first creating paralog clusters among the ingroup members, with ingroup genes with greater similarity to each other than to their most similar outgroup gene being placed in the paralogous cluster.
Fourfold degenerate transversion rates (4DTV, one measure of the divergence time between genes based on the number of neutral substitutions) are used to preclude clustering of paralogs with low degrees of similarity that are drawn together simply by the lack of a high scoring hit to an outgroup gene. After paralog clustering along both branches, clusters are merged across branches to capture orthologous relationships.
This merging uses a combination of Mutual Best Hit (MBH, based on blastp sequence similarity) and synteny. Synteny is the hypothesis that disjoint genomic regions that contain a fair number of genes with high quality hits to each other likely share a common origin in a speciation, whole genome duplication, or segmental duplication event, though only the former is directly relevant in orthology identification. This process is repeated as the clustering moves from the crown nodes, which are individual organisms, to the root.
By construction, for any given node and gene family, the members of that family remain together in any larger gene family defined at a more ancient node closer to the root. The clustering is rigorously hierarchical. Every gene from an organism present at a particular node is in one and only one cluster at that node, so the clusters are hard rather than fuzzy.
Some clusters may contain only one extant gene. These singleton clusters can come from fast evolution leading to so much sequence divergence that sequence similarity based clustering is confounded, from gene loss, or from gene calling errors.
Clustering statistics
Release 2.0.4 reports the following gene family and singleton counts by node. The largest cluster size and defline are shown for each node.
| Node | Gene families | Singletons | Largest cluster size and defline |
|---|---|---|---|
| Eumetazoa | 31,902 | 203,761 | 496, Notch 1 related |
| Bilateria | 30,367 | 192,975 | 489, Helix loop helix DNA binding domain containing protein |
| Protostome | 16,314 | 75,876 | 141, UDP glucuronosyltransferase related |
| Deuterostome | 21,031 | 126,047 | 230, Major Facilitator Superfamily domain containing protein |
| Arthropod | 8,663 | 32,389 | 49, Odorant binding protein AgamOBP3 |
| Chordate | 17,527 | 91,781 | 189, RAB5A |
| Vertebrate | 18,382 | 87,487 | 121, G protein signaling regulator |
| Teleost | 16,213 | 33,235 | 26, Glicacolin related and C1q domain containing protein |
| Tetrapod | 15,654 | 62,031 | 30, ETS domain containing predicted transcription factor |
Metazome team
| Role | Contributors |
|---|---|
| Software | David M. Goodstein, Rusty Howson, Rochak Neupane, Shengqiang Shu |
| Analysis | Bill Dirks, Uffe Hellsten, Therese Mitros, Dan Rokhsar |
Data sources
Metazome integrates public annotations and sequence resources. The links below identify the sources referenced in this release record and the public resources used by the comparative framework.