Early evolution of the T-box transcription factor family

Never-ending story of Brachyury: From short-tailed mice to tailless primates

The history of developmental biology starts from the almost simultaneous discoveries of the Organizer of axial structures in amphibians by Spemann and Mangold in Freiburg and of the Brachyury mutant in mammals by the Dobrovolskaya-Zavadskaya laboratory at the Curie Institute and its follow-up studies in the Leslie Dunn laboratory at Columbia University. Following the Organizer's discovery, the inductive activity of several other embryonic tissues was found, including that of the ear primordium by Boris Balinsky in Kiev. Initially, the experimental embryological and genetic lines of research existed independently of each other, but after they met at the bench of Salome Gluecksohn, they strengthened and cross-fertilized each other, eventually leading to developmental genetics, which later became known as developmental biology. It appears that the regulatory activities of Brachyury and related T-box proteins in general are at the heart of the development of all vertebrates. These activities are fundamental and have been discovered in several model organisms subjected to mutagenesis, exemplified by the story of George Streisinger's discovery of the no tail mutant in zebrafish. This essay describes the history of Brachyury studies, their connection to an idea of embryonic induction by Organizer, and an impact of Brachyury and related genes on various fields of research from embryology and cell biology to medical genetics and evolutionary theory.

New insights into mesoderm and endoderm development, and the nature of the onychophoran blastopore

BACKGROUND

Early during onychophoran development and prior to the formation of the germ band, a posterior tissue thickening forms the posterior pit. Anterior to this thickening forms a groove, the embryonic slit, that marks the anterior-posterior orientation of the developing embryo. This slit is by some authors considered the blastopore, and thus the origin of the endoderm, while others argue that the posterior pit represents the blastopore. This controversy is of evolutionary significance because if the slit represents the blastopore, then this would support the amphistomy hypothesis that suggests that a slit-like blastopore in the bilaterian ancestor evolved into protostomy and deuterostomy.

RESULTS

In this paper, we summarize our current knowledge about endoderm and mesoderm development in onychophorans and provide additional data on early endoderm- and mesoderm-determining marker genes such as Blimp, Mox, and the T-box genes.

CONCLUSION

We come to the conclusion that the endoderm of onychophorans forms prior to the development of the embryonic slit, and thus that the slit is not the primary origin of the endoderm. It is thus unlikely that the embryonic slit represents the blastopore. We suggest instead that the posterior pit indeed represents the lips of the blastopore, and that the embryonic slit (and surrounding tissue) represents a morphologically superficial archenteron-like structure. We conclude further that both endoderm and mesoderm development are under control of conserved gene regulatory networks, and that many of the features found in arthropods including the model Drosophila melanogaster are likely derived.

Conserved enhancers control notochord expression of vertebrate Brachyury

The cell type-specific expression of key transcription factors is central to development and disease. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three conserved Brachyury-controlling notochord enhancers, T3, C, and I, in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, in cis deletion of all three enhancers in mouse abolishes Brachyury/T/Tbxt expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. The three Brachyury-driving notochord enhancers are conserved beyond mammals in the brachyury/tbxtb loci of fishes, dating their origin to the last common ancestor of jawed vertebrates. Our data define the vertebrate enhancers for Brachyury/T/TBXTB notochord expression through an auto-regulatory mechanism that conveys robustness and adaptability as ancient basis for axis development.

The Origin of Metazoan Multicellularity: A Potential Microbial Black Swan Event

The emergence of animals from their unicellular ancestors is a major evolutionary event. Thanks to the study of diverse close unicellular relatives of animals, we now have a better grasp of what the unicellular ancestor of animals was like. However, it is unclear how that unicellular ancestor of animals became the first animals. To explain this transition, two popular theories, the choanoblastaea and the synzoospore, have been proposed. We will revise and expose the flaws in these two theories while showing that, due to the limits of our current knowledge, the origin of animals is a biological black swan event. As such, the origin of animals defies retrospective explanations. Therefore, we should be extra careful not to fall for confirmation biases based on few data and, instead, embrace this uncertainty and be open to alternative scenarios. With the aim to broaden the potential explanations on how animals emerged, we here propose two novel and alternative scenarios. In any case, to find the answer to how animals evolved, additional data will be required, as will the hunt for microscopic creatures that are closely related to animals but have not yet been sampled and studied.

The evolutionary history of Brachyury genes in Hydrozoa involves duplications, divergence, and neofunctionalization

Brachyury, a member of T-box gene family, is widely known for its major role in mesoderm specification in bilaterians. It is also present in non-bilaterian metazoans, such as cnidarians, where it acts as a component of an axial patterning system. In this study, we present a phylogenetic analysis of Brachyury genes within phylum Cnidaria, investigate differential expression and address a functional framework of Brachyury paralogs in hydrozoan Dynamena pumila. Our analysis indicates two duplication events of Brachyury within the cnidarian lineage. The first duplication likely appeared in the medusozoan ancestor, resulting in two copies in medusozoans, while the second duplication arose in the hydrozoan ancestor, resulting in three copies in hydrozoans. Brachyury1 and 2 display a conservative expression pattern marking the oral pole of the body axis in D. pumila. On the contrary, Brachyury3 expression was detected in scattered presumably nerve cells of the D. pumila larva. Pharmacological modulations indicated that Brachyury3 is not under regulation of cWnt signaling in contrast to the other two Brachyury genes. Divergence in expression patterns and regulation suggest neofunctionalization of Brachyury3 in hydrozoans.

Conserved enhancer logic controls the notochord expression of vertebrate Brachyury

The cell type-specific expression of key transcription factors is central to development. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three Brachyury-controlling notochord enhancers T3, C, and I in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, deletion of all three enhancers in mouse abolishes Brachyury/T expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. Sequence and functional conservation of Brachyury-driving notochord enhancers with the brachyury/tbxtb loci from diverse lineages of fishes dates their origin to the last common ancestor of jawed vertebrates. Our data define the enhancers for Brachyury/T/TBXTB notochord expression as ancient mechanism in axis development.

The transcription factor optomotor-blind restricts apterous expression through TrxG and PcG genes

