Identification of novel genes affecting mesoderm formation and morphogenesis through an enhanced large scale functional screen in Xenopus

Atypical RhoUV GTPases in development and disease

RhoU and RhoV are members of the Rho family of small GTPases that comprise their own subfamily. RhoUV GTPases are classified as atypical due to the kinetics of their GTP/GDP binding cycles. They also possess unique N- and C-termini that regulate their subcellular localization and activity. RhoU and RhoV have been linked to cytoskeletal regulation, cell adhesion, and cell migration. They each exhibit distinct expression patterns during embryonic development and diseases such as cancer metastasis, suggesting they have specialized functions. In this review, we will discuss the known functions of RhoU and RhoV, with a focus on their roles in early development, organogenesis, and disease.

Spatial transcriptomics reveals a cnidarian segment polarity program in Nematostella vectensis

During early animal evolution, the emergence of axially polarized segments was central to the diversification of complex bilaterian body plans. Nevertheless, precisely how and when segment polarity pathways arose remains obscure. Here, we demonstrate the molecular basis for segment polarization in developing larvae of the sea anemone Nematostella vectensis. Utilizing spatial transcriptomics, we first constructed a 3D gene expression atlas of developing larval segments. Capitalizing on accurate in silico predictions, we identified Lbx and Uncx, conserved homeodomain-containing genes that occupy opposing subsegmental domains under the control of both bone morphogenetic protein (BMP) signaling and the Hox-Gbx cascade. Functionally, Lbx mutagenesis eliminated all molecular evidence of segment polarization at the larval stage and caused an aberrant mirror-symmetric pattern of retractor muscles (RMs) in primary polyps. These results demonstrate the molecular basis for segment polarity in a non-bilaterian animal, suggesting that polarized metameric structures were present in the Cnidaria-Bilateria common ancestor over 600 million years ago.

Spatial transcriptomics reveals a conserved segment polarity program that governs muscle patterning in Nematostella vectensis

During early animal evolution, the emergence of axially-polarized segments was central to the diversification of complex bilaterian body plans. Nevertheless, precisely how and when segment polarity pathways arose remains obscure. Here we demonstrate the molecular basis for segment polarization in developing larvae of the pre-bilaterian sea anemone Nematostella vectensis . Utilizing spatial transcriptomics, we first constructed a 3-D gene expression atlas of developing larval segments. Capitalizing on accurate in silico predictions, we identified Lbx and Uncx, conserved homeodomain-containing genes that occupy opposing subsegmental domains under the control of both BMP signaling and the Hox-Gbx cascade. Functionally, Lbx mutagenesis eliminated all molecular evidence of segment polarization at larval stage and caused an aberrant mirror-symmetric pattern of retractor muscles in primary polyps. These results demonstrate the molecular basis for segment polarity in a pre-bilaterian animal, suggesting that polarized metameric structures were present in the Cnidaria-Bilateria common ancestor over 600 million years ago.

Highlights

Nematostella endomesodermal tissue forms metameric segments and displays a transcriptomic profile similar to that observed in bilaterian mesoderm Construction of a comprehensive 3-D gene expression atlas enables systematic dissection of segmental identity in endomesoderm Lbx and Uncx , two conserved homeobox-containing genes, establish segment polarity in Nematostella The Cnidarian-Bilaterian common ancestor likely possessed the genetic toolkit to generate polarized metameric structures.

NUCKS1 is a highly modified, chromatin-associated protein involved in a diverse set of biological and pathophysiological processes

The Nuclear Casein and Cyclin-dependent Kinase Substrate 1 (NUCKS1) protein is highly conserved in vertebrates, predominantly localized to the nucleus and one of the most heavily modified proteins in the human proteome. NUCKS1 expression is high in stem cells and the brain, developmentally regulated in mice and associated with several diverse malignancies in humans, including cancer, metabolic syndrome and Parkinson's disease. NUCKS1 function has been linked to modulating chromatin architecture and transcription, DNA repair and cell cycle regulation. In this review, we summarize and discuss the published information on NUCKS1 and highlight the questions that remain to be addressed to better understand the complex biology of this multifaceted protein.

The Xenopus laevis teratogenesis assay for developmental toxicity of phthalate plasticizers and alternatives

