Hox genes and the phylogeny of the arthropods

Evolution of the insect Hox gene cluster: Comparative analysis across 243 species

The Hox gene cluster is an iconic example of evolutionary conservation between divergent animal lineages, providing evidence for ancient similarities in the genetic control of embryonic development. However, there are differences between taxa in gene order, gene number and genomic organisation implying conservation is not absolute. There are also examples of radical functional change of Hox genes; for example, the ftz, zen and bcd genes in insects play roles in segmentation, extraembryonic membrane formation and body polarity, rather than specification of anteroposterior position. There have been detailed descriptions of Hox genes and Hox gene clusters in several insect species, including important model systems, but a large-scale overview has been lacking. Here we extend these studies using the publicly-available complete genome sequences of 243 insect species from 13 orders. We show that the insect Hox cluster is characterised by large intergenic distances, consistently extreme in Odonata, Orthoptera, Hemiptera and Trichoptera, and always larger between the 'posterior' Hox genes. We find duplications of ftz and zen in many species and multiple independent cluster breaks, although certain modules of neighbouring genes are rarely broken apart suggesting some organisational constraints. As more high-quality genomes are obtained, a challenge will be to relate structural genomic changes to phenotypic change across insect phylogeny.

Expression of posterior Hox genes and opisthosomal appendage development in a mygalomorph spider

Spiders represent an evolutionary successful group of chelicerate arthropods. The body of spiders is subdivided into two regions (tagmata). The anterior tagma, the prosoma, bears the head appendages and four pairs of walking legs. The segments of the posterior tagma, the opisthosoma, either lost their appendages during the course of evolution or their appendages were substantially modified to fulfill new tasks such as reproduction, gas exchange, and silk production. Previous work has shown that the homeotic Hox genes are involved in shaping the posterior appendages of spiders. In this paper, we investigate the expression of the posterior Hox genes in a tarantula that possesses some key differences of posterior appendages compared to true spiders, such as the lack of the anterior pair of spinnerets and a second set of book lungs instead of trachea. Based on the observed differences in posterior Hox gene expression in true spiders and tarantulas, we argue that subtle changes in the Hox gene expression of the Hox genes abdA and AbdB are possibly responsible for at least some of the morphological differences seen in true spiders versus tarantulas.

Diversity, duplication, and genomic organization of homeobox genes in Lepidoptera

Homeobox genes encode transcription factors with essential roles in patterning and cell fate in developing animal embryos. Many homeobox genes, including Hox and NK genes, are arranged in gene clusters, a feature likely related to transcriptional control. Sparse taxon sampling and fragmentary genome assemblies mean that little is known about the dynamics of homeobox gene evolution across Lepidoptera or about how changes in homeobox gene number and organization relate to diversity in this large order of insects. Here we analyze an extensive data set of high-quality genomes to characterize the number and organization of all homeobox genes in 123 species of Lepidoptera from 23 taxonomic families. We find most Lepidoptera have around 100 homeobox loci, including an unusual Hox gene cluster in which the lab gene is repositioned and the ro gene is next to pb A topologically associating domain spans much of the gene cluster, suggesting deep regulatory conservation of the Hox cluster arrangement in this insect order. Most Lepidoptera have four Shx genes, divergent zen-derived loci, but these loci underwent dramatic duplication in several lineages, with some moths having over 165 homeobox loci in the Hox gene cluster; this expansion is associated with local LINE element density. In contrast, the NK gene cluster content is more stable, although there are differences in organization compared with other insects, as well as major rearrangements within butterflies. Our analysis represents the first description of homeobox gene content across the order Lepidoptera, exemplifying the potential of newly generated genome assemblies for understanding genome and gene family evolution.

Functional Identification and Characterization of the Diuretic Hormone 31 (DH31) Signaling System in the Green Shore Crab, Carcinus maenas

The functional characterization of crustacean neuropeptides and their cognate receptors has not kept pace with the recent advances in sequence determination and, therefore, our understanding of the physiological roles of neuropeptides in this important arthropod sub-phylum is rather limited. We identified a candidate receptor-ligand pairing for diuretic hormone 31 (DH31) in a neural transcriptome of the crab, Carcinus maenas. In insects, DH31 plays species -specific but central roles in many facets of physiology, including fluid secretion, myoactivity, and gut peristalsis but little is known concerning its functions in crustaceans. The C. maenas DH31 transcript codes for a 147 amino acid prepropeptide, and a single receptor transcript translates to a secretin-like (Class B1) G protein-coupled receptor (GPCR). We used an in vitro aequorin luminescence Ca2+ mobilization assay to demonstrate that this candidate DH31R is activated byCarcinus and insect DH31s in a dose-dependent manner (EC50 15-30 nM). Whole mount immunohistochemical and in situ hybridization localization revealed extensive DH31 expressing neurons throughout the central nervous system, most notably in the abdominal ganglion where large, unpaired cells give rise to medial nerves, which terminate in extensive DH31 immunopositive dendritic fields intimately associated with oesophageal musculature. This system constitutes a large and hitherto undescribed neurohemal area adjacent to key muscle groups associated with the gastric system. DH31 expressing neurons were also seen in the cardiac, commissural, oesophageal, and stomatogastric ganglia and intense labeling was seen in dendrites innervating fore- and hindgut musculature but not with limb muscles. These labeling patterns, together with measurement of DH31R mRNA in the heart and hindgut, prompted us test the effects of DH31 on semi-isolated heart preparations. Cardiac superfusion with peptide evoked increased heart rates (10-100 nM). The neuroanatomical distribution of DH31 and its receptor transcripts, particularly that associated with gastric and cardiac musculature, coupled with the cardio- acceleratory effects of the peptide implicate this peptide in key myoactive roles, likely related to rhythmic coordination.

Characterization of NlHox3, an essential gene for embryonic development in Nilaparvata lugens