The establishment of body pattern is a fundamental process in developmental biology. In Drosophila, the wing disc is subdivided into dorsal (D) and ventral (V) compartments by the D/V boundary. The dorsal fate is adopted by expressing the selector gene apterous (ap). ap expression is regulated by three combinational cis-regulatory modules which are activated by EGFR pathway, Ap-Vg auto-regulatory and epigenetic mechanisms. Here, we found that the Tbx family transcription factor Optomotor-blind (Omb) restricted ap expression in the ventral compartment. Loss of omb induced autonomous initiation of ap expression in the middle third instar larvae in the ventral compartment. Oppositely, over-activation of omb inhibited ap in the medial pouch. All three enhancers apE, apDV and apP were upregulated in omb null mutants, indicating a combinational regulation of ap modulators. However, Omb affected ap expression neither by directly regulating EGFR signaling, nor via Vg regulation. Therefore, a genetic screen of epigenetic regulators, including the Trithorax group (TrxG) and Polycomb group (PcG) genes was performed. We found that knocking down the TrxG gene kohtalo (kto), domino (dom) or expressing the PcG gene grainy head (grh), the ectopic ap in omb mutants was repressed. The inhibition of apDV by kto knockdown and grh activation could contribute to ap repression. Moreover, Omb and the EGFR pathway are genetically parallel in ap regulation in the ventral compartment. Collectively, Omb is a repressive signal for ap expression in the ventral compartment, which requires TrxG and PcG genes.

An ancestral Wnt-Brachyury feedback loop in axial patterning and recruitment of mesoderm-determining target genes

Transcription factors are crucial drivers of cellular differentiation during animal development and often share ancient evolutionary origins. The T-box transcription factor Brachyury plays a pivotal role as an early mesoderm determinant and neural repressor in vertebrates; yet, the ancestral function and key evolutionary transitions of the role of this transcription factor remain obscure. Here, we present a genome-wide target-gene screen using chromatin immunoprecipitation sequencing in the sea anemone Nematostella vectensis, an early branching non-bilaterian, and the sea urchin Strongylocentrotus purpuratus, a representative of the sister lineage of chordates. Our analysis reveals an ancestral gene regulatory feedback loop connecting Brachyury, FoxA and canonical Wnt signalling involved in axial patterning that predates the cnidarian-bilaterian split about 700 million years ago. Surprisingly, we also found that part of the gene regulatory network controlling the fate of neuromesodermal progenitors in vertebrates was already present in the common ancestor of cnidarians and bilaterians. However, while several endodermal and neuronal Brachyury target genes are ancestrally shared, hardly any of the key mesodermal downstream targets in vertebrates are found in the sea anemone or the sea urchin. Our study suggests that a limited number of target genes involved in mesoderm formation were newly acquired in the vertebrate lineage, leading to a dramatic shift in the function of this ancestral developmental regulator.

Shared cis-regulatory modules control expression of the tandem paralogs midline and H15 in the follicular epithelium

The posterior end of the follicular epithelium is patterned by midline (MID) and its paralog H15, the Drosophila homologs of the mammalian Tbx20 transcription factor. We have previously identified two cis-regulatory modules (CRMs) that recapitulate the endogenous pattern of mid in the follicular epithelium. Here, using CRISPR/Cas9 genome editing, we demonstrate redundant activity of these mid CRMs. Although the deletion of either CRM alone generated marginal change in mid expression, the deletion of both CRMs reduced expression by 60%. Unexpectedly, the deletion of the 5' proximal CRM of mid eliminated H15 expression. Interestingly, expression of these paralogs in other tissues remained unaffected in the CRM deletion backgrounds. These results suggest that the paralogs are regulated by a shared CRM that coordinates gene expression during posterior fate determination. The consistent overlapping expression of mid and H15 in various tissues may indicate that the paralogs could also be under shared regulation by other CRMs in these tissues.

Eomes function is conserved between zebrafish and mouse and controls left-right organiser progenitor gene expression via interlocking feedforward loops

The T-box family transcription factor Eomesodermin (Eomes) is present in all vertebrates, with many key roles in the developing mammalian embryo and immune system. Homozygous Eomes mutant mouse embryos exhibit early lethality due to defects in both the embryonic mesendoderm and the extraembryonic trophoblast cell lineage. In contrast, zebrafish lacking the predominant Eomes homologue A (Eomesa) do not suffer complete lethality and can be maintained. This suggests fundamental differences in either the molecular function of Eomes orthologues or the molecular configuration of processes in which they participate. To explore these hypotheses we initially analysed the expression of distinct Eomes isoforms in various mouse cell types. Next we compared the functional capabilities of these murine isoforms to zebrafish Eomesa. These experiments provided no evidence for functional divergence. Next we examined the functions of zebrafish Eomesa and other T-box family members expressed in early development, as well as its paralogue Eomesb. Though Eomes is a member of the Tbr1 subfamily we found evidence for functional redundancy with the Tbx6 subfamily member Tbx16, known to be absent from eutherians. However, Tbx16 does not appear to synergise with Eomesa cofactors Mixl1 and Gata5. Finally, we analysed the ability of Eomesa and other T-box factors to induce zebrafish left-right organiser progenitors (known as dorsal forerunner cells) known to be positively regulated by vgll4l, a gene we had previously shown to be repressed by Eomesa. Here we demonstrate that Eomesa indirectly upregulates vgll4l expression via interlocking feedforward loops, suggesting a role in establishment of left-right asymmetry. Conversely, other T-box factors could not similarly induce left-right organiser progenitors. Overall these findings demonstrate conservation of Eomes molecular function and participation in similar processes, but differential requirements across evolution due to additional co-expressed T-box factors in teleosts, albeit with markedly different molecular capabilities. Our analyses also provide insights into the role of Eomesa in left-right organiser formation in zebrafish.

A mosaic of conserved and novel modes of gene expression and morphogenesis in mesoderm and muscle formation of a larval bivalve