Contamination of phthalate ester plasticizers threatens the wildlife as well as human health. To evaluate the developmental toxicity of commonly used phthalate esters and emerging alternatives, the frog embryo teratogenesis assay-Xenopus (FETAX) was conducted for dibutyl-phthalate (DBP), benzyl-butyl-phthalate (BBP), dioctyl-terephthalate (DOTP), di(2-propylheptyl)-phthalate (DPHP), diisononyl-phthalate (DINP), diisodecyl-phthalate (DIDP), diethyl hexyl cyclohexane (DEHCH), and diisononyl-cyclohexane-1,2-dicarboxylate (DINCH). The 96-hrs LC₅₀ for DBP, BBP, DOTP, DIDP, DINCH, DINP, DPHP, and DEHCH were 18.3, 20.1, 588.7, 718.0, 837.5, 859.3, 899.0, and 899.0 mg/L, respectively. The 96-hrs EC₅₀ of developmental abnormality of DBP, BBP, DPHP, DOTP, DINP, DEHCH, DINCH, and DIDP were 7.5, 18.2, 645.1, 653.6, 664.4, 745.6, 813.7, and 944.5 mg/L, respectively. The lowest observed effective concentration for embryonic survival, malformation, and growth was DINP, DBP, BBP, DIDP, DPHP, DINCH, DEHCH, and DOTP in increasing order. In tadpoles, DBP, BBP, DEHCH, DINP, and DIDP caused inositol-requiring enzyme 1 or protein kinase R-like endoplasmic reticulum kinase pathway endoplasmic reticulum stress (ERS) in order, and BBP, DBP, DOTP, DPHP, DINP, and DIDP caused long term ERS-related apoptosis or mitochondrial apoptosis in order. Together, in Xenopus embryos, the developmental toxicity and the cellular stress-inducing potential of tested plasticizers were DEHCH, DINCH, DPHP, DIDP, DINP, DOTP, BBP, and DBP in increasing order. In consideration of public as well as environmental health this information would be helpful for industrial choice of phthalate ester plasticizers and their alternatives.

Comprehensive expression analysis for the core cell cycle regulators in the chicken embryo reveals novel tissue-specific synexpression groups and similarities and differences with expression in mouse, frog and zebrafish

The core cell cycle machinery is conserved from yeast to humans, and hence it is assumed that all vertebrates share the same set of players. Yet during vertebrate evolution, the genome was duplicated twice, followed by a further genome duplication in teleost fish. Thereafter, distinct genes were retained in different vertebrate lineages; some individual gene duplications also occurred. To which extent these diversifying tendencies were compensated by retaining the same expression patterns across homologous genes is not known. This study for the first time undertook a comprehensive expression analysis for the core cell cycle regulators in the chicken, focusing in on early neurula and pharyngula stages of development, with the latter representing the vertebrate phylotypic stage. We also compared our data with published data for the mouse, Xenopus and zebrafish, the other established vertebrate models. Our work shows that, while many genes are expressed widely, some are upregulated or specifically expressed in defined tissues of the chicken embryo, forming novel synexpression groups with markers for distinct developmental pathways. Moreover, we found that in the neural tube and in the somite, mRNAs of some of the genes investigated accumulate in a specific subcellular localisation, pointing at a novel link between the site of mRNA translation, cell cycle control and interkinetic nuclear movements. Finally, we show that expression patterns of orthologous genes may differ in the four vertebrate models. Thus, for any study investigating cell proliferation, cell differentiation, tissue regeneration, stem cell behaviour and cancer/cancer therapy, it has to be carefully examined which of the observed effects are due to the specific model organism used, and which can be generalised.

Tadpole tail regeneration in Xenopus

Some organisms have a remarkable ability to heal wounds without scars and to regenerate complex tissues following injury. By gaining a more complete understanding of the biological mechanisms that promote scar-free healing and tissue regeneration, it is hoped that novel treatments that can enhance the healing and regenerative capacity of human patients can be found. In the present article, we briefly examine the genetic, molecular and cellular mechanisms underlying the regeneration of the Xenopus tadpole tail.

Engineering Xenopus embryos for phenotypic drug discovery screening

