Assembly of continuous high-resolution draft genome sequence of Hemicentrotus pulcherrimus using long-read sequencing.

The update of the draft genome assembly of sea urchin, Hemicentrotus pulcherrimus, which is widely studied in East Asia as a model organism of early development, was performed using Oxford nanopore long-read sequencing. The updated assembly provided ~600-Mb genome sequences divided into 2,163 contigs with N50 = 516 kb. BUSCO completeness score and transcriptome model mapping ratio (TMMR) of the present assembly were obtained as 96.5% and 77.8%, respectively. These results were more continuous with higher resolution than those by the previous version of H. pulcherrimus draft genome, HpulGenome_v1, where the number of scaffolds = 16,251 with a total of ~100 Mb, N50 = 143 kb, BUSCO completeness score = 86.1%, and TMMR = 55.4%. The obtained genome contained 36,055 gene models that were consistent with those in other echinoderms. Additionally, two tandem repeat sequences of early histone gene locus containing 47 copies and 34 copies of all histone genes, and 185 of the homologous sequences of the interspecifically conserved region of the Ars insulator, ArsInsC, were obtained. These results provide further advance for genome-wide research of development, gene regulation, and intranuclear structural dynamics of multicellular organisms using H. pulcherrimus.

TrBase: A genome and transcriptome database of Temnopleurus reevesii.

