- Open Access
Lim homeobox genes in the Ctenophore Mnemiopsis leidyi: the evolution of neural cell type specification
EvoDevo volume 3, Article number: 2 (2012)
Nervous systems are thought to be important to the evolutionary success and diversification of metazoans, yet little is known about the origin of simple nervous systems at the base of the animal tree. Recent data suggest that ctenophores, a group of macroscopic pelagic marine invertebrates, are the most ancient group of animals that possess a definitive nervous system consisting of a distributed nerve net and an apical statocyst. This study reports on details of the evolution of the neural cell type specifying transcription factor family of LIM homeobox containing genes (Lhx), which have highly conserved functions in neural specification in bilaterian animals.
Using next generation sequencing, the first draft of the genome of the ctenophore Mnemiopsis leidyi has been generated. The Lhx genes in all animals are represented by seven subfamilies (Lhx1/5, Lhx3/4, Lmx, Islet, Lhx2/9, Lhx6/8, and LMO) of which four were found to be represented in the ctenophore lineage (Lhx1/5, Lhx3/4, Lmx, and Islet). Interestingly, the ctenophore Lhx gene complement is more similar to the sponge complement (sponges do not possess neurons) than to either the cnidarian-bilaterian or placozoan Lhx complements. Using whole mount in situ hybridization, the Lhx gene expression patterns were examined and found to be expressed around the blastopore and in cells that give rise to the apical organ and putative neural sensory cells.
This research gives us a first look at neural cell type specification in the ctenophore M. leidyi. Within M. leidyi, Lhx genes are expressed in overlapping domains within proposed neural cellular and sensory cell territories. These data suggest that Lhx genes likely played a conserved role in the patterning of sensory cells in the ancestor of sponges and ctenophores, and may provide a link to the expression of Lhx orthologs in sponge larval photoreceptive cells. Lhx genes were later co-opted into patterning more diversified complements of neural and non-neural cell types in later evolving animals.
LIM Homeobox (Lhx) genes were first isolated from the nematode, Caenorhabditis elegans, where the Lhx homolog, MEC-3, was shown to be required for the proper differentiation of touch receptor neurons . Subsequent studies in C. elegans and rat isolated LIN-11 and Islet1, respectively, which together with MEC-3 are the founding members for the LIM family acronym [2, 3]. Phylogenetically, Lhx genes were originally subdivided into six subfamilies, Lhx1/5, Lhx2/9, Lhx3/4, Lhx6/8, Islet, and Lmx . Lhx proteins are composed of tandem zinc-finger LIM domains at the N-terminus, which function by binding specific co-factors that mediate their function, while the helix-turn-helix homeodomain (HD) interacts with DNA in a sequence-specific manner . While LIM domains and homeodomains are found in non-metazoan eukaryotes, the specific combination of LIM-LIM-HD is only found in animals .
Lhx genes have roles in cell specification, tissue differentiation and neural patterning. In both vertebrate and invertebrate taxa, Lhx genes have conserved roles in the patterning of sensory neurons, interneurons, and motor neurons (reviewed in [4, 7]). It has been suggested that cells expressing different combinations of Lhx genes form a "LIM code" that is important in specifying cell types within a tissue or organ , and motor neuron axon pathway finding .
The C. elegans Lhx1/5 genes, LIN-11 and MEC-3, are required for the terminal differentiation of non-overlapping sensory, motor neurons and interneurons [1, 9, 10]. Mouse null mutants of the Lhx1 family die during mid-gastrulation due to massive head defects , while the Lhx5 homologue is expressed in the anterior neural plate and parts of the developing diencephalon . At later stages, the expression pattern of Lhx5 extends to parts of the midbrain, hindbrain and spinal cord . A study (reviewed in  on the LIM code for axon path finding in motor neuron subtypes revealed that depending on the combination of Lhx genes expressed in motor neurons, different motor neuron subtypes are produced. Lhx1 when expressed with Islet2 produces motor neurons that project into ventral limb bud musculature, while Islet1, Islet2, and Lhx3 expressing motor neurons project into the medial motor column, and different motor neurons are produced when Isl1 and Isl2 are expressed either together or separately. In addition to neural patterning, Lhx1 is also known for its role in blastoporal organizer activity during gastrulation in Xenopus, ascidians, amphioxus, and cnidarians  as well as endoderm specification in ascidians , amphioxus  and mice .
The Lhx2/9 group, also known as the apterous group (named after the Drosophila gene), has diverse roles in patterning the nervous system, wing development, muscle development, axon guidance, and neurotransmitter choice [18, 19]. The C. elegans Lhx2/9 homologue ttx-3 is expressed in a pair of interneurons and anterior muscle cells that project into the nerve ring . The vertebrate Lhx2/9 homologues are expressed in the nerve cords, eyes, olfactory organs and limbs .
The Lhx3/4 subgroup expression patterns were found to be restricted to post mitotic neurons in Drosophila. The function of Lhx3/4 genes in Drosophila is in axon guidance of specific motor neurons when combined with the expression of Lhx gene islet . In the developing chick embryo, Lhx3/4 has been shown to specify the formation of interneurons in the presence of one of its binding partners Ldb1 (LIM domain-binding protein 1), while additional expression within the same cells with the Lhx gene Islet1 specifies the formation of motor neurons , which show that Lhx genes function in a combinatorial fashion called the 'LIM code' to specify specific neural identities [4, 8, 23].
The Drosophila Lhx6/8 homologue arrowhead is expressed in the nervous system and has yet to be functionally analyzed . The C. elegans homologue is expressed in sensory, motor and interneurons in the brain . The vertebrate members of this subfamily are expressed in the forebrain and first branchial arch .
The islet gene subfamily is expressed in mesodermal derived cells and in subsets of motor and interneurons of the central nervous system (CNS) of Drosophila . In vertebrates, Islet paralogs are expressed in many tissues, including heart, liver, pancreas, brain, and eyes [28–30].
The C. elegans homologue to the group Lmx is expressed in post-mitotic neurons controlling axon guidance and the synthesis of the neurotransmitter GABA. It is also expressed in endothelial cells of the uterus and in the excretory system . The vertebrate chick homologue Lmx-1b is involved in dorsal ventral patterning of limbs and in patterning the otic vesicle .
Novel neural cell types arose at some point during the evolution of metazoans, as nervous systems evolved and diversified throughout the bilaterian radiation. Nested between single-celled choanoflagellate-like ancestors (Figure 1) and the Bilaterian radiation of metazoans are four groups: Porifera, Ctenophora, Cnidaria, and Placozoa. While ctenophores and cnidarians have distinct nerve cells and other neural structures, poriferans and placozoans do not [33, 34]. Interestingly, the genomes of the sponge, Amphimedon queenslandica, and the placozoan, Trichoplax adhaerens, have revealed the presence of most of the genes involved in forming the post-synaptic scaffold, as well as neurotransmitter biosynthesis [33–35]. Even choanoflagellates, the closest extant sister group of the Metazoa, possess genes for many of the post-synaptic scaffolding proteins, showing that many of the components to build a neuron likely predated its evolution . Thus, it is of some interest to determine how the molecular components were assembled in early branching metazoans to form functionally integrated nervous systems.
Cnidarians are now well accepted to be the sister group to all bilaterians [37–41]. The anthozoan cnidarian, Nematostella vectensis, possesses members of all six Lhx subfamilies, as does the placozoan, Trichoplax, suggesting that complete diversification already occurred in the ParaHoxozoan common ancestor . The sponge, Amphimedon, only has members of three subfamilies: Lhx1/5, Lhx3/4, and Islet. Amphimedon Lhx1/5 and Lhx3/4 were found to be expressed within putative photosensory cells surrounding larval pigment cells, while Islet is ubiquitously expressed . In Trichoplax, all Lhx genes were shown to be expressed in adult animals through the use of RT-PCR; however, spatio-temporal patterns were not examined . Expression patterns of NvLhx6/8, NvLhx1/5, and NvLmx in Nematostella overlap with previously described neural territories . Interestingly, while the Amphimedon and Nematostella Lhx genes are all on separate genomic scaffolds, a single Trichoplax scaffold contains a cluster of Lmx, Lhx3/4, Lhx6/8, a LIM-only gene, and Lhx2/9 further downstream, suggesting that the Lhx family arose via tandem duplications .
Ctenophores occupy a highly debated phylogenetic position, once grouped with cnidarians forming the group Coelenterata [42, 43]. Ctenophores are now thought to have diverged prior to cnidarians and may actually be the earliest extant metazoan phyla that possess definitive neurons [37, 38]. Comparisons of gene complements between the four basal metazoans will be helpful in giving support for one or more of the various competing phylogenetic hypotheses (reviewed in ). Recent gene content studies [45–47] give consistent support for Porifera and Ctenophora diverging prior to the ParaHoxozoa (Placozoa, Cnidaria, Bilateria), however identifying the earliest branching taxon remains problematic.
For this study the genomic complement of Lhx genes was examined in the lobate ctenophore M. leidyi. Using the recently sequenced genome of M. leidyi, four Lhx genes were predicted by preliminary searches . We further analyzed these genes by obtaining their full sequences through RACE PCR. Using whole mount in situ hybridization, we examined the expression patterns of the four Lhx genes throughout development. These genes are expressed in discreet sensory cell types, and in an overlapping fashion within the apical sensory organ, a highly innervated nervous structure.
Gene identification, structure and genomic organization
Partial sequence information from the four M. leidyi Lhx homeobox genes previously reported  were used to design non-degenerate primers. Full-length cDNA transcripts for all four Mnemiopsis Lhx genes were isolated by 5' and 3' RACE RT-PCR from mixed-stage embryonic cDNA. The four sequences have been submitted to GenBank: MlLhx1/5 [Genbank: JF912807], MlLhx3/4 [Genbank: JF912808], MlIslet [Genbank: JF912806], and MlLmx [Genbank: JF912809]. All of the Lhx genes isolated contained two tandem N-terminal LIM domains followed by a homeodomain sequence as identified by the SMART domain prediction  (Figure 2).
The genomic organization of the Mnemiopsis Lhx genes showed no apparent linkage groups and they were found to be located on separate genomic scaffolds. The size of the scaffolds is as follows: MlLhx1/5 is found on scaffold ML1325 with a length of 127,187 bases. MlLmx is found on scaffold ML0569 with a length of 734,417 bases. MlIslet is found on scaffold ML0530 with a length of 126,527 bases. MlLhx3/4 is found on scaffold ML0681 with a length of 490,590 bases.
In bilaterians, there are two conserved introns in the Lmx genes, within the homeodomain sequence. Both conserved intronic breaks are present in Trichoplax  a non-bilaterian species; however, the first intron position is slightly translocated in Mnemiopsis (Figure 3), and the second intron, while conserved in many other taxa, is not conserved in Nematostella. The homeobox sequences of Lhx1/5, Lhx3/4, and Islet were also surveyed for conserved intron positions (Additional file 1). Although introns were found within the homeodomain sequences of Lhx1/5, Lhx3/4, and Islet, they did not appear to have any conservation with any of the other taxa sampled.
Phylogenetic analyses of the isolated ctenophore Lim homeodomain-containing genes were performed using sequence data from the Lhx complement of the cnidarian N. vectensis , the placozoan Trichoplax adherens , the sponge. A. queenslandica , and from the published genomes of Homo sapiens, Danio rerio, Gallus gallus and Drosophila melanogaster. Maximum likelihood and Bayesian analysis were performed individually on the aligned tandem LIM domains and the homeodomain (Figure 4). Similar to analyses using only the homeodomains , our analyses show relatively high support for the four Mnemiopsis Lhx genes within the Lhx1/5, Lhx3/4, Islet, and Lmx subclasses. The sponge A. queenslandica lacks an Lmx gene homologue, but due to low branching support between subclasses, we cannot determine whether the demosponge lineage lost the Lmx gene, or if this lineage diverged before Lmx arose. These predictions need to be tested with other sponge genomes.
Included in this phylogenetic study are three LIM domain-containing gene families that lack a homeodomain and are the most closely related to Lhx genes, which were used for outgrouping: nuclear LIM only (LMO), Actin binding LIM (Ablim), and Paxilin. The most closely related family to the Lhx genes are found within the nuclear LIM only (LMO) gene family . Both M. leidyi and A. queenslandica do not possess genes within this family, and neither do the unicellular non-metazoans Monosiga brevicollis and Capsaspora owczarzaki surveyed in this study. The absence of LMO genes in these taxa in addition to the current genomic and phylogenetic data available and their internal branching within the Lhx hierarchy in our phylogenetic tree (although with low support values) suggests that these genes are likely a subfamily of Lhx genes and arose through a gene duplication shortly thereafter losing their homeodomains after the split of ctenophores and sponges and before the rise of ParaHoxozoa.
The four M. leidyi Lhx genes correspond closely to the A. queenslandica complement in gene content. In addition we also found that there are no homologues to the LIM-domain binding factors (Ldb) in the currently available sponge and ctenophore genomes, although this gene is present in all other animals. These data suggest that the ancestor of sponges and ctenophores contained three or four of the seven subclasses of Lhx genes, with Parahoxozoans acquiring the two closely related subclasses (Lhx 2/9, Lhx 6/8), Ldb factors, and the LMO gene family (derived from an already existing Lhx gene family) after their divergence from sponges and ctenophores (Table 1).
Developmental expression of M. leidyi Lhx genes
To determine potential roles for Lhx genes in M. leidyi development, we looked at the spatial and temporal expression of all four Lhx genes through whole mount in situ hybridization. MlIslet is the earliest gene detected, at gastrulation or approximately four hours post-fertilization (hpf), in two to three rows of cells at the aboral pole along the sagittal axis (Figure 5A). This expression continues through development (Figure 5B-E), eventually forming the polar fields, as well as the most sagittal regions of the apical organ floor. The other Lhx genes are detected just after gastrulation, at approximately 5 hpf. MlLhx1/5 is expressed in a subset of mesodermal micromeres born at the oral pole (site of gastrulation) that enter the blastocoel, as well as in cells surrounding the blastopore (Figure 5I-J). These mesodermal oral micromere derivatives proliferate and migrate to areas underlying the comb plates (Figure 5K-N) to form the presumptive photocytes, or light-producing cells. The blastoporal expression continues through development, with expression later confined to individual cells of the pharynx. At the cydippid stage, there is an additional expression domain in the floor of the apical organ (Figure 5O-P). These four groups of cells possibly correspond to the putative photoreceptor cells .
MlLhx3/4 is expressed in cells at the aboral pole, eventually being confined to the central part of the apical organ floor (Figure 5Q-X). MlLmx is expressed in two groups of cells at the aboral pole just after gastrulation, along the tentacular plane (Figure 5Y). These cells eventually form part of the tentacle bulbs (Figure 5C'-F'). At 12 to 14 hpf, there is an additional expression domain in the ectoderm of the aboral pole in four groups of cells (Figure 5C', E'), which later appear to overlap with MlLhx1/5 in the region of the putative photoreceptors in the floor of the apical organ (Figure 5O).
In summary, all four Lhx genes are expressed in the apical organ at the cydippid stage, in addition to other regions (Figure 6). MlLhx1/5, MlLmx, and MlIslet are expressed in an overlapping region, while MlLhx3/4 is expressed in other gene-specific regions of the apical organ. The other regions of expression include the polar fields (MlIslet), photocytes and pharynx (MlLhx1/5), and tentacle bulb apparatus (MlLmx).
Lhx gene complements
By probing the full genome of M. leidyi, we find that it has a genetic complement more similar to that of the sponge lineage than to the complete repertoire of Lhx gene subfamilies found in placozoans, cnidarians, and bilaterians. These data provide further evidence of the basal position of ctenophore and sponge lineages relative to Parahoxazoans. While multiple independent gene loss in both ctenophore and sponge clades may be the cause for this apparent synapomorphy, growing arguments based on genomic content, are suggesting a close link between sponges and ctenophores, including the existence of homeodomain containing transcription factor complement , nuclear receptor repertoire , Wnt signaling pathway , and TGF beta pathway .
Mnemiopsis expression data
Lhx genes are expressed in regions which contain the highest concentration of neural elements in the cnidarian Nematostella, including the apical tuft and nerve rings around the mouth and in the pharynx . In the sponge Amphimedon, there are broader expression domains; however, there are overlapping domains associated with the photosensory pigment ring, suggesting that Lhx genes are playing a combinatorial role in neural or sensory cell fate specification across the Metazoa, including the basal metazoans . The expression data here suggest Lhx genes are deployed in a similar fashion in the ctenophore Mnemiopsis. There are both overlapping and gene-specific expression domains associated with the apical sensory organ, a highly innervated structure involved in gravity, light, and pressure detection . MlLhx1/5, MlLmx and MlIslet have overlapping expression in the putative photoreceptor cells (Figure 6), while MlLhx3/4 and MlIslet are also expressed in other parts of the apical organ. It is possible that a 'LIM code' is involved in specifying different neural territories or cell types of the apical organ. Although the molecular interactions of only Lhx3/4 and Islet have been studied in depth , it is not unlikely that other Lhx genes may be able to form complexes to differentially regulate gene expression patterns. The expression of Lhx genes also forms distinct boundaries within the apical organ and associated polar fields, which may contribute to regional specificity, axonal projection boundaries, or neural transmitter phenotypes of these neural structures as is seen in bilaterians .
Previous studies have described other genes in different regions of the apical organ of Mnemiopsis, including the T-box genes Tbx2/3 and Brachyury , the homeobox genes Prd1 and Prd3 , the Wnt genes MlWnt6 and MlWntX , as well as the zinc fingers MlGli and MlGlis . Histological and electron microscopic studies of the ctenophore apical organ have suggested it is a highly innervated sensory structure [58–71]; however, gene expression studies suggest that it is much more complex than previously thought. These data are supported by recent work in the ctenophore, Pleurobrachia pileus; immunohistochemical studies have shown that there is a highly complex nerve net underlying the apical organ, as well as distinct groups of neurons within the apical organ floor that belie its morphological simplicity . Using this information about the neural anatomy of ctenophores, we suggest that ctenophore Lhx genes are expressed in a combinatorial fashion in the developing sensory cells associated with the apical organ.
In addition to the apical organ, Lhx expression domains are also seen in putative neurosensory regions such as the polar fields expression of MlIslet. While these ciliated structures have been identified morphologically and histologically [59, 72], their function has remained largely unknown. They have been proposed to be olfactory organs [58, 73–75], however, there is little evidence to confirm this. The recent immunohistological study of P. pileus by Jager et al.  identified a special nerve net that extended throughout the polar fields, as well as extensive neuro-sensory structures, termed 'Z bodies.' The gene expression of MlIslet overlaps with the Z bodies, but is not specific to them; instead the expression domain encompasses the entirety of the polar fields.
The expression of MlLhx1/5 in the presumptive photocytes is also quite intriguing. Ctenophores are capable of bioluminescence via calcium-activated photoproteins , but the mechanism or function has yet to be elucidated. The photocytes are located around the meridional canals underlying the comb rows . In addition to photoprotein genes, these cells also express an opsin gene , which suggests that these cells may represent an ancestral neuroeffector cell type. They express a gene that senses the environment (light-sensing via opsin), and an effector (photoprotein, generating light production) preceded by a transcription factor (MlLhx1/5), such that these photocytes appear to be capable of both sensing as well as responding to stimuli.
MlLhx1/5 expression is also found around the blastopore and is regionally localized to a portion of the pharynx in later stages. Considering that Lhx1/5 is expressed around the blastopore in many bilaterians, as well as the cnidarian Nematostella, we can speculate that besides neural patterning, Lhx1/5 was an ancestral 'blastoporal' gene. However, whether it functions in organizing activity in ctenophores, as it does in bilaterians, remains to be determined.
Recent gene content studies [45–47] give support for Porifera and Ctenophora diverging prior to the ParaHoxozoa (Placozoa, Cnidaria, Bilaterian) during early animal evolution. This is further supported by the present study, in which we have further classified the genomic complement of M. leidyi Lhx genes. Comprehensive phylogenetic analyses of LIM genes  in addition to our Lhx gene phylogenetic analyses of cnidarians, ctenophores, sponges and placozoans, indicate that the ancestor of ctenophores and sponges likely had three or four Lhx gene subfamilies: Lhx1/5, Islet, Lhx 3/4, and Lmx. The common ancestor of ParaHoxozoa contained three more Lhx gene subclasses, the LIM domain binding proteins, and the gene family LMO which was formed by duplications from a pre-existing Lhx gene subclass. Mnemiopsis Lhx genes are expressed in both overlapping and non-overlapping domains within proposed neural and sensory cell territories. These data suggest that Lhx genes first likely played a role in the patterning of sensory cells, as is seen in sponge larval photoreceptive cells . In other animal lineages with more complex nervous systems, Lhx genes appear to have a conserved role in patterning neural and sensory cell types. This apparent link shows support for the hypothesis that neural cell types may have evolved from non-neural sensory cells. An alternate hypothesis may be that sponges lost components of a simple nervous system after the split with ctenophores, or the ctenophore nervous system arose independently as a case of convergent evolution. These competing hypotheses still need to be considered due to the unexpected complex repertoire of structural and patterning nervous system genes found within sponges [34, 35]. Future genomic, expression, and functional studies are needed to characterize the molecular nature of the ctenophore nervous system.
Animal collection, RNA extraction and cDNA synthesis
M. leidyi adults were collected from Eel Pond or the NOAA Rock Jetty in Woods Hole, MA, during the months of June and July and spawned as previously described . RNA was extracted from embryos at regular intervals from fertilization to 36 hpf using TRI Reagent (Molecular Research Center, Cincinnati, OH, USA) . RNA was reverse transcribed to generate cDNA using the SMARTer RACE cDNA Amplification Kit (BD Biosciences, San Jose, CA, USA).
Identification of LIM homeobox genes in ctenophores
The Mnemiopsis genome was scanned in silico for genes of interest using a reciprocal BLAST approach. Human, Drosophila, and Nematostella orthologs were used as queries for TBLASTN searches. Candidate matches were then used in BLASTP searches of the human genome to find the closest hit. If the closest match was not the original ortholog or the E-value was greater than 0.001, then it was coded as being absent from the genome. For all genes of interest, gene-specific primers were designed for RACE PCR (MacVector, Cary, NC, USA). RACE PCR fragments were then conceptually spliced and aligned back to genomic contigs for comparison of intron-exon boundaries using Sequencher (Gene Codes, Ann Arbor, MI, USA). The following genes were isolated and fully sequenced: MlLhx1/5 [Genbank: JF912807], MlLhx3/4 [Genbank: JF912808], MlIslet [Genbank: JF912806], and MlLmx [Genbank: JF912809].
Whole mount in-situ hybridization
Embryos were fixed at various stages from freshly collected uncleaved embryos (0 hpf) to cydippids (24 to 36 hpf) . They were stored in methanol at -20°C until ready to use. Digoxygenin-labeled M. leidyi Lhx gene riboprobes (Ambion/Applied Biosystems, Inc., Foster City, CA, USA of the following sizes: MlIslet 1200 bp, MlLmx 1100 bp, MlLhx1/5 850 bp, and MlLhx3/4 1500 bp, were hybridized for 48 hours at 60°C at 0.1 ng/ul and detected using an alkaline phosphatase conjugated antibody (Roche Applied Science Inc, Indianapolis, IN, USA) and the colorimetric substrate NBT/BCIP (Roche Applied Science Inc, Indianapolis, IN, USA) Following detection, specimens were washed with phosphate-buffered saline and transferred through a glycerol series up to 70% glycerol. They were then mounted on glass slides, viewed under an Axioskop 2 compound microscope, and imaged using an AxioCam HRc with Axiovision software (Zeiss Inc, Jena, Germany). Color balance and brightness were adjusted using Adobe Photoshop CS3. The only modification to the in situ protocol is a change in acetic anhydride treatment (treated in 0.1 M triethanolamine rather than 1% w/v). For the most recently updated protocols, contact the authors. All in situ images presented here, as well as additional developmental stages and/or views, are available online via the comparative gene expression database, Kahikai http://www.kahikai.com.
M. leidyi Lhx protein sequences were aligned to the complements of Lhx proteins of the various species: H. sapiens, Drosophila melanogaster, N. vectensis, T. adhaerens, and A. queenslandica, and outgroup sequences of LIM-only domain containing gene 4 (LMO4) in addition to the above species when present, G. gallus, D. rerio, and the related unicellular choanoflagellates M. brevicollis, and C. owczarzaki LIM domain-containing genes paxilin and ablim (full sequence data provided in Additional file 2, using MUSCLE (default parameters) and trimmed by eye in Jalview  to include the conserved LIM domains and Homeodomains. Construction of the Lhx gene orthology was conducted in RaxML v7.0.0 , using maximum likelihood analysis under the JTT I+G model of evolution determined by ProtTest v1.4 analysis . A total of 1,000 searches was performed and 500 bootstrap replicates were applied to the tree with the best likelihood score used to generate branch support values. Bayesian phylogenetic analyses were also performed with MrBayes 3.1  using the JTT I+G model with two runs of 2,000,000 generations sampled every 100 generations. The first 5,000 trees were disregarded as burn-in.
actin binding LIM
Basic Local Alignment Search Tool
central nervous system
hours post fertilization
- Ldb :
LIM domain-binding protein
LIM homeobox containing genes
LIM-only domain containing gene
LIM-only domain containing gene 4
Placozoa: Cnidaria: Bilateria
polymerase chain reaction
randomized accelerated ML.
Way JC, Chalfie M: mec-3, a homeobox-containing gene that specifies differentiation of the touch receptor neurons in C. elegans. Cell. 1988, 54: 5-16. 10.1016/0092-8674(88)90174-2.
Freyd G, Kim SK, Horvitz HR: Novel cysteine-rich motif and homeodomain in the product of the Caenorhabditis elegans cell lineage gene lin-11. Nature. 1990, 344: 876-879. 10.1038/344876a0.
Karlsson O, Thor S, Norberg T, Ohlsson H, Edlund T: Insulin gene enhancer binding protein Isl-1 is a member of a novel class of proteins containing both a homeo- and a Cys-His domain. Nature. 1990, 344: 879-882. 10.1038/344879a0.
Hobert O, Westphal H: Functions of LIM-homeobox genes. Trends Genet. 2000, 16: 75-83. 10.1016/S0168-9525(99)01883-1.
Shirasaki R, Pfaff SL: Transcriptional codes and the control of neuronal identity. Annu Rev Neurosci. 2002, 25: 251-281. 10.1146/annurev.neuro.25.112701.142916.
Srivastava M, Larroux C, Lu DR, Mohanty K, Chapman J, Degnan BM, Rokhsar DS: Early evolution of the LIM homeobox gene family. BMC Biol. 2010, 8: 4-10.1186/1741-7007-8-4.
Kadrmas JL, Beckerle MC: The LIM domain: from the cytoskeleton to the nucleus. Nat Rev Mol Cell Biol. 2004, 5: 920-931. 10.1038/nrm1499.
Gill GN: Decoding the LIM development code. Trans Am Clin Climatol Assoc. 2003, 114: 179-189.
Way JC, Chalfie M: The mec-3 gene of Caenorhabditis elegans requires its own product for maintained expression and is expressed in three neuronal cell types. Genes Dev. 1989, 3: 1823-1833. 10.1101/gad.3.12a.1823.
Hobert O, D'Alberti T, Liu Y, Ruvkun G: Control of neural development and function in a thermoregulatory network by the LIM homeobox gene lin-11. J Neurosci. 1998, 18: 2084-2096.
Shawlot W, Behringer RR: Requirement for Lim1 in head-organizer function. Nature. 1995, 374: 425-430. 10.1038/374425a0.
Sheng HZ, Bertuzzi S, Chiang C, Shawlot W, Taira M, Dawid I, Westphal H: Expression of murine Lhx5 suggests a role in specifying the forebrain. Dev Dyn. 1997, 208: 266-277. 10.1002/(SICI)1097-0177(199702)208:2<266::AID-AJA13>3.0.CO;2-1.
Zhao Y, Sheng HZ, Amini R, Grinberg A, Lee E, Huang S, Taira M, Westphal H: Control of hippocampal morphogenesis and neuronal differentiation by the LIM homeobox gene Lhx5. Science. 1999, 284: 1155-1158. 10.1126/science.284.5417.1155.
Yasuoka Y, Kobayashi M, Kurokawa D, Akasaka K, Saiga H, Taira M: Evolutionary origins of blastoporal expression and organizer activity of the vertebrate gastrula organizer gene lhx1 and its ancient metazoan paralog lhx3. Development. 2009, 136: 2005-14. 10.1242/dev.028530.
Satou Y, Imai KS, Satoh N: Early embryonic expression of a LIM-homeobox gene Cs-lhx3 is downstream of beta-catenin and responsible for the endoderm differentiation in Ciona savignyi embryos. Development. 2001, 128: 3559-3570.
Wang Y, Zhang PJ, Yasui K, Saiga H: Expression of Bblhx3, a LIM-homeobox gene, in the development of amphioxus Branchiostoma belcheri tsingtauense. Mech Dev. 2002, 117: 315-319. 10.1016/S0925-4773(02)00197-1.
Perea-Gomez A, Shawlot W, Sasaki H, Behringer RR, Ang S: HNF3beta and Lim1 interact in the visceral endoderm to regulate primitive streak formation and anterior-posterior polarity in the mouse embryo. Development. 1999, 126: 4499-4511.
Benveniste RJ, Thor S, Thomas JB, Taghert PH: Cell type-specific regulation of the Drosophila FMRF-NH2 neuropeptide gene by Apterous, a LIM homeodomain transcription factor. Development. 1998, 125: 4757-4765.
Lundgren SE, Callahan CA, Thor S, Thomas JB: Control of neuronal pathway selection by the Drosophila LIM homeodomain gene apterous. Development. 1995, 121: 1769-1773.
Hobert O, Mori I, Yamashita Y, Honda H, Ohshima Y, Liu Y, Ruvkun G: Regulation of interneuron function in the C. elegans thermoregulatory pathway by the ttx-3 LIM homeobox gene. Neuron. 1997, 19: 345-357. 10.1016/S0896-6273(00)80944-7.
Thor S, Andersson SG, Tomlinson A, Thomas JB: A LIM-homeodomain combinatorial code for motor-neuron pathway selection. Nature. 1999, 397: 76-80. 10.1038/16275.
Bhati M, Lee C, Nancarrow AL, Lee M, Craig VJ, Bach I, Guss JM, Mackay JP, Matthews JM: Implementing the LIM code: the structural basis for cell type-specific assembly of LIM-homeodomain complexes. EMBO J. 2008, 27: 2018-29. 10.1038/emboj.2008.123.
Tosney KW, Hotary KB, Lance-Jones C: Specifying the target identity of motoneurons. Bioessays. 1995, 17: 379-382. 10.1002/bies.950170503.
Curtiss J, Heilig JS: Establishment of Drosophila imaginal precursor cells is controlled by the Arrowhead gene. Development. 1995, 121: 3819-3828.
Sagasti A, Hobert O, Troemel ER, Ruvkun G, Bargmann CI: Alternative olfactory neuron fates are specified by the LIM homeobox gene lim-4. Genes Dev. 1999, 13: 1794-1806. 10.1101/gad.13.14.1794.
Grigoriou M, Tucker AS, Sharpe PT, Pachnis V: Expression and regulation of Lhx6 and Lhx7, a novel subfamily of LIM homeodomain encoding genes, suggests a role in mammalian head development. Development. 1998, 125: 2063-2074.
Thor S, Thomas JB: The Drosophila islet gene governs axon pathfinding and neurotransmitter identity. Neuron. 1997, 18: 397-409. 10.1016/S0896-6273(00)81241-6.
Li Y, Zhang Y, He B, Wang Y, Yuan Z, Yuan W, Liao P, Deng Y, Xiao J, Zhu C, Wang Y, Wu X, Liu M: Cloning and expression of a novel human gene, Isl-2, encoded a LIM-homeodomain protein. Mol Biol Rep. 2007, 34: 19-26. 10.1007/s11033-006-9003-0.
Pfaff SL, Mendelsohn M, Stewart CL, Edlund T, Jessell TM: Requirement for LIM homeobox gene Isl1 in motor neuron generation reveals a motor neuron-dependent step in interneuron differentiation. Cell. 1996, 84: 309-320. 10.1016/S0092-8674(00)80985-X.
Kikuchi Y, Segawa H, Tokumoto M, Tsubokawa T, Hotta Y, Uyemura K, Okamoto H: Ocular and cerebellar defects in zebrafish induced by overexpression of the LIM domains of the islet-3 LIM/homeodomain protein. Neuron. 1997, 18: 369-382. 10.1016/S0896-6273(00)81239-8.
Hobert O, Tessmar K, Ruvkun G: The Caenorhabditis elegans lim-6 LIM homeobox gene regulates neurite outgrowth and function of particular GABAergic neurons. Development. 1999, 126: 1547-1562.
Giraldez F: Regionalized organizing activity of the neural tube revealed by the regulation of lmx1 in the otic vesicle. Dev Biol. 1998, 203: 189-200. 10.1006/dbio.1998.9023.
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, Grimwood 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-60. 10.1038/nature07191.
Srivastava M, Simakov O, Chapman J, Fahey B, Gauthier ME, 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-6. 10.1038/nature09201.
Sakarya O, Armstrong K, Adamska M, Adamski M, Wang I-F, Tidor B, Degnan BM, Oakley TH, Kosik KS: A post-synaptic scaffold at the origin of the animal kingdom. PloS one. 2007, 2: e506-10.1371/journal.pone.0000506.
Alié A, Manuel M: The backbone of the post-synaptic density originated in a unicellular ancestor of choanoflagellates and metazoans. BMC Evol Biol. 2010, 10: 34-10.1186/1471-2148-10-34.
Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, Sørensen MV, Haddock SHD, Schmidt-Rhaesa A, Okusu A, Kristensen RM, Wheeler WC, Martindale MQ, Giribet G: Broad phylogenomic sampling improves resolution of the animal tree of life. Nature. 2008, 452: 745-9. 10.1038/nature06614.
Hejnol A, Obst M, Stamatakis A, Ott M, Rouse GW, Edgecombe GD, Martinez P, Baguna J, Bailly X, Jondelius U, Wiens M, Muller WE, Seaver E, Wheeler WC, Martindale MQ, Giribet G, Dunn CW: Assessing the root of bilaterian animals with scalable phylogenomic methods. Proc Biol Sci. 2009, 276: 4261-4270. 10.1098/rspb.2009.0896.
Medina M, Collins AG, Silberman JD, Sogin ML: Evaluating hypotheses of basal animal phylogeny using complete sequences of large and small subunit rRNA. Proc Natl Acad Sci USA. 2001, 98: 9707-9712. 10.1073/pnas.171316998.
Collins AG: Evaluating multiple alternative hypotheses for the origin of Bilateria: an analysis of 18S rRNA molecular evidence. Proc Nat Acad Sci USA. 1998, 95: 15458-15463. 10.1073/pnas.95.26.15458.
Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, Salamov A, Terry A, Shapiro H, Lindquist E, Kapitonov VV, Jurka J, Genikhovich G, Grigoriev IV, Lucas SM, Steele RE, Finnerty JR, Technau U, Martindale MQ, Rokhsar DS: Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science. 2007, 317: 86-94. 10.1126/science.1139158.
Hyman LH: The invertebrates: Protozoa through Ctenophora. 1940, New York; McGraw-Hill Book Company, Inc
Philippe H, Derelle R, Lopez P, Pick K, Borchiellini C, Boury-Esnault N, Vacelet J, Renard E, Houliston E, Quéinnec E, Da Silva C, Wincker P, Le Guyader H, Leys S, Jackson DJ, Schreiber F, Erpenbeck D, Morgenstern B, Wörheide G, Manuel M: Phylogenomics revives traditional views on deep animal relationships. Curr Biol. 2009, 19: 706-712. 10.1016/j.cub.2009.02.052.
Edgecombe GD, Giribet G, Dunn CW, Hejnol A, Kristensen RM, Neves RC, Rouse GW, Worsaae K, Sørensen MV: Higher-level metazoan relationships: recent progress and remaining questions. Organisms Diversity & Evolution. 2011
Ryan JF, Pang K, Mullikin JC, Martindale MQ, Baxevanis AD: The homeodomain complement of the ctenophore Mnemiopsis leidyi suggests that Ctenophora and Porifera diverged prior to the ParaHoxozoa. EvoDevo. 2010, 1: 9-10.1186/2041-9139-1-9.
Pang K, Ryan JF, Mullikin JC, Baxevanis AD, Martindale MQ: Genomic insights into Wnt signaling in an early diverging metazoan, the ctenophore Mnemiopsis leidyi. EvoDevo. 2010, 1: 10-10.1186/2041-9139-1-10.
Reitzel AM, Pang K, Ryan JF, Mullikin JC, Martindale MQ, Baxevanis AD, Tarrant AM: Nuclear receptors from the ctenophore Mnemiopsis leidyi lack a zinc-finger DNA-binding domain: lineage-specific loss or ancestral condition in the emergence of the nuclear receptor superfamily?. EvoDevo. 2011, 2: 3-10.1186/2041-9139-2-3.
Letunic I, Doerks T, Bork P: SMART 6: recent updates and new developments. Nucleic Acids Res. 2009, 37: D229-232. 10.1093/nar/gkn808.
Ryan JF, Burton PM, Mazza ME, Kwong GK, Mullikin JC, Finnerty JR: The cnidarian-bilaterian ancestor possessed at least 56 homeoboxes : evidence from the starlet sea anemone, Nematostella vectensis. Genome Biol. 2006, 1-20.
Larroux C, Luke GN, Koopman P, Rokhsar DS, Shimeld SM, Degnan BM: Genesis and expansion of metazoan transcription factor gene classes. Mol Biol Evol. 2008, 25: 980-996. 10.1093/molbev/msn047.
Horridge GA: Presumed photoreceptive cilia in a ctenophore. Quart J Microsc Sci. 1964, 105: 311-317.
Pang K, Ryan JF, Mullikin JC, Baxevanis AD, Martindale MQ: TGF beta pathway. EvoDevo.
Tamm SL: Ctenophora. Electrical conduction and behaviour in "simple" invertebrates. Edited by: Shelton GAB. 1982, Oxford: Oxford University Press, 266-358.
Bachy I, Vernier P, Retaux S: The LIM-homeodomain gene family in the developing Xenopus brain: conservation and divergences with the mouse related to the evolution of the forebrain. J Neurosci. 2001, 21: 7620-7629.
Yamada A, Pang K, Martindale MQ, Tochinai S: Surprisingly complex T-box gene complement in diploblastic metazoans. Evol Dev. 2007, 9: 220-230. 10.1111/j.1525-142X.2007.00154.x.
Pang K, Martindale MQ: Developmental expression of homeobox genes in the ctenophore Mnemiopsis leidyi. Dev Genes Evol. 2008, 218: 307-319. 10.1007/s00427-008-0222-3.
Layden MJ, Meyer NP, Pang K, Seaver EC, Martindale MQ: Expression and phylogenetic analysis of the zic gene family in the evolution and development of metazoans. EvoDevo. 2010, 1: 12-10.1186/2041-9139-1-12.
Chun C: Die Ctenophoren des Golfes von Neapel und der angrenzenden Meeres-Abschnitte. Fauna und Flora des Golfes Neapel. 1880, Leipzig: Engelmann, 1: 1-311.
Hertwig R: Ueber den Bau der Ctenophoren. 1880, Jena: G. Fischer
Samassa P: Zur Histologie der Ctenophoren. Arch mikr Anat. 1892, 40: 157-242. 10.1007/BF02954492.
Bethe A: Der subepitheliale Nervenplexus der Ctenophoren. Biol Zbl. 1895, 15: 140-145.
Heider K: Über das Nervensystem von Beroe ovata. Nachrichten der Gesellschaft der Wissenschaften zu Göttingen, Math.-phys. KI. Göttingen: Vandenhoeck & Ruprecht. 1927
Horridge GA, Mackay B: Naked axons and symmetrical synapses in coelenterates. J Cell Sci. 1962, 3: 531-541.
Horridge GA, Mackay B: Neurociliary synapses in Pleurobrachia (Ctenophora). J Cell Sci. 1964, 3: 163-174.
Horridge GA: Relations between nerves and cilia in Ctenophores. Am Zool. 1965, 5: 357-375.
Hernandez-Nicaise ML: The nervous system of ctenophores. I. Structure and ultrastructure of the epithelial nerve-nets. Z Zellforsch Mikrosk Anat. 1973, 137: 223-250. 10.1007/BF00307432.
Hernandez-Nicaise ML: The nervous system of ctenophores. II. The nervous elements of the mesoglea of beroids and cydippids. Z Zellforsch Mikrosk Anat. 1973, 143: 117-133. 10.1007/BF00307455.
Hernandez-Nicaise ML: The nervous system of ctenophores. III. Ultrastructure of synapses. J Neurocytol. 1973, 2: 249-263. 10.1007/BF01104029.
Hernandez-Nicaise ML: Système nerveux et intégration chez. les ctenaires; etude ultrastructurale et comportementale. 1974, Lyon, France: Univ. Claude Bernard
Tamm S, Tamm SL: A giant nerve net with multi-effector synapses underlying epithelial adhesive strips in the mouth of Beroe (Ctenophora). J Neurocytol. 1995, 24: 711-723. 10.1007/BF01179820.
Tamm SL, Tamm S: Novel bridge of axon-like processes of epithelial cells in the aboral sense organ of ctenophores. J Morphol. 2002, 254: 99-120. 10.1002/jmor.10019.
Jager M, Chiori R, Alié A, Dayraud C, Quéinnec E, Manuel M: New insights on ctenophore neural anatomy: immunofluorescence study in Pleurobrachia pileus (Müller, 1776). J Exp Zool B Mol Dev Evol. 2011, 316B: 171-187. 10.1002/jez.b.21386.
Fol H: Ein Beitrag zur Anatomie und Entwickelungsgeschichte. einiger Rippenquallen. 1869
Horridge GA: Recent studies on the Ctenophora. Coelenterate Biology. Edited by: Muscatine L, Lenhoff HM. 1974, NewYork: Academic Press, 439-468.
Hernandez-Nicaise ML: Ctenophora. Biology of the Iintegument, Volume 1: Invertebrates. Edited by: Bereiter-Hahn J, Matoltsy AG, Richards KS. 1984, Berlin Heidelberg, New York: Springer, 96-111.
Ward WW: Extraction of renilla-type luciferin from the calcium-activated photoproteins aequorin, mnemiopsin, and berovin. Proc Nat Acad Sci USA. 1975, 72: 2530-2534. 10.1073/pnas.72.7.2530.
Anctil M: Ultrastructure of the luminescent system of the ctenophore Mnemiopsis leidyi. Cell Tissue Res. 1985, 242: 333-340.
Schnitzler CE, Pang K, Powers ML, Reitzel AM, Ryan JF, Simmons DK, Park M, Gupta J, Brooks SY, Blakesley RW, Haddock HSD, Martindale MQ, Baxevanis AD: Bioluminescence and the evolution of photoproteins: A ctenophore genome lights the way.
Koch JB, Ryan JF, Baxevanis AD: The diversification of the LIM superclass at the base of the metazoa increased subcellular complexity and promoted multicellular specialization. Plos One.
Pang K, Martindale MQ: Comb jellies (Ctenophora): a model for basal metazoan evolution and development. Emerging Model Organisms. Edited by: Crotty DA, Gann A. 2009, Cold Spring Harbor: Cold Spring Harbor Laboratory Press, 1: 167-195.
Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ: Jalview Version 2--a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009, 25: 1189-1191. 10.1093/bioinformatics/btp033.
Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.
Abascal F, Zardoya R, Posada D: ProtTest: Selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
This work was supported by NASA and NSF grants to MQM, and the Intramural Research Program of the National Human Genome Research Institute (National Institutes of Health). We would like to thank Bernard J Koch, Joseph F Ryan and Andreas D Baxevanis for continued motivation and inspiration. In addition we would like to thank Alice Young, Brian Schmidt, Natalie Gurson, Richelle Legaspi, and Betsy Novotny, who were involved with the Mnemiopsis genomic sequencing at NISC.
The authors declare that they have no competing interests.
DS drafted the manuscript, performed the sequence alignment, and performed the phylogenetic analyses. DS and KP performed the animal collection, gene isolation, sequencing and in situ hybridizations. KP helped in drafting the manuscript and in the conception of the study. MQM participated in the design and conception of the study and helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Alignments of Lhx gene homeodomain sequences. Alignment of Lhx1/5, Lhx3/4, and Islet genes homeodomain regions. Introns positions indicated by (0). Amq, Amphimedon queenslandica; Hs, Homo sapiens; Ml,Mnemiopsis leidyi; Nv, Nematostella vectensis; Ta, Trichoplax adhaerens. A. Human Lhx1 and Lhx5 have a shared intron not found within any other species. The Mnemiopsis Lhx1/5 gene contains an intron that is also not shared with any other species in this study and appears to be species specific. B. The Lhx3/4 gene in Mnemiopsis has two introns not found in any other species examined and appears to be species specific. Humans also have an intron position not shared with the other species. C. The Islet gene in Mnemiopsis has a species specific intron position interrupting its homeodomain sequence. Nematostellla has an intron interrupting the homeodomain sequence in a different location than the ctenophore sequence. Both the Mnemiopsis and Nematostellla intron positions do not overlap and seem to not be related. Humans, placozoans, and sponges do not have introns interrupting the homeodomain sequences. (PNG 341 KB)
Additional file 2: Full sequence data. The full sequence data of all of the LIM genes used in this study. Gene accession numbers, corresponding gene names used in the phylogeny, and species name are listed in table format. (XLS 47 KB)
Authors’ original submitted files for images
About this article
Cite this article
Simmons, D.K., Pang, K. & Martindale, M.Q. Lim homeobox genes in the Ctenophore Mnemiopsis leidyi: the evolution of neural cell type specification. EvoDevo 3, 2 (2012). https://doi.org/10.1186/2041-9139-3-2
- Polar Field
- Apical Organ
- Neural Cell Type
- Aboral Pole
- Homeodomain Sequence