Hox genes encode transcriptional regulatory proteins that control axial patterning in all bilaterians. The brown planthopper (BPH), Nilaparvata lugens (Hemiptera: Delphacidae), is a destructive insect pest of rice plants in Asian countries. During analysis of the N. lugens transcriptome, we identified a Hox3-like gene (NlHox3) that was highly and specifically expressed in the embryonic stage. We performed functional analysis on the gene to identify its roles in embryonic development and its potential use as a target in RNA interference (RNAi) based pest control. The sequence analysis showed that NlHox3 was homologous to the Hox3 gene and was most closely related with zen of Drosophila. There were no significant differences in oviposition between the treated and control females after injecting double-stranded RNA of NlHox3 (dsNlHox3) into newly emerged female adult BPHs; however, there was a significant difference in the hatchability of those eggs laid, which no egg from the treated group hatched normally. Injecting female adult BPHs with dsNlHox3 led to necrosis of these offspring embryos, with eye reversal and undeveloped organs, suggesting that NlHox3 was an essential gene for embryonic development and might be a potential target for RNAi-based control of this insect pest.

Use of Continuous Traits Can Improve Morphological Phylogenetics

The recent surge in enthusiasm for simultaneously inferring relationships from extinct and extant species has reinvigorated interest in statistical approaches for modeling morphological evolution. Current statistical methods use the Mk model to describe substitutions between discrete character states. Although representing a significant step forward, the Mk model presents challenges in biological interpretation, and its adequacy in modeling morphological evolution has not been well explored. Another major hurdle in morphological phylogenetics concerns the process of character coding of discrete characters. The often subjective nature of discrete character coding can generate discordant results that are rooted in individual researchers' subjective interpretations. Employing continuous measurements to infer phylogenies may alleviate some of these issues. Although not widely used in the inference of topology, models describing the evolution of continuous characters have been well examined, and their statistical behavior is well understood. Also, continuous measurements avoid the substantial ambiguity often associated with the assignment of discrete characters to states. I present a set of simulations to determine whether use of continuous characters is a feasible alternative or supplement to discrete characters for inferring phylogeny. I compare relative reconstruction accuracy by inferring phylogenies from simulated continuous and discrete characters. These tests demonstrate significant promise for continuous traits by demonstrating their higher overall accuracy as compared to reconstruction from discrete characters under Mk when simulated under unbounded Brownian motion, and equal performance when simulated under an Ornstein-Uhlenbeck model. Continuous characters also perform reasonably well in the presence of covariance between sites. I argue that inferring phylogenies directly from continuous traits may be benefit efforts to maximize phylogenetic information in morphological data sets by preserving larger variation in state space compared to many discretization schemes. I also suggest that the use of continuous trait models in phylogenetic reconstruction may alleviate potential concerns of discrete character model adequacy, while identifying areas that require further study in this area. This study provides an initial controlled demonstration of the efficacy of continuous characters in phylogenetic inference.

Functional Characterization and Signaling Systems of Corazonin and Red Pigment Concentrating Hormone in the Green Shore Crab, Carcinus maenas

Neuropeptides play a central role as neurotransmitters, neuromodulators and hormones in orchestrating arthropod physiology. The post-genomic surge in identified neuropeptides and their putative receptors has not been matched by functional characterization of ligand-receptor pairs. Indeed, until very recently no G protein-coupled receptors (GPCRs) had been functionally defined in any crustacean. Here we explore the structurally-related, functionally-diverse gonadotropin-releasing hormone paralogs, corazonin (CRZ) and red-pigment concentrating hormone (RPCH) and their G-protein coupled receptors (GPCRs) in the crab, Carcinus maenas. Using aequorin luminescence to measure in vitro Ca2+ mobilization we demonstrated receptor-ligand pairings of CRZ and RPCH. CRZR-activated cell signaling in a dose-dependent manner (EC50 0.75 nM) and comparative studies with insect CRZ peptides suggest that the C-terminus of this peptide is important in receptor-ligand interaction. RPCH interacted with RPCHR with extremely high sensitivity (EC50 20 pM). Neither receptor bound GnRH, nor the AKH/CRZ-related peptide. Transcript distributions of both receptors indicate that CRZR expression was, unexpectedly, restricted to the Y-organs (YO). Application of CRZ peptide to YO had no effect on ecdysteroid biosynthesis, excepting a modest stimulation in early post-molt. CRZ had no effect on heart activity, blood glucose levels, lipid mobilization or pigment distribution in chromatophores, a scenario that reflected the distribution of its mRNA. Apart from the well-known activity of RPCH as a chromatophorotropin, it also indirectly elicited hyperglycemia (which was eyestalk-dependent). RPCHR mRNA was also expressed in the ovary, indicating possible roles in reproduction. The anatomy of CRZ and RPCH neurons in the nervous system is described in detail by immunohistochemistry and in situ hybridization. Each peptide has extensive but non-overlapping distribution in the CNS, and neuroanatomy suggests that both are possibly released from the post-commissural organs. This study is one of the first to deorphanize a GPCR in a crustacean and to provide evidence for hitherto unknown and diverse functions of these evolutionarily-related neuropeptides.

Progress, pitfalls and parallel universes: a history of insect phylogenetics

The phylogeny of insects has been both extensively studied and vigorously debated for over a century. A relatively accurate deep phylogeny had been produced by 1904. It was not substantially improved in topology until recently when phylogenomics settled many long-standing controversies. Intervening advances came instead through methodological improvement. Early molecular phylogenetic studies (1985-2005), dominated by a few genes, provided datasets that were too small to resolve controversial phylogenetic problems. Adding to the lack of consensus, this period was characterized by a polarization of philosophies, with individuals belonging to either parsimony or maximum-likelihood camps; each largely ignoring the insights of the other. The result was an unfortunate detour in which the few perceived phylogenetic revolutions published by both sides of the philosophical divide were probably erroneous. The size of datasets has been growing exponentially since the mid-1980s accompanied by a wave of confidence that all relationships will soon be known. However, large datasets create new challenges, and a large number of genes does not guarantee reliable results. If history is a guide, then the quality of conclusions will be determined by an improved understanding of both molecular and morphological evolution, and not simply the number of genes analysed.

A molecular palaeobiological exploration of arthropod terrestrialization