Sea urchins have a long history as model organisms in biology, but their use in genetics is limited because of their long breeding cycle. In sea urchin genetics, genome editing technology was first established in Hemicentrotus pulcherrimus, whose genome has already been published. However, because this species also has a long breeding cycle, new model sea urchins that are more suitable for genetics have been sought. Here, we report a draft genome of another Western Pacific species, Temnopleurus reevesii, which we established as a new model sea urchin recently since this species has a comparable developmental process to other model sea urchins but a short breeding cycle of approximately half a year. The genome of T. reevesii was assembled into 28,742 scaffold sequences with an N50 length of 67.6 kb and an estimated genome size of 905.9 Mb. In the assembled genome, 27,064 genes were identified, 23,624 of which were expressed in at least one of the seven developmental stages. To provide genetic information, we constructed the genome database TrBase (https://cell-innovation.nig.ac.jp/Tree/). We also constructed the Western Pacific Sea Urchin Genome Database (WestPac-SUGDB) (https://cell-innovation.nig.ac.jp/WPAC/) with the aim of establishing a portal site for genetic information on sea urchins in the West Pacific. This site contains genomic information on two species, T. reevesii and H. pulcherrimus, and is equipped with homology search programs for comparing the two datasets. Therefore, TrBase and WestPac-SUGDB are expected to contribute not only to genetic research using sea urchins but also to comparative genomics and evolutionary research.

Human disease-associated extracellular matrix orthologs ECM3 and QBRICK regulate primary mesenchymal cell migration in sea urchin embryos.

Sea urchin embryos have been one of model organisms to investigate cellular behaviors because of their simple cell composition and transparent body. They also give us an opportunity to investigate molecular functions of human proteins of interest that are conserved in sea urchin. Here we report that human disease-associated extracellular matrix orthologues ECM3 and QBRICK are necessary for mesenchymal cell migration during sea urchin embryogenesis. Immunofluorescence has visualized the colocalization of QBRICK and ECM3 on both apical and basal surface of ectoderm. On the basal surface, QBRICK and ECM3 constitute together a mesh-like fibrillar structure along the blastocoel wall. When the expression of ECM3 was knocked down by antisense-morpholino oligonucleotides, the ECM3-QBRICK fibrillar structure completely disappeared. When QBRICK was knocked down, the ECM3 was still present, but the basally localized fibers became fragmented. The ingression and migration of primary mesenchymal cells were not critically affected, but their migration at later stages was severely affected in both knock-down embryos. As a consequence of impaired primary mesenchymal cell migration, improper spicule formation was observed. These results indicate that ECM3 and QBRICK are components of extracellular matrix, which play important role in primary mesenchymal cell migration, and that sea urchin is a useful experimental animal model to investigate human disease-associated extracellular matrix proteins.

Usage of the Sea Urchin Hemicentrotus pulcherrimus Database, HpBase.

HpBase ( http://cell-innovation.nig.ac.jp/Hpul/ ) is a database that provides genome and transcriptome resources of the sea urchin Hemicentrotus pulcherrimus. In addition to downloading the bulk data, several analysis tools for resource use are available: gene search, homology search, and genome browsing. HpBase also discloses the protocols for biological experiments using H. pulcherrimus that have been accumulated so far. Therefore, HpBase can assist efficient use of genome resources for researchers from various fields-evolutionary, developmental, and cell biology. In this chapter we present an overview and usage of tools in HpBase.

Meis transcription factor maintains the neurogenic ectoderm and regulates the anterior-posterior patterning in embryos of a sea urchin, Hemicentrotus pulcherrimus.

Precise body axis formation is an essential step in the development of multicellular organisms, for most of which the molecular gradient and/or specifically biased localization of cell-fate determinants in eggs play important roles. In sea urchins, however, any biased proteins and mRNAs have not yet been identified in the egg except for vegetal cortex molecules, suggesting that sea urchin development is mostly regulated by uniformly distributed maternal molecules with contributions to axis formation that are not well characterized. Here, we describe that the maternal Meis transcription factor regulates anterior-posterior axis formation through maintenance of the most anterior territory in embryos of a sea urchin, Hemicentrotus pulcherrimus. Loss-of-function experiments revealed that Meis is intrinsically required for maintenance of the anterior neuroectoderm specifier foxQ2 after hatching and, consequently, the morphant lost anterior neuroectoderm characteristics. In addition, the expression patterns of univin and VEGF, the lateral ectoderm markers, and the mesenchyme-cell pattern shifted toward the anterior side in Meis morphants more than they did in control embryos, indicating that Meis contributes to the precise anteroposterior patterning by regulating the anterior neuroectodermal fate.

HpBase: A genome database of a sea urchin, Hemicentrotus pulcherrimus.

To understand the mystery of life, it is important to accumulate genomic information for various organisms because the whole genome encodes the commands for all the genes. Since the genome of Strongylocentrotus purpratus was sequenced in 2006 as the first sequenced genome in echinoderms, the genomic resources of other North American sea urchins have gradually been accumulated, but no sea urchin genomes are available in other areas, where many scientists have used the local species and reported important results. In this manuscript, we report a draft genome of the sea urchin Hemincentrotus pulcherrimus because this species has a long history as the target of developmental and cell biology in East Asia. The genome of H. pulcherrimus was assembled into 16,251 scaffold sequences with an N50 length of 143 kbp, and approximately 25,000 genes were identified in the genome. The size of the genome and the sequencing coverage were estimated to be approximately 800 Mbp and 100×, respectively. To provide these data and information of annotation, we constructed a database, HpBase (http://cell-innovation.nig.ac.jp/Hpul/). In HpBase, gene searches, genome browsing, and blast searches are available. In addition, HpBase includes the "recipes" for experiments from each lab using H. pulcherrimus. These recipes will continue to be updated according to the circumstances of individual scientists and can be powerful tools for experimental biologists and for the community. HpBase is a suitable dataset for evolutionary, developmental, and cell biologists to compare H. pulcherrimus genomic information with that of other species and to isolate gene information.

Complete mitochondrial genome of sea urchin: Hemicentrotus pulcherrimus (Camarodonta, Strongylocentrotidae).

The complete mitochondrial genome of sea urchin Hemicentrotus pulcherrimus (A. Agassiz, 1863) is determined. The genome is a circular molecule of 15,705 bp in length and the gene content and order are identical to those of other strongylocentrotid species. Lengths of all protein coding genes (PCGs) are identical to those of M. nudus. S. droebachiensis and S. pallidus but ND6 has a deletion of 3 amino acids at the C-terminal. The nucleotide composition of the genome shows a bias toward A + T (59.8%) and the A + T bias in PCGs (59.4%) is slightly higher than the other strongylocentrotid species. The mitogenome of H. pulcherrimus provides genomic information of the single species genus, Hemicentrotus and helps understand the genetic relationship and the evolution of the genus in relation to the other genera of the family Strongylocentrotidae.

Unc-5/netrin-mediated axonal projection during larval serotonergic nervous system formation in the sea urchin, Hemicentrotus pulcherrimus.

The molecular structure and role of two splice-isoforms of Unc-5 (Hp-Unc-5v1 and v2) in Unc-5/netrin interaction during serotonergic axonal projection were elucidated in this study. Hp-Unc-5v1 was found to be comprised of two immunoglobulin-like domains, two thrombospondin domains in the extracellular region, and ZU-5, DB, and Death domains in the cytoplasmic region, whereas Hp-Unc-5v2 lacked one thrombospondin domain, the transmembrane domain, and all cytoplasmic domains. Hp-Unc-5v1 was transcribed in unfertilized eggs, which continued until the 3-day post-fertilization (-dpf) 2-arm pluteus stage, but was suspended at the mesenchyme blastula stage (mB1), whereas Hp-Unc-5v2 was not transcribed in unfertilized eggs, but was from after fertilization to the same developmental stage of mB1 as Hp-Unc-5v1. Relative accumulation of transcripts of both splice-isoforms peaked at the prism stage and declined thereafter, and they were localized at the vegetal pole region of early gastrulae, around the blastopore in mid- to late gastrulae, at fore- and mid-gut regions and on the basal side of dorsal ectoderm in 28-hour post-fertilization prism larvae, and within axons at and after the 2-dpf pluteus stage. Hp-Unc-5v2:GFP was detected in the entire serotonergic cell body and extracellularly on the basal surface of oral ectoderm in 2-dpf plutei and exclusively within axons in 4-dpf plutei. Overexpression of Hp-Unc-5v2 resulted in decreased axonal projection in plutei. Knockdown of Hp-Unc-5v1 by morpholino antisense oligonucleotide resulted in severe deficiency of axonal projection. Interference of Unc-5/netrin interaction using an exogenous synthetic SQDFGKTW peptide from the VI domain in Hp-netrin inhibited axonal projection and larval swimming.

The 3'UTR of nanos2 directs enrichment in the germ cell lineage of the sea urchin.

Nanos is a translational regulator required for the survival and maintenance of primordial germ cells during embryogenesis. Three nanos homologs are present in the genome of the sea urchin Strongylocentrotus purpuratus (Sp), and each nanos mRNA accumulates specifically in the small micromere (sMic) lineage. We found that a highly conserved element in the 3' UTR of nanos2 is sufficient for reporter expression selectively in the sMic lineage: microinjection into a Sp fertilized egg of an RNA that contains the GFP open reading frame followed by Sp nanos2 3'UTR leads to selective reporter enrichment in the small micromeres in blastulae. The same result was seen with nanos2 from the sea urchin Hemicentrotus pulcherrimus (Hp). In both species, the 5'UTR alone is not sufficient for the sMic localization but it always increased the sMic reporter enrichment when present with the 3'UTR. We defined an element conserved between Hp and Sp in the nanos2 3'UTR which is necessary and sufficient for protein enrichment in the sMic, and refer to it as GNARLE (Global Nanos Associated RNA Lability Element). We also found that the nanos2 3'UTR is essential for the selective RNA retention in the small micromeres; GNARLE is required but not sufficient for this process. These results show that a combination of selective RNA retention and translational control mechanisms instills nanos accumulation uniquely in the sMic lineage.

Zinc finger homeobox is required for the differentiation of serotonergic neurons in the sea urchin embryo.

Serotonergic neurons differentiate in the neurogenic animal plate ectoderm of the sea urchin embryo. The regulatory mechanisms that control the specification or differentiation of these neurons in the sea urchin embryo are not yet understood, although, after the genome was sequenced, many genes encoding transcription factors expressed in this region were identified. Here, we report that zinc finger homeobox (zfhx1/z81) is expressed in serotonergic neural precursor cells, using double in situ hybridization screening with a serotonergic neural marker, tryptophan 5-hydroxylase (tph) encoding a serotonin synthase that is required for the differentiation of serotonergic neurons. zfhx1/z81 begins to be expressed at gastrula stage in individual cells in the anterior neuroectoderm, some of which also express delta. zfhx1/z81 expression gradually disappears as neural differentiation begins with tph expression. When the translation of Zfhx1/Z81 is blocked by morpholino injection, embryos express neither tph nor the neural marker synaptotagminB in cells of the animal plate, and serotonergic neurons do not differentiate. In contrast, Zfhx1/Z81 morphants do express fez, another neural precursor marker, which appears to function in the initial phase of specification/differentiation of serotonergic neurons. In addition, zfhx1/z81 is one of the targets suppressed in the animal plate by anti-neural signals such as Nodal as well as Delta-Notch. We conclude that Zfhx1/Z81 functions during the specification of individual anterior neural precursors and promotes the expression of tph and synaptotagminB, required for the differentiation of serotonergic neurons.

Nucleosome exclusion from the interspecies-conserved central AT-rich region of the Ars insulator.

The Ars insulator is a boundary element identified in the upstream region of the arylsulfatase (HpArs) gene in the sea urchin, Hemicentrotus pulcherrimus, and possesses the ability to both block enhancer-promoter communications and protect transgenes from silent chromatin. To understand the molecular mechanism of the Ars insulator, we investigated the correlation between chromatin structure, DNA structure and insulator activity. Nuclease digestion of nuclei isolated from sea urchin embryos revealed the presence of a nuclease-hypersensitive site within the Ars insulator. Analysis of micrococcal nuclease-sensitive sites in the Ars insulator, reconstituted with nucleosomes, showed the exclusion of nucleosomes from the central AT-rich region. Furthermore, the central AT-rich region in naked DNA was sensitive to nucleotide base modification by diethylpyrocarbonate (DEPC). These observations suggest that non-B-DNA structures in the central AT-rich region may inhibit nucleosomal formation, which leads to nuclease hypersensitivity. Furthermore, comparison of nucleotide sequences between the HpArs gene and its ortholog in Strongylocentrotus purpuratus revealed that the central AT-rich region of the Ars insulator is conserved, and this conserved region showed significant enhancer blocking activity. These results suggest that the central AT-rich nucleosome-free region plays an important role in the function of the Ars insulator.

HpSumf1 is involved in the activation of sulfatases responsible for regulation of skeletogenesis during sea urchin development.

Sulfatases such as arylsulfatase and heparan sulfate 6-O-endosulfatase play important roles in morphogenesis during sea urchin development. For the activation of these sulfatases, Cα-formylglycine formation by sulfatase modifying factor (Sumf) is required. In this study, to clarify the regulatory mechanisms for the activation of sulfatases during sea urchin development, we examined the expression and function of the Hemicentrotus pulcherrimus homologs of Sumf1 and Sumf2 (HpSumf1 and HpSumf2, respectively). Expression of HpSumf1 but not HpSumf2 mRNA was dynamically changed during early development. Functional analyses of recombinant HpSumf1 and HpSumf2 using HEK293T cells expressing mouse arylsulfatase A (ArsA) indicated that HpSumf1 and HpSumf2 were both able to activate mammalian ArsA. Knockdown of HpSumf1 using morpholino antisense oligonucleotides caused abnormal spicule formation in the sea urchin embryo. Injection of HpSumf2 mRNA had no effect on skeletogenesis, while injection of HpSumf1 mRNA induced severe supernumerary spicule formation. Taken together, these findings suggest that HpSumf1 is involved in the activation of sulfatases required for control of skeletogenesis.

ankAT-1 is a novel gene mediating the apical tuft formation in the sea urchin embryo.

In sea urchin embryos, the apical tuft forms within the neurogenic animal plate. When FoxQ2, one of the earliest factors expressed specifically in the animal plate by early blastula stage, is knocked down, the structure of the apical tuft is altered. To determine the basis of this phenotype, we identified FoxQ2-dependent genes using microarray analysis. The most strongly down-regulated gene in FoxQ2 morphants encodes a protein with ankyrin repeats region in its N-terminal domain. We named this gene ankAT-1, Ankyrin-containing gene specific for Apical Tuft. Initially its expression in the animal pole region of very early blastula stage embryos is FoxQ2-independent but becomes FoxQ2-dependent beginning at mesenchyme blastula stage and continuing in the animal plate of 3-day larvae. Furthermore, like FoxQ2, this gene is expressed throughout the expanded apical tuft region that forms in embryos lacking nuclear ß-catenin. When AnkAT-1 is knocked-down by injecting a morpholino, the cilia at the animal plate in the resulting embryos are much shorter and their motility is less than that of motile cilia in other ectoderm cells, and remains similar to that of long apical tuft cilia. We conclude that AnkAT-1 is involved in regulating the length of apical tuft cilia.

Targeted mutagenesis in the sea urchin embryo using zinc-finger nucleases.

We showed that engineered zinc-finger nucleases (ZFNs), which consist of a zinc-finger DNA-binding array and a nuclease domain of the restriction enzyme FokI, can introduce mutations at a specific genomic site in the sea urchin embryo. Using bacterial one-hybrid screening with zinc-finger randomized libraries and a single-strand annealing assay in cultured cells, ZFNs targeting the sea urchin Hemicentrotus pulcherrimus homologue of HesC (HpHesC) were efficiently selected. Consistent with the phenotype observed in embryos injected with an antisense morpholino oligonucleotide against HpHesC, an increase in the primary mesenchyme cell population was observed in embryos injected with a pair of HpHesC ZFN mRNAs. In addition, sequence analysis of the mutations showed that deletions and insertions occurred at the HpHesC target site in the embryos injected with the HpHesC ZFN mRNAs. These results suggest that targeted gene disruption using ZFNs is feasible for the sea urchin embryo.

Excision and transposition activity of Tc1/mariner superfamily transposons in sea urchin embryos.

Tc1/mariner superfamily transposons are used as transformation vectors in various model organisms. The utility of this transposon family is evidenced by the fact that Tc1/mariner transposons have loose host specificity. However, the activity of these transposons has been observed in only a few organisms, and a recent study in the ascidian Ciona intestinalis suggests that not all Tc1/ mariner transposons show loose host specificity. To understand host specificity, we used sea urchins, since they have a long history as materials of embryology and developmental biology. Transposon techniques have not been reported in this organism, despite the likelihood that these techniques would open up many experimental possibilities. Here we tested the activity of three Tc1/ mariner transposons (Minos, Sleeping Beauty, and Frog Prince) in the sea urchin Hemicentrotus pulcherrimus. Minos has both excision and transposition activity in H. pulcherrimus embryos, whereas no excision activity was detected for Sleeping Beauty or Frog Prince. This study suggests that Minos is active in a broad range of non-host organisms and can be used as a transformation tool in sea urchin embryos.

Spatiotemporal expression pattern of an encephalopsin orthologue of the sea urchin Hemicentrotus pulcherrimus during early development, and its potential role in larval vertical migration.

We have cloned and studied Hp-ECPN, an encephalopsin orthologue of the sea urchin Hemicentrotus pulcherrimus. Hp-ecpn cDNA was produced and found to contain a 1461-bp open reading frame that encodes 486 amino acids. Accumulation of Hp-ecpn mRNA and protein expression occurred at the 14 h postfertilization (hpf) swimming blastula stage and thereafter. The Hp-ECPN protein was N-glycosylated, and the amino acid sequence was similar to that of vertebrate encephalopsins. Whole-mount immunohistochemistry revealed the presence of Hp-ECPN in cells (ECPN cells) that appeared initially around the tip of the archenteron in 20 hpf early gastrulae. By the 54 hpf pluteus stage, ECPN cells had spread through the aboral ectoderm, and, by the eight-arm pluteus stage, were restricted to the tips of the larval arms and the posterior end of the body. The number of ECPN cells increased under conditions of continuous light, but decreased under continuous dark. Knockdown of Hp-ecpn mRNA using morpholino antisense oligonucleotides decreased the number of ECPN cells considerably, and inhibited the vertical swimming of the larvae. This suggested that Hp-ECPN plays a role in photosensitive larval swimming vertical migration. In adult tissues, the ECPN cells were detected exclusively in tube feet.

Spatio-temporal expression of a Netrin homolog in the sea urchin Hemicentrotus pulcherrimus (HpNetrin) during serotonergic axon extension.

A netrin homolog of the sea urchin, Hemicentrotus pulcherrimus (HpNetrin) was characterized in terms of its expression behavior. The predicted amino acid sequence was an ortholog of hemichordate netrin-1. Reverse transcriptase-PCR, immunoblotting, and whole mount immunohistochemistry showed that gene transcription and protein expression occurred from 15-hour post-fertilization (hpf) swimming blastula stage to, at least 53-hpf pluteus stage. HpNetrin was detected on the entire basal surface of the ectoderm in swimming and 16-hpf mesenchyme blastulae. However, by 24-hpf prism stage, immunoreaction on the oral ectoderm decreased, whereas it increased in the aboral ectoderm including near the animal plate ectoderm (area-I). By 48-hpf pluteus stage, the HpNetrin-rich area-I comprised a 40mm wide dorsal midline belt (DMB) that stretched from the dorsal posterior edge of the animal plate to the posterior end of the larval body. Serotonergic cells were first detected in the HpNetrin-moderate area between the anterior DMB and the HpNetrin-poor oral ectoderm near the animal plate in 24-hpf prism larvae. By 48-hpf pluteus stage, these cells extended axons toward the middle-ridge to form a neural plexus of the apical ganglion. At 53-hpf pluteus stage, these axons extended further away from the apical ganglion directly or through the DMB toward the HpNetrin-poor contralateral ectoderm. The protein expression and axon extension pattern were reproduced in embryonic cell-aggregates formed from artificially dissociated 20-hpf gastrulae and resembled small blastula. In Hpnetrin morpholino anti-sense oligonucleotide-injected plutei, serotonergic axon extension was severely inhibited. Thus, HpNetrin functions as a serotonergic axon guidance cue in this basal deuterostome.

cDNA cloning, expression, and characterization of methyl-CpG-binding domain type 2/3 proteins from starfish and sea urchin.

Two kinds of cDNAs that are highly homologous to mammalian MBD2 and MBD3 cDNAs were cloned from ovary of the starfish Asterina pectinifera. They are splicing variants and designated sMBD2/3a and sMBD2/3b cDNAs. sMBD2/3a cDNA spans 1378 bp and consists of a 48-bp upstream untranslated region, a 807-bp open reading frame encoding sMBD2/3a, and a 523-bp downstream untranslated region. sMBD2/3a and sMBD2/3b cDNAs encode proteins with predicted molecular weights of 30,724 and 29,635 consisting of 268 and 260 amino acid residues, respectively. The deduced amino acid sequences of these two are identical from residues 1 to 255, but different from residues 256 to the C-terminal ends. sMBD2/3a is expressed in all the tissues of starfish, whereas sMBD2/3b is highly expressed in ovary and oocytes, slightly in testis, but not in somatic cells. As suggested from the whole-genome sequence of Strongylocentrotus purpuratus, a sea urchin MBD2/3 cDNA was cloned from eggs of Hemicentrotus pulcherrimus and designated suMBD2/3 cDNA. It encodes a protein with predicted molecular weight of 30,778 consisting of 274 amino acid residues. All the three echinodermal MBD2/3 proteins consist of a methy-CpG-binding domain (MBD) and a coiled-coil domain, and only sMBD2/3a contains a glutamate-rich C-terminal region, a key mark in vertebrate MBD3. The three MBD2/3 proteins expressed in Escherichia coli and purified to homogeneity were capable to bind specifically to methylated DNA. It was shown that sMBD2/3a exists as dimer or in the monomer-dimer equilibrium, whereas sMBD2/3b and suMBD2/3 exist as monomer and dimer, respectively.

DNA variations within the sea urchin Otx gene enhancer.

We performed both intra- and interspecific comparisons of the Otx gene in the sea urchin to investigate DNA variations within the enhancer elements. Intraspecific comparisons within Hemicentrotus pulcherrimus showed that indel variations were rare within the Otx enhancer, whereas SNP variations were found uniformly within the whole test region. A similar pattern of DNA variation was observed in comparisons between closely related species. On the other hand, both nucleotide substitution and indel variations were at high levels between distant species. Additionally, the regions corresponding to the Otx enhancer in two related species showed substantial activities during development. Our results suggest the possibility that a stabilizing selection has occurred during the evolution of the Otx gene enhancer in the sea urchin that maintains its expression pattern.

Real-time monitoring of functional interactions between upstream and core promoter sequences in living cells of sea urchin embryos.

There are some functional compatibilities between upstream and core promoter sequences for transcriptional activation in yeast, Drosophila and mammalian cells. Here we examined whether similar compatibilities exist in sea urchin embryos, and if so, whether they are dynamically regulated during early development. Two reporter plasmids, each containing a test promoter conjugated to either CFP or YFP, were concurrently introduced into embryos, and their expression patterns were studied by fluorescence microscopy. The upstream sequence of the Hemicentrotus pulcherrimus (Hp) OtxE promoter drives the expression of its own core promoter and that of Strongylocentrotus purpuratus (Sp) Spec2a in different embryonic regions, especially at the late gastrula stage. Interestingly, when the four putative transcription factor binding sites of this upstream sequence were individually mutated, the resulting sequences directed different spatiotemporal expression from the same set of two core promoters, indicating that combinations of upstream factors may determine core promoter usage in sea urchin embryos. In addition, the insertion or deletion of consensus or nonconsensus TATA sequences changed the expression profile significantly, irrespective of whether the upstream sequence was intact or mutated. Thus, the TATA sequence may serve as a primary determinant for core promoter selection in these cells.