The Dlx gene complement of the leopard shark, Triakis semifasciata, resembles that of mammals: implications for genomic and morphological evolution of jawed vertebrates

Ontogenetic morphometry of the brown smoothhound shark Mustelus henlei with implications for ecology and evolution

The central tenet of ecomorphology links ecological and morphological variation through the process of selection. Traditionally used to rationalise morphological differences between taxa, an ecomorphological approach is increasingly being utilised to study morphological differences expressed through ontogeny. Elasmobranchii (sharks, rays and skates) is one clade in which such ontogenetic shifts in body form have been reported. Such studies are limited to a relatively small proportion of total elasmobranch ecological and morphological diversity, and questions remain regarding the extent to which ecological selection are driving observed morphometric trends. In this study, we report ontogenetic growth trajectories obtained via traditional linear morphometrics from a large data set of the brown smoothhound shark (Mustelus henlei). We consider various morphological structures including the caudal, dorsal and pectoral fins, as well as several girth measurements. We use an ecomorphological approach to infer the broad ecological characteristics of this population and refine understanding of the selective forces underlying the evolution of specific morphological structures. We suggest that observed scaling trends in M. henlei are inconsistent with migratory behaviour, but do not contradict a putative trophic niche shift. We also highlight the role of predation pressure and sex-based ecological differences in driving observed trends in morphometry, a factor which has previously been neglected when considering the evolution of body form in sharks.

Direct evidence of allele equivalency at the Dlx5/6 locus

The retention of paralogous regulatory genes is a vertebrate hallmark and likely underpinned vertebrate origins. Dlx genes belong to a family of paralogous transcription factors whose evolutionary history of gene expansion and divergence is apparent from the gene synteny, shared exon-intron structure, and coding sequence homology found in extant vertebrate genomes. Dlx genes are expressed in a nested combination within the first pharyngeal arch and knockout studies in mice clearly point to a "Dlx code" that operates to define maxillary and mandibular position in the first arch. The nature of that code is not yet clear; an important goal for understanding Dlx gene function in both patterning and differentiation lies in distinguishing functional inputs that are paralog-specific (a qualitative model) versus Dlx family-generic (a quantitative model) and, in the latter case, the relative contribution made by each paralog. Here, multiple developmental deficiencies were identified in derivatives of the first pharyngeal arch in neonatal Dlx5/6(+/-) mice that resembled those seen in either paralog-specific null mutants. These data clearly demonstrate a substantial degree of allele equivalency and support a quantitative model of Dlx function during craniofacial morphogenesis. genesis 54:272-276, 2016. © 2016 Wiley Periodicals, Inc.

Molecular phylogeny of four homeobox genes from the purple sea star Pisaster ochraceus

Homeobox genes cloned from the purple sea star Pisaster ochraceus (Phylum Echinodermata/Class Asteroidea) were used along with related sequences available from members of other representative animal phyla to generate molecular phylogenies for Distal-less/Dlx, Hox5, Hox7, and Hox9/10 homeobox genes. Phylogenetic relationships were inferred based on the predicted 60 amino acid homeodomain, using amino acid (AA) and nucleotide (NT) models as well as the recently developed codon substitution models of sequence evolution. The resulting phylogenetic trees were mostly congruent with the consensus species-tree, grouping these newly identified genes with those isolated from other Asteroidea. This analysis also allowed a preliminary comparison of the performance of codon models with that of NT and AA evolutionary models in the inference of homeobox phylogeny. We found that, overall, the NT models displayed low reliability in recovering major clades at the Superphylum/Phylum level, and that codon models were slightly more dependable than AA models. Remarkably, in the majority of cases, codon substitution models seemed to outperform both AA and NT models at both the Class level and homeobox paralogy-group level of classification.

Heterogeneous conservation of Dlx paralog co-expression in jawed vertebrates

Background

The Dlx gene family encodes transcription factors involved in the development of a wide variety of morphological innovations that first evolved at the origins of vertebrates or of the jawed vertebrates. This gene family expanded with the two rounds of genome duplications that occurred before jawed vertebrates diversified. It includes at least three bigene pairs sharing conserved regulatory sequences in tetrapods and teleost fish, but has been only partially characterized in chondrichthyans, the third major group of jawed vertebrates. Here we take advantage of developmental and molecular tools applied to the shark Scyliorhinus canicula to fill in the gap and provide an overview of the evolution of the Dlx family in the jawed vertebrates. These results are analyzed in the theoretical framework of the DDC (Duplication-Degeneration-Complementation) model.

Results

The genomic organisation of the catshark Dlx genes is similar to that previously described for tetrapods. Conserved non-coding elements identified in bony fish were also identified in catshark Dlx clusters and showed regulatory activity in transgenic zebrafish. Gene expression patterns in the catshark showed that there are some expression sites with high conservation of the expressed paralog(s) and other expression sites with events of paralog sub-functionalization during jawed vertebrate diversification, resulting in a wide variety of evolutionary scenarios within this gene family.