Understanding animal terrestrialization, the process through which animals colonized the land, is crucial to clarify extant biodiversity and biological adaptation. Arthropoda (insects, spiders, centipedes and their allies) represent the largest majority of terrestrial biodiversity. Here we implemented a molecular palaeobiological approach, merging molecular and fossil evidence, to elucidate the deepest history of the terrestrial arthropods. We focused on the three independent, Palaeozoic arthropod terrestrialization events (those of Myriapoda, Hexapoda and Arachnida) and showed that a marine route to the colonization of land is the most likely scenario. Molecular clock analyses confirmed an origin for the three terrestrial lineages bracketed between the Cambrian and the Silurian. While molecular divergence times for Arachnida are consistent with the fossil record, Myriapoda are inferred to have colonized land earlier, substantially predating trace or body fossil evidence. An estimated origin of myriapods by the Early Cambrian precedes the appearance of embryophytes and perhaps even terrestrial fungi, raising the possibility that terrestrialization had independent origins in crown-group myriapod lineages, consistent with morphological arguments for convergence in tracheal systems.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.

Composition and genomic organization of arthropod Hox clusters

Background

The ancestral arthropod is believed to have had a clustered arrangement of ten Hox genes. Within arthropods, Hox gene mutations result in transformation of segment identities. Despite the fact that variation in segment number/character was common in the diversification of arthropods, few examples of Hox gene gains/losses have been correlated with morphological evolution. Furthermore, a full appreciation of the variation in the genomic arrangement of Hox genes in extant arthropods has not been recognized, as genome sequences from each major arthropod clade have not been reported until recently. Initial genomic analysis of the chelicerate Tetranychusurticae suggested that loss of Hox genes and Hox gene clustering might be more common than previously assumed. To further characterize the genomic evolution of arthropod Hox genes, we compared the genomic arrangement and general characteristics of Hox genes from representative taxa from each arthropod subphylum.

Results

In agreement with others, we find arthropods generally contain ten Hox genes arranged in a common orientation in the genome, with an increasing number of sampled species missing either Hox3 or abdominal-A orthologs. The genomic clustering of Hox genes in species we surveyed varies significantly, ranging from 0.3 to 13.6 Mb. In all species sampled, arthropod Hox genes are dispersed in the genome relative to the vertebrate Mus musculus. Differences in Hox cluster size arise from variation in the number of intervening genes, intergenic spacing, and the size of introns and UTRs. In the arthropods surveyed, Hox gene duplications are rare and four microRNAs are, in general, conserved in similar genomic positions relative to the Hox genes.

Conclusions

The tightly clustered Hox complexes found in the vertebrates are not evident within arthropods, and differential patterns of Hox gene dispersion are found throughout the arthropods. The comparative genomic data continue to support an ancestral arthropod Hox cluster of ten genes with a shared orientation, with four Hox gene-associated miRNAs, although the degree of dispersion between genes in an ancestral cluster remains uncertain. Hox3 and abdominal-A orthologs have been lost in multiple, independent lineages, and current data support a model in which inversions of the Abdominal-B locus that result in the loss of abdominal-A correlate with reduced trunk segmentation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0048-4) contains supplementary material, which is available to authorized users.

Hox genes, evo-devo, and the case of the ftz gene

The discovery of the broad conservation of embryonic regulatory genes across animal phyla, launched by the cloning of homeotic genes in the 1980s, was a founding event in the field of evolutionary developmental biology (evo-devo). While it had long been known that fundamental cellular processes, commonly referred to as housekeeping functions, are shared by animals and plants across the planet-processes such as the storage of information in genomic DNA, transcription, translation and the machinery for these processes, universal codon usage, and metabolic enzymes-Hox genes were different: mutations in these genes caused "bizarre" homeotic transformations of insect body parts that were certainly interesting but were expected to be idiosyncratic. The isolation of the genes responsible for these bizarre phenotypes turned out to be highly conserved Hox genes that play roles in embryonic patterning throughout Metazoa. How Hox genes have changed to promote the development of diverse body plans remains a central issue of the field of evo-devo today. For this Memorial article series, I review events around the discovery of the broad evolutionary conservation of Hox genes and the impact of this discovery on the field of developmental biology. I highlight studies carried out in Walter Gehring's lab and by former lab members that have continued to push the field forward, raising new questions and forging new approaches to understand the evolution of developmental mechanisms.

Evolution of epithelial morphogenesis: phenotypic integration across multiple levels of biological organization

Morphogenesis involves the dynamic reorganization of cell and tissue shapes to create the three-dimensional body. Intriguingly, different species have evolved different morphogenetic processes to achieve the same general outcomes during embryonic development. How are meaningful comparisons between species made, and where do the differences lie? In this Perspective, we argue that examining the evolution of embryonic morphogenesis requires the simultaneous consideration of different levels of biological organization: (1) genes, (2) cells, (3) tissues, and (4) the entire egg, or other gestational context. To illustrate the importance of integrating these levels, we use the extraembryonic epithelia of insects—a lineage-specific innovation and evolutionary hotspot—as an exemplary case study. We discuss how recent functional data, primarily from RNAi experiments targeting the Hox3/Zen and U-shaped group transcription factors, provide insights into developmental processes at all four levels. Comparisons of these data from several species both challenge and inform our understanding of homology, in assessing how the process of epithelial morphogenesis has itself evolved.

Hox gene duplications correlate with posterior heteronomy in scorpions

The evolutionary success of the largest animal phylum, Arthropoda, has been attributed to tagmatization, the coordinated evolution of adjacent metameres to form morphologically and functionally distinct segmental regions called tagmata. Specification of regional identity is regulated by the Hox genes, of which 10 are inferred to be present in the ancestor of arthropods. With six different posterior segmental identities divided into two tagmata, the bauplan of scorpions is the most heteronomous within Chelicerata. Expression domains of the anterior eight Hox genes are conserved in previously surveyed chelicerates, but it is unknown how Hox genes regionalize the three tagmata of scorpions. Here, we show that the scorpion Centruroides sculpturatus has two paralogues of all Hox genes except Hox3, suggesting cluster and/or whole genome duplication in this arachnid order. Embryonic anterior expression domain boundaries of each of the last four pairs of Hox genes (two paralogues each of Antp, Ubx, abd-A and Abd-B) are unique and distinguish segmental groups, such as pectines, book lungs and the characteristic tail, while maintaining spatial collinearity. These distinct expression domains suggest neofunctionalization of Hox gene paralogues subsequent to duplication. Our data reconcile previous understanding of Hox gene function across arthropods with the extreme heteronomy of scorpions.

Somatic and Germline Diversification of a Putative Immunoreceptor within One Phylum: Dscam in Arthropods