Many rare human inherited diseases remain untreatable despite the fact that the disease causing genes are known and adequate mouse disease models have been developed. In vivo phenotypic drug screening relies on isolating drug candidates by their ability to produce a desired therapeutic phenotype in whole organisms. Embryos of zebrafish and Xenopus frogs are abundant, small and free-living. They can be easily arrayed in multi-well dishes and treated with small organic molecules. With the development of novel genome modification tools, such a zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and CRISPR/Cas, it is now possible to efficiently engineer non-mammalian models of inherited human diseases. Here, we will review the rapid progress made in adapting these novel genome editing tools to Xenopus. The advantages of Xenopus embryos as in vivo models to study human inherited diseases will be presented and their utility for drug discovery screening will be discussed. Being a tetrapod, Xenopus complements zebrafish as an indispensable non-mammalian animal model for the study of human disease pathologies and the discovery of novel therapeutics for inherited diseases.

A functional genome-wide in vivo screen identifies new regulators of signalling pathways during early Xenopus embryogenesis

Embryonic development requires exquisite regulation of several essential processes, such as patterning of tissues and organs, cell fate decisions, and morphogenesis. Intriguingly, these diverse processes are controlled by only a handful of signalling pathways, and mis-regulation in one or more of these pathways may result in a variety of congenital defects and diseases. Consequently, investigating how these signalling pathways are regulated at the molecular level is essential to understanding the mechanisms underlying vertebrate embryogenesis, as well as developing treatments for human diseases. Here, we designed and performed a large-scale gain-of-function screen in Xenopus embryos aimed at identifying new regulators of MAPK/Erk, PI3K/Akt, BMP, and TGF-β/Nodal signalling pathways. Our gain-of-function screen is based on the identification of gene products that alter the phosphorylation state of key signalling molecules, which report the activation state of the pathways. In total, we have identified 20 new molecules that regulate the activity of one or more signalling pathways during early Xenopus development. This is the first time that such a functional screen has been performed, and the findings pave the way toward a more comprehensive understanding of the molecular mechanisms regulating the activity of important signalling pathways under normal and pathological conditions.

Apoptosis of limb innervating motor neurons and erosion of motor pool identity upon lineage specific dicer inactivation

Diversification of mammalian spinal motor neurons into hundreds of subtypes is critical for the maintenance of body posture and coordination of complex movements. Motor neuron differentiation is controlled by extrinsic signals that regulate intrinsic genetic programs specifying and consolidating motor neuron subtype identity. While transcription factors have been recognized as principal regulators of the intrinsic program, the role of posttranscriptional regulations has not been systematically tested. MicroRNAs produced by Dicer mediated cleavage of RNA hairpins contribute to gene regulation by posttranscriptional silencing. Here we used Olig2-cre conditional deletion of Dicer gene in motor neuron progenitors to examine effects of miRNA biogenesis disruption on postmitotic spinal motor neurons. We report that despite the initial increase in the number of motor neuron progenitors, disruption of Dicer function results in a loss of many limb- and sympathetic ganglia-innervating spinal motor neurons. Furthermore, it leads to defects in motor pool identity specification. Thus, our results indicate that miRNAs are an integral part of the genetic program controlling motor neuron survival and acquisition of subtype specific properties.

From expression cloning to gene modeling: the development of Xenopus gene sequence resources

The Xenopus community has made concerted efforts over the last 10–12 years systematically to improve the available sequence information for this amphibian model organism ideally suited to the study of early development in vertebrates. Here I review progress in the collection of both sequence data and physical clone reagents for protein coding genes. I conclude that we have cDNA sequences for around 50% and full-length clones for about 35% of the genes in Xenopus tropicalis, and similar numbers but a smaller proportion for Xenopus laevis. In addition, I demonstrate that the gaps in the current genome assembly create problems for the computational elucidation of gene sequences, and suggest some ways to ameliorate the effects of this. genesis 50:143–154, 2012. © 2012 Wiley Periodicals, Inc.

Xenopus research: metamorphosed by genetics and genomics

Research using Xenopus takes advantage of large, abundant eggs and readily manipulated embryos in addition to conserved cellular, developmental and genomic organization with mammals. Research on Xenopus has defined key principles of gene regulation and signal transduction, embryonic induction, morphogenesis and patterning as well as cell cycle regulation. Genomic and genetic advances in this system, including the development of Xenopus tropicalis as a genetically tractable complement to the widely used Xenopus laevis, capitalize on the classical strengths and wealth of achievements. These attributes provide the tools to tackle the complex biological problems of the new century, including cellular reprogramming, organogenesis, regeneration, gene regulatory networks and protein interactions controlling growth and development, all of which provide insights into a multitude of human diseases and their potential treatments.