Conclusion

Dlx gene expression patterns in the catshark show that there has been little neo-functionalization in Dlx genes over gnathostome evolution. In most cases, one tandem duplication and two rounds of vertebrate genome duplication have led to at least six Dlx coding sequences with redundant expression patterns followed by some instances of paralog sub-functionalization. Regulatory constraints such as shared enhancers, and functional constraints including gene pleiotropy, may have contributed to the evolutionary inertia leading to high redundancy between gene expression patterns.

Developmental evidence for serial homology of the vertebrate jaw and gill arch skeleton

Gegenbaur’s classical hypothesis of jaw-gill arch serial homology is widely cited, but remains unsupported by either paleontological evidence (e.g. a series of fossils reflecting the stepwise transformation of a gill arch into a jaw) or developmental genetic data (e.g. shared molecular mechanisms underlying segment identity in the mandibular, hyoid and gill arch endoskeletons). Here we show that nested expression of Dlx genes – the “Dlx code” that specifies upper and lower jaw identity in mammals and teleosts – is a primitive feature of the mandibular, hyoid and gill arches of jawed vertebrates. Using fate-mapping techniques, we demonstrate that the principal dorsal and ventral endoskeletal segments of the jaw, hyoid and gill arches of the skate Leucoraja erinacea derive from molecularly equivalent mesenchymal domains of combinatorial Dlx gene expression. Our data suggest that vertebrate jaw, hyoid and gill arch cartilages are serially homologous, and were primitively patterned dorsoventrally by a common Dlx blueprint.

Pattern and polarity in the development and evolution of the gnathostome jaw: both conservation and heterotopy in the branchial arches of the shark, Scyliorhinus canicula

The acquisition of jaws constitutes a landmark event in vertebrate evolution, one that in large part potentiated their success and diversification. Jaw development and patterning involves an intricate spatiotemporal series of reciprocal inductive and responsive interactions between the cephalic epithelia and the cranial neural crest (CNC) and cephalic mesodermal mesenchyme. The coordinated regulation of these interactions is critical for both the ontogenetic registration of the jaws and the evolutionary elaboration of variable jaw morphologies and designs. Current models of jaw development and evolution have been built on molecular and cellular evidence gathered mostly in amniotes such as mice, chicks and humans, and augmented by a much smaller body of work on the zebrafish. These have been partnered by essential work attempting to understand the origins of jaws that has focused on the jawless lamprey. Chondrichthyans (cartilaginous fish) are the most distant group to amniotes within extant gnathostomes, and comprise the crucial clade uniting amniotes and agnathans; yet despite their critical phylogenetic position, evidence of the molecular and cellular underpinnings of jaw development in chondrichthyans is still lacking. Recent advances in genome and molecular developmental biology of the lesser spotted dogfish shark, Scyliorhinus canicula, make it ideal for the molecular study of chondrichthyan jaw development. Here, following the 'Hinge and Caps' model of jaw development, we have investigated evidence of heterotopic (relative changes in position) and heterochronic (relative changes in timing) shifts in gene expression, relative to amniotes, in the jaw primordia of S. canicula embryos. We demonstrate the presence of clear proximo-distal polarity in gene expression patterns in the shark embryo, thus establishing a baseline molecular baüplan for branchial arch-derived jaw development and further validating the utility of the 'Hinge and Caps' model in comparative studies of jaw development and evolution. Moreover, we correlate gene expression patterns with the absence of a lambdoidal junction (formed where the maxillary first arch meets the frontonasal processes) in chondrichthyans, further highlighting the importance of this region for the development and evolution of jaw structure in advanced gnathostomes.

Non-parsimonious evolution of hagfish Dlx genes

Background

The number of members of the Dlx gene family increased during the two rounds of whole-genome duplication that occurred in the common ancestor of the vertebrates. Because the Dlx genes are involved in the development of the cranial skeleton, brain, and sensory organs, their expression patterns have been analysed in various organisms in the context of evolutionary developmental biology. Six Dlx genes have been isolated in the lampreys, a group of living jawless vertebrates (cyclostomes), and their expression patterns analysed. However, little is known about the Dlx genes in the hagfish, the other cyclostome group, mainly because the embryological analysis of this animal is difficult.

Results

To identify the hagfish Dlx genes and describe their expression patterns, we cloned the cDNA from embryos of the Japanese inshore hagfish Eptatretus burgeri. Our results show that the hagfish has at least six Dlx genes and one pseudogene. In a phylogenetic analysis, the hagfish Dlx genes and those of the lampreys tended to be excluded from the clade of the gnathostome Dlx genes. In several cases, the lamprey Dlx genes clustered with the clade consisting of two hagfish genes, suggesting that independent gene duplications have occurred in the hagfish lineage. Analysis of the expression of these genes showed distinctive overlapping expression patterns in the cranial mesenchymal cells and the inner ear.

Conclusions