Arthropod Dscam, the homologue of the human Down Syndrome cell adhesion molecule, is a receptor used by the nervous and immune systems. Unlike in vertebrates, evolutionary pressure has selected and maintained a vast Dscam diversity of isoforms, known to specifying neuronal identity during the nervous system differentiation. This chapter examines the different modes of Dscam diversification in the context of arthropods' evolution and that of their immune system, where its role is controversial. In the single Dscam gene of insects and crustaceans, mutually exclusive alternative splicing affects three clusters of duplicated exons encoding the variable parts of the receptor. The Dscam gene produces over 10,000 isoforms. In the more basal arthropods such as centipedes, Dscam diversity results from a combination of many germline genes (over 80) with, in about half of those, the possibility of alternative splicing affecting only one exon cluster. In the even more basal arthropods, such as chelicerates, no splicing possibility is detected, but there exist dozens of germline Dscam genes. Compared to controlling the expression of multiple germline genes, the somatic mutually alternative splicing within a single gene may offer a simplified way of expressing a large Dscam repertoire. Expressed by hemocytes, Dscam is considered a phagocytic receptor but is also encountered in solution. More information is necessary about its binding to pathogens, its role in phagocytosis, its possible role in specifying hemocyte identity, its kinetics of expression, and the regulation of its RNA splicing to understand how its diversity is linked to immunity.

The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

Arthropods are the most abundant animals on earth. Among them, insects clearly dominate on land, whereas crustaceans hold the title for the most diverse invertebrates in the oceans. Much is known about the biology of these groups, not least because of genomic studies of the fruit fly Drosophila, the water flea Daphnia, and other species used in research. Here we report the first genome sequence from a species belonging to a lineage that has previously received very little attention—the myriapods. Myriapods were among the first arthropods to invade the land over 400 million years ago, and survive today as the herbivorous millipedes and venomous centipedes, one of which—Strigamia maritima—we have sequenced here. We find that the genome of this centipede retains more characteristics of the presumed arthropod ancestor than other sequenced insect genomes. The genome provides access to many aspects of myriapod biology that have not been studied before, suggesting, for example, that they have diversified receptors for smell that are quite different from those used by insects. In addition, it shows specific consequences of the largely subterranean life of this particular species, which seems to have lost the genes for all known light-sensing molecules, even though it still avoids light.

The ovary structure and oogenesis in the basal crustaceans and hexapods. Possible phylogenetic significance

Recent large-scale phylogenetic analyses of exclusively molecular or combined molecular and morphological characters support a close relationship between Crustacea and Hexapoda. The growing consensus on this phylogenetic link is reflected in uniting both taxa under the name Pancrustacea or Tetraconata. Several recent molecular phylogenies have also indicated that the monophyletic hexapods should be nested within paraphyletic crustaceans. However, it is still contentious exactly which crustacean taxon is the sister group to Hexapoda. Among the favored candidates are Branchiopoda, Malacostraca, Remipedia and Xenocarida (Remipedia + Cephalocarida). In this context, we review morphological and ultrastructural features of the ovary architecture and oogenesis in these crustacean groups in search of traits potentially suitable for phylogenetic considerations. We have identified a suite of morphological characters which may prove useful in further comparative studies.

Phylogenetic position of Myriapoda revealed by 454 transcriptome sequencing

Myriapods had been considered closely allied to hexapods (insects and relatives). However, analyses of molecular sequence data have consistently placed Myriapoda either as a sister group of Pancrustacea, comprising crustaceans and hexapods, and thereby supporting the monophyly of Mandibulata, or retrieved Myriapoda as a sister group of Chelicerata (spiders, ticks, mites and allies). In addition, the relationships among the four myriapod groups (Pauropoda, Symphyla, Diplopoda, Chilopoda) are unclear. To resolve the phylogeny of myriapods and their relationship to other main arthropod groups, we collected transcriptome data from the symphylan Symphylella vulgaris, the centipedes Lithobius forficatus and Scolopendra dehaani, and the millipedes Polyxenus lagurus, Glomeris pustulata and Polydesmus angustus by 454 sequencing. We concatenated a multiple sequence alignment that contained 1550 orthologous single copy genes (1,109,847 amino acid positions) from 55 euarthropod and 14 outgroup taxa. The final selected alignment included 181 genes and 37,425 amino acid positions from 55 taxa, with eight myriapods and 33 other euarthropods. Bayesian analyses robustly recovered monophyletic Mandibulata, Pancrustacea and Myriapoda. Most analyses support a sister group relationship of Symphyla in respect to a clade comprising Chilopoda and Diplopoda. Inclusion of additional sequence data from nine myriapod species resulted in an alignment with poor data density, but broader taxon average. With this dataset we inferred Diplopoda+Pauropoda as closest relatives (i.e., Dignatha) and recovered monophyletic Helminthomorpha. Molecular clock calculations suggest an early Cambrian emergence of Myriapoda ∼513 million years ago and a late Cambrian divergence of myriapod classes. This implies a marine origin of the myriapods and independent terrestrialization events during myriapod evolution.

Genomic sequence and experimental tractability of a new decapod shrimp model, Neocaridina denticulata

The speciose Crustacea is the largest subphylum of arthropods on the planet after the Insecta. To date, however, the only publically available sequenced crustacean genome is that of the water flea, Daphnia pulex, a member of the Branchiopoda. While Daphnia is a well-established ecotoxicological model, previous study showed that one-third of genes contained in its genome are lineage-specific and could not be identified in any other metazoan genomes. To better understand the genomic evolution of crustaceans and arthropods, we have sequenced the genome of a novel shrimp model, Neocaridina denticulata, and tested its experimental malleability. A library of 170-bp nominal fragment size was constructed from DNA of a starved single adult and sequenced using the Illumina HiSeq2000 platform. Core eukaryotic genes, the mitochondrial genome, developmental patterning genes (such as Hox) and microRNA processing pathway genes are all present in this animal, suggesting it has not undergone massive genomic loss. Comparison with the published genome of Daphnia pulex has allowed us to reveal 3750 genes that are indeed specific to the lineage containing malacostracans and branchiopods, rather than Daphnia-specific (E-value: 10⁻⁶). We also show the experimental tractability of N. denticulata, which, together with the genomic resources presented here, make it an ideal model for a wide range of further aquacultural, developmental, ecotoxicological, food safety, genetic, hormonal, physiological and reproductive research, allowing better understanding of the evolution of crustaceans and other arthropods.