A functional screen for genes involved in Xenopus pronephros development

In Xenopus, the pronephros is the functional larval kidney and consists of two identifiable components; the glomus, the pronephric tubules, which can be divided into four separate segments, based on marker gene expression. The simplicity of this organ, coupled with the fact that it displays the same basic organization and function as more complex mesonephros and metanephros, makes this an attractive model to study vertebrate kidney formation. In this study, we have performed a functional screen specifically to identify genes involved in pronephros development in Xenopus. Gain-of-function screens are performed by injecting mRNA pools made from a non-redundant X. tropicalis full-length plasmid cDNA library into X. laevis eggs, followed by sib-selection to identify the single clone that caused abnormal phenotypes in the pronephros. Out of 768 egg and gastrula stage cDNA clones, 31 genes, approximately 4% of the screened clones, affected pronephric marker expression examined by whole mount in situ hybridization or antibody staining. Most of the positive clones had clear expression patterns in pronephros and predicted/established functions highly likely to be involved in developmental processes. In order to carry out a more detailed study, we selected Sox7, Cpeb3, P53csv, Mecr and Dnajc15, which had highly specific expression patterns in the pronephric region. The over-expression of these five selected clones indicated that they caused pronephric abnormalities with different temporal and spatial effects. These results suggest that our strategy to identify novel genes involved in pronephros development was highly successful, and that this strategy is effective for the identification of novel genes involved in late developmental events.

Studying Wnt signaling in Xenopus

Xenopus is an established and powerful model system for the study of Wnt signaling in vertebrates. Above all, the relatively large size of the embryos enables microinjection experiments, which have led to key discoveries not only about the functional role of Wnt signaling in vertebrate embryos, but also about the molecular mechanisms of Wnt signaling in vertebrate cells. A major advantage of the Xenopus model is the ability to obtain large numbers of embryos, which develop relatively rapidly and which can be studied in natural separation from sentient adult parental animals. In order to obtain Xenopus embryos, ovulation in females is induced with a simple hormone injection, the eggs collected and fertilized with sperm from males. The Xenopus model system has been further strengthened by recent advances such as morpholino technology and efficient transgenic methods, as well as the development of Xenopus tropicalis as a diploid genetic model system with a shorter generation time and a genome similar to higher vertebrates.

Novel genes involved in Ciona intestinalis embryogenesis: characterization of gene knockdown embryos

The sequenced genome of the urochordate ascidian Ciona intestinalis contains nearly 2,500 genes that have vertebrate homologues, but their functions are as yet unknown. To identify novel genes involved in early chordates embryogenesis, we previously screened 200 Ciona genes by knockdown experiments using specific morpholino oligonucleotides and found that suppression of the translation of 40 genes caused embryonic defects (Yamada et al. [2003] Development 130:6485-6495). We have since examined an additional 304 genes, that is, screening 504 genes overall, and a total of 111 genes showed morphological defects when gene function was suppressed. We further examined the role of these genes in the differentiation of six major tissues of the embryo: endoderm, muscle, epidermis, neural tissue, mesenchyme, and notochord. Based on the similarity of phenotypes of gene knockdown embryos, genes were categorized into several groups, with the suggestion that the genes within a given group are involved in similar developmental processes. For example, five were shown to be novel genes that are likely involved in beta-catenin-mediated endoderm formation. The type of large-scale screening used is, therefore, a powerful approach to identify novel genes with significant developmental functions, the details of which will be determined in future studies.

Maintenance of motor neuron progenitors in Xenopus requires a novel localized cyclin

The ventral spinal cord contains a pool of motor neuron progenitors (pMNs), which sequentially generate motor neurons and oligodendrocytes in the embryo. The mechanisms responsible for the maintenance of pMNs are not clearly understood. We have identified a novel cyclin, cyclin Dx (ccndx), which is specifically expressed in pMNs in Xenopus. Here, we show that inhibition of ccndx causes paralysis in embryos. Furthermore, we show that maintenance of pMNs requires ccndx function. In addition, inhibition of ccndx results in the specific loss of differentiated motor neurons. However, the expression of interneuron or sensory neuron markers is unaffected in these embryos, suggesting that the role of ccndx is specifically to maintain pMNs. Thus, we have identified, for the first time, a tissue-specific cell-cycle regulator that is essential for the maintenance of a pool of neural progenitors in the vertebrate spinal cord.