Independent duplication, pseudogenization, and loss of the Dlx genes probably occurred in the hagfish lineage after its split from the other vertebrate lineages. This pattern is reminiscent of the non-parsimonious evolution of its morphological traits, including its inner ear and vertebrae, which indicate that this group is an early-branching lineage that diverged before those characters evolved.

The Dlx genes as clues to vertebrate genomics and craniofacial evolution

The group of Dlx genes belongs to the homeobox-containing superfamily, and its members are involved in various morphogenetic processes. In vertebrate genomes, Dlx genes exist as multiple paralogues generated by tandem duplication followed by whole genome duplications. In this review, we provide an overview of the Dlx gene phylogeny with an emphasis on the chordate lineage. Referring to the Dlx gene repertoire, we discuss the establishment and conservation of the nested expression patterns of the Dlx genes in craniofacial development. Despite the accumulating genomic sequence resources in diverse vertebrates, embryological analyses of Dlx gene expression and function remain limited in terms of species diversity. By supplementing our original analysis of shark embryos with previous data from other osteichthyans, such as mice and zebrafish, we support the previous speculation that the nested Dlx expression in the pharyngeal arch is likely a shared feature among all the extant jawed vertebrates. Here, we highlight several hitherto unaddressed issues regarding the evolution and function of Dlx genes, with special reference to the craniofacial development of vertebrates.

The homology of odontodes in gnathostomes: insights from Dlx gene expression in the dogfish, Scyliorhinus canicula

BACKGROUND

Teeth and tooth-like structures, together named odontodes, are repeated organs thought to share a common evolutionary origin. These structures can be found in gnathostomes at different locations along the body: oral teeth in the jaws, teeth and denticles in the oral-pharyngeal cavity, and dermal denticles on elasmobranch skin. We, and other colleagues, had previously shown that teeth in any location were serially homologous because: i) pharyngeal and oral teeth develop through a common developmental module; and ii) the expression patterns of the Dlx genes during odontogenesis were highly divergent between species but almost identical between oral and pharyngeal dentitions within the same species. Here we examine Dlx gene expression in oral teeth and dermal denticles in order to test the hypothesis of serial homology between these odontodes.

RESULTS

We present a detailed comparison of the first developing teeth and dermal denticles (caudal primary scales) of the dogfish (Scyliorhinus canicula) and show that both odontodes develop through identical stages that correspond to the common stages of oral and pharyngeal odontogenesis. We identified six Dlx paralogs in the dogfish and found that three showed strong transcription in teeth and dermal denticles (Dlx3, Dlx4 and Dlx5) whereas a weak expression was detected for Dlx1 in dermal denticles and teeth, and for Dlx2 in dermal denticles. Very few differences in Dlx expression patterns could be detected between tooth and dermal denticle development, except for the absence of Dlx2 expression in teeth.

CONCLUSIONS

Taken together, our histological and expression data strongly suggest that teeth and dermal denticles develop from the same developmental module and under the control of the same set of Dlx genes. Teeth and dermal denticles should therefore be considered as serial homologs developing through the initiation of a common gene regulatory network (GRN) at several body locations. This mechanism of heterotopy supports the 'inside and out' model that has been recently proposed for odontode evolution.

Annelid Distal-less/Dlx duplications reveal varied post-duplication fates

Background

Dlx (Distal-less) genes have various developmental roles and are widespread throughout the animal kingdom, usually occurring as single copy genes in non-chordates and as multiple copies in most chordate genomes. While the genomic arrangement and function of these genes is well known in vertebrates and arthropods, information about Dlx genes in other organisms is scarce. We investigate the presence of Dlx genes in several annelid species and examine Dlx gene expression in the polychaete Pomatoceros lamarckii.

Results

Two Dlx genes are present in P. lamarckii, Capitella teleta and Helobdella robusta. The C. teleta Dlx genes are closely linked in an inverted tail-to-tail orientation, reminiscent of the arrangement of vertebrate Dlx pairs, and gene conversion appears to have had a role in their evolution. The H. robusta Dlx genes, however, are not on the same genomic scaffold and display divergent sequences, while, if the P. lamarckii genes are linked in a tail-to-tail orientation they are a minimum of 41 kilobases apart and show no sign of gene conversion. No expression in P. lamarckii appendage development has been observed, which conflicts with the supposed conserved role of these genes in animal appendage development. These Dlx duplications do not appear to be annelid-wide, as the polychaete Platynereis dumerilii likely possesses only one Dlx gene.

Conclusions

On the basis of the currently accepted annelid phylogeny, we hypothesise that one Dlx duplication occurred in the annelid lineage after the divergence of P. dumerilii from the other lineages and these duplicates then had varied evolutionary fates in different species. We also propose that the ancestral role of Dlx genes is not related to appendage development.

Ancestral and derived attributes of the dlx gene repertoire, cluster structure and expression patterns in an African cichlid fish

Background