The mesoderm gives rise to several key morphological features of bilaterian animals including endoskeletal elements and the musculature. A number of regulatory genes involved in mesoderm and/or muscle formation (e.g., Brachyury (Bra), even-skipped (eve), Mox, myosin II heavy chain (mhc)) have been identified chiefly from chordates and the ecdysozoans Drosophila and Caenorhabditis elegans, but data for non-model protostomes, especially those belonging to the ecdysozoan sister clade, Lophotrochozoa (e.g., flatworms, annelids, mollusks), are only beginning to emerge. Within the lophotrochozoans, Mollusca constitutes the most speciose and diverse phylum. Interestingly, however, information on the morphological and molecular underpinnings of key ontogenetic processes such as mesoderm formation and myogenesis remains scarce even for prominent molluscan sublineages such as the bivalves. Here, we investigated myogenesis and developmental expression of Bra, eve, Mox, and mhc in the quagga mussel Dreissena rostriformis, an invasive freshwater bivalve and an emerging model in invertebrate evodevo. We found that all four genes are expressed during mesoderm formation, but some show additional, individual sites of expression during ontogeny. While Mox and mhc are involved in early myogenesis, eve is also expressed in the embryonic shell field and Bra is additionally present in the foregut. Comparative analysis suggests that Mox has an ancestral role in mesoderm and possibly muscle formation in bilaterians, while Bra and eve are conserved regulators of mesoderm development of nephrozoans (protostomes and deuterostomes). The fully developed Dreissena veliger larva shows a highly complex muscular architecture, supporting a muscular ground pattern of autobranch bivalve larvae that includes at least a velum muscle ring, three or four pairs of velum retractors, one or two pairs of larval retractors, two pairs of foot retractors, a pedal plexus, possibly two pairs of mantle retractors, and the muscles of the pallial line, as well as an anterior and a posterior adductor. As is typical for their molluscan kin, remodelling and loss of prominent larval features such as the velum musculature and various retractor systems appear to be also common in bivalves.

Supplementary information

The online version contains supplementary material available at 10.1007/s13127-022-00569-5.

The Non-continuum Nature of Eukaryotic Transcriptional Regulation

Eukaryotic transcription factors are versatile mediators of specificity in gene regulation. This versatility is achieved through mutual specification by context-specific DNA binding on the one hand, and identity-specific protein-protein partnerships on the other. This interactivity, known as combinatorial control, enables a repertoire of complex transcriptional outputs that are qualitatively disjoint, or non-continuum, with respect to binding affinity. This feature contrasts starkly with prokaryotic gene regulators, whose activities in general vary quantitatively in step with binding affinity. Biophysical studies on prokaryotic model systems and more recent investigations on transcription factors highlight an important role for folded state dynamics and molecular hydration in protein/DNA recognition. Analysis of molecular models of combinatorial control and recent literature in low-affinity gene regulation suggest that transcription factors harbor unique conformational dynamics that are inaccessible or unused by prokaryotic DNA-binding proteins. Thus, understanding the intrinsic dynamics involved in DNA binding and co-regulator recruitment appears to be a key to understanding how transcription factors mediate non-continuum outcomes in eukaryotic gene expression, and how such capability might have evolved from ancient, structurally conserved counterparts.

Role of T-box genes in cancer, epithelial-mesenchymal transition, and cancer stem cells

Sharing a common DNA binding motif called T-box, transcription factor T-box gene family controls embryonic development and is also involved in cancer progression and metastasis. Cancer metastasis shows therapy resistance and involves complex processes. Among them, epithelial-mesenchymal transition (EMT) triggers cancer cell invasiveness and the acquisition of stemness of cancer cells, called cancer stem cells (CSCs). CSCs are a small fraction of tumor bulk and are capable of self-renewal and tumorsphere formation. Recent progress has highlighted the critical roles of T-box genes in cancer progression, EMT, and CSC function, and such regulatory functions of T-box genes have emerged as potential therapeutic candidates for cancer. Herein we summarize the current understanding of the regulatory mechanisms of T-box genes in cancer, EMT, and CSCs, and discuss the implications of targeting T-box genes as anticancer therapeutics.

Intracellular and Intercellular Gene Regulatory Network Inference From Time-Course Individual RNA-Seq

Gene regulatory network (GRN) inference is an effective approach to understand the molecular mechanisms underlying biological events. Generally, GRN inference mainly targets intracellular regulatory relationships such as transcription factors and their associated targets. In multicellular organisms, there are both intracellular and intercellular regulatory mechanisms. Thus, we hypothesize that GRNs inferred from time-course individual (whole embryo) RNA-Seq during development can reveal intercellular regulatory relationships (signaling pathways) underlying the development. Here, we conducted time-course bulk RNA-Seq of individual mouse embryos during early development, followed by pseudo-time analysis and GRN inference. The results demonstrated that GRN inference from RNA-Seq with pseudo-time can be applied for individual bulk RNA-Seq similar to scRNA-Seq. Validation using an experimental-source-based database showed that our approach could significantly infer GRN for all transcription factors in the database. Furthermore, the inferred ligand-related and receptor-related downstream genes were significantly overlapped. Thus, the inferred GRN based on whole organism could include intercellular regulatory relationships, which cannot be inferred from scRNA-Seq based only on gene expression data. Overall, inferring GRN from time-course bulk RNA-Seq is an effective approach to understand the regulatory relationships underlying biological events in multicellular organisms.

Studying evolution of the primary body axis in vivo and in vitro

The metazoan body plan is established during early embryogenesis via collective cell rearrangements and evolutionarily conserved gene networks, as part of a process commonly referred to as gastrulation. While substantial progress has been achieved in terms of characterizing the embryonic development of several model organisms, underlying principles of many early patterning processes nevertheless remain enigmatic. Despite the diversity of (pre-)gastrulating embryo and adult body shapes across the animal kingdom, the body axes, which are arguably the most fundamental features, generally remain identical between phyla. Recently there has been a renewed appreciation of ex vivo and in vitro embryo-like systems to model early embryonic patterning events. Here, we briefly review key examples and propose that similarities in morphogenesis and associated gene expression dynamics may reveal an evolutionarily conserved developmental mode as well as provide further insights into the role of external or extraembryonic cues in shaping the early embryo. In summary, we argue that embryo-like systems can be employed to inform previously uncharted aspects of animal body plan evolution as well as associated patterning rules.

Tracing the evolutionary history of Ca2+-signaling modulation by human Bcl-2: Insights from the Capsaspora owczarzaki IP3 receptor ortholog

Recently, a functional IP₃R ortholog (CO.IP₃R-A) capable of IP₃-induced Ca²⁺ release has been discovered in Capsaspora owczarzaki, a close unicellular relative to Metazoa. In contrast to mammalian IP₃Rs, CO.IP₃R-A is not modulated by Ca²⁺, ATP or PKA. Protein-sequence analysis revealed that CO.IP₃R-A contained a putative binding site for anti-apoptotic Bcl-2, although Bcl-2 was not detected in Capsaspora owczarzaki and only appeared in Metazoa. Here, we examined whether human Bcl-2 could form a complex with CO.IP₃R-A channels and modulate their Ca²⁺-flux properties using ectopic expression approaches in a HEK293 cell model in which all three IP₃R isoforms were knocked out. We demonstrate that human Bcl-2 via its BH4 domain could functionally interact with CO.IP₃R-A, thereby suppressing Ca²⁺ flux through CO.IP₃R-A channels. The BH4 domain of Bcl-2 was sufficient for interaction with CO.IP₃R-A channels. Moreover, mutating the Lys17 of Bcl-2's BH4 domain, the residue critical for Bcl-2-dependent modulation of mammalian IP₃Rs, abrogated Bcl-2's ability to bind and inhibit CO.IP₃R-A channels. Hence, this raises the possibility that a unicellular ancestor of animals already had an IP₃R that harbored a Bcl-2-binding site. Bcl-2 proteins may have evolved as controllers of IP₃R function by exploiting this pre-existing site, thereby counteracting Ca²⁺-dependent apoptosis.