Defining synphenotype groups in Xenopus tropicalis by use of antisense morpholino oligonucleotides

To identify novel genes involved in early development, and as proof-of-principle of a large-scale reverse genetics approach in a vertebrate embryo, we have carried out an antisense morpholino oligonucleotide (MO) screen in Xenopus tropicalis, in the course of which we have targeted 202 genes expressed during gastrula stages. MOs were designed to complement sequence between −80 and +25 bases of the initiating AUG codons of the target mRNAs, and the specificities of many were tested by (i) designing different non-overlapping MOs directed against the same mRNA, (ii) injecting MOs differing in five bases, and (iii) performing “rescue” experiments. About 65% of the MOs caused X. tropicalis embryos to develop abnormally (59% of those targeted against novel genes), and we have divided the genes into “synphenotype groups,” members of which cause similar loss-of-function phenotypes and that may function in the same developmental pathways. Analysis of the expression patterns of the 202 genes indicates that members of a synphenotype group are not necessarily members of the same synexpression group. This screen provides new insights into early vertebrate development and paves the way for a more comprehensive MO-based analysis of gene function in X. tropicalis.

Genome sequencing projects have provided remarkable insights into the expression and regulation of many genes. For some species, such as the invertebrates Caenorhabditis elegans and Drosophila melanogaster, it has been possible to assign functions to these genes on a genome-wide scale. For the vertebrates, similar efforts are being made in mouse and zebrafish, but work in the former species is expensive and slow, and the zebrafish experienced a whole genome duplication event, so that some genes may have retained redundant functions. Here, this study uses antisense morpholino oligonucleotides (MOs) to show that the diploid amphibian Xenopus tropicalis provides a powerful alternative species. The authors have designed MOs to target sequences around the initiating AUG codons of 202 genes expressed during early development and confirmed that these function in a specific manner. About 65% of the MOs caused embryos to develop abnormally, and the authors have divided the genes into “synphenotype groups,” members of which cause similar loss-of-function phenotypes. Expression pattern analysis indicates that members of a synphenotype group are not necessarily members of the same synexpression group. This screen provides new insights into vertebrate development and paves the way for a comprehensive MO-based analysis of gene function in X. tropicalis.

Mesoderm induction: from caps to chips

Vertebrate mesoderm induction is one of the classical problems in developmental biology. Various developmental biology approaches, particularly in Xenopus and zebrafish, have identified many of the key factors that are involved in this process and have provided major insights into how these factors interact as part of a signalling and transcription-factor network. These data are beginning to be refined by high-throughput approaches such as microarray assays. Future challenges include understanding how the prospective mesodermal cells integrate the various signals they receive and how they resolve this information to regulate their morphogenetic behaviours and cell-fate decisions.

Genetic and genomic prospects for Xenopus tropicalis research

Research using Xenopus laevis has made enormous contributions to our understanding of vertebrate development, control of the eukaryotic cell cycle and the cytoskeleton. One limitation, however, has been the lack of systematic genetic studies in Xenopus to complement molecular and cell biological investigations. Work with the closely related diploid frog Xenopus tropicalis is beginning to address this limitation. Here, we review the resources that will make genetic studies using X. tropicalis a reality.

Xenomics

Xenopus genomics, or Xenomics for short, is coming of age. Indeed, biological insight into processes such as growth factor signaling and patterning of the early embryo is now being gained by combining the value of Xenopus as a model organism for cell and developmental biology with genomic approaches. In this review I address these recent advances and explore future possibilities gained from combining this powerful experimental system with genomic approaches, as well as how our quest to understand basic biological principles will be greatly facilitated though the marriage of Xenopus and genomics.

Global gene expression profiling and cluster analysis in Xenopus laevis