Cichlid fishes have undergone rapid, expansive evolutionary radiations that are manifested in the diversification of their trophic morphologies, tooth patterning and coloration. Understanding the molecular mechanisms that underlie the cichlids' unique patterns of evolution requires a thorough examination of genes that pattern the neural crest, from which these diverse phenotypes are derived. Among those genes, the homeobox-containing Dlx gene family is of particular interest since it is involved in the patterning of the brain, jaws and teeth.

Results

In this study, we characterized the dlx genes of an African cichlid fish, Astatotilapia burtoni, to provide a baseline to later allow cross-species comparison within Cichlidae. We identified seven dlx paralogs (dlx1a, -2a, -4a, -3b, -4b, -5a and -6a), whose orthologies were validated with molecular phylogenetic trees. The intergenic regions of three dlx gene clusters (dlx1a-2a, dlx3b-4b, and dlx5a-6a) were amplified with long PCR. Intensive cross-species comparison revealed a number of conserved non-coding elements (CNEs) that are shared with other percomorph fishes. This analysis highlighted additional lineage-specific gains/losses of CNEs in different teleost fish lineages and a novel CNE that had previously not been identified. Our gene expression analyses revealed overlapping but distinct expression of dlx orthologs in the developing brain and pharyngeal arches. Notably, four of the seven A. burtoni dlx genes, dlx2a, dlx3b, dlx4a and dlx5a, were expressed in the developing pharyngeal teeth.

Conclusion

This comparative study of the dlx genes of A. burtoni has deepened our knowledge of the diversity of the Dlx gene family, in terms of gene repertoire, expression patterns and non-coding elements. We have identified possible cichlid lineage-specific changes, including losses of a subset of dlx expression domains in the pharyngeal teeth, which will be the targets of future functional studies.

Distal-less-like protein distribution in the larval lamprey forebrain

A polyclonal antibody against the Drosophila distal-less (DLL) protein, cross-reactive with cognate vertebrate proteins, was employed to map DLL-like expression in the midlarval lamprey forebrain. This work aimed to characterize in detail the separate diencephalic and telencephalic DLL expression domains, in order to test our previous modified definition of the lamprey prethalamus [Pombal and Puelles (1999) J Comp Neurol 414:391-422], adapt our earlier schema of prosomeric subdivisions in the lamprey forebrain to more recent versions of this model [Pombal et al. (2009) Brain Behav Evol 74:7-19] and reexamine the pallio-subpallial regionalization of the lamprey telencephalon. We observed a large-scale conservation of the topologic distribution of the DLL protein, in consonance with patterns of Dlx expression present in other vertebrates studied. Moreover, evidence was obtained of substantial numbers of DLL-positive neurons in the olfactory bulb and the cerebral hemispheres, in a pattern consistent with possible tangential migration out of the subpallium into the overlying pallium, as occurs in mammals, birds, frogs and teleost fishes.

A conserved cluster of three PRD-class homeobox genes (homeobrain, rx and orthopedia) in the Cnidaria and Protostomia

Background

Homeobox genes are a superclass of transcription factors with diverse developmental regulatory functions, which are found in plants, fungi and animals. In animals, several Antennapedia (ANTP)-class homeobox genes reside in extremely ancient gene clusters (for example, the Hox, ParaHox, and NKL clusters) and the evolution of these clusters has been implicated in the morphological diversification of animal bodyplans. By contrast, similarly ancient gene clusters have not been reported among the other classes of homeobox genes (that is, the LIM, POU, PRD and SIX classes).

Results

Using a combination of in silico queries and phylogenetic analyses, we found that a cluster of three PRD-class homeobox genes (Homeobrain (hbn), Rax (rx) and Orthopedia (otp)) is present in cnidarians, insects and mollusks (a partial cluster comprising hbn and rx is present in the placozoan Trichoplax adhaerens). We failed to identify this 'HRO' cluster in deuterostomes; in fact, the Homeobrain gene appears to be missing from the chordate genomes we examined, although it is present in hemichordates and echinoderms. To illuminate the ancestral organization and function of this ancient cluster, we mapped the constituent genes against the assembled genome of a model cnidarian, the sea anemone Nematostella vectensis, and characterized their spatiotemporal expression using in situ hybridization. In N. vectensis, these genes reside in a span of 33 kb with the same gene order as previously reported in insects. Comparisons of genomic sequences and expressed sequence tags revealed the presence of alternative transcripts of Nv-otp and two highly unusual protein-coding polymorphisms in the terminal helix of the Nv-rx homeodomain. A population genetic survey revealed the Rx polymorphisms to be widespread in natural populations. During larval development, all three genes are expressed in the ectoderm, in non-overlapping territories along the oral-aboral axis, with distinct temporal expression.

Conclusion

