Expression and phylogenetic analysis of the zic gene family in the evolution and development of metazoans
© Layden et al; licensee BioMed Central Ltd. 2010
Received: 31 May 2010
Accepted: 5 November 2010
Published: 5 November 2010
zic genes are members of the gli/glis/nkl/zic super-family of C2H2 zinc finger (ZF) transcription factors. Homologs of the zic family have been implicated in patterning neural and mesodermal tissues in bilaterians. Prior to this study, the origin of the metazoan zic gene family was unknown and expression of zic gene homologs during the development of early branching metazoans had not been investigated.
Phylogenetic analyses of novel zic candidate genes identified a definitive zic homolog in the placozoan Trichoplax adhaerens, two gli/glis/nkl- like genes in the ctenophore Mnemiopsis leidyi, confirmed the presence of three gli/glis/nkl-like genes in Porifera, and confirmed the five previously identified zic genes in the cnidarian Nematostella vectensis. In the cnidarian N. vectensis, zic homologs are expressed in ectoderm and the gastrodermis (a bifunctional endomesoderm), in presumptive and developing tentacles, and in oral and sensory apical tuft ectoderm. The Capitella teleta zic homolog (Ct-zic) is detectable in a subset of the developing nervous system, the foregut, and the mesoderm associated with the segmentally repeated chaetae. Lastly, expression of gli and glis homologs in Mnemiopsis. leidyi is detected exclusively in neural cells in floor of the apical organ.
Based on our analyses, we propose that the zic gene family arose in the common ancestor of the Placozoa, Cnidaria and Bilateria from a gli/glis/nkl-like gene and that both ZOC and ZF-NC domains evolved prior to cnidarian-bilaterian divergence. We also conclude that zic expression in neural ectoderm and developing neurons is pervasive throughout the Metazoa and likely evolved from neural expression of an ancestral gli/glis/nkl/zic gene. zic expression in bilaterian mesoderm may be related to the expression in the gastrodermis of a cnidarian-bilaterian common ancestor.
The cnidarians are the only basal metazoan known to possess zic gene homologs. A single zic gene was identified in polymerase chain reaction (PCR) surveys of Hydra vulgaris  and the jellyfish Scolionema suvaense . Nematostella vectensis possesses five zic genes, of which four are arranged tandemly on the same genomic scaffold , suggesting tandem duplication from a single ancestral gene . Recent sequencing of the Amphimedon queenslandica sponge genome identified one putative ancestral gli/glis/nkl gene, one glis-like gene and one gli-like gene . However, no zic homolog was identified . zic homologs have not been verified in placozoans and ctenophores.
Expression of zic homologs during bilaterian development suggests that the urbilaterian zic gene was expressed in neural and mesodermal domains. zic expression in protostomes (Ecdysozoa and Lophotrochozoa) has been characterized in Drosophila melanogaster, Caenorhabditis elegans and Tubifex tubifex. The fruit fly D. melanogaster zic (odd-paired; opa) is expressed in all ectodermal and mesodermal precursors in the presumptive segmented region of the embryo [7, 8]. At later stages opa is expressed in the neural ectoderm  and in a subset of visceral mesoderm around the midgut . In the nematode C. elegans, the zic homolog, ref-2, is expressed in a subset of neural precursors [10, 11]. In the lophotrochozoan oligochaete T. tubifex, Ttuzic expression is observed in the mesodermal m-blast cells and germ bands, mesoderm associated with chaetal sacs, anterior ectoderm and the developing brain . In summary, protostome zic genes are expressed in a sub-set of neural-ectoderm and a sub-set of mesoderm. However, the structures expressing zic within a particular germ layer vary between species (subset of mesodermal precursors in T. tubifex versus broad mesodermal precursor in D. melanogaster). Moreover, expression even within similar tissues (for example, neural) varies in both timing and location (anterior brain in T. tubifex versus ventral nerve cord (VNC) in D. melanogaster).
In deuterostomes, zic expression and function have been best characterized in the chordate lineage. The five vertebrate zic gene expression patterns show extensive overlap during development. In the mouse, zebrafish and Xenopus embryos, zic is expressed in presumptive ectoderm at or just prior to gastrulation . The ectodermal expression resolves to the neural plate and is ultimately restricted to the dorsal neural tube and the lateral boundary of the neural plate where neural crest cells will form [2, 3]. In addition, vertebrate zic genes are expressed in the dorsal somitic mesoderm and MmuZic1 and MmuZic3 are expressed in the developing retina . In ascidians, two zic genes have been identified (macho-1 and zicL/N), which display similar patterns among the three species examined [14–18]. zicL/N is expressed in presumptive muscle and notochord mesodermal lineages and in neural precursors [14–16]. macho-1 mRNA is localized to mesodermal precursor blastomeres [17, 18]. The amphioxus zic gene (amphizic) is detected in the dorsal ectoderm and underlying mesoderm in gastrula, in the mesodermally-derived roof of the anterior archenteron, in the dorsal neural tube and in dorsal-lateral somatic mesoderm in the developing somites .
Previously, the only zic gene expression pattern to be characterized in early branching metazoans was in the cnidarian H. vulgaris. In these animals zic (hyzic) expression is detected in adult polyps within the dividing i-cells , which are stem cells that generate cnidocytes, neurons and other cell types . However, embryonic expression of zic gene homologs has yet to be characterized in cnidarians.
In order to better understand the diverse roles of zic genes and their potential role in the evolution and development of metazoans, we have investigated the origins and expression of the zic gene family in a diverse group of previously unsampled metazoans. Our findings indicate that the zic gene family probably arose from a gli/glis/nkl-like precursor gene in the common ancestor of Placozoa/Cnidaria/Bilateria, as supported by a definitive zic ortholog in the Trichoplax adhaerens genome and a lack of zic homologs in either of the more basally branching Ctenophora or Porifera. Additionally, we present expression patterns for zic genes during the development of the cnidarian N. vectensis, the lophotrochozoan C. teleta and for gli and glis genes in the ctenophore M. leidyi.
Zic genes evolved prior to the split between Placozoa and Cnidaria
List of genes and abbreviations used in this study.
Genes used: Accession No./JGI model ID/Sponge ID
MIGlis [HM265718], MIGLi [Genbank submission number 1388697]
AmqGli2/3b [Aqu1.219964], AmqGli2/3a [Aqu1.217717], AmqGlis1/3 [Aqu1.213405]
TadZic[XP_002108473], Tad188.8.131.52 [XP_002118240], Tad1.1.692.1 [XP_002107781]
NvZicA [AB231867], NvZicB [AB231868], NvZicC [AB231868], NvZicD [AB231868], NvZicE [AB231868], NvGli [EU162649], NvGlis [XP_001637240], NvNkl [XP_001636498] Nvpredprot [XP_001636147.1]
SsuZic [AB231883 ]
HmaZic [XP_002166449 ]
DjaZicA [AB231880], DjaZicB [AB231881]
DacZicA [AB266039 ], DacZicB [AB266038 ]
TtuZic [AB231870 ]
Ct-Zic [HM235719], Ct-Glis [e_gw184.108.40.206], Ct-Gli [e_gw220.127.116.11]
opa [NP_524228], ci [NP_524617], Sug [AAS65032], imd [AAF57692]
AgaSug [XP_309008] ]
SpuZic [XP_792929], SpuGlis [XP_798511]
ApeZic [AB231872 ]
SkoZic [NP_001158430], SkoGlis2 [XP_002738784]
CinZicL [NP_001071853], Cinmacho1 [NP_001027958], CinGli [NP_001071951], CinGlis [NP_001071922]
MteZicL [BAE54350], Mtemacho1 [BAE54349]
BflZic (aka amphiZic) [CAB96573], BflGlis [XP_002611119], Bflpredprot [XP_002610135.1]
XlaZic4 [BAF36750], XlaZic2 [AAH82436], XlaGli2 [AAD28180], XlaGlis2 [NP_001082092], XlaGli3 [NP_001081440.1] XlaZic3 [NP_001081088.1], XlaZic5 [BAA95699.1]
MmuZic5 [NP_075363], MmuZic2 [NP_033600], MmuZic1 [AAH60247], MmuGlis2 [AF325913], MmuGli1 [BAA85004.1], MmuGli2 [NP_001074594], MmuGli3 [NP_032156], MmuGlis1 [AAM93156], MmuGlis3 [ZB131654.1] MmuZic3 [CAM20754.1], MmuZic4 [NP_033602.2], MmuZF143 [EDL16976.1]
RraZic2 [NP_001101862], RraGli2 [NP_001100639],
HsaZic2 [AAC96325], HsaGli3 [CAB59315], HsaGlis3 [CAH70655], HsaGlis1 [NP_671726], HsaGli1 [AAM13391.1] HsaGlis2 [NP_115964], HsaZic5 [NP_149123], HsaZic3 [NP_003404], HsaZic1 [NP_003403.2], HsaZic4 [AAH29507.1] HsaZF143 [CAC17610.1]
CfaGlis3 [XP_541295], CfaGlis1 [XP_546702.2]
BtaGlis1 [XP_615122], BtaGlis2 [DAA15722.1], BtaGlis3 [XP_002689658.1]
DreZic4 [NP_001070080], DreGli1 [NP_840081], DreGli2b [NP_001015069], DreZic5 [AAQ67349.1]
DreZic1 [NP571008.1], DreZic3 [NP001001950.1], DreZic2 [AAG35717.2], DreZf143 [AA124736.1]
GgaZic1 [BAB92091], GgaGlis1 [XP_422485], GgaGli1 [P55878], GgaZF143 [XP_426401.2],
Phylogenetic relationships among candidate genes were determined by performing Bayesian  and Maximum Likelihood  (Figure 2) analyses using sequences beginning with the first cysteine residue of the first ZF domain (ZF1) through the last histidine residue in the fifth ZF domain (ZF5) (Additional File 1). Our analyses confirmed previously identified cnidarian zic genes (NvzicA-E)  and suggested that the two novel Nematostella candidate genes were not zic genes but, rather, were members of the glis and nkl clades. Bayesian and Maximum Likelihood analyses grouped NvGlis with the Glis family (Bayesian posterior probability (PP) = 97, Maximum Likelihood bootstrap value (BS) = 0.75; Figure 2), and grouped NvNkl with the Nkl family (PP = 1, BS = 0.99; Figure 2). Two of the three placozoan T. adhaerens candidate genes failed to group within the gli/glis/nkl/zic superfamily (see Tad1.1.1692.1 and Tad18.104.22.168; Figure 2), but reciprocal BLAST searches suggested that they were related to ZF143 genes rather than members of the gli/glis/nkl/zic gene family. This is supported by phylogenetic analysis for Tad1.1.1692.1 (PP = 1, BS = 0.96; Figure 2). Tad22.214.171.124 groups strongly with the ZF143 (PP = 1, BS = 0.99; Figure 2) genes rather than the gli/glis/nkl/zic genes suggesting that it is also not a gli/glis/nkl/zic related gene. However the third T. adhaerens gene, Tadzic, grouped within the zic gene sub-family with strong support (PP = 1, BS = 0.99; Figure 2). In the polychaete C. teleta, phylogenetic analyses strongly support the assignment of the three candidate genes, Ct-zic, Ct-gli and Ct-glis, into the zic, gli and glis sub-families, respectively (PP = 1, 1, 0.97 and BS = 0.99, 0.79 and 0.75, respectively; Figure 2). In addition, we cloned a partial sequence out of the acoel C. longifissura, Clozic, which grouped within the zic sub-family (PP = 1, BS = 0.97; Figure 2). The identification of Tadzic in Placozoa is the putative earliest metazoan zic ortholog identified to date. We conclude that zic genes were present in the common ancestor of Placozoa and later branching metazoans.
Two gli/glis/nkl/zic sequences were identified in M. leidyi despite extensive searches of the M. leidyi sequenced genome (50X sequence coverage, Ryan et al., manuscript in preparation). One M. leidyi sequence, Mlglis, grouped with the glis sub-family (PP = 0.97, BS = 0.75; Figure 2). The second M. leidyi sequence, Mlgli, grouped with the gli sub-family (PP = 1, BS = 0.79; Figure 2). However, this grouping is based on the presence of only four of the five ZF domains, ZF2-ZF5 (Additional File 1). Attempts to identify the ZF1 domain via RACE PCR and sequence predictions were unsuccessful (data not shown). Thus, we cannot confirm that Mlgli contains all five ZFs. However, because it groups so strongly with the gli sub-family, we included it in our analysis. We confirmed the reported orthology of the three previously described A. queenslandica sequences. We identified an ancestral/sister gli/glis/nkl sequence (Amqgli2/3b, PP = 0.99, BS = 0.82; Figure 2), one gli-like sequence (Amqgli2/3a, PP = 1, BS = 0.79), and one glis-like sequences (Amqglis1/3, PP = 0.95, BS = 0.65). Examination of the spacing between the cysteine residues in the ZF1 peptide sequence is also consistent with Mlglis, Amqgli2/3a, Amqgli2/3b and Amqglis1/3 sequences being more closely related to the gli/glis/nkl sub-families than the zic gene sub-family (Additional File 1). Degenerate PCR in a second sponge species (Ephydatia muelleri) failed to identify a poriferan zic sequence (data not shown). We conclude that no definitive zic ortholog exists in the ctenophore or poriferan animals analysed, but that both groups contain gli/glis/nkl-like homologs in their genomes.
Clozic and Ct-zic, but not Tadzic, contain ZF-NC domains
As previous studies  have suggested that zic genes, but not gli or glis genes, contain a ZF-NC domain, we attempted to identify ZF-NC domains in the newly described zic genes Tadzic, Ct-zic and Clozic as well as in the gli/glis/nkl genes. We identified no ZF-NC or ZOC domains outside of the zic homologs. The ZF-NC domain is located immediately N-terminal to the ZF1 domain, so we aligned all the sequences N-terminal to the ZF1 domains to identify candidate ZF-NC domains. We identified clear ZF-NC domains in CloZic, Ct-Zic and a partially conserved ZF-NC in TadZic (Figure 1B). The ZF-NC consensus sequence has been previously described as GAF(F/L)RYMRQP-(0-7AA)-IKQE . Ct-Zic and CloZic each have 87% and 93% conserved similarity with this consensus sequence across their respective ZF-NC domains (Figure 1B). However, the TadZic sequence is poorly conserved except for the presence of the four C-terminal amino acids, a conserved proline (Figure 1B), and a conservative arginine to histidine change at position 5 (Figure 1B). While CloZic and Ct-Zic contain definitive ZF-NC domains, it is not clear if the TadZic sequence represents a functional ZF-NC domain.
A single zic ortholog containing a ZOC domain was shared by Cnidaria and Bilateria
Using non-biased alignment of the N-terminal portions of all Zic proteins assayed, we identified a conserved RDFL-(1-2AA)-RR ZOC consensus sequence. The ZOC domain previously described consensus sequence was (S/T)RDFLxxxR [1, 4]. However, our alignment suggests that the RDFL-(1-2AA)-RR consensus sequence is more prevalent throughout the metazoans than the previously described consensus. Using the novel consensus sequence, we identified a ZOC domain in NvZicA,C,D,E, Ct-Zic, and CloZic (Figure 1C). The NvZicA ZOC sequence is 100% identical to the new ZOC consensus sequence. This is particularly important, because ZOC domains have not been previously described in cnidarians . Our finding pushes back the emergence of the ZOC domain to at least the cnidarian-bilaterian common ancestor. The NvzicC-E genes possess more weakly conserved putative ZOC domains. However, in the case of NvzicD, the amino acid content of the domain HDYIHH, is 100% similar to the RDFL-RR sequence (Figure 1C), suggesting that the NvzicD ZOC domain may still be functional. Altering the definition of the ZOC domain results in the inclusion of a number of additional genes and taxa into the ZOC possessing zic group. They are the cnidarians (see Nvzic A,C,D,E; Figure 1C), the vertebrate Zic4 and Zic5 sequences (see MmuZic5, XlaZic4; Figure 1C), nematodes (see CelRef2; Figure 1C), platyhelminthes (see DjaZicA; Figure 1C) and urochodates (see amphiZic; Figure 1C). We did not find a ZOC-like sequence for TadZic. However, due to sequence and annotation errors, it is possible that we have not identified a full length N-terminal sequence and, thus, a ZOC-like domain may be encoded in the full Tadzic coding sequence. Using sequence similarity, it appears that putative ZOC domains can be identified in cnidarians, acoels, platyhelminthes, nematodes, annelids, vertebrates, insects and urochordates (Figure 1C). Thus, the earliest described ZOC sequence appears to have arisen in the cnidarian-bilaterian ancestor, rather than the urbilaterian ancestor.
Nvzic expression during development of the cnidarian sea anemone N. vectensis
We investigated the expression pattern of the five Nvzic genes during development of the cnidarian N. vectensis. NvzicA and NvzicB appear to be expressed at very low levels, as we were unable to obtain reproducible expression for both genes, although they are detected by RT-PCR in 24, 48 and 72 hours post fertilization hpf cDNA samples suggesting that they are expressed at low levels beginning just after gastrulation and maintain expression into at least the planula larval stages (data not shown).
NvzicE expression is first detected during embryonic stages, which is earlier than either NvzicC or NvzicD. It is weakly expressed throughout the embryonic ectoderm shortly after gastrulation (Figure 3K). At planula stages, weak ubiquitous expression of NvzicE is detected throughout the ectoderm but enrichment is observed in oral ectoderm and presumptive tentacle bud domain and at the aboral pole in the apical tuft (Figure 3L; arrow and arrowhead, respectively). By late planula stages NvzicE expression is down regulated in the body wall ectoderm but remains detectable in the apical tuft ectoderm (Figure 3M; arrow) and the presumptive tentacle domain (Figure 3M). By the tentacle bud stage, aboral expression of NvzicE is no longer detectable (Figure 3N), due to the loss of the apical tuft after settlement. However, expression persists in tentacle endoderm (Figure 3O) and, occasionally, single NvzicE- positive cells appear in the tentacle ectoderm in juvenile polyps (Figure 3R). In summary, NvzicC-E display similar expression patterns in the developing tentacle endoderm and ectoderm, the endodermal component of the mesenteries and in the oral ectoderm. In addition, NvzicE is expressed in the apical tuft ectoderm of the late embryonic stages and in individual cells throughout the body of the developing animal.
Ct-zic expression during development of the bilaterian lophotrochozoan C. teleta
Prior to this study zic gene expression had only been investigated in a single lophotrochozoan, the oligochaete T. tubifex. Most protostome zic expression has been described in the ecdysozoans. Thus, investigating additional lophotrochozoan animals will provide a more complete comparison for understanding the relationship among bilaterian zic gene expression patterns.
Neural ectodermal expression of Ct-zic is also detected along the lateral edge of the ventral neural ectoderm at stage 4 (Figure 4C and 4O, arrows). The lateral ectodermal Ct-zic domain expands toward the posterior of the trunk (Figure 4F), and is then down-regulated in an anterior to posterior wave, visible by late stage 6 (Figure 4G; white lines mark anterior border of expression). From stages 7 - 9, Ct-zic expression in the lateral ectoderm is maintained in ventral-lateral cells in the growth zone (Figure 4I, arrow), where new segments are forming. Although the exact lateral border of the ventral neural ectoderm, which forms the ventral nerve chord (VNC), has not been mapped, we think that lateral ectodermal Ct-zic expression is within the ventral neural ectoderm. In support of this hypothesis, early neurogenic markers (Ct-ash and Ct-neurogenin; unpublished observations) positionally overlap with Ct-zic expression, although expression of the neurogenic markers extends more medially than Ct-zic at similar stages. Furthermore, we observe Ct-zic in the VNC beginning at stage 6 (Figure 4G, arrowhead). Ct-zic positive cells in the VNC are first detected in anterior segments where Ct-zic expression is no longer apparent in the ventral-lateral ectoderm (Figure 4G). By stage 8, Ct-zic positive cells are positioned medially and laterally within the VNC (Figures 4K and 4Q). Thus, we conclude that Ct-zic is expressed at the lateral edge of the presumptive ventral neural ectoderm, in the VNC, and in a ventral-lateral ectodermal domain in the growth zone where new segments are forming.
Mesodermal expression of Ct-zic is detectable at stages 3 to 9. Beginning at stage 3, expression is in one or a few cells in an anterior subset of the bilateral mesodermal bands (Figure 4A and 4B). At stages 4 and 5, mesodermal expression is detected in domains lateral to the foregut (Figure 4C and 4E, arrowheads). The expression around the foregut expands during stages 6 - 7 (Figure 4G and 4I) and is in a sub-set of the Ct-twist positive visceral mesoderm surrounding the foregut  (Figure 4I, notched arrowhead). In addition, Ct-zic expression is observed within a subset of Ct-twist-positive anterior mesoderm in the head (Figure 4I, arrowhead).
There is additional Ct-zic mesodermal expression in four longitudinal stripes (two on each side of the larva) of segmentally-reiterated cell clusters. This expression domain is first visible in a few cells at stage 5 (Figure 4E, 4F and 4O, arrowheads). The number of Ct-zic-positive cell clusters increases along the anterior-posterior (A-P) axis as the animal elongates and new segments form (Figure 4J and 4L, arrowheads). The A-P expansion of this pattern is reminiscent of the expansion of Ttu-zic expression in chaetal sac associated mesoderm . By stages 7 and 8, when chaetae are present, Ct-zic expression can be clearly observed in cells wrapping the chaetal sacs (Figure 4J, 4L and 4P) but not in the chaetal sacs themselves (Figure 4P, arrows), indicating that the zic-positive cells are chaetal-sac associated mesoderm. Ct-zic expression is maintained in these structures through stage 9. In summary, mesodermal expression of Ct-zic is found in the anterior portion of mesoderm bands at stage 3, in the visceral mesoderm surrounding the foregut, in the head mesoderm at the anterior tip of the animal and in mesoderm associated with forming chaetae.
Mlglis and Mlgli expression during development of the ctenophore M. leidyi
Mlglis expression during M. leidyi development is exclusively associated with the apical organ and is not detected in any other tissues. Expression is first detected upon the completion of gastrulation in two opposing groups of cells in the presumptive apical organ ectoderm at the aboral pole of the developing embryo (Figure 5A). Over the next hour of development, expression expands into two arches around the floor plate ectoderm of the apical organ (Figure 5B). The number of cells in each domain expands from approximately four to twenty, though it is not clear if the expansion is due to nascent expression of Mlglis in adjacent cells or via cell proliferation of the Mlglis expressing cells. Over the next three hours, Mlglis domains assume a dumbbell shape and eventually bifurcate (Figures 5C-E), resulting in four distinct domains of expression in the floor plate of the apical organ (Figure 5F). The apical organ contains a statocyst consisting of four groups of cells that give rise to balancing cilia that support the mineral containing lithocytes. The statocyst functions as a balance organ allowing recognition of proper body orientation within the water column. This is achieved by sensing the gravitational pull on the lithocytes perched on the four groups of balancing cilia . Mlglis expression is detected in the statocyst cells associated with the balancing cilia (Figure 5G and 5H, arrows) which, based on function, are putative neuronal or neuronal-like cells. In conclusion, Mlglis is expressed in the presumptive balancing cilia of the apical organ statocyst beginning shortly after gastrulation and expression is maintained throughout development.
Mlgli expression is first detected in 24-hour-old larvae in four groups of putative sensory cells in the floor of the apical organ. Mlgli- expressing cells are symmetrically arranged within the floor of the apical organ in the 24 hpf cydippid larva (Figure 5I and 5J). Based on their position, the Mlgli- expressing cells and the Mlglis expressing cells represent two distinct populations of cells (compare Figure 5J with 5H). The Mlgli expressing cells are not associated with the balancing cilia and are positioned further from the centre along the adesophageal plane of the apical organ than the Mlglis expressing balancing cilia cells. The domains exhibiting Mlgli expression are also associated with putative sensory structures. Mlgli expression in the apical organ is reminiscent of previously identified structures, which morphologically resemble photoreceptors by electron microscopy . We conclude that Mlgli is expressed in the putative photoreceptive cells in the apical organ and has a distinct pattern from that observed for Mlglis.
The origin of the metazoan zic gene
Regardless of when zic genes evolved, they probably diverged from a gli/glis/nkl-like gene. The basally branching ctenophore and poriferans likely possessed at least two gli/glis-like genes, arguing that a gli/glis/nkl-like gene, which is lacking in non-metazoans, was present and already duplicated at least once at, or near the time of, the emergence of the earliest metazoans. The presence of the poriferan Amqgli2/3b sequence suggests that a third gli/glis/nkl-like gene may have existed prior to the divergence of poriferans and ctenophores. However, Amqgli2/3b may also represent a lineage specific duplication. As the gli/glis/nkl clades represent a single sister group to the zic clade, it is unclear whether a gli, glis or nkl gene duplication/divergence would have given rise to the zic family. No clear nkl genes are present in the basal metazoans, except for the cnidarians, arguing that duplication of a gli or glis-like gene resulted in the divergence of the zic family. However, the existence and phylogenetic position of the Amqgli2/3b gene raises concerns about this model. Duplication or divergence of an Amqgli2/3b-like gene could have also given rise to the zic family. Further resolution of the branching pattern of basal metazoans and sampling for gli/glis/nkl/zic genes in ctenophore and poriferan species is required before more a finite hypotheses can be drawn about the evolutionary origin and precursor of the zic gene family.
Two distinct gli similar (glis) sub-families exist
Though the glis and nkl gene sub-family names have been used interchangeably, nkl should uniquely identify the glis2/sug gene family. glis2 genes were originally described as neurogenic kruppel like (nkl) genes , while the glis1/3 genes had been identified as gli similar . Our data suggest that the nkl and glis genes represent two distinct clades. The glis clade consists of glis1, glis3 and the D. melanogaster lameduck (lmd) homologs (Figure 2), and the nkl clade consists of glis2 and D. melanogaster sugarbabe (sug) homologs (Figure 2). This distinction is also supported by a recent investigation of C2H2 ZF genes in amphioxus. The authors of that study identified two distinct strongly supported glis clades for glis2 and glis1/3 genes . We propose that the nkl gene sub-family name should be used to uniquely identify the sug and glis2 genes as a distinct gene sub-family.
ZOC and ZF-NC domains pre-date emergence of Bilateria
The cnidarian-bilaterian zic ortholog contained a ZF-NC-like domain. The first eight amino acids of the published ZF-NC conserved sequence (GAFFRYMR(0-7aa)IKQE)  are 86% identical and 100% similar in the putative ZF-NC domain of NvZicA (GAFFRFMR; Figure 1B). In most ZF-NC sequences characterized to date, the arginine residue at position 8 is followed by a QPIKQE sequence (Figure 1B). The corresponding NvZicA sequence (SPAKDN) is 50% similar to the C-terminal portion of the ZF-NC domain (Figure 1B). Lastly, although the first part of the consensus domain is missing in the TadZic sequence, the last four amino acids (AKQE) are 75% identical and 100% similar to the consensus sequence. As no functional mapping has been conducted on the ZF-NC domain and there are strong sequence similarities between the ZF-NC regions of early branching metazoans and bilaterians, it seems likely that presence of a ZF-NC domain may predate the urbilaterian zic gene (Figure 6C, green square).
NvZicA-E each contain all, or some of, the reported, and/or proposed, ZOC consensus sequence or have amino acids in similar functional groups to the consensus sequence (Figure 1C). In the case of NvZicE, the actual sequence is not identical to the consensus sequence, but 100% similarity is observed in the five N-terminal amino acids (Figure 1C). In fact, the NvZicE sequence is arguably more similar to the consensus sequence than the oligochaete Ttu-Zic sequence (compare NvZicE and TtuZic; Figure 1C), which has a ZOC domain . Thus, we suggest that the conserved ZOC and ZOC-like domains in the Nematostella Nvzic genes supports the inclusion of the ZOC domain in the cnidarian-bilaterian ancestral zic gene (Figure 6C, red square). Our results suggest that previously undescribed ZOC domains are present in platyhelminthes, nematode and urochordate zic genes. Thus, all sampled bilaterian clades to date now have at least one lineage with a zic gene containing a ZOC domain.
ZOC domains have been shown to have two functions in in vitro transcription assays. The MmuZic2 ZOC domain has been demonstrated to increase Zic transcription factor activity and to bind the transcription factor Mfa-1 . Currently it is not clear if the ZOC domain increases Zic transcription factor activity via intramolecular interactions, intermolecular interactions or both in vivo. Future studies focused on characterizing the activity of cnidarian ZOC sequences are required in order to verify the existence of functional ZOC domains in the cnidarians. Identification and comparative analysis of ZOC interacting proteins in both cnidarian and bilaterian lineages may be useful for determining the molecular relationships of zic gene function in metazoan development and evolution.
Evolution of mesodermal zic expression
Mesodermal expression of zic genes is present throughout the Bilateria. In D. melanogaster, opa is expressed broadly throughout the mesoderm primordia of the segmented region of the embryo . In both lophotrochozoans studied thus far (T. tubifex and C. teleta) zic is expressed in a segmentally reiterated pattern in mesoderm associated with the chaetal sacs  (Figure 4) and, in chordates, zic expression is observed in a segmented pattern in dorsal domains of forming somites. In addition, within protostomes, zic homologs are expressed in visceral mesoderm [surrounding foregut in C. teleta (Figure 4) and around a portion of the forming midgut in D. melanogaster] . In Cnidaria, we observe Nvzic expression in the gastrodermis, a bifunctional endomesodermal tissue, in the forming tentacles, around the pharynx and in the endomesodermal component of the directive mesenteries (Figure 4).
One hypothesis raised by our data is that bilaterian mesodermal zic expression is derived, at least in part, from the endomesodermal expression that existed in the common ancestor of the Bilateria and Cnidaria. The expression of NvzicD and NvzicE in tentacle endomesoderm and NvzicD in endomesoderm surrounding the pharynx supports this hypothesis. The tentacles are capable of multiple complex movements, including capture and transport of food to the oral opening due to the high density of myoepithelial cells. It is possible that there are similarities between tentacle endomesoderm and bilaterian musculature. Future work characterizing cell types expressing NvzicD and NvzicE and comparing the zic function in N. vectensis endomesoderm will provide an insight to the overall relationship of cnidarian endomesoderm to bilaterian endoderm and mesoderm.
zic expression in neural development
Our data support an ancestral role for zic genes as regulators of neural sub-type development. Expression and function of zic homologs in neural development is present throughout the Bilateria. However, only in the chordate lineage does zic expression and function include the entire neurogenic region in the early embryo. This suggests that, though zic genes have a conserved neural role in bilaterian animals, the ancestral role was probably sub-type specification. Nvzic genes are expressed in the presumptive and developing tentacle buds in both the endoderm and ectoderm in the cnidarian N. vectensis (Figure 3J and 3O and 4E, respectively). The tentacles are highly neuralized structures , suggesting that zic may be regulating tentacular neurogenesis. In support of this hypothesis, a number of neural genes (Nvelav, NvanthoRF, Nvmushashi, Nvgcm) and cross-reactive antibodies against the serotonin and GABA neuropeptides are detected at high levels in the developing tentacle endoderm and ectoderm . In addition, NvzicE expression in the forming apical tuft sensory organ (Figure 3M) is consistent with a role in neural development. zic expression in the apical tuft of echinoderms has also been described . We describe expression of Mlglis and Mlgli in the ctenophore M. leidyi in neural structures of the developing apical organ (Figure 5). Thus, we suggest that the role of zic genes in neural development may have been shared with a gli/glis/nkl/zic ancestral gene. This is also supported by the observations that gli, glis and nkl genes all have described roles in metazoan neurogenesis [2, 3]. Previous studies and our data, taken together, support the hypothesis that one of the ancestral roles of the zic gene family was to specify neural sub-types.
Zic homologs may have been co-opted to pattern segmentally reiterated structures along the A-P axis in Bilateria
In the lophotrochozoans T. tubifex and C. teleta, zic is expressed in mesoderm associated with segmentally reiterated chaetae. In amphioxus and vertebrates, zic homologs are expressed in the dorsal region of somites. Both somitic expression and mesodermal expression associated with chaetae indicate a role in segmental musculature formation for zic homologs in lophotrochozoans and chordates. C. teleta zic homologs are also expressed in a reiterated pattern in the ventral nerve cord, reflecting the segmental organization of its body plan (Figure 4). Likewise, opa has a role promoting the proneural gene achaete expression in each segment during D. melanogaster neurogenesis . opa mutants were originally identified for their pair-rule phenotype (every other segment disrupted) . However, opa's ubiquitous expression throughout the presumptive segmented region is distinct from the typical seven-stripe pair-rule pattern. In addition, opa is required for proper expression of the wg (wnt) segment polarity gene in all parasegments . Together, these data suggest that zic homologs function to pattern segmentally reiterated structures and molecular domains in bilaterians, but are unlikely to function directly in the segmentation process. Rather, the lack of similarity between zic expressing reiterated structures across bilaterian lineages suggests that zic genes, due to their ancestral roles in neural and mesoderm development in a cnidarian-like ancestor, have been co-opted multiple times downstream of the segmentation programme to pattern various segmentally iterated structures.
Based on our analysis, we propose that the first metazoan zic arose from a gli/glis/nkl- like gene prior to the emergence of the Placozoa, but after the divergence of both Ctenophora and Porifera. ZOC and ZF-NC domains were probably present in the ancestral zic gene shared by the cnidarian-bilaterian common ancestor. We hypothesize that the metazoan zic neural expression, which is present in all animals assayed thus far, may be derived from the expression of an ancestral gene in sensory cells. We also propose that the mesodermal expression of bilaterian zic genes may be derived from gastrodermal expression of zic homologs in a cnidarian-bilaterian ancestor. Lastly, we suggest that zic genes have been co-opted multiple times down stream of segmentation programmes to pattern segmentally reiterated structures in some bilaterian clades. In summary, our results raise an interesting hypothesis about the molecular link between cnidarian endomesoderm and bilaterian mesoderm, and suggest that a role for zic genes in neurogenesis may be an ancient metazoan feature.
Identification of candidate gli/glis/nkl/zic genes
In order to identify candidate zic genes we searched published genomic resources, Trichoplax adhaerens http://genome.jgi-psf.org/Triad1/Triad1.home.html, Nematostella vectensis http://genome.jgi-psf.org/Nemve1/Nemve1.home.html, Capitella teleta http://genome.jgi-psf.org/Capca1/Capca1.home.html, the unpublished genome of the ctenophore Mnemiopsis leidyi (NIH), the Amphimedon queenslandica http://spongezome.metazome.net/cgi-bin/gbrowse/amphimedon/ and using keyword or TBLASTN and BLASTX searches http://www.ncbi.nlm.nih.gov/BLAST/. We also used previously identified degenerate PCR primers  in order to attempt to clone zic fragments from the mixed stage embryonic cDNA acoel Convolutriloba longifissura, sponge Ephydatia muelleri and ctenophore Mnemiopsis leidyi using previously described zic specific primers . Degenerate fragments were cloned by TA cloning into the pGEM-T Easy vector (Promega, CA, USA) and sequenced (Macrogen, Seoul, Korea). Sequences were used to design nested 5' and 3' RACE primers for each candidate sequence. Full-length sequences were obtained by compiling 5' and 3' RACE sequences. Finally, sequences were confirmed by amplifying the full length sequence from mixed stage embryonic cDNA.
Full length sequences for Gli, Glis and Zic proteins were aligned using the Muscle Alignment  web server http://www.ebi.ac.uk/Tools/muscle/index.html. The protein sequences were trimmed using the N-terminal cysteine residue of the ZF1 domain and C-terminal histidine residue of the ZF5 domain to demarcate the five C2H2 ZF domains (Additional File 1). A Bayesian phylogenetic analysis was carried out using MrBayes 3.1.2  with a mixed protein model for 1,000,000 generations sampled every 100 generations with four chains. The Jones model was chosen with a posterior probability of 1. A summary 'consensus' tree was generated in MrBayes using the last 7500 trees (30,000 total trees) representing 750,000 stationary generations. Posterior probabilities were calculated from the 'consensus' tree. A maximum likelihood analysis (using PhyML 3.0 ) was also conducted. The Jones model is not present in PhyML 3.0 so the JTT model of evolution (selected via ProtTest ) was used with 1000 bootstrap replicates to find the most likely tree: the JTT model is not present in MrBayes. Thus, identical models could not be used for each analysis.
Alignment of ZOC and ZF-NC domains
The C-terminal domain of full-length alignments derived from Muscle alignment (above) was trimmed from the third amino acid after the first cysteine residue of ZF1. The remaining N-terminal portions of the proteins were then realigned using Muscle and viewed using MacVector version 11.1.1 (MacVector, Inc, NC, USA).
RNA in situ hybridization and image analysis
We performed in situ hybridizations, as previously described, for Mnemiopsis leidyi , N. vectensis  and C. teleta . All images for M. leidyi and N. vectensis were acquired using a Zeiss Axioskop 2 in conjunction with the Axiocam HRc and Axio vision 4.7 software (Zeiss), Jena, Germany. C. teleta images were acquired using a Zeiss Axioskop 2 microscope and a SPOT Flex digital camera (Diagnostic Instruments, Inc, MI, USA) in conjunction with the Spot Advanced version 4.6 software (Diagnostic Instruments, Inc). Additional image processing was done with Helicon Focus software (Helicon Soft Ltd. Kharkov, Ukraine) for some images as indicated in figure legends. Additional images for N. vectensis and M. leidyi genes shown here are available in the Kahi Kai image database http://www.kahikai.org/.
- Ct-zic Capitella teleta zic:
hours post fertilization
polymerase chain reaction
ventral nerve cord
Zic 1-3 oddpaired conserved.
We would like to thank Sally Leys and Pamela Windsor for providing Ephydatia muelleri cDNA and RACE cDNA, and Andreas Hejnol for providing Convolutriloba longifissura RACE cDNA and help in obtaining Convolutriloba longifissura embryonic material. We would also like to thank Joseph F Ryan, Andreas D Baxevanis and the NIH/NISC Comparative Sequencing Program for allowing access to the draft sequences of the ctenophore Mnemiopsis leidyi genome. This work was made possible by funding from The National Institutes of Health, Eunice Kennedy Shriver National Institute of Child Health and Human Development F32HD055000 and the National Center for Research Resources G12RR003061.
- Aruga J, Kamiya A, Takahashi H, Fujimi TJ, Shimizu Y, Ohkawa K, Yazawa S, Umesono Y, Noguchi H, Shimizu T, Saitou N, Mikoshiba K, Sakaki Y, Agata K, Toyoda A: A wide-range phylogenetic analysis of Zic proteins: implications for correlations between protein structure conservations and body plan complexity. Genomics. 2006, 87: 783-792. 10.1016/j.ygeno.2006.02.011.View ArticlePubMedGoogle Scholar
- Aruga J: The role of Zic genes in neural development. Molec Cell Neurosci. 2004, 26: 205-221. 10.1016/j.mcn.2004.01.004.View ArticlePubMedGoogle Scholar
- Merzdorf CS: Emerging roles for zic g enes in early development. Devl Dynamics. 2007, 236: 922-940. 10.1002/dvdy.21098.View ArticleGoogle Scholar
- Mizugishi K, Hatayama M, Tohmonda T, Ogawa M, Inoue T, Mikoshiba K, Aruga J: Myogenic repressor I-mfa interferes with the function of of Zic family proteins. Biochem Biophysl Res Commun. 2004, 320: 233-240. 10.1016/j.bbrc.2004.05.158.View ArticleGoogle Scholar
- Lindgens D, Holstein TW, Technau U: Hyzic, the Hydra homolog of the zic/odd-paired gene, is involved in the early specification of the sensory nematocytes. Development. 2004, 131: 191-201. 10.1242/dev.00903.View ArticlePubMedGoogle Scholar
- Srivastava M, Simakov O, Chapman J, Fahey B, Gauthier MEA, Mitros T, Richards GS, Conaco C, Dacre M, Hellsten U, Larroux C, Putnam NH, Stanke M, Adamska M, Darling A, Degnan SM, Oakley TH, Plachetzki DC, Zhai Y, Adamski M, Calcino A, Cummins SF, Goodstein DM, Harris C, Jackson DJ, Leys SP, Shu S, Woodcroft BJ, Vervoort M, Kosik KS, Manning G, Degnan BM, Rokhsar DS: The Amphimedon queenslandica genome and the evolution of animal complexity. Nature. 2010, 466: 720-727. 10.1038/nature09201.PubMed CentralView ArticlePubMedGoogle Scholar
- Benedyk MJ, Mullen JR, DiNardo S: Odd-paired: a zinc finger pair-rule protein required for the timely activation of engrailed and wingless in Drosophila embryos. Genes Dev. 1994, 8: 105-117. 10.1101/gad.8.1.105.View ArticlePubMedGoogle Scholar
- Cimbora DM, Sakonju S: Drosophila midgut morphogenesis requires the function of the segmentation gene odd-paired. Dev Biol. 1995, 169: 580-595. 10.1006/dbio.1995.1171.View ArticlePubMedGoogle Scholar
- Skeath JB, Panganiban G, Selegue J, Carroll SB: Gene regulation in two dimensions: the proneural achaete and scute gnees are controlled by combinations of axis-patterning genes through a common intergenic control region. Genes Dev. 1992, 6: 2606-2619. 10.1101/gad.6.12b.2606.View ArticlePubMedGoogle Scholar
- Alper S, Kenyon C: The zinc finger protein REF-2 functions with the Hox genes to inhibit cell fusion in the ventral epidermis of C. elegans. Development. 2002, 129: 3335-3348.PubMedGoogle Scholar
- Bertrand V, Holber O: Linking asymmetric cell division to the terminal differentiation program of postmitotic neurons in C. elegans. Dev Cell. 2009, 16 (4): 563-575. 10.1016/j.devcel.2009.02.011.PubMed CentralView ArticlePubMedGoogle Scholar
- Takahashi H, Shimizu T, Aruga J: Expression pattern of the annelid Zic in embryonic development of the oligochaete Tubifex tubifex. Dev Genes Evol. 2008, 218: 553-560. 10.1007/s00427-008-0252-x.View ArticlePubMedGoogle Scholar
- Nagai T, Aruga J, Takada S, Gunther T, Sporle R, Schughard K, Mikoshiba K: The expression of the mouse Zic1, Zic2, and Zic3 gene suggests an essential role for Zic genes in body pattern formation. Dev Biol. 1997, 182: 299-313. 10.1006/dbio.1996.8449.View ArticlePubMedGoogle Scholar
- Gostling NJ, Shimeld SM: Protochordate Zic genes define primitive somite compartments and highlight molecular changes underlying neural crest evolution. Evol Dev. 2003, 5 (2): 136-144. 10.1046/j.1525-142X.2003.03020.x.View ArticlePubMedGoogle Scholar
- Wada S, Saiga H: Hrzic N, a new Zic family gene of ascidians, plays essential roles in the neural tube and notochord development. Development. 2002, 129: 5597-5608. 10.1242/dev.00156.View ArticlePubMedGoogle Scholar
- Imai KS, Satou Y, Satoh N: Multiple functions of a Zic-like gene in the differentiation of notochord, central nervous system and muscle in Ciona savignyi embryos. Development. 2002, 129: 2723-2732.PubMedGoogle Scholar
- Nishida H, Sawada K: macho-1 encodes a localized mRNA in ascidian eggs that specifies muscle fate during embryogenesis. Nature. 2001, 409: 724-729. 10.1038/35055568.View ArticlePubMedGoogle Scholar
- Gyoja F: Expression of a muscle determinant gene, macho-1, in the anural ascidian Molgula tectiformis. Dev Genes Evol. 2006, 216: 285-289. 10.1007/s00427-005-0056-1.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Meyer NP, Seaver EC: Neurogenesis in an annelid: characterization of brain neural precursors in the polychaete Capitella sp. I. Dev Biol. 2009, 335: 237-252. 10.1016/j.ydbio.2009.06.017.View ArticlePubMedGoogle Scholar
- Meyer NP, Boyle MJ, Martindale MQ, Seaver EC: A comprehensive fate map by intracellular injection of identified blastomeres in the marine polychaete Capitella teleta. EvoDevo. 2010, 1: 8-10.1186/2041-9139-1-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Dill KK, Thamm K, Seaver EC: Characterization of twist and snail gene expression during mesoderm and nervous system development in the polychaete annelid Capitella sp. I. Dev Genes Evol. 2007, 217: 435-447. 10.1007/s00427-007-0153-4.View ArticlePubMedGoogle Scholar
- Hernandez-Nicaise ML: Ctenophora. Microscopic Anatomy of Invertebrates, Placaozoa, Porifera, Cnidaria and Ctenophora. Edited by: Herrison FW, Westfall JA. 1991, New York: Wiley-Liss, 2: 359-418.Google Scholar
- Horridge GA: Presumed photoreceptive cilia in a ctenophore. Q J Micr Sci. 1964, 105: 311-317.Google Scholar
- Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, Sorensen MV, Haddock SHD, Schmidt Rhaesa A, Okusu A, Kristensen RM, Wheeler WC, Martindale MQ, Giribet G: Broad phylogenetic sampling improves resolution of the animal tree of life. Nature. 2008, 452: 745-749. 10.1038/nature06614.View ArticlePubMedGoogle Scholar
- Hejnol A, Obst M, Stamatakis A, Ott M, Rouse GW, Edgecombe GD, Martinez P, Baguna J, Bailly X, Jondelius U, Wiens M, Werner E, Muller G, Seaver E, Wheeler WC, Martindale MQ, Giribet G, Dunn CW: Assessing the root of the bilaterian animals with scalable phylogenomic methods. Proc Royal Soc Bioll Sci. 2009, 476: 4261-4270. 10.1098/rspb.2009.0896.View ArticleGoogle Scholar
- Pick KS, Philippe H, Schreiber F, Erpenbeck D, Jackson DJ, Wrede P, Wiens M, Alie A, Morgenstern B, Manuel M, Worheide G: Improved phylogenomic taxon sampling noticeably affects non-bilaterian relationships. Molec Biol Evol. 2010, 27: 1983-87. 10.1093/molbev/msq089.PubMed CentralView ArticlePubMedGoogle Scholar
- Collins AG: Evaluating multiple alternative hypotheses for the origin of Bilateria: An analysis of 18 S rRNA molecular evidence. PNAS. 1998, 95: 15458-15463. 10.1073/pnas.95.26.15458.PubMed CentralView ArticlePubMedGoogle Scholar
- Philippe H, Derelle R, Lopez P, Pick K, Borchiellini C, Boury-Esnault N, Vacelet J, Renard E, Houliston E, Queinnec E, Da Silva C, Winicker P, Le Guayader H, Leys S, Jackson DJ, Schreiber F, Erpenbeck D, Morgenstern B, Woerheide G: Phylogeomics revives traditional views on deep animal relationships. Curr Biol. 2009, 19: 706-712. 10.1016/j.cub.2009.02.052.View ArticlePubMedGoogle Scholar
- Srivastava M, Begovic E, Chapman J, Putnam NH, Hellsten U, Kawashima T, Kuo A, Mitros T, Salamov A, Carpenter ML, Signorovitch AY, Moreno MA, Kamm K, Grimmwood J, Schmutz J, Shapiro H, Grigoriev IV, Buss LW, Schierwater B, Dellaporta SL, Rokhsar DS: The Trichoplax genome and the nature of placozoans. Nature. 2008, 454: 955-960. 10.1038/nature07191.View ArticlePubMedGoogle Scholar
- Lamar E, Kinter C, Goulding M: Identification of NKL, a novel Gli-Kruppel zinc-finger protein that promotes neural differentiation. Development. 2001, 128: 1335-1346.PubMedGoogle Scholar
- Kim YS, Lewandoski M, Perantoni AO, Kurebay S, Nankanishi G, Jetten AM: Identification of Glis1, a novel Gli-related, Kruppel-like zinc finger protein containing transactivation and repressor functions. J Biol Chem. 2002, 277: 30901-30913. 10.1074/jbc.M203563200.View ArticlePubMedGoogle Scholar
- Shimeld SM: C2H2 zinc finger genes of the Gli, Zic, KLF, SP, Wilms' tumour, Huckebein, Smail, Ovo, Spalt, Odd, Blimp-1, Fez and related families from Branchiostoma floridae. Dev Genes Evol. 2008, 218: 639-649. 10.1007/s00427-008-0248-6.View ArticlePubMedGoogle Scholar
- Marlow HQ, Srivastava M, Matus DQ, Rokhsar D, Martindale MQ: Anatomy and development of the nervous system of Nematostella vectensis, and Anthozoan Cnidarian. Dev Neurobiol. 2009, 69: 235-254. 10.1002/dneu.20698.View ArticlePubMedGoogle Scholar
- Materna SC, Howard-Ashby M, Gray RF, Davidson EH: The C2H2 zinc finger genes of Strongylocentrotus purpuratus and their expression in embryonic development. Dev Biol. 2006, 300: 108-120. 10.1016/j.ydbio.2006.08.032.View ArticlePubMedGoogle Scholar
- Nusslein-Volhard C, Wieschaus E: Mutations affecting segment number and polarity in Drosophila. Nature. 1980, 287: 795-801. 10.1038/287795a0.View ArticlePubMedGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.PubMed CentralView ArticlePubMedGoogle Scholar
- Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.View ArticlePubMedGoogle Scholar
- Pang K, Martindale MQ: Developmental expression of homeobox genes in the ctenophore Mnemiopsis leidyi. Dev Genes Evol. 2008, 18: 307-319. 10.1007/s00427-008-0222-3.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.