We have undertaken a large-scale microarray gene expression analysis using cDNAs corresponding to 21,000 Xenopus laevis ESTs. mRNAs from 37 samples, including embryos and adult organs, were profiled. Cluster analysis of embryos of different stages was carried out and revealed expected affinities between gastrulae and neurulae, as well as between advanced neurulae and tadpoles, while egg and feeding larvae were clearly separated. Cluster analysis of adult organs showed some unexpected tissue-relatedness, e.g. kidney is more related to endodermal than to mesodermal tissues and the brain is separated from other neuroectodermal derivatives. Cluster analysis of genes revealed major phases of co-ordinate gene expression between egg and adult stages. During the maternal-early embryonic phase, genes maintaining a rapidly dividing cell state are predominantly expressed (cell cycle regulators, chromatin proteins). Genes involved in protein biosynthesis are progressively induced from mid-embryogenesis onwards. The larval-adult phase is characterised by expression of genes involved in metabolism and terminal differentiation. Thirteen potential synexpression groups were identified, which encompass components of diverse molecular processes or supra-molecular structures, including chromatin, RNA processing and nucleolar function, cell cycle, respiratory chain/Krebs cycle, protein biosynthesis, endoplasmic reticulum, vesicle transport, synaptic vesicle, microtubule, intermediate filament, epithelial proteins and collagen. Data filtering identified genes with potential stage-, region- and organ-specific expression. The dataset was assembled in the iChip microarray database, , which allows user-defined queries. The study provides insights into the higher order of vertebrate gene expression, identifies synexpression groups and marker genes, and makes predictions for the biological role of numerous uncharacterized genes.

Expression cloning screening of a unique and full-length set of cDNA clones is an efficient method for identifying genes involved in Xenopus neurogenesis

Functional screens, where a large numbers of cDNA clones are assayed for certain biological activity, are a useful tool in elucidating gene function. In Xenopus, gain of function screens are performed by pool screening, whereby RNA transcribed in vitro from groups of cDNA clones, ranging from thousands to a hundred, are injected into early embryos. Once an activity is detected in a pool, the active clone is identified by sib-selection. Such screens are intrinsically biased towards potent genes, whose RNA is active at low quantities. To improve the sensitivity and efficiency of a gain of function screen we have bioinformatically processed an arrayed and EST sequenced set of 100,000 gastrula and neurula cDNA clones, to create a unique and full-length set of approximately 2500 clones. Reducing the redundancy and excluding truncated clones from the starting clone set reduced the total number of clones to be screened, in turn allowing us to reduce the pool size to just eight clones per pool. We report that the efficiency of screening this clone set is five-fold higher compared to a redundant set derived from the same libraries. We have screened 960 cDNA clones from this set, for genes that are involved in neurogenesis. We describe the overexpression phenotypes of 18 single clones, the majority of which show a previously uncharacterised phenotype and some of which are completely novel. In situ hybridisation analysis shows that a large number of these genes are specifically expressed in neural tissue. These results demonstrate the effectiveness of a unique full-length set of cDNA clones for uncovering players in a developmental pathway.

An atlas of differential gene expression during early Xenopus embryogenesis

We have carried out a large-scale, semi-automated whole-mount in situ hybridization screen of 8369 cDNA clones in Xenopus laevis embryos. We confirm that differential gene expression is prevalent during embryogenesis since 24% of the clones are expressed non-ubiquitously and 8% are organ or cell type specific marker genes. Sequence analysis and clustering yielded 723 unique genes displaying a differential expression pattern. Of these, 18% were already described in Xenopus, 47% have homologs and 35% are lacking significant sequence similarity in databases. Many of them encode known developmental regulators. We classified 363 of the 723 genes for which a Gene Ontology annotation for molecular function could be attributed and found 'DNA binding' and 'enzyme' the most represented terms. The most common protein domains encoded in these embryonic, differentially expressed genes are the homeobox and RNA Recognition Motif (RRM). Fifty-nine putative orthologs of human disease genes, and 254 organ or cell specific marker genes were identified. Markers were found for nasal placode and archenteron roof, organs for which a specific marker was previously unavailable. Markers were also found for novel subdomains of various other organs. The tissues for which most markers were found are muscle and epidermis. Expression of cell cycle regulators fell in two classes, containing proliferation-promoting and anti-proliferative genes, respectively. We identified 66 new members of the BMP4, chromatin, endoplasmic reticulum, and karyopherin synexpression groups, thus providing a first glimpse of their probable cellular roles. Cluster analysis of tissues to measure tissue relatedness yielded some unorthodox affinities besides expectable lineage relationships. In conclusion, this study represents an atlas of gene expression patterns, which reveals embryonic regionalization, provides novel marker genes, and makes predictions about the functional role of unknown genes.