We report the first evidence for a PRD-class homeobox cluster that appears to have been conserved since the time of the cnidarian-bilaterian ancestor, and possibly even earlier, given the presence of a partial cluster in the placozoan Trichoplax. Very similar clusters comprising these three genes exist in Nematostella and diverse protostomes. Interestingly, in chordates, one member of the ancestral cluster (homeobrain) has apparently been lost, and there is no linkage between rx and orthopedia in any of the vertebrates. In Nematostella, the spatial expression of these three genes along the body column is not colinear with their physical order in the cluster but the temporal expression is, therefore, using the terminology that has been applied to the Hox cluster genes, the HRO cluster would appear to exhibit temporal but not spatial colinearity. It remains to be seen whether the mechanisms responsible for the evolutionary conservation of the HRO cluster are the same mechanisms responsible for cohesion of the Hox cluster and other ANTP-class homeobox clusters that have been widely conserved throughout animal evolution.

Functional conservation of a forebrain enhancer from the elephant shark (Callorhinchus milii ) in zebrafish and mice

Background

The phylogenetic position of the elephant shark (Callorhinchus milii ) is particularly relevant to study the evolution of genes and gene regulation in vertebrates. Here we examine the evolution of Dlx homeobox gene regulation during vertebrate embryonic development with a particular focus on the forebrain. We first identified the elephant shark sequence orthologous to the URE2 cis -regulatory element of the mouse Dlx1/Dlx2 locus (herein named CmURE2). We then conducted a comparative study of the sequence and enhancer activity of CmURE2 with that of orthologous regulatory sequences from zebrafish and mouse.

Results

The CmURE2 sequence shows a high percentage of identity with its mouse and zebrafish counterparts but is overall more similar to mouse URE2 (MmURE2) than to zebrafish URE2 (DrURE2). In transgenic zebrafish and mouse embryos, CmURE2 displayed enhancer activity in the forebrain that overlapped with that of DrURE2 and MmURE2. However, we detected notable differences in the activity of the three sequences in the diencephalon. Outside of the forebrain, CmURE2 shows enhancer activity in areas such as the pharyngeal arches and dorsal root ganglia where its' counterparts are also active.

Conclusions

Our transgenic assays show that part of the URE2 enhancer activity is conserved throughout jawed vertebrates but also that new characteristics have evolved in the different groups. Our study demonstrates that the elephant shark is a useful outgroup to study the evolution of regulatory mechanisms in vertebrates and to address how changes in the sequence of cis -regulatory elements translate into changes in their regulatory activity.

The evolution of gnathostome development: Insight from chondrichthyan embryology

Chondrichthyans (cartilaginous fishes) represent one of the two lineages of gnathostomes, the other being the osteicthyans (bony fishes). Classical studies on chondrichthyan embryology have strongly impacted our views of vertebrate body plan evolution, while recent studies highlight oviparous chondrichthyans as emerging vertebrate model systems that are amenable to experimental embryological manipulation. Here, we review three particular areas of interest in the field of chondrichthyan developmental biology-gastrulation, neural development, and appendage patterning-and we discuss recent findings within a broader chondrichthyan-osteichthyan comparative framework. In some cases, comparative studies of chondrichthyan and osteichthyan development reveal conserved patterns of gene expression in common developmental contexts. Studies of chondrichthyan gastrulation reveal conserved patterns of developmental gene expression, despite highly divergent modes of mesendoderm internalization, while molecular characterization of chondrichthyan neurogenic placodes indicates a conservation of placode transcription factor expression across gnathostome phylogeny. In other cases, comparative studies of chondrichthyan and osteichthyan development yield evidence of shared patterning mechanisms functioning in different developmental contexts. This is exemplified by studies on the development of chondrichthyan appendages-paired fins, median fins, and gill rays. These have demonstrated that a retinoic acid-responsive Shh-expressing signaling center functions to pattern the endoskeleton of gnathostome paired fins and chondrichthyan gill rays, while expression patterns of Tbx18 and HoxD family members are shared by gnathostome paired fins and chondrichthyan median fins. These findings fuel novel hypotheses of developmental genetic homology, and demonstrate how comparative studies of gnathostome development can provide insight into the evolutionary processes that underlie morphological diversity.

Gene classification based on amino acid motifs and residues: the DLX (distal-less) test case

Background

Comparative studies using hundreds of sequences can give a detailed picture of the evolution of a given gene family. Nevertheless, retrieving only the sequences of interest from public databases can be difficult, in particular, when working with highly divergent sequences. The difficulty increases substantially when one wants to include in the study sequences from many (or less well studied) species whose genomes are non-annotated or incompletely annotated.

Methodology/Principal Findings

In this work we evaluate the usefulness of different approaches of gene retrieval and classification, using the distal-less (DLX) gene family as a test case. Furthermore, we evaluate whether the use of a large number of gene sequences from a wide range of animal species, the use of multiple alternative alignments, and the use of amino acids aligned with high confidence only, is enough to recover the accepted DLX evolutionary history.

Conclusions/Significance