Onychophoran Hox genes and the evolution of arthropod Hox gene expression

Introduction

Onychophora is a relatively small phylum within Ecdysozoa, and is considered to be the sister group to Arthropoda. Compared to the arthropods, that have radiated into countless divergent forms, the onychophoran body plan is overall comparably simple and does not display much in-phylum variation. An important component of arthropod morphological diversity consists of variation of tagmosis, i.e. the grouping of segments into functional units (tagmata), and this in turn is correlated with differences in expression patterns of the Hox genes. How these genes are expressed in the simpler onychophorans, the subject of this paper, would therefore be of interest in understanding their subsequent evolution in the arthropods, especially if an argument can be made for the onychophoran system broadly reflecting the ancestral state in the arthropods.

Results

The sequences and embryonic expression patterns of the complete set of ten Hox genes of an onychophoran (Euperipatoides kanangrensis) are described for the first time. We find that they are all expressed in characteristic patterns that suggest a function as classical Hox genes. The onychophoran Hox genes obey spatial colinearity, and with the exception of Ultrabithorax (Ubx), they all have different and distinct anterior expression borders. Notably, Ubx transcripts form a posterior to anterior gradient in the onychophoran trunk. Expression of all onychophoran Hox genes extends continuously from their anterior border to the rear end of the embryo.

Conclusions

The spatial expression pattern of the onychophoran Hox genes may contribute to a combinatorial Hox code that is involved in giving each segment its identity. This patterning of segments in the uniform trunk, however, apparently predates the evolution of distinct segmental differences in external morphology seen in arthropods. The gradient-like expression of Ubx may give posterior segments their specific identity, even though they otherwise express the same set of Hox genes. We suggest that the confined domains of Hox gene expression seen in arthropods evolved from an ancestral onychophoran-like Hox gene pattern. Reconstruction of the ancestral arthropod Hox pattern and comparison with the patterns in the different arthropod classes reveals phylogenetic support for Mandibulata and Tetraconata, but not Myriochelata and Atelocerata.

Embryonic neurogenesis in Pseudopallene sp. (Arthropoda, Pycnogonida) includes two subsequent phases with similarities to different arthropod groups

BACKGROUND

Studies on early neurogenesis have had considerable impact on the discussion of the phylogenetic relationships of arthropods, having revealed striking similarities and differences between the major lineages. In Hexapoda and crustaceans, neurogenesis involves the neuroblast, a type of neural stem cell. In each hemi-segment, a set of neuroblasts produces neural cells by repeated asymmetrical and interiorly directed divisions. In Euchelicerata and Myriapoda, neurogenesis lacks neural stem cells, featuring instead direct immigration of neural cell groups from fixed sites in the neuroectoderm. Accordingly, neural stem cells were hitherto assumed to be an evolutionary novelty of the Tetraconata (Hexapoda + crustaceans). To further test this hypothesis, we investigated neurogenesis in Pycnogonida, or sea spiders, a group of marine arthropods with close affinities to euchelicerates.

RESULTS

We studied neurogenesis during embryonic development of Pseudopallene sp. (Callipallenidae), using fluorescent histochemical staining and immunolabelling. Embryonic neurogenesis has two phases. The first phase shows notable similarities to euchelicerates and myriapods. These include i) the lack of morphologically different cell types in the neuroectoderm; ii) the formation of transiently identifiable, stereotypically arranged cell internalization sites; iii) immigration of predominantly post-mitotic ganglion cells; and iv) restriction of tangentially oriented cell proliferation to the apical cell layer. However, in the second phase, the formation of a central invagination in each hemi-neuromere is accompanied by the differentiation of apical neural stem cells. The latter grow in size, show high mitotic activity and an asymmetrical division mode. A marked increase of ganglion cell numbers follows their differentiation. Directly basal to the neural stem cells, an additional type of intermediate neural precursor is found.

CONCLUSIONS

Embryonic neurogenesis of Pseudopallene sp. combines features of central nervous system development that have been hitherto described separately in different arthropod taxa. The two-phase character of pycnogonid neurogenesis calls for a thorough reinvestigation of other non-model arthropods over the entire course of neurogenesis. With the currently available data, a common origin of pycnogonid neural stem cells and tetraconate neuroblasts remains unresolved. To acknowledge this, we present two possible scenarios on the evolution of arthropod neurogenesis, whereby Myriapoda play a key role in the resolution of this issue.

Molecular phylogenetic analyses support the monophyly of Hexapoda and suggest the paraphyly of Entognatha

BACKGROUND

Molecular phylogenetic analyses have revealed that Hexapoda and Crustacea form a common clade (the Pancrustacea), which is now widely accepted among zoologists; however, the origin of Hexapoda remains unresolved. The main problems are the unclear relationships among the basal hexapod lineages, Protura (proturans), Collembola (springtails), Diplura (diplurans), and Ectognatha (bristletails, silverfishes, and all winged insects). Mitogenomic analyses have challenged hexapod monophyly and suggested the reciprocal paraphyly of Hexapoda and Crustacea, whereas studies based on nuclear molecular data support the monophyletic origin of hexapods. Additionally, there are significant discrepancies with respect to these issues between the results of morphological and molecular studies. To investigate these problems, we performed phylogenetic analyses of Pancrustacea based on the protein sequences of three orthologous nuclear genes encoding the catalytic subunit of DNA polymerase delta and the largest and second largest subunits of RNA polymerase II from 64 species of arthropods, including representatives of all hexapod orders.

RESULTS

Phylogenetic analyses were conducted based on the inferred amino acid (aa) sequences (~3400 aa in total) of the three genes using the maximum likelihood (ML) method and Bayesian inference. Analyses were also performed with additional datasets generated by excluding long-branch taxa or by using different outgroups. These analyses all yielded essentially the same results. All hexapods were clustered into a common clade, with Branchiopoda as its sister lineage, whereas Crustacea was paraphyletic. Within Hexapoda, the lineages Ectognatha, Palaeoptera, Neoptera, Polyneoptera, and Holometabola were each confirmed to be monophyletic with robust support, but monophyly was not supported for Entognatha (Protura + Collembola + Diplura), Ellipura (Protura + Collembola), or Nonoculata (Protura + Diplura). Instead, our results showed that Protura is the sister lineage to all other hexapods and that Diplura or Diplura + Collembola is closely related to Ectognatha.