Evolutionary Cell Type Mapping with Single-Cell Genomics

A fundamental characteristic of animal multicellularity is the spatial coexistence of functionally specialized cell types that are all encoded by a single genome sequence. Cell type transcriptional programs are deployed and maintained by regulatory mechanisms that control the asymmetric, differential access to genomic information in each cell. This genome regulation ultimately results in specific cellular phenotypes. However, the emergence, diversity, and evolutionary dynamics of animal cell types remain almost completely unexplored beyond a few species. Single-cell genomics is emerging as a powerful tool to build comprehensive catalogs of cell types and their associated gene regulatory programs in non-traditional model species. We review the current state of sampling efforts across the animal tree of life and challenges ahead for the comparative study of cell type programs. We also discuss how the phylogenetic integration of cell atlases can lead to the development of models of cell type evolution and a phylogenetic taxonomy of cells.

The origin of animals: an ancestral reconstruction of the unicellular-to-multicellular transition

How animals evolved from a single-celled ancestor, transitioning from a unicellular lifestyle to a coordinated multicellular entity, remains a fascinating question. Key events in this transition involved the emergence of processes related to cell adhesion, cell–cell communication and gene regulation. To understand how these capacities evolved, we need to reconstruct the features of both the last common multicellular ancestor of animals and the last unicellular ancestor of animals. In this review, we summarize recent advances in the characterization of these ancestors, inferred by comparative genomic analyses between the earliest branching animals and those radiating later, and between animals and their closest unicellular relatives. We also provide an updated hypothesis regarding the transition to animal multicellularity, which was likely gradual and involved the use of gene regulatory mechanisms in the emergence of early developmental and morphogenetic plans. Finally, we discuss some new avenues of research that will complement these studies in the coming years.

Targeted brachyury degradation disrupts a highly specific autoregulatory program controlling chordoma cell identity

Chordomas are rare spinal tumors addicted to expression of the developmental transcription factor brachyury. In chordomas, brachyury is super-enhancer associated and preferentially downregulated by pharmacologic transcriptional CDK inhibition, leading to cell death. To understand the underlying basis of this sensitivity, we dissect the brachyury transcription regulatory network and compare the consequences of brachyury degradation with transcriptional CDK inhibition. Brachyury defines the chordoma super-enhancer landscape and autoregulates through binding its super-enhancer, and its locus forms a transcriptional condensate. Transcriptional CDK inhibition and brachyury degradation disrupt brachyury autoregulation, leading to loss of its transcriptional condensate and transcriptional program. Compared with transcriptional CDK inhibition, which globally downregulates transcription, leading to cell death, brachyury degradation is much more selective, inducing senescence and sensitizing cells to anti-apoptotic inhibition. These data suggest that brachyury downregulation is a core tenet of transcriptional CDK inhibition and motivates developing strategies to target brachyury and its autoregulatory feedback loop.

Brachyury defines the chordoma super-enhancer landscape

Brachyury autoregulates through a transcriptional condensate

CDK7/12/13i and brachyury degradation target the brachyury transcriptional condensate

Brachyury degradation inhibits chordoma identity genes and induces senescence

T-Box Transcription Factor 22 Is an Immune Microenvironment-Related Biomarker Associated With the BRAF V600E Mutation in Papillary Thyroid Carcinoma

Papillary thyroid cancer (PTC) is the most common malignant disease in endocrine systems. T-box transcription factor 22 (TBX22) is a phylogenetically conserved family member that has not been widely characterized in cancers. In this study, we explored the potential clinical significance and biological functions of TBX22 in PTC. Comprehensive analyses of TBX22 were based on the public databases and our local qRT-PCR cohort. We observed that TBX22 was significantly downregulated in PTC compared with normal tissues. TBX22 was associated with several clinicopathological factors in PTC. Low TBX22 expression correlated with BRAF ^V600E and TERT mutation. Functional enrichment analysis revealed that cancer-related pathways and immune progress were closely associated with TBX22 in PTC. In TBX22-low PTC, high immune infiltration levels with increased CD8⁺ T cells, natural killer, M1 macrophages, and T-regulatory cells were observed. TBX22 was negatively correlated with the activity of different steps of the anticancer immunity cycle. Functionally, overexpression of TBX22 inhibited the proliferation, invasion, and migration in PTC cells, while knocking down of TBX22 showed the opposite effects. The present findings disclose that TBX22, as an immune microenvironment-related biomarker, could be an important tumor suppresser gene and might inform the management of PTC patients better.

Towards understanding the origin of animal development

Almost all animals undergo embryonic development, going from a single-celled zygote to a complex multicellular adult. We know that the patterning and morphogenetic processes involved in development are deeply conserved within the animal kingdom. However, the origins of these developmental processes are just beginning to be unveiled. Here, we focus on how the protist lineages sister to animals are reshaping our view of animal development. Most intriguingly, many of these protistan lineages display transient multicellular structures, which are governed by similar morphogenetic and gene regulatory processes as animal development. We discuss here two potential alternative scenarios to explain the origin of animal embryonic development: either it originated concomitantly at the onset of animals or it evolved from morphogenetic processes already present in their unicellular ancestors. We propose that an integrative study of several unicellular taxa closely related to animals will allow a more refined picture of how the last common ancestor of animals underwent embryonic development.

The genetic factors of bilaterian evolution

The Cambrian explosion was a unique animal radiation ~540 million years ago that produced the full range of body plans across bilaterians. The genetic mechanisms underlying these events are unknown, leaving a fundamental question in evolutionary biology unanswered. Using large-scale comparative genomics and advanced orthology evaluation techniques, we identified 157 bilaterian-specific genes. They include the entire Nodal pathway, a key regulator of mesoderm development and left-right axis specification; components for nervous system development, including a suite of G-protein-coupled receptors that control physiology and behaviour, the Robo-Slit midline repulsion system, and the neurotrophin signalling system; a high number of zinc finger transcription factors; and novel factors that previously escaped attention. Contradicting the current view, our study reveals that genes with bilaterian origin are robustly associated with key features in extant bilaterians, suggesting a causal relationship.