The canonical DLX homeobox gene sequence here derived, together with the characteristic amino acid variants here identified in the DLX homeodomain region, can be used to retrieve and classify DLX genes in a simple and efficient way. A program is made available that allows the easy retrieval of synteny information that can be used to classify gene sequences. Maximum likelihood trees using hundreds of sequences can be used for gene identification. Nevertheless, for the DLX case, the proposed DLX evolutionary is not recovered even when multiple alignment algorithms are used.

Multiple chromosomal rearrangements structured the ancestral vertebrate Hox-bearing protochromosomes

While the proposal that large-scale genome expansions occurred early in vertebrate evolution is widely accepted, the exact mechanisms of the expansion--such as a single or multiple rounds of whole genome duplication, bloc chromosome duplications, large-scale individual gene duplications, or some combination of these--is unclear. Gene families with a single invertebrate member but four vertebrate members, such as the Hox clusters, provided early support for Ohno's hypothesis that two rounds of genome duplication (the 2R-model) occurred in the stem lineage of extant vertebrates. However, despite extensive study, the duplication history of the Hox clusters has remained unclear, calling into question its usefulness in resolving the role of large-scale gene or genome duplications in early vertebrates. Here, we present a phylogenetic analysis of the vertebrate Hox clusters and several linked genes (the Hox "paralogon") and show that different phylogenies are obtained for Dlx and Col genes than for Hox and ErbB genes. We show that these results are robust to errors in phylogenetic inference and suggest that these competing phylogenies can be resolved if two chromosomal crossover events occurred in the ancestral vertebrate. These results resolve conflicting data on the order of Hox gene duplications and the role of genome duplication in vertebrate evolution and suggest that a period of genome reorganization occurred after genome duplications in early vertebrates.

Timing of genome duplications relative to the origin of the vertebrates: did cyclostomes diverge before or after?

Two rounds of whole-genome duplications are thought to have played an important role in the establishment of gene repertoires in vertebrates. These events occurred during chordate evolution after the split of the urochordate and cephalochordate lineages but before the radiation of extant gnathostomes (jawed vertebrates). During this interval, diverse agnathans (jawless fishes), including cyclostomes (hagfishes and lampreys), diverged. However, there is no solid evidence for the timing of these genome duplications in relation to the divergence of cyclostomes from the gnathostome lineage. We conducted cDNA sequencing in diverse early vertebrates for members of homeobox-containing (Dlx and ParaHox) and other gene families that would serve as landmarks for genome duplications. Including these new sequences, we performed a molecular phylogenetic census using the maximum likelihood method for 55 gene families. In most of these gene families, we detected many more gene duplications before the cyclostome-gnathostome split, than after. Many of these gene families (e.g., visual opsins, RAR, Notch) have multiple paralogs in conserved, syntenic genomic regions that must have been generated by large-scale duplication events. Taken together, this indicates that the genome duplications occurred before the cyclostome-gnathostome split. We propose that the redundancy in gene repertoires possessed by all vertebrates, including hagfishes and lampreys, was introduced primarily by genome duplications. Apart from subsequent lineage-specific modifications, these ancient genome duplication events might serve generally to distinguish vertebrates from invertebrates at the genomic level.

Phylogenetic and chromosomal analyses of multiple gene families syntenic with vertebrate Hox clusters

Background

Ever since the theory about two rounds of genome duplication (2R) in the vertebrate lineage was proposed, the Hox gene clusters have served as the prime example of quadruplicate paralogy in mammalian genomes. In teleost fishes, the observation of additional Hox clusters absent in other vertebrate lineages suggested a third tetraploidization (3R). Because the Hox clusters occupy a quite limited part of each chromosome, and are special in having position-dependent regulation within the multi-gene cluster, studies of syntenic gene families are needed to determine the extent of the duplicated chromosome segments. We have analyzed in detail 14 gene families that are syntenic with the Hox clusters to see if their phylogenies are compatible with the Hox duplications and the 2R/3R scenario. Our starting point was the gene family for the NPY family of peptides located near the Hox clusters in the pufferfish Takifugu rubripes, the zebrafish Danio rerio, and human.

Results

Seven of the gene families have members on at least three of the human Hox chromosomes and two families are present on all four. Using both neighbor-joining and quartet-puzzling maximum likelihood methods we found that 13 families have a phylogeny that supports duplications coinciding with the Hox cluster duplications. One additional family also has a topology consistent with 2R but due to lack of urochordate or cephalocordate sequences the time window when these duplications could have occurred is wider. All but two gene families also show teleost-specific duplicates.

Conclusion

Based on this analysis we conclude that the Hox cluster duplications involved a large number of adjacent gene families, supporting expansion of these families in the 2R, as well as in the teleost 3R tetraploidization. The gene duplicates presumably provided raw material in early vertebrate evolution for neofunctionalization and subfunctionalization.

Evolution and development of complex biomechanical systems: 300 million years of fish jaws

