|Home | About | Journals | Submit | Contact Us | Français|
Choanoflagellates are the closest known relatives of metazoans. To discover potential molecular mechanisms underlying the evolution of metazoan multicellularity, we sequenced and analysed the genome of the unicellular choanoflagellate Monosiga brevicollis. The genome contains approximately 9,200 intron-rich genes, including a number that encode cell adhesion and signalling protein domains that are otherwise restricted to metazoans. Here we show that the physical linkages among protein domains often differ between M. brevicollis and metazoans, suggesting that abundant domain shuffling followed the separation of the choanoflagellate and metazoan lineages. The completion of the M. brevicollis genome allows us to reconstruct with increasing resolution the genomic changes that accompanied the origin of metazoans.
Choanoflagellates have long fascinated evolutionary biologists for their marked similarity to the ‘feeding cells’ (choanocytes) of sponges and the possibility that they might represent the closest living relatives of metazoans1,2. Over the past decade or so, evidence supporting this relationship has accumulated from phylogenetic analyses of nuclear and mitochondrial genes3–6, comparative genomics between the mitochondrial genomes of choanoflagellates, sponges and other metazoans7,8, and the finding that choanoflagellates express homologues of metazoan signalling and adhesion genes9–12. Furthermore, species-rich phylogenetic analyses demonstrate that choanoflagellates are not derived from metazoans, but instead represent a distinct lineage that evolved before the origin and diversification of metazoans (Fig. 1a, Supplementary Fig. 1 and Supplementary Note 3.1)8,13. By virtue of their position on the tree of life, studies of choanoflagellates provide an unparallelled window into the nature of the unicellular and colonial progenitors of metazoans14.
Choanoflagellates are abundant and globally distributed microbial eukaryotes found in marine and freshwater environments15,16. Like sponge choanocytes, each cell bears an apical flagellum surrounded by a distinctive collar of actin-filled microvilli, with which choanoflagellates trap bacteria and detritus (Fig. 1b). Using this highly effective means of prey capture, choanoflagellates link bacteria to higher trophic levels and thus have critical roles in oceanic carbon cycling and in the microbial food web17,18.
More than 125 choanoflagellate species have been identified, and all species have a unicellular life-history stage. Some can also form simple colonies of equipotent cells, although these differ substantially from the obligate associations of differentiated cells in metazoans19. Studies of basal metazoans indicate that the ancestral metazoan was multicellular and had differentiated cell types, an epithelium, a body plan and regulated development including gastrulation. In contrast, the last common ancestor of choanoflagellates and metazoans was unicellular or possibly capable of forming simple colonies, underscoring the abundant biological innovation that accompanied metazoan origins.
Despite their evolutionary and ecological importance, little is known about the genetics and cell biology of choanoflagellates. To gain insight into the biology of choanoflagellates and to reconstruct the genomic changes attendant on the early evolution of metazoans, we sequenced the genome of the choanoflagellate M. brevicollis and compared it with genomes from metazoans and other eukaryotes.
The ~41.6 megabase (Mb) M. brevicollis genome contains approximately 9,200 genes (Supplementary Note 1 and Supplementary Note 2) and is comparable in size to the genomes of filamentous fungi (~30–40 Mb) and other free-living unicellular eukaryotes (for example, small diatoms at ~20–35 Mb20 and ichthyosporeans at ~20–25 Mb21). Metazoan genomes are typically larger, with few exceptions22.
M. brevicollis genes have several distinguishing structural features (Table 1). Whereas the M. brevicollis genome is compact, its genes are almost as intron-rich as human genes (6.6 introns per M. brevicollis gene versus 7.7 introns per human gene). M. brevicollis introns are short (averaging 174 bp) relative to metazoan introns, and with few exceptions do not include the extremely long introns found in some metazoan genes (Supplementary Fig. 2 and Supplementary Note 3.3).
Comparisons of intron positions in a set of conserved genes from M. brevicollis, diverse metazoans and representative intron-rich fungi, plants and a ciliate reveal that the last common ancestor of choanoflagellates and metazoans had genes at least as intron-rich as those of modern choanoflagellates (Fig. 2, Supplementary Fig 3 and Supplementary Fig 4, and Supplementary Note 3.3). Notably, these analyses reveal that the eumetazoan ancestor contained a substantially higher density of introns than the last common ancestor of choanoflagellates and metazoans. This is consistent with a proliferation of introns during the early evolution of the Metazoa23.
The M. brevicollis genome provides unprecedented insight into the early evolution of metazoan genes. Annotation of the M. brevicollis genome using Pfam and SMART (two protein domain prediction databases) reveals 78 protein domains that are exclusive to choanoflagellates and metazoans, only two of which have been reported previously in choanoflagellates (Supplementary Table 4)10. Because genomic features shared by M. brevicollis and metazoans were probably present in their last common ancestor, this study extends the evolutionary history of a cohort of important protein domains to the premetazoan era. Many of these domains are central to cell signalling and adhesion processes in metazoans, suggesting a role in the origin of multicellularity. In contrast, metazoan genomic features that are missing from the M. brevicollis genome may have evolved within the metazoan lineage, or may have existed in the last common ancestor with choanoflagellates and were subsequently lost on the stem leading to M. brevicollis. Presumably, there are many genomic features that evolved in the metazoan lineage, and the M. brevicollis genome provides our first glimpse at the complement of genes and protein domains that predate metazoan origins.
To investigate further the extent to which molecular components required for metazoan multicellularity evolved before the origin of metazoans, we performed targeted searches in the M. brevicollis genome and representative metazoan, fungal and plant genomes for homologues of critical metazoan cell adhesion, cell signalling and transcription factor protein families.
A critical step in the transition to multicellularity was the evolution of mechanisms for stable cell adhesion. M. brevicollis encodes a diverse array of cell adhesion and extracellular matrix (ECM) protein domains previously thought to be restricted to metazoans (Fig. 3). At least 23 M. brevicollis genes encode one or more cadherin domains, homologues of which are required for cell sorting and adhesion during metazoan embryogenesis24, and 12 genes encode C-type lectins, 2 of which are transmembrane proteins. Whereas soluble C-type lectins have functions ranging from pathogen recognition to ECM organization, transmembrane C-type lectins mediate specific adhesive activities such as contact between leukocytes and vascular endothelial cells, cell recognition, and molecular uptake by endocytosis25–27.
The genome of M. brevicollis also contains integrin-α and immunoglobulin domains—cell adhesion domains formerly thought to be restricted to metazoans. In metazoans, integrin-α- and integrin-β-domain-containing proteins heterodimerize before binding to ECM proteins such as collagen28. We find that M. brevicollis has at least 17 integrin-α-domain-containing proteins, but no integrin-β domains. Metazoan immunoglobulin-domain-containing proteins have both adhesive and immune functions. The M. brevicollis genome encodes a total of five immunoglobulin domains that show affinity for the I-set, V-set or C2-set subfamilies, but not the vertebrate-specific C1-set subfamily. In contrast to M. brevicollis, metazoan genomes possess from ~150 to ~1,500 immunoglobulin domains (Supplementary Table 7), suggesting that the radiation of the immunoglobulin superfamily occurred after the divergence of choanoflagellates and metazoans.
The finding in M. brevicollis of cell adhesion domains that were previously known only in metazoans has two important implications. First, the common ancestor of metazoans and choanoflagellates possessed several of the critical structural components used for multicellularity in modern metazoans. Second, given the absence of evidence for stable cell adhesion in M. brevicollis, this also suggests that homologues of metazoan cell adhesion domains may act to mediate interactions between M. brevicollis and its extracellular environment.
As the targets of many adhesion receptors, the question of whether metazoan-type ECM proteins and domains evolved before or after the transition to multicellularity is of great interest. In metazoans, collagens are ECM proteins that polymerize to form a major component of the basement membrane of epithelia and have been invoked as a potential ‘key innovation’ during the transition to multicellularity29. We find five collagen-domain-encoding genes in the M. brevicollis genome, two of which encode the diagnostic Gly-X-Y repetitive sequence motif (where X and Y are frequently proline and hydroxyproline, respectively) in an arrangement similar to metazoan collagens30. Other ECM-associated domains previously known only from metazoans that occur in M. brevicollis include laminin domains (an important class that contributes to the basement membrane), the reeler domain (found in the neuronal ECM protein reelin31) and the ependymin domain (an extracellular glycoprotein found in cerebrospinal fluid32; Fig. 3 and Supplementary Table 4).
The discovery of putatively secreted ECM proteins in a free-living choanoflagellate suggests that elements of the metazoan ECM evolved in contact with the external environment before being sequestered within an epithelium. Although some choanoflagellates secrete extracellular structures or adhere to form colonial assemblages 19,33,34, M. brevicollis is not known to do so. Instead, these ECM protein homologues in M. brevicollis may mediate an analogous process such as substrate attachment.
Against the backdrop of abundant conservation of cell adhesion and ECM protein domains among the genomes of M. brevicollis and metazoans, it is important to note the differences. Individual cell adhesion and ECM-associated domains in the M. brevicollis genome often occur in unique arrangements, and clear orthologues of specific metazoan adhesion proteins are rarely found. Although the domains associated with metazoan adhesion and ECM proteins were present in the ancestor of choanoflagellates and metazoans, the canonical metazoan adhesion protein architectures35 probably evolved after the divergence of the two lineages.
Our analysis of the M. brevicollis genome reveals little evidence that metazoan-specific signalling pathways were present in the last common ancestor of choanoflagellates and metazoans. Many pathways are missing entirely, and M. brevicollis genes with some similarity to metazoan signalling machinery are largely found to share conserved domains without aligning across the full span of what are often complex multidomain proteins (for example, epidermal growth factor (EGF) repeats are common to Notch, and also to many other proteins; Supplementary Table 8). Specifically, no receptors or ligands were identified from the NHR (nuclear hormone receptor), WNT and TGF-β signalling pathways. The only evidence of the JAK (Janus kinase)/STAT (signal transducer and activator of transcription) pathway is an apparent STAT-like gene that encodes a STAT DNA-binding domain and a partial SH2 domain. Convincing evidence is also lacking for the Toll signalling pathway—a signalling system important both for development and for innate immunity in metazoans.
Nonetheless, the genome of M. brevicollis does provide insights into the evolution of Notch and hedgehog (Hh) signalling pathways. Cassettes of protein domains found in metazoan Notch receptors (EGF, NL and ANK (ankyrin repeats)) are encoded on separate M. brevicollis genes in arrangements that differ from metazoan Notch proteins, and definitive domains, such as the NOD domain and the MNNL region, are absent (Fig. 4a).
Homologues of hedgehog, dispatched and patched genes are also present; however, there is no evidence for smoothened nor for its defining frizzled domain. In metazoans, hedgehog consists of an amino-terminal signalling domain and carboxy-terminal hedgehog/ intein (Hint) domain responsible for autocatalytically cleaving the protein. In one M. brevicollis hedgehog-like protein, a hedgehog N-terminal signalling domain is found at the N terminus of a large transmembrane protein that, instead of a Hint domain, includes von Willebrand A, cadherin, TNFR (tumor necrosis factor receptor), furin and EGF domains. Similar proteins are found in the sponge Amphimedon queenslandica and in the cnidarian Nematostella vectensis36, revealing that the M. brevicollis genome captures an ancestral arrangement of protein domains rather than representing a lineage-specific domain-shuffling event. Another M. brevicollis hedgehog-like protein contains a Hint domain—a key region involved in the autocatalytic processing of hedgehog (Fig. 4b). The identification of a hedgehog-like gene in a choanoflagellate is not without precedent. A distinct Hint-domain-containing protein, named Hoglet, was identified in the distantly related Monosiga ovata12, supporting the idea that isolated signalling components were present in the last common ancestor of choanoflagellates and metazoans.
Phosphotyrosine (pTyr)-based signalling was considered unique to metazoans until its recent discovery in choanoflagellates9,11. The key domains involved in pTyr signalling are found in abundance in the M. brevicollis genome: tyrosine kinase domains that phosphorylate tyrosine (~120 occurrences), pTyr-specific phosphatases (PTP) that remove the phosphate modification (~30) and SH2 domains that bind pTyr-containing peptides (~80) (Supplementary Fig. 7). In contrast, these domains are rare in non-metazoans; for example, S. cerevisiae has no tyrosine kinase domains, only three PTP domains and a single SH2 domain. These findings support a model in which the full set of pTyr signalling machinery evolved before the separation of the choanoflagellate and metazoan lineages.
Although pTyr signalling machinery is present in metazoans and choanoflagellates, the mode of usage in M. brevicollis may be distinct from that in metazoans. A simple metric for the use of a particular domain is the range of domain types with which it is found in combination37. In the M. brevicollis genome, more than half of the observed pairwise domain combinations involving tyrosine kinase, PTP and SH2 domains are distinct from those seen in any metazoan genome (Fig. 5 and Supplementary Note 3.7). In contrast, for other sets of common signalling domains (those involved in phosphoserine/ threonine (pSer/Thr), Ras–GTP and Rho–GTP signalling), most observed combinations are shared between M. brevicollis and metazoans. These observations are consistent with a simple model in which pSer/Thr, Ras–GTP and Rho–GTP signalling were fully elaborated before the branching of the choanoflagellate and metazoan lineages (consistent with the presence of these systems in other eukaryotes, including fungi, Dictyostelium and plants). In contrast, simple pTyr signalling may have emerged in the common ancestor and diverged radically between choanoflagellates and metazoans.
The core transcriptional apparatus of M. brevicollis is, in many ways, typical of most eukaryotes examined to date (Supplementary Table 10) including, for example, all 12 RNA polymerase II subunits and most of the transcription elongation factors (TFIIS, NELF, PAF, DSIF and P-TEFb, but not elongin). However, homologues of the largest subunit of TFIIF and several subunits of TFIIH are apparently lacking from the genome and the expressed-sequence-tag collection (Supplementary Fig. 8), reminiscent of the absence of several basal factors from the Giardia lamblia genome, suggesting alternative strategies for interacting with core promoter elements38. Similarly, only a limited number of general co-activators are identifiable in M. brevicollis, including the components of several chromatin-remodelling complexes (Supplementary Fig. 9 and Supplementary Note 3.8).
Perhaps not surprisingly, M. brevicollis possesses members from most of the ubiquitous families of eukaryotic transcription factors (Supplementary Fig. 10). Most of the predicted transcription factors are zinc-coordinating; approximately 44% are C2H2-type zinc fingers. Eight proteins (5% of a total of 155 predicted transcription factors) are forkhead transcription factors, otherwise known only from metazoans and fungi.
The homeodomain transcription factors are an ancient protein family found in all known eukaryotes. At least two major superclasses of homeodomain proteins evolved before the origin of metazoans: those containing homeodomains of 60 amino acids (the ‘typical’, or non-TALE superclass), and those containing homeodomains of more than 63 amino acids (the TALE superclass)39. The M. brevicollis genome encodes only two homeodomain proteins, both of which group with the MEIS sub-class of TALE homeodomains (Supplementary Fig. 12). Apparently, genes encoding non-TALE homeodomain proteins have been lost in the lineage leading to M. brevicollis. Bona fide HOX class homeobox genes—a subclass of the non-TALE superclass—are absent from both M. brevicollis and the Amphimedon queenslandica (demosponge) genome sequence reads, indicating that this characteristic metazoan gene family probably emerged along the stem leading to eumetazoans40.
M. brevicollis contains a subset of the transcription factor families previously thought to be specific to metazoans. Members of the metazoan p53, Myc and Sox/TCF families were identified, whereas many transcription factor families associated with metazoan patterning and development (ETS, HOX, NHR, POU and T-box) seem to be absent (Fig. 3).
Choanoflagellates, sponges and other metazoans last shared a unicellular common ancestor in the late Precambrian period, more than 600 million years ago41,42. Although the origin of metazoans was a pivotal event in life’s history, little is known about the genetic underpinnings of the requisite transition to multicellularity. Comparisons of modern genomes provide our most direct insights into the ancient genomic conditions from which metazoans emerged. By comparing choanoflagellate and metazoan genomes, we infer that their common ancestor had intron-rich genes, some of which encoded protein domains characteristically associated with cell adhesion and the ECM in animals.
In addition to containing protein domains associated with metazoan cell adhesion, M. brevicollis possesses a surprising abundance of tyrosine kinases and their downstream signalling targets. In contrast, components of most other intercellular signalling pathways, as well as many of the diverse transcription factors that comprise the developmental toolkit of modern animals, are absent. These presumably reached their modern form on the metazoan stem, although it is formally possible that they were in place much earlier and degenerated in the M. brevicollis lineage. Likewise, it is possible that the last common ancestor of choanoflagellates and metazoans had an early form of multicellularity that became more robust in metazoans and was lost in the choanoflagellate lineage. In any event, the evolutionary distance between choanoflagellates and metazoans is substantial, and evidently few, if any, intermediate lineages survived. There are, for example, no other known microbial eukaryotes that possess any of the eight developmental signalling pathways characteristic of metazoans.
The mechanism of invention of new genes on the metazoan stem, and their integration to create the emergent network of cell signalling and transcriptional regulation fundamental to metazoan biology, remains mysterious. Domain shuffling, which has frequently been proposed as an important mechanism for the evolution of metazoan multidomain proteins43,44, is implicated by the presence of essential metazoan signalling domains in M. brevicollis that appear in unique combinations relative to animals. For pTyr-based signalling in particular, the marked divergence of domain combinations suggests that this mode of cellular interaction existed in a nascent form in the common choanoflagellate–metazoan ancestor, and was subsequently specialized and elaborated on in each lineage.
Given the limited transcription factor diversity in M. brevicollis, it is notable that the genome encodes representatives of the otherwise metazoan-specific p53, Mycand Sox/TCF transcription factor families. These transcription factors may have had early and critical roles in the evolution of metazoan ancestors by regulating the differential expression of genes to allow multiple cell types to exist in a single organism, and their study in choanoflagellates is a promising future direction.
The M. brevicollis sequence opens the door to genome-enabled studies of choanoflagellates, a diverse group of microbial eukaryotes that are important in their own right as bacterial predators in both marine and freshwater ecosystems. Although M. brevicollis is strictly unicellular, other choanoflagellates facultatively form colonies, and the modulation of these associations by cell signalling, adhesion, transcriptional regulation and environmental influences is poorly understood. An integrative approach that unites studies of choanoflagellate genomes, cell biology and ecology with the biogeochemistry of the Precambrian promises to reveal the intrinsic and extrinsic factors underlying metazoan origins.
All analyses described were performed on Version 1.0 of the genome sequence. Details can be found in the Supplementary Information.
Using physical separation techniques combined with antibiotic treatments, a culture line with only a single bacterial food source, Flavobacterium sp., was developed. The GC content of Flavobacterium DNA (33%) is sufficiently different from that of M. brevicollis (55%) to allow separation over a CsCl gradient. M. brevicollis genomic DNA was used to construct replicate libraries containing inserts of 2–3 kilobases (kb), 6–8 kb and 35–40 kb.
The draft sequence of the M. brevicollis genome was generated from ~8.5-fold redundant paired-end whole-genome shotgun sequence coverage (Supplementary Information). Sequence data derived from six whole-genome shotgun libraries were assembled using release 2.9.2 of the whole-genome shotgun assembler Jazz. Completeness of the draft genome was assessed by capturing ~98.5% of sequenced expressed sequence tags.
Gene models (9,196) were predicted and annotated using the Joint Genome Institute (JGI) Annotation Pipeline (Supplementary Information).
Homologues of 473 highly conserved genes from M. brevicollis and representative eukaryotes were aligned to reveal the position and phylogenetic distribution of 1,989 highly reliable intron splice sites at 1,054 conserved positions. The evolutionary history of introns in orthologous genes was inferred using Dollo parsimony, Roy-Gilbert maximum likelihood and Csuros maximum likelihood45–47.
Gene models containing metazoan signalling, adhesion and transcription factor domains were identified by text and protein domain ID searches of the JGI M. brevicollis genome portal, local BLAST searches within the M.brevicollis genome scaffolds, the online Pfam and SMART tools, and reciprocal BLAST searches in the NCBI non-redundant protein database (Supplementary Information).
This work was performed under the auspices of the US Department of Energy's Office of Science, Biological and Environmental Research Program, and by the University of California, Lawrence Livermore National Laboratory, Lawrence Berkeley National Laboratory, and Los Alamos National Laboratory. Work in the King laboratory is supported by funding from the Gordon and Betty Moore Foundation, the Pew Scholars program, and R. Melmon. The Rokhsar group is supported by the Gordon and Betty Moore Foundation and R. Melmon. We thank J. Stajich, P. Johnson and R. Lusk for discussions, E. Hare, E. Meltzer and K. Osoegawa for technical advice, E. Begovic for assistance with Fig. 1, M. Dayel and N. Patel for critical reading of the manuscript, and S. Carroll for early support of this project. N.K. is a Scholar in the Integrated Microbial Biodiversity Program of the Canadian Institute for Advanced Research.
Supplementary Information is linked to the online version of the paper at www.nature.com/nature.