Brachyury in the gastrula of basal vertebrates

The Brachyury gene encodes a transcription factor that is conserved across all animals. In non-chordate metazoans, brachyury is primarily expressed in ectoderm regions that are added to the endodermal gut during development, and often form a ring around the site of endoderm internalization in the gastrula, the blastopore. In chordates, this brachyury ring is conserved, but the gene has taken on a new role in the formation of the mesoderm. In this phylum, a novel type of mesoderm that develops into notochord and somites has been added to the ancestral lateral plate mesoderm. Brachyury contributes to a shift in cell fate from neural ectoderm to posterior notochord and somites during a major lineage segregation event that in Xenopus and in the zebrafish takes place in the early gastrula. In the absence of this brachyury function, impaired formation of posterior mesoderm indirectly affects the gastrulation movements of peak involution and convergent extension. These movements are confined to specific regions and stages, leaving open the question why brachyury expression in an extensive, coherent ring, before, during and after gastrulation, is conserved in the two species whose gastrulation modes differ considerably, and also in many other metazoan gastrulae of diverse structure.

Insights into the origin of metazoan multicellularity from predatory unicellular relatives of animals

BACKGROUND

The origin of animals from their unicellular ancestor was one of the most important events in evolutionary history, but the nature and the order of events leading up to the emergence of multicellular animals are still highly uncertain. The diversity and biology of unicellular relatives of animals have strongly informed our understanding of the transition from single-celled organisms to the multicellular Metazoa. Here, we analyze the cellular structures and complex life cycles of the novel unicellular holozoans Pigoraptor and Syssomonas (Opisthokonta), and their implications for the origin of animals.

RESULTS

Syssomonas and Pigoraptor are characterized by complex life cycles with a variety of cell types including flagellates, amoeboflagellates, amoeboid non-flagellar cells, and spherical cysts. The life cycles also include the formation of multicellular aggregations and syncytium-like structures, and an unusual diet for single-celled opisthokonts (partial cell fusion and joint sucking of a large eukaryotic prey), all of which provide new insights into the origin of multicellularity in Metazoa. Several existing models explaining the origin of multicellular animals have been put forward, but these data are interestingly consistent with one, the "synzoospore hypothesis."

CONCLUSIONS

The feeding modes of the ancestral metazoan may have been more complex than previously thought, including not only bacterial prey, but also larger eukaryotic cells and organic structures. The ability to feed on large eukaryotic prey could have been a powerful trigger in the formation and development of both aggregative (e.g., joint feeding, which also implies signaling) and clonal (e.g., hypertrophic growth followed by palintomy) multicellular stages that played important roles in the emergence of multicellular animals.

Origin and evolution of eukaryotic transcription factors

Transcription factors (TFs) have a central role in genome regulation directing gene transcription through binding specific DNA sequences. Eukaryotic genomes encode a large diversity of TF classes, each defined by unique DNA-interaction domains. Recent advances in genome sequencing and phylogenetic placement of diverse eukaryotic and archaeal species are re-defining the evolutionary history of eukaryotic TFs. The emerging view from a comparative genomics perspective is that the Last Eukaryotic Common Ancestor (LECA) had an extensive repertoire of TFs, most of which represent eukaryotic evolutionary novelties. This burst of TF innovation coincides with the emergence of genomic nuclear segregation and complex chromatin organization.

Penaeid shrimp brachyury: sequence analysis and expression during gastrulation

Gastrulation occurs by a variety of morphogenetic movements, often correlated with diverse expression of the T-box transcription factor Brachyury (Bra). Bra may be expressed in ectoderm, mesoderm, or endoderm, but its role in cell fate specification or regulation of gastrulation movements has not been studied in the development of crustaceans. Penaeid shrimp (Decapoda: Dendrobranchiata: Penaeidae) develop by complete cleavage and gastrulation by invagination to a free-swimming nauplius larva. Penaeid gastrulation diverges from other decapods and from insects, occurring early at a low cell number with the formation of a radial invagination. Toward a better understanding of gastrulation movements in penaeid shrimp, bra was identified from newly available penaeid shrimp genomes and transcriptomes of Litopenaeus vannamei, Marsupenaeus japonicus, and Penaeus monodon. Additional bra homologs were obtained from the outgroups Sicyonia ingentis (Decapoda: Dendrobranchiata: Sicyoniidae) and the caridean shrimp Caridina multidentata (Decapoda: Pleocymata). The genes encoded penaeid shrimp Bra proteins of 551-552 amino acids, containing the highly conserved T-box DNA-binding region. The N-terminal Smad1-binding domain, conserved in most animals, was absent in shrimp Bra. The R1 repressor domain was the best conserved of the C-terminal regulatory domains, which were widely divergent compared to other species. The penaeid shrimp bra gene consisted of six exons, with splice sites conserved with other phyla across the animal kingdom. Real-time qPCR and FPKM analysis showed that shrimp bra mRNA was strongly expressed during gastrulation. These findings begin to address the evolution of gastrulation in shrimp at the molecular level.

Embracing Uncertainty in Reconstructing Early Animal Evolution

The origin of animals, one of the major transitions in evolution, remains mysterious. Many key aspects of animal evolution can be reconstructed by comparing living species within a robust phylogenetic framework. However, uncertainty remains regarding the evolutionary relationships between two ancient animal lineages - sponges and ctenophores - and the remaining animal phyla. Comparative morphology and some phylogenomic analyses support the view that sponges represent the sister lineage to the rest of the animals, while other phylogenomic analyses support ctenophores, a phylum of carnivorous, gelatinous marine organisms, as the sister lineage. Here, we explore why different studies yield different answers and discuss the implications of the two alternative hypotheses for understanding the origin of animals. Reconstruction of ancient evolutionary radiations is devilishly difficult and will likely require broader sampling of sponge and ctenophore genomes, improved analytical strategies and critical analyses of the phylogenetic distribution and molecular mechanisms underlying apparently conserved traits. Rather than staking out positions in favor of the ctenophores-sister or the sponges-sister hypothesis, we submit that research programs aimed at understanding the biology of the first animals should instead embrace the uncertainty surrounding early animal evolution in their experimental designs.

Gene family innovation, conservation and loss on the animal stem lineage