CONCLUSION

This is the first study to include all hexapod orders in a phylogenetic analysis using multiple nuclear protein-coding genes to investigate the phylogeny of Hexapoda, with an emphasis on Entognatha. The results strongly support the monophyletic origin of hexapods but reject the monophyly of Entognatha, Ellipura, and Nonoculata. Our results provided the first molecular evidence in support of Protura as the sister group to other hexapods. These findings are expected to provide additional insights into the origin of hexapods and the processes involved in the adaptation of insects to life on land.

Posterior Hox gene reduction in an arthropod: Ultrabithorax and Abdominal-B are expressed in a single segment in the mite Archegozetes longisetosus

Background

Hox genes encode transcription factors that have an ancestral role in all bilaterian animals in specifying regions along the antero-posterior axis. In arthropods (insects, crustaceans, myriapods and chelicerates), Hox genes function to specify segmental identity, and changes in Hox gene expression domains in different segments have been causal to the evolution of novel arthropod morphologies. Despite this, the roles of Hox genes in arthropods that have secondarily lost or reduced their segmental composition have been relatively unexplored. Recent data suggest that acariform mites have a reduced segmental component of their posterior body tagma, the opisthosoma, in that only two segments are patterned during embryogenesis. This is in contrast to the observation that in many extinct and extant chelicerates (that is, horseshoe crabs, scorpions, spiders and harvestmen) the opisthosoma is comprised of ten or more segments. To explore the role of Hox genes in this reduced body region, we followed the expression of the posterior-patterning Hox genes Ultrabithorax (Ubx) and Abdominal-B (Abd-B), as well as the segment polarity genes patched (ptc) and engrailed (en), in the oribatid mite Archegozetes longisetosus.

Results

We find that the expression patterns of ptc are in agreement with previous reports of a reduced mite opisthosoma. In comparison to the ptc and en expression patterns, we find that Ubx and Abd-B are expressed in a single segment in A. longisetosus, the second opisthosomal segment. Abd-B is initially expressed more posteriorly than Ubx, that is, into the unsegmented telson; however, this domain clears in subsequent stages where it remains in the second opisthosomal segment.

Conclusions

Our findings suggest that Ubx and Abd-B are expressed in a single segment in the opisthosoma. This is a novel observation, in that these genes are expressed in several segments in all studied arthropods. These data imply that a reduction in opisthosomal segmentation may be tied to a dramatically reduced Hox gene input in the opisthosoma.

Isolation of Hox cluster genes from insects reveals an accelerated sequence evolution rate

Among gene families it is the Hox genes and among metazoan animals it is the insects (Hexapoda) that have attracted particular attention for studying the evolution of development. Surprisingly though, no Hox genes have been isolated from 26 out of 35 insect orders yet, and the existing sequences derive mainly from only two orders (61% from Hymenoptera and 22% from Diptera). We have designed insect specific primers and isolated 37 new partial homeobox sequences of Hox cluster genes (lab, pb, Hox3, ftz, Antp, Scr, abd-a, Abd-B, Dfd, and Ubx) from six insect orders, which are crucial to insect phylogenetics. These new gene sequences provide a first step towards comparative Hox gene studies in insects. Furthermore, comparative distance analyses of homeobox sequences reveal a correlation between gene divergence rate and species radiation success with insects showing the highest rate of homeobox sequence evolution.

The homeodomain complement of the ctenophore Mnemiopsis leidyi suggests that Ctenophora and Porifera diverged prior to the ParaHoxozoa

Background

The much-debated phylogenetic relationships of the five early branching metazoan lineages (Bilateria, Cnidaria, Ctenophora, Placozoa and Porifera) are of fundamental importance in piecing together events that occurred early in animal evolution. Comparisons of gene content between organismal lineages have been identified as a potentially useful methodology for phylogenetic reconstruction. However, these comparisons require complete genomes that, until now, did not exist for the ctenophore lineage. The homeobox superfamily of genes is particularly suited for these kinds of gene content comparisons, since it is large, diverse, and features a highly conserved domain.

Results

We have used a next-generation sequencing approach to generate a high-quality rough draft of the genome of the ctenophore Mnemiopsis leidyi and subsequently identified a set of 76 homeobox-containing genes from this draft. We phylogenetically categorized this set into established gene families and classes and then compared this set to the homeodomain repertoire of species from the other four early branching metazoan lineages. We have identified several important classes and subclasses of homeodomains that appear to be absent from Mnemiopsis and from the poriferan Amphimedon queenslandica. We have also determined that, based on lineage-specific paralog retention and average branch lengths, it is unlikely that these missing classes and subclasses are due to extensive gene loss or unusually high rates of evolution in Mnemiopsis.

Conclusions

This paper provides a first glimpse of the first sequenced ctenophore genome. We have characterized the full complement of Mnemiopsis homeodomains from this species and have compared them to species from other early branching lineages. Our results suggest that Porifera and Ctenophora were the first two extant lineages to diverge from the rest of animals. Based on this analysis, we also propose a new name - ParaHoxozoa - for the remaining group that includes Placozoa, Cnidaria and Bilateria.

A congruent solution to arthropod phylogeny: phylogenomics, microRNAs and morphology support monophyletic Mandibulata

While a unique origin of the euarthropods is well established, relationships between the four euarthropod classes-chelicerates, myriapods, crustaceans and hexapods-are less clear. Unsolved questions include the position of myriapods, the monophyletic origin of chelicerates, and the validity of the close relationship of euarthropods to tardigrades and onychophorans. Morphology predicts that myriapods, insects and crustaceans form a monophyletic group, the Mandibulata, which has been contradicted by many molecular studies that support an alternative Myriochelata hypothesis (Myriapoda plus Chelicerata). Because of the conflicting insights from published molecular datasets, evidence from nuclear-coding genes needs corroboration from independent data to define the relationships among major nodes in the euarthropod tree. Here, we address this issue by analysing two independent molecular datasets: a phylogenomic dataset of 198 protein-coding genes including new sequences for myriapods, and novel microRNA complements sampled from all major arthropod lineages. Our phylogenomic analyses strongly support Mandibulata, and show that Myriochelata is a tree-reconstruction artefact caused by saturation and long-branch attraction. The analysis of the microRNA dataset corroborates the Mandibulata, showing that the microRNAs miR-965 and miR-282 are present and expressed in all mandibulate species sampled, but not in the chelicerates. Mandibulata is further supported by the phylogenetic analysis of a comprehensive morphological dataset covering living and fossil arthropods, and including recently proposed, putative apomorphies of Myriochelata. Our phylogenomic analyses also provide strong support for the inclusion of pycnogonids in a monophyletic Chelicerata, a paraphyletic Cycloneuralia, and a common origin of Arthropoda (tardigrades, onychophorans and arthropods), suggesting that previous phylogenies grouping tardigrades and nematodes may also have been subject to tree-reconstruction artefacts.