The jaws of teleost fishes are diverse and complex musculoskeletal systems. The focus in this review is on the major biomechanical systems in the teleost head, and the range and interplay of functional, developmental, and genetic influences that shape the modular and integrated evolution of elements. Insights possible from comparative studies are discussed in the context of traditional and new models for studies of craniofacial evolution and development.

Expression and function of Dlx genes in the osteoblast lineage

Our laboratory and others have shown that overexpression of Dlx5 stimulates osteoblast differentiation. Dlx5(-/-)/Dlx6(-/-) mice have more severe craniofacial and limb defects than Dlx5(-/-), some of which are potentially due to defects in osteoblast maturation. We wished to investigate the degree to which other Dlx genes compensate for the lack of Dlx5, thus allowing normal development of the majority of skeletal elements in Dlx5(-/-) mice. Dlx gene expression in cells from different stages of the osteoblast lineage isolated by FACS sorting showed that Dlx2, Dlx5 and Dlx6 are expressed most strongly in less mature osteoblasts, whereas Dlx3 is very highly expressed in differentiated osteoblasts and osteocytes. In situ hybridization and Northern blot analysis demonstrated the presence of endogenous Dlx3 mRNA within osteoblasts and osteocytes. Dlx3 strongly upregulates osteoblastic markers with a potency comparable to Dlx5. Cloned chick or mouse Dlx6 showed stimulatory effects on osteoblast differentiation. Our results suggest that Dlx2 and Dlx6 have the potential to stimulate osteoblastic differentiation and may compensate for the absence of Dlx5 to produce relatively normal osteoblastic differentiation in Dlx5 knockout mice, while Dlx3 may play a distinct role in late stage osteoblast differentiation and osteocyte function.

Evolutionary genetics of vertebrate tissue mineralization: the origin and evolution of the secretory calcium-binding phosphoprotein family

Three principal mineralized tissues are present in teeth; a highly mineralized surface layer (enamel or enameloid), body dentin, and basal bone. Similar tissues have been identified in the dermal skeleton of Paleozoic jawless vertebrates, suggesting their ancient origin. These dental tissues form on protein matrix and their mineralization is controlled by distinctive proteins. We have shown that many secretory calcium-binding phosphoproteins (SCPPs) are involved in tetrapod tissue mineralization. These SCPPs all originated from the common ancestral gene SPARCL1 (secreted protein, acidic, cysteine-rich like 1) that initially arose from SPARC. The SCPP family also includes a bird eggshell matrix protein, mammalian milk casein, and salivary proteins. The eggshell SCPP plays crucial roles in rigid eggshell production, milk SCPPs in efficient lactation and in the evolution of complex dentition, and salivary SCPPs in maintaining tooth integrity. A comparative analysis of the mammalian, avian, and amphibian genomes revealed a tandem duplication history of the SCPP genes in tetrapods. Although these tetrapod SCPP genes are fewer in teleost genomes, independent parallel duplication has created distinct SCPP genes in this lineage. These teleost SCPPs are also used for enameloid and dentin mineralization, implying essential roles of SCPPs for dental tissue mineralization in osteichthyans. However, the SCPPs used for tetrapod enamel and teleost enameloid, as well as tetrapod dentin and teleost dentin, are all different. Thus, the evolution of vertebrate mineralized tissues seems to be explained by phenogenetic drift: while mineralized tissues are retained during vertebrate evolution, the underlying genetic basis has extensively drifted.

Non-overlapping expression patterns of the clustered Dll-A/B genes in the ascidian Ciona intestinalis

The Ci-Dll-A and Ci-Dll-B genes of Ciona intestinalis are arranged in a convergently transcribed gene cluster. This genomic arrangement is similar to that of the multiple bigene clusters of the Dlx homologs in vertebrates. Analysis of whole genome sequences showed that linkage to the Hox cluster is conserved with the vertebrate clusters. Phylogenetic analysis supports gene trees consistent with homology of the ascidian and vertebrate Dlx clusters, and in combination with the apparent conservation of genomic arrangement, it is concluded that the ascidian cluster is most likely homologous with the vertebrate clusters. Using whole-mount in situ hybridization, Ci-Dll-B transcripts were detected in all ectodermal lineages through gastrulation. Expression is radically downregulated in the neurula with detectable expression disappearing around the time that Ci-Dll-A expression appears in the anterior ectoderm. By the late tailbud stage Ci-Dll-Atranscripts were detected in the bilateral atrial primordia and persisted in the atrial rudiments to the larval stage, suggesting a role in development of these neural placode-like structures. This non-overlapping expression contradicts a common pattern seen in clustered genes, where as adjacent paralogs have largely overlapping expression domains. Enhancer sharing is often proposed as an explanation for the overlapping expression of gene cluster members. For this case of non-overlapping expression a model is proposed in which repressors acting at different stages override one or more shared enhancers. The enhancer sharing prevents breakup of the cluster while the independent temporal suppressors hide the presence of the shared enhancers.