Choanoflagellates, the closest living relatives of animals, can provide unique insights into the changes in gene content that preceded the origin of animals. However, only two choanoflagellate genomes are currently available, providing poor coverage of their diversity. We sequenced transcriptomes of 19 additional choanoflagellate species to produce a comprehensive reconstruction of the gains and losses that shaped the ancestral animal gene repertoire. We identified ~1944 gene families that originated on the animal stem lineage, of which only 39 are conserved across all animals in our study. In addition, ~372 gene families previously thought to be animal-specific, including Notch, Delta, and homologs of the animal Toll-like receptor genes, instead evolved prior to the animal-choanoflagellate divergence. Our findings contribute to an increasingly detailed portrait of the gene families that defined the biology of the Urmetazoan and that may underpin core features of extant animals.

All animals, from sea sponges and reef-building corals to elephants and humans, share a single common ancestor that lived over half a billion years ago. This single-celled predecessor evolved the ability to develop into a creature made up of many cells with specialized jobs. Reconstructing the steps in this evolutionary process has been difficult because the earliest animals were soft-bodied and microscopic and did not leave behind fossils that scientists can study.

Though their bodies have since disintegrated, many of the instructions for building the first animals live on in genes that were passed on to life forms that still exist. Scientists are trying to retrace those genes back to the first animal by comparing the genomes of living animals with their closest relatives, the choanoflagellates. Choanoflagellates are single-celled, colony-forming organisms that live in waters around the world. Comparisons with choanoflagellates may help scientists identify which genes were necessary to help animals evolve and diversify into so many different species. So far, 1,000 animal and two choanoflagellate genomes have been sequenced. But the gene repertoires of most species of choanoflagellates have yet to be analyzed.

Now, Richter et al. have cataloged the genes of 19 more species of choanoflagellates. This added information allowed them to recreate the likely gene set of the first animal and to identify genetic changes that occurred during animal evolution. The analyses showed that modern animals lost about a quarter of the genes present in their last common ancestor with choanoflagellates and gained an equal number of new genes. Richter et al. identified several dozen core animal genes that were gained and subsequently preserved throughout animal evolution. Many of these are necessary so that an embryo can develop properly, but the precise roles of some core genes remain a mystery. Most other genes that emerged in the first animals have been lost in at least one living animal.

The study of Richter et al. also showed that some very important genes in animals, including genes essential for early development and genes that help the immune system detect pathogens, predate animals. These key genes trace back to animals’ last common ancestor with choanoflagellates and may have evolved new roles in animals.

Transfection of Capsaspora owczarzaki, a close unicellular relative of animals

How animals emerged from their unicellular ancestor remains a major evolutionary question. New genome data from the closest unicellular relatives of animals have provided important insights into the evolution of animal multicellularity. We know that the unicellular ancestor of animals had an unexpectedly complex genetic repertoire, including many genes that are key to animal development and multicellularity. Thus, assessing the function of these genes among unicellular relatives of animals is key to understanding how they were co-opted at the onset of the Metazoa. However, such analyses have been hampered by the lack of genetic tools. Progress has been made in choanoflagellates and teretosporeans, two of the three lineages closely related to animals, whereas no tools are yet available for functional analysis in the third lineage: the filastereans. Importantly, filastereans have a striking repertoire of genes involved in transcriptional regulation and other developmental processes. Here, we describe a reliable transfection method for the filasterean Capsaspora owczarzaki. We also provide a set of constructs for visualising subcellular structures in live cells. These tools convert Capsaspora into a unique experimentally tractable organism to use to investigate the origin and evolution of animal multicellularity.

Summary: Development of genetic tools in a close unicellular relative of animals, the filasterean Capsaspora owczarzaki, will open new opportunities for functional studies to understand the transition to animal multicellularity.

Characterization of a group I Nme protein of Capsaspora owczarzaki-a close unicellular relative of animals

Nucleoside diphosphate kinases are enzymes present in all domains of life. In animals, they are called Nme or Nm23 proteins, and are divided into group I and II. Human Nme1 was the first protein identified as a metastasis suppressor. Because of its medical importance, it has been extensively studied. In spite of the large research effort, the exact mechanism of metastasis suppression remains unclear. It is unknown which of the biochemical properties or biological functions are responsible for the antimetastatic role of the mammalian Nme1. Furthermore, it is not clear at which point in the evolution of life group I Nme proteins acquired the potential to suppress metastasis, a process that is usually associated with complex animals. In this study we performed a series of tests and assays on a group I Nme protein from filasterean Capsaspora owczarzaki, a close unicellular relative of animals. The aim was to compare the protein to the well-known human Nme1 and Nme2 homologs, as well as with the homolog from a simple animal-sponge (Porifera), in order to see how the proteins changed with the transition to multicellularity, and subsequently in the evolution of complex animals. We found that premetazoan-type protein is highly similar to the homologs from sponge and human, in terms of biochemical characteristics and potential biological functions. Like the human Nme1 and Nme2, it is able to diminish the migratory potential of human cancer cells in culture.

Cas9-mediated excision of Nematostella brachyury disrupts endoderm development, pharynx formation and oral-aboral patterning

The mesoderm is a key novelty in animal evolution, although we understand little of how the mesoderm arose. brachyury, the founding member of the T-box gene family, is a key gene in chordate mesoderm development. However, the brachyury gene was present in the common ancestor of fungi and animals long before mesoderm appeared. To explore ancestral roles of brachyury prior to the evolution of definitive mesoderm, we excised the gene using CRISPR/Cas9 in the diploblastic cnidarian Nematostella vectensis Nvbrachyury is normally expressed in precursors of the pharynx, which separates endoderm from ectoderm. In knockout embryos, the pharynx does not form, embryos fail to elongate, and endoderm organization, ectodermal cell polarity and patterning along the oral-aboral axis are disrupted. Expression of many genes both inside and outside the Nvbrachyury expression domain is affected, including downregulation of Wnt genes at the oral pole. Our results point to an ancient role for brachyury in morphogenesis, cell polarity and the patterning of both ectodermal and endodermal derivatives along the primary body axis.

Novel Predators Reshape Holozoan Phylogeny and Reveal the Presence of a Two-Component Signaling System in the Ancestor of Animals