Gene expression suggests conserved aspects of Hox gene regulation in arthropods and provides additional support for monophyletic Myriapoda

Antisense transcripts of Ultrabithorax (aUbx) in the millipede Glomeris and the centipede Lithobius are expressed in patterns complementary to that of the Ubx sense transcripts. A similar complementary expression pattern has been described for non-coding RNAs (ncRNAs) of the bithoraxoid (bxd) locus in Drosophila, in which the transcription of bxd ncRNAs represses Ubx via transcriptional interference. We discuss our findings in the context of possibly conserved mechanisms of Ubx regulation in myriapods and the fly.

Bicistronic transcription of Ubx and Antennapedia (Antp) has been reported previously for a myriapod and a number of crustaceans. In this paper, we show that Ubx/Antp bicistronic transcripts also occur in Glomeris and an onychophoran, suggesting further conserved mechanisms of Hox gene regulation in arthropods.

Myriapod monophyly is supported by the expression of aUbx in all investigated myriapods, whereas in other arthropod classes, including the Onychophora, aUbx is not expressed. Of the two splice variants of Ubx/Antp only one could be isolated from myriapods, representing a possible further synapomorphy of the Myriapoda.

Evolution of a genomic regulatory domain: the role of gene co-option and gene duplication in the Enhancer of split complex

The Drosophila Enhancer of split complex [E(spl)-C] is a remarkable complex of genes many of which are effectors or modulators of Notch signaling. The complex contains different classes of genes including four bearded genes and seven basic helix-loop-helix (bHLH) genes. We examined the evolution of this unusual complex by identifying bearded and bHLH genes in the genome sequences of Arthropods. We find that a four-gene E(spl)-C, containing three bHLH genes and one bearded gene, is an ancient component of the genomes of Crustacea and Insects. The complex is well conserved in insects but is highly modified in Drosophila, where two of the ancestral genes of the complex are missing, and the remaining two have been duplicated multiple times. Through examining the expression of E(spl)-C genes in honeybees, aphids, and Drosophila, we determined that the complex ancestrally had a role in Notch signaling. The expression patterns of genes found inserted into the complex in some insects, or that of ancestral E(spl)-C genes that have moved out of the complex, imply that the E(spl)-C is a genomic domain regulated as a whole by Notch signaling. We hypothesize that the E(spl)-C is a Notch-regulated genomic domain conserved in Arthropod genomes for around 420 million years. We discuss the consequence of this conserved domain for the recruitment of novel genes into the Notch signaling cascade.

Arthropod phylogeny: an overview from the perspectives of morphology, molecular data and the fossil record

Monophyly of Arthropoda is emphatically supported from both morphological and molecular perspectives. Recent work finds Onychophora rather than Tardigrada to be the closest relatives of arthropods. The status of tardigrades as panarthropods (rather than cycloneuralians) is contentious from the perspective of phylogenomic data. A grade of Cambrian taxa in the arthropod stem group includes gilled lobopodians, dinocaridids (e.g., anomalocaridids), fuxianhuiids and canadaspidids that inform on character acquisition between Onychophora and the arthropod crown group. A sister group relationship between Crustacea (itself likely paraphyletic) and Hexapoda is retrieved by diverse kinds of molecular data and is well supported by neuroanatomy. This clade, Tetraconata, can be dated to the early Cambrian by crown group-type mandibles. The rival Atelocerata hypothesis (Myriapoda+Hexapoda) has no molecular support. The basal node in the arthropod crown group is embroiled in a controversy over whether myriapods unite with chelicerates (Paradoxopoda or Myriochelata) or with crustaceans and hexapods (Mandibulata). Both groups find some molecular and morphological support, though Mandibulata is presently the stronger morphological hypothesis. Either hypothesis forces an unsampled ghost lineage for Myriapoda from the Cambrian to the mid Silurian.

Velvet worm development links myriapods with chelicerates

Despite the advent of modern molecular and computational methods, the phylogeny of the four major arthropod groups (Chelicerata, Myriapoda, Crustacea and Hexapoda, including the insects) remains enigmatic. One particular challenge is the position of myriapods as either the closest relatives to chelicerates (Paradoxopoda/Myriochelata hypothesis), or to crustaceans and hexapods (Mandibulata hypothesis). While neither hypothesis receives conclusive support from molecular analyses, most morphological studies favour the Mandibulata concept, with the mandible being the most prominent feature of this group. Although no morphological evidence was initially available to support the Paradoxopoda hypothesis, a putative synapomorphy of chelicerates and myriapods has recently been put forward based on studies of neurogenesis. However, this and other morphological characters remain of limited use for phylogenetic systematics owing to the lack of data from an appropriate outgroup. Here, we show that several embryonic characters are synapomorphies uniting the chelicerates and myriapods, as revealed by an outgroup comparison with the Onychophora or velvet worms. Our findings, thus provide, to our knowledge, first morphological/embryological support for the monophyly of the Paradoxopoda and suggest that the mandible might have evolved twice within the arthropods.

A 454 sequencing approach for large scale phylogenomic analysis of the common emperor scorpion (Pandinus imperator)