The HOX genes are expressed, in vivo, in human tooth germs: in vitro cAMP exposure of dental pulp cells results in parallel HOX network activation and neuronal differentiation

Homeobox-containing genes play a crucial role in odontogenesis. After the detection of Dlx and Msx genes in overlapping domains along maxillary and mandibular processes, a homeobox odontogenic code has been proposed to explain the interaction between different homeobox genes during dental lamina patterning. No role has so far been assigned to the Hox gene network in the homeobox odontogenic code due to studies on specific Hox genes and evolutionary considerations. Despite its involvement in early patterning during embryonal development, the HOX gene network, the most repeat-poor regions of the human genome, controls the phenotype identity of adult eukaryotic cells. Here, according to our results, the HOX gene network appears to be active in human tooth germs between 18 and 24 weeks of development. The immunohistochemical localization of specific HOX proteins mostly concerns the epithelial tooth germ compartment. Furthermore, only a few genes of the network are active in embryonal retromolar tissues, as well as in ectomesenchymal dental pulp cells (DPC) grown in vitro from adult human molar. Exposure of DPCs to cAMP induces the expression of from three to nine total HOX genes of the network in parallel with phenotype modifications with traits of neuronal differentiation. Our observations suggest that: (i) by combining its component genes, the HOX gene network determines the phenotype identity of epithelial and ectomesenchymal cells interacting in the generation of human tooth germ; (ii) cAMP treatment activates the HOX network and induces, in parallel, a neuronal-like phenotype in human primary ectomesenchymal dental pulp cells.

Expression of Dlx genes during the development of the zebrafish pharyngeal dentition: evolutionary implications

In order to investigate similarities and differences in genetic control of development among teeth within and between species, we determined the expression pattern of all eight Dlx genes of the zebrafish during development of the pharyngeal dentition and compared these data with that reported for mouse molar tooth development. We found that (i) dlx1a and dlx6a are not expressed in teeth, in contrast to their murine orthologs, Dlx1 and Dlx6; (ii) the expression of the six other zebrafish Dlx genes overlaps in time and space, particularly during early morphogenesis; (iii) teeth in different locations and generations within the zebrafish dentition differ in the number of genes expressed; (iv) expression similarities and differences between zebrafish Dlx genes do not clearly follow phylogenetic and linkage relationships; and (v) similarities and differences exist in the expression of zebrafish and mouse Dlx orthologs. Taken together, these results indicate that the Dlx gene family, despite having been involved in vertebrate tooth development for over 400 million years, has undergone extensive diversification of expression of individual genes both within and between dentitions. The latter type of difference may reflect the highly specialized dentition of the mouse relative to that of the zebrafish, and/or genome duplication in the zebrafish lineage facilitating a redistribution of Dlx gene function during odontogenesis.

Reassessing the Dlx code: the genetic regulation of branchial arch skeletal pattern and development

The branchial arches are meristic vertebrate structures, being metameric both between each other within the rostrocaudal series along the ventrocephalic surface of the embryonic head and within each individual arch: thus, just as each branchial arch must acquire a unique identity along the rostrocaudal axis, each structure within the proximodistal axis of an arch must also acquire a unique identity. It is believed that regional specification of metameric structures is controlled by the nested expression of related genes resulting in a regional code, a principal that is though to be demonstrated by the regulation of rostrocaudal axis development in animals exerted by the nested HOM-C/Hox homeobox genes. The nested expression pattern of the Dlx genes within the murine branchial arch ectomesenchyme has more recently led to the proposal of a Dlx code for the regional specification along the proximodistal axis of the branchial arches (i.e. it establishes intra-arch identity). This review re-examines this hypothesis, and presents new work on an allelic series of Dlx loss-of-function mouse mutants that includes various combinations of Dlx1, Dlx2, Dlx3, Dlx5 and Dlx6. Although we confirm fundamental aspects of the hypothesis, we further report a number of novel findings. First, contrary to initial reports, Dlx1, Dlx2 and Dlx1/2 heterozygotes exhibit alterations of branchial arch structures and Dlx2-/- and Dlx1/2-/- mutants have slight alterations of structures derived from the distal portions of their branchial arches. Second, we present evidence for a role for murine Dlx3 in the development of the branchial arches. Third, analysis of compound Dlx mutants reveals four grades of mandibular arch transformations and that the genetic interactions of cis first-order (e.g. Dlx5 and Dlx6), trans second-order (e.g. Dlx5 and Dlx2) and trans third-order paralogues (e.g. Dlx5 and Dlx1) result in significant and distinct morphological differences in mandibular arch development. We conclude by integrating functions of the Dlx genes within the context of a hypothesized general mechanism for the establishment of pattern and polarity in the first branchial arch of gnathostomes that includes regionally secreted growth factors such as Fgf8 and Bmp and other transcription factors such as Msx1, and is consistent both with the structure of the conserved gnathostome jaw bauplan and the elaboration of this bauplan to meet organismal end-point designs.