Our understanding of the origin of animals has been transformed by characterizing their most closely related, unicellular sisters: the choanoflagellates, filastereans, and ichthyosporeans. Together with animals, these lineages make up the Holozoa [1, 2]. Many traits previously considered "animal specific" were subsequently found in other holozoans [3, 4], showing that they evolved before animals, although exactly when is currently uncertain because several key relationships remain unresolved [2, 5]. Here we report the morphology and transcriptome sequencing from three novel unicellular holozoans: Pigoraptor vietnamica and Pigoraptor chileana, which are related to filastereans, and Syssomonas multiformis, which forms a new lineage with Corallochytrium in phylogenomic analyses. All three species are predatory flagellates that feed on large eukaryotic prey, and all three also appear to exhibit complex life histories with several distinct stages, including multicellular clusters. Examination of genes associated with multicellularity in animals showed that the new filastereans contain a cell-adhesion gene repertoire similar to those of other species in this group. Syssomonas multiformis possessed a smaller complement overall but does encode genes absent from the earlier-branching ichthyosporeans. Analysis of the T-box transcription factor domain showed expansion of T-box transcription factors based on combination with a non-T-box domain (a receiver domain), which has not been described outside of vertebrates. This domain and other domains we identified in all unicellular holozoans are part of the two-component signaling system that has been lost in animals, suggesting the continued use of this system in the closest relatives of animals and emphasizing the importance of studying loss of function as well as gain in major evolutionary transitions.

The origin of Metazoa: a unicellular perspective

The first animals evolved from an unknown single-celled ancestor in the Precambrian period. Recently, the identification and characterization of the genomic and cellular traits of the protists most closely related to animals have shed light on the origin of animals. Comparisons of animals with these unicellular relatives allow us to reconstruct the first evolutionary steps towards animal multicellularity. Here, we review the results of these investigations and discuss their implications for understanding the earliest stages of animal evolution, including the origin of metazoan genes and genome function.

How complexity increases in development: An analysis of the spatial-temporal dynamics of Gene expression in Ciona intestinalis

The increase in complexity in an embryo over developmental time is perhaps one of the most intuitive processes of animal development. It is also intuitive that the embryo becomes progressively compartmentalized over time and space. In spite of this intuitiveness, there are no systematic attempts to quantify how this occurs. Here, we present a quantitative analysis of the compartmentalization and spatial complexity of Ciona intestinalis over developmental time by analyzing thousands of gene expression spatial patterns from the ANISEED database. We measure compartmentalization in two ways: as the relative volume of expression of genes and as the disparity in gene expression between body parts. We also use a measure of the curvature of each gene expression pattern in 3D space. These measures show a similar increase over time, with the most dramatic change occurring from the 112-cell stage to the early tailbud stage. Combined, these measures point to a global pattern of increase in complexity in the Ciona embryo. Finally, we cluster the different regions of the embryo depending on their gene expression similarity, within and between stages. Results from this clustering analysis, which partially correspond to known fate maps, provide a global quantitative overview about differentiation and compartmentalization between body parts at each developmental stage.

Conservation and divergence of bHLH genes in the calcisponge Sycon ciliatum

Background

Basic Helix-Loop-Helix (bHLH) genes encode a large family of eukaryotic transcription factors, categorized into six high-order groups: pan-eukaryotic group B involved in regulation of cell cycle, metabolism, and development; holozoan-specific groups C and F involved in development and maintenance of homeostasis; and metazoan-specific groups A, D and E including well-studied genes, such as Atonal, Twist and Hairy, with diverse developmental roles including control of morphogenesis and specification of neurons. Current scenarios of bHLH evolution in animals are mainly based on the bHLH gene set found in the genome of demosponge Amphimedon queenslandica. In this species, the majority of the 21 identified bHLH genes belong to group B, and the single group A gene is orthologous to several neurogenic bilaterian subfamilies, including atonal and neurogenin.

Results

Given recently discovered differences in developmental toolkit components between siliceous and calcareous sponges, we have carried out genome-wide analysis of bHLH genes in Sycon ciliatum, an emerging calcisponge model. We identified 30 bHLH genes in this species, representing 12 individual families, including four group A families not found in Amphimedon, and two larger family groupings. Notably, the families represented in Sycon are only partially overlapping with those represented in Amphimedon. Developmental expression analysis of a subset of the identified genes revealed patterns consistent with deeply conserved roles, such as specification of sensory cells by Atona-related and stem cells by Myc genes.

Conclusions

Our results demonstrate independent gene loss events in demosponges and calcisponges, implying a complex bHLH toolkit in the last common metazoan ancestor.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0060-8) contains supplementary material, which is available to authorized users.

The Remarkably Diverse Family of T-Box Factors in Caenorhabditis elegans

The nematode Caenorhabditis elegans is a simple metazoan animal that is widely used as a model to understand the genetic control of development. The completely sequenced C. elegans genome contains 22 T-box genes, and they encode factors that show remarkable diversity in sequence, DNA-binding specificity, and function. Only three of the C. elegans T-box factors can be grouped into the conserved subfamilies found in other organisms, while the remaining factors are significantly diverged and unlike those in most other animals. While some of the C. elegans factors can bind canonical T-box binding elements, others bind and regulate target gene expression through distinct sequences. The nine genetically characterized T-box factors have varied functions in development and morphogenesis of muscle, hypodermal tissues, and neurons, as well as in early blastomere fate specification, cell migration, apoptosis, and sex determination, but the functions of most of the C. elegans T-box factors have not yet been extensively characterized. Like T-box factors in other animals, interaction with a Groucho-family corepressor and posttranslational SUMOylation have been shown to affect C. elegans T-box factor activity, and it is likely that additional mechanisms affecting T-box factor activity will be discovered using the effective genetic approaches in this organism.

T-Box Genes in the Kidney and Urinary Tract

T-box (Tbx) genes encode an ancient group of transcription factors that play important roles in patterning, specification, proliferation, and differentiation programs in vertebrate organogenesis. This is testified by severe organ malformation syndromes in mice homozygous for engineered null alleles of specific T-box genes and by the large number of human inherited organ-specific diseases that have been linked to mutations in these genes. One of the organ systems that has not been associated with loss of specific T-box gene function in human disease for long is the excretory system. However, this has changed with the finding that mutations in TBX18, a member of a vertebrate-specific subgroup within the Tbx1-subfamily of T-box transcription factor genes, cause congenital anomalies of the kidney and urinary tract, predominantly hydroureter and ureteropelvic junction obstruction. Gene expression analyses, loss-of-function studies, and lineage tracing in the mouse suggest a primary role for this transcription factor in specifying the ureteric mesenchyme in the common anlage of the kidney, the ureter, and the bladder. We review the function of Tbx18 in ureterogenesis and discuss the body of evidence that Tbx18 and other members of the T-box gene family, namely, Tbx1, Tbx2, Tbx3, and Tbx20, play additional roles in development and homeostasis of other components of the excretory system in vertebrates.