In recent years, phylogenetic tree reconstructions that rely on multiple gene alignments that had been deduced from expressed sequence tags (ESTs) have become a popular method in molecular systematics. Here, we present a 454 pyrosequencing approach to infer the transcriptome of the Emperor scorpion Pandinus imperator. We obtained 428,844 high-quality reads (mean length=223+/-50 b) from total cDNA, which were assembled into 8334 contigs (mean length 422+/-313 bp) and 26,147 singletons. About 1200 contigs were successfully annotated by BLAST and orthology search. Specific analyses of eight distinct hemocyanin sequences provided further proof for the quality of the 454 reads and the assembly process. The P. imperator sequences were included in a concatenated alignment of 149 orthologous genes of 67 metazoan taxa that covers 39,842 amino acids. After removal of low-quality regions, 11,168 positions were employed for phylogenetic reconstructions. Using Bayesian and maximum likelihood methods, we obtained strongly supported monophyletic Ecdysozoa, Arthropoda (excluding Tardigrada), Euarthropoda, Pancrustacea and Hexapoda. We also recovered the Myriochelata (Chelicerata+Myriapoda). Within the chelicerates, Pycnogonida form the sister group of Euchelicerata. However, Arachnida were found paraphyletic because the Acari (mites and ticks) were recovered as sister group of a clade comprising Xiphosura, Scorpiones and Araneae. In summary, we have shown that 454 pyrosequencing is a cost-effective method that provides sufficient data and coverage depth for gene detection and multigene-based phylogenetic analyses.

The involvement of engrailed and wingless during segmentation in the onychophoran Euperipatoides kanangrensis (Peripatopsidae: Onychophora) (Reid 1996)

As the putative sister group to the arthropods, onychophorans can provide insight into ancestral developmental mechanisms in the panarthropod clade. Here, we examine the expression during segmentation of orthologues of wingless (Wnt1) and engrailed, two genes that play a key role in defining segment boundaries in Drosophila and that appear to play a role in segmentation in many other arthropods. Both are expressed in segmentally reiterated stripes in all forming segments except the first (brain) segment, which only shows an engrailed stripe. Engrailed is expressed before segments are morphologically visible and is expressed in both mesoderm and ectoderm. Segmental wingless expression is not detectable until after mesodermal somites are clearly distinct. Early engrailed expression lies in and extends to both sides of the furrow that first demarcates segments in the ectoderm, but is largely restricted to the posterior part of somites. Wingless expression lies immediately anterior to engrailed expression, as it does in many arthropods, but there is no precise cellular boundary between the two expression domains analogous to the overt parasegment boundary seen in Drosophila. Engrailed stripes extend along the posterior part of each limb bud, including the antenna, while wingless is restricted to the distal tip of the limbs and the neurectoderm basal to the limbs.

Isolation and expression of Pax6 and atonal homologues in the American horseshoe crab, Limulus polyphemus

Pax6 regulates eye development in many animals. In addition, Pax6 activates atonal transcription factors in both invertebrate and vertebrate eyes. Here, we investigate the roles of Pax6 and atonal during embryonic development of Limulus polyphemus rudimentary lateral, medial and ventral eyes, and the initiation of lateral ommatidial eye and medial ocelli formation. Limulus eye development is of particular interest because these animals hold a unique position in arthropod phylogeny and possess multiple eye types. Furthermore, the molecular underpinnings of eye development have yet to be investigated in chelicerates. We characterized a Limulus Pax6 gene, with multiple splice products and predicted protein isoforms, and one atonal homologue. Unexpectedly, neither gene is expressed in the developing eye types examined, although both genes are present in the lateral sense organ, a structure of unknown function.

The complete mitochondrial genome of Symphylella sp. (Myriapoda: Symphyla): Extensive gene order rearrangement and evidence in favor of Progoneata

We determined the complete 14,667bp mitochondrial DNA sequence of Symphylella sp., the first representative of the Scolopendrellidae (Arthropoda: Myriapoda: Symphyla). With respect to the ancestral arthropod mitochondrial gene order, two protein-coding genes, the rRNAs and 10 of the tRNAs appear to be rearranged. This rearrangement is novel in the arthropods and genes with identical transcriptional polarity are clustered except for trnE, trnN and putative control region (CR), resembling two previously reported diplopod genomes. A duplication/loss (random and non-random)-recombination model was proposed to account for the generation of the gene order in Symphylella sp. All phylogenetic analysis yielded strong support for a clade of Symphyla plus Diplopoda, i.e., Progoneata. However, the phylogenetic position of Myriapoda within Arthropoda remains unclear. The amino acid dataset gives strong support for an affinity to Pancrustacea, while the nucleotide dataset weakly supports Myriapoda grouped with Chelicerata.

Phylogenetic relationships and the evolution of regulatory gene sequences in the parrotfishes

Regulatory genes control the expression of other genes and are key components of developmental processes such as segmentation and embryonic construction of the skull in vertebrates. Here we examine the variability and evolution of three vertebrate regulatory genes, addressing issues of their utility for phylogenetics and comparing the rates of genetic change seen in regulatory loci to the rates seen in other genes in the parrotfishes. The parrotfishes are a diverse group of colorful fishes from coral reefs and seagrasses worldwide and have been placed phylogenetically within the family Labridae. We tested phylogenetic hypotheses among the parrotfishes, with a focus on the genera Chlorurus and Scarus, by analyzing eight gene fragments for 42 parrotfishes and eight outgroup species. We sequenced mitochondrial 12s rRNA (967 bp), 16s rRNA (577 bp), and cytochrome b (477 bp). From the nuclear genome, we sequenced part of the protein-coding genes rag2 (715 bp), tmo4c4 (485 bp), and the developmental regulatory genes otx1 (672 bp), bmp4 (488bp), and dlx2 (522 bp). Bayesian, likelihood, and parsimony analyses of the resulting 4903 bp of DNA sequence produced similar topologies that confirm the monophyly of the scarines and provide a phylogeny at the species level for portions of the genera Scarus and Chlorurus. Four major clades of Scarus were recovered, with three distributed in the Indo-Pacific and one containing Caribbean/Atlantic taxa. Molecular rates suggest a Miocene origin of the parrotfishes (22 mya) and a recent divergence of species within Scarus and Chlorurus, within the past 5 million years. Developmentally important genes made a significant contribution to phylogenetic structure, and rates of genetic evolution were high in bmp4, similar to other coding nuclear genes, but low in otx1 and the dlx2 exons. Synonymous and non-synonymous substitution patterns in developmental regulatory genes support the hypothesis of stabilizing selection during the history of these genes, with several phylogenetic regions of accelerated non-synonymous change detected in the phylogeny.