The Unicellular Ancestry of Groucho-Mediated Repression and the Origins of Metazoan Transcription Factors

Groucho is a co-repressor that interacts with many transcription factors playing a crucial role in animal development. The evolutionary origins of Groucho are not clear. It is generally regarded as being a distinct animal-specific protein, although with similarities to the yeast Tup-like proteins. Here, it is shown that Groucho has true orthologs in unicellular relatives of animals. Based on their phylogenetic distribution, and an analysis of ligand-binding residues, these genes are unlikely to be orthologs of the fungal Tup-like genes. By identifying conserved candidate Groucho interaction motifs (GIMs) in nonmetazoan transcription factors, it is demonstrated that the details of molecular interactions between Groucho and transcription factors are likely to have been established prior to the origin of animals, but that the association of GIMs with many transcription factor types can be regarded as a metazoan innovation.

Sponges as models to study emergence of complex animals

The emergence of complex animal life forms remains poorly understood despite substantial interest and research in this area. To be informative, the ideal models to study transitions from single-cell organisms to the first animals and then to mammalian-level complexity should be phylogenetically strategically placed and retain ancestral characters. Sponges (Porifera) are likely to be the earliest branching animal phylum. When analysed from morphological, genomic and developmental perspectives, sponges appear to combine features of single-cell eukaryotic organisms and the complex multicellular animals (Eumetazoa). Intriguingly, homologues of components of the eumetazoan regulatory networks specifying the endoderm, the germ-cells and stem cells and (neuro) sensory cells are expressed in sponge choanocytes, archaeocytes and larval sensory cells. Studies using sponges as model systems are already bringing insights into animal evolution, and have opened avenues to further research benefitting from the recent spectacular expansion of genomic technologies.

The Dynamic Regulatory Genome of Capsaspora and the Origin of Animal Multicellularity

The unicellular ancestor of animals had a complex repertoire of genes linked to multicellular processes. This suggests that changes in the regulatory genome, rather than in gene innovation, were key to the origin of animals. Here, we carry out multiple functional genomic assays in Capsaspora owczarzaki, the unicellular relative of animals with the largest known gene repertoire for transcriptional regulation. We show that changing chromatin states, differential lincRNA expression, and dynamic cis-regulatory sites are associated with life cycle transitions in Capsaspora. Moreover, we demonstrate conservation of animal developmental transcription-factor networks and extensive network interconnection in this premetazoan organism. In contrast, however, Capsaspora lacks animal promoter types, and its regulatory sites are small, proximal, and lack signatures of animal enhancers. Overall, our results indicate that the emergence of animal multicellularity was linked to a major shift in genome cis-regulatory complexity, most notably the appearance of distal enhancer regulation.

T-bet as a key regulator of mucosal immunity

Initially understood to be a key regulator of interferon-γ-producing helper T cells, our knowledge of T-bet's functional roles has expanded to encompass a growing range of cellular lineages. In addition to regulating other interferon-γ-producing adaptive immune cells, it is now clear that T-bet plays a fundamental role in the regulation of innate immune responses across mucosal surfaces. This homeostatic role is demonstrated by the spontaneous colitis that occurs when T-bet is deleted from innate immune cells in RAG(-/-) mice. Using this model as a focal point, we review our understanding of T-bet's regulation of adaptive and innate immune systems, focusing particularly on mucosal populations including innate lymphoid cells, dendritic cells and intraepithelial lymphocytes. With the increasingly diverse effects of T-bet on different lineages, the classical binding-centric paradigm of T-bet's molecular functionality has increasingly struggled to account for the versatility of T-bet's biological effects. Recent recognition of the synergistic interactions between T-bet and other canonical transcription factors has led to a co-operative paradigm that has provided greater explanatory power. Synthesizing insights from ChIP-seq and comparative biology, we expand the co-operative paradigm further and suggest a network approach as a powerful way to understand and model T-bet's diverse functionality.

Bilaterian-like promoters in the highly compact Amphimedon queenslandica genome

The regulatory systems underlying animal development must have evolved prior to the emergence of eumetazoans (cnidarians and bilaterians). Although representatives of earlier-branching animals - sponges ctenophores and placozoans - possess most of the developmental transcription factor families present in eumetazoans, the DNA regulatory elements that these transcription factors target remain uncharted. Here we characterise the core promoter sequences, U1 snRNP-binding sites (5' splice sites; 5'SSs) and polyadenylation sites (PASs) in the sponge Amphimedon queenslandica. Similar to unicellular opisthokonts, Amphimedon's genes are tightly packed in the genome and have small introns. In contrast, its genes possess metazoan-like core promoters populated with binding motifs previously deemed to be specific to vertebrates, including Nrf-1 and Krüppel-like elements. Also as in vertebrates, Amphimedon's PASs and 5'SSs are depleted downstream and upstream of transcription start sites, respectively, consistent with non-elongating transcripts being short-lived; PASs and 5'SSs are more evenly distributed in bidirectional promoters in Amphimedon. The presence of bilaterian-like regulatory DNAs in sponges is consistent with these being early and essential innovations of the metazoan gene regulatory repertoire.

Evolution of Prdm Genes in Animals: Insights from Comparative Genomics

Prdm genes encode transcription factors with a subtype of SET domain known as the PRDF1-RIZ (PR) homology domain and a variable number of zinc finger motifs. These genes are involved in a wide variety of functions during animal development. As most Prdm genes have been studied in vertebrates, especially in mice, little is known about the evolution of this gene family. We searched for Prdm genes in the fully sequenced genomes of 93 different species representative of all the main metazoan lineages. A total of 976 Prdm genes were identified in these species. The number of Prdm genes per species ranges from 2 to 19. To better understand how the Prdm gene family has evolved in metazoans, we performed phylogenetic analyses using this large set of identified Prdm genes. These analyses allowed us to define 14 different subfamilies of Prdm genes and to establish, through ancestral state reconstruction, that 11 of them are ancestral to bilaterian animals. Three additional subfamilies were acquired during early vertebrate evolution (Prdm5, Prdm11, and Prdm17). Several gene duplication and gene loss events were identified and mapped onto the metazoan phylogenetic tree. By studying a large number of nonmetazoan genomes, we confirmed that Prdm genes likely constitute a metazoan-specific gene family. Our data also suggest that Prdm genes originated before the diversification of animals through the association of a single ancestral SET domain encoding gene with one or several zinc finger encoding genes.