Organization of the Hox gene cluster of the silkworm, Bombyx mori: a split of the Hox cluster in a non-Drosophila insect

Evolutionary plasticity in nematode Hox gene complements and genomic loci arrangement

Hox genes are central to metazoan body plan formation, patterning and evolution, playing a critical role in cell fate decisions early in embryonic development in invertebrates and vertebrates. While the archetypical Hox gene cluster consists of members of nine ortholog groups (HOX1-HOX9), arrayed in close linkage in the order in which they have their anterior-posterior patterning effects, nematode Hox gene sets do not fit this model. The Caenorhabditis elegans Hox gene set is not clustered and contains only six Hox genes from four of the ancestral groups. The pattern observed in C. elegans is not typical of the phylum, and variation in orthologue set presence and absence and in genomic organisation has been reported. Recent advances in genome sequencing have resulted in the availability of many novel genome assemblies in Nematoda, especially from taxonomic groups that had not been analysed previously. Here, we explored Hox gene complements in high-quality genomes of 80 species from all major clades of Nematoda to understand the evolution of this key set of body pattern genes and especially to probe the origins of the "dispersed" cluster observed in C. elegans. We also included the recently available high-quality genomes of some Nematomorpha as an outgroup. We find that nematodes can have Hox genes from up to six orthology groups. While nematode Hox "clusters" are often interrupted by unrelated genes we identify species in which the cluster is intact and not dispersed.

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.

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.

Characterizing Hox genes in mayflies (Ephemeroptera), with Hexagenia limbata as a new mayfly model

BACKGROUND

Hox genes are key regulators of appendage development in the insect body plan. The body plan of mayfly (Ephemeroptera) nymphs differs due to the presence of abdominal appendages called gills. Despite mayflies' phylogenetic position in Paleoptera and novel morphology amongst insects, little is known of their developmental genetics, such as the appendage-regulating Hox genes. To address this issue we present an annotated, early instar transcriptome and embryonic expression profiles for Antennapedia, Ultrabithorax, and Abdominal A proteins in the mayfly Hexagenia limbata, identify putative Hox protein sequences in the mayflies H. limbata, Cloeon dipterum, and Ephemera danica, and describe the genomic organization of the Hox gene cluster in E. danica.

RESULTS

Transcriptomic sequencing of early instar H. limbata nymphs yielded a high-quality assembly of 83,795 contigs, of which 22,975 were annotated against Folsomia candida, Nilaparvata lugens, Zootermopsis nevadensis and UniRef90 protein databases. Homeodomain protein phylogeny and peptide annotations identified coding sequences for eight of the ten canonical Hox genes (excluding zerknüllt/Hox3 and fushi tarazu) in H. limbata and C. dipterum, and all ten in E. danica. Mayfly Hox protein sequences and embryonic expression patterns of Antp, Ubx, and Abd-A appear highly conserved with those seen in other non-holometabolan insects. Similarly, the genomic organization of the Hox cluster in E. danica resembles that seen in most insects.

CONCLUSIONS

We present evidence that mayfly Hox peptide sequences and the embryonic expression patterns for Antp, Ubx, and Abd-A are extensively conserved with other insects, as is organization of the mayfly Hox gene cluster. The protein data suggest mayfly Antp, Ubx, and Abd-A play appendage promoting and repressing roles during embryogenesis in the thorax and abdomen, respectively, as in other insects. The identified expression of eight Hox genes, including Ubx and abd-A, in early instar nymphs further indicates a post-embryonic role, possibly in gill development. These data provide a basis for H. limbata as a complementary Ephemeridae model to the growing repertoire of mayfly model species and molecular techniques.

High-Quality Genome Assemblies Reveal Long Non-coding RNAs Expressed in Ant Brains

Ants are an emerging model system for neuroepigenetics, as embryos with virtually identical genomes develop into different adult castes that display diverse physiology, morphology, and behavior. Although a number of ant genomes have been sequenced to date, their draft quality is an obstacle to sophisticated analyses of epigenetic gene regulation. We reassembled de novo high-quality genomes for two ant species, Camponotus floridanus and Harpegnathos saltator. Using long reads enabled us to span large repetitive regions and improve genome contiguity, leading to comprehensive and accurate protein-coding annotations that facilitated the identification of a Gp-9-like gene as differentially expressed in Harpegnathos castes. The new assemblies also enabled us to annotate long non-coding RNAs in ants, revealing caste-, brain-, and developmental-stage-specific long non-coding RNAs (lncRNAs) in Harpegnathos. These upgraded genomes, along with the new gene annotations, will aid future efforts to identify epigenetic mechanisms of phenotypic and behavioral plasticity in ants.

Analysis of the genome of the New Zealand giant collembolan (Holacanthella duospinosa) sheds light on hexapod evolution

BACKGROUND

The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events.

RESULTS

Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180×. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes.

CONCLUSIONS

The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda.

Ultrabithorax and abdominal-A specify the abdominal appendage in a dosage-dependent manner in silkworm, Bombyx mori

In insects, there is a considerable diversity in leg distribution on the body, including number, segmental arrangement, morphological identity and consequent function, but the genetic basis for these differences is not well understood. Here by positional cloning, we showed that a ~355 kb region, including Bombyx mori Ultrabithorax (BmUbx) and abdominal-A (Bmabd-A), was responsible for the silkworm mutant Kh-extra-crescents-like (E^Kh-l) that displayed additional thoracic limb-like legs on the first abdominal segment (A1) and occasionally on the second abdominal segment (A2). We found that BmUbx gene was downregulated at both messenger RNA level and protein level in E^Kh-l embryo, while its expression domain in the E^Kh-l embryo was almost the same as that in the wild type. Whereas Bmabd-A was upregulated at both levels and was ectopically overexpressed on the supernumerary leg-bearing segments in E^Kh-l. Compared with the previously reported E^cs-l mutant in which increased expression of both BmUbx and Bmabd-A gave rise to ectopic proleg-like appendages on the same segments, we propose that overexpressed Bmabd-A gene is capable to promote the outgrowth of extra leg appendages on A1 and A2 segments, whereas BmUbx gene is required to specify accurate morphologies of the ectopic legs in a dosage-dependent manner in silkworm. These results provide insights into how these hox genes regulate the leg morphologic diversity on the same segments.

Fine Mapping of a Degenerated Abdominal Legs Mutant (Edl) in Silkworm, Bombyx mori

In insects, abdominal appendages, also called prolegs, vary due to adaptive evolution. Mutations on prolegs within species provide insights to better understand the mechanisms underlying appendage development and diversity. In silkworm Bombyx mori, extra-crescents and degenerated abdominal legs (E^dl) mutant, belonging to the E pseudoallele group, is a spontaneous mutation that adds crescents and degenerates prolegs on the third abdominal segment (A3). This mutation may be a homeotic transformation of A3 to A2. In this study, the E^dl locus was mapped within approximately a 211 Kb region that is 10 Kb upstream of Bmabdominal-A (Bmabd-A). RT-quantitative PCR (RT-qPCR) and Western blot analysis of Bmabd-A expression showed a slight but significant decrease, while the expression of BmUltrabithorax (BmUbx) was up-regulated in the E^dl mutant compared to wildtype (Dazao). Moreover, we also found that BmDistal-less (BmDll), which regulated the development of distal proleg structures, was missing at the tips of the A3 prolegs in the E^dl mutant compared to BmDll expression in normally developed prolegs in both the wildtype and mutant. Collectively, we identified approximately a 211 Kb region in the E^dl locus that regulates BmUbx and Bmabd-A expression and found that changes in BmUbx and Bmabd-A expression may lead to the loss of distal proleg structures in B. mori.

Regulation of Silk Genes by Hox and Homeodomain Proteins in the Terminal Differentiated Silk Gland of the Silkworm Bombyx mori

The silk gland of the silkworm Bombyx mori is a long tubular organ that is divided into several subparts along its anteroposterior (AP) axis. As a trait of terminal differentiation of the silk gland, several silk protein genes are expressed with unique regional specificities. Most of the Hox and some of the homeobox genes are also expressed in the differentiated silk gland with regional specificities. The expression patterns of Hox genes in the silk gland roughly correspond to those in embryogenesis showing "colinearity". The central Hox class protein Antennapedia (Antp) directly regulates the expression of several middle silk gland-specific silk genes, whereas the Lin-1/Isl-1/Mec3 (LIM)-homeodomain transcriptional factor Arrowhead (Awh) regulates the expression of posterior silk gland-specific genes for silk fiber proteins. We summarize our results and discuss the usefulness of the silk gland of Bombyx mori for analyzing the function of Hox genes. Further analyses of the regulatory mechanisms underlying the region-specific expression of silk genes will provide novel insights into the molecular bases for target-gene selection and regulation by Hox and homeodomain proteins.

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.

Comprehensive analysis of Hox gene expression in the amphipod crustacean Parhyale hawaiensis

Hox genes play crucial roles in establishing regional identity along the anterior-posterior axis in bilaterian animals, and have been implicated in generating morphological diversity throughout evolution. Here we report the identification, expression, and initial genomic characterization of the complete set of Hox genes from the amphipod crustacean Parhyale hawaiensis. Parhyale is an emerging model system that is amenable to experimental manipulations and evolutionary comparisons among the arthropods. Our analyses indicate that the Parhyale genome contains a single copy of each canonical Hox gene with the exception of fushi tarazu, and preliminary mapping suggests that at least some of these genes are clustered together in the genome. With few exceptions, Parhyale Hox genes exhibit both temporal and spatial colinearity, and expression boundaries correlate with morphological differences between segments and their associated appendages. This work represents the most comprehensive analysis of Hox gene expression in a crustacean to date, and provides a foundation for functional studies aimed at elucidating the role of Hox genes in arthropod development and evolution.

Conservation and variation in Hox genes: how insect models pioneered the evo-devo field

Evolutionary developmental biology, or evo-devo, broadly investigates how body plan diversity and morphological novelties have arisen and persisted in nature. The discovery of Hox genes in Drosophila, and their subsequent identification in most other metazoans, led biologists to try to understand how embryonic genes crucial for proper development have changed to promote the vast morphological variation seen in nature. Insects are ideal model systems for studying this diversity and the mechanisms underlying it because phylogenetic relationships are well established, powerful genetic tools have been developed, and there are many examples of evolutionary specializations that have arisen in nature in different insect lineages, such as the jumping leg of orthopterans and the helmet structures of treehoppers. Here, we briefly introduce the field of evo-devo and Hox genes, discuss functional tools available to study early developmental genes in insects, and provide examples in which changes in Hox genes have contributed to changes in body plan or morphology.

Conserved boundary elements from the Hox complex of mosquito, Anopheles gambiae

The conservation of hox genes as well as their genomic organization across the phyla suggests that this system of anterior–posterior axis formation arose early during evolution and has come under strong selection pressure. Studies in the split Hox cluster of Drosophila have shown that proper expression of hox genes is dependent on chromatin domain boundaries that prevent inappropriate interactions among different types of cis-regulatory elements. To investigate whether boundary function and their role in regulation of hox genes is conserved in insects with intact Hox clusters, we used an algorithm to locate potential boundary elements in the Hox complex of mosquito, Anopheles gambiae. Several potential boundary elements were identified that could be tested for their functional conservation. Comparative analysis revealed that like Drosophila, the bithorax region in A. gambiae contains an extensive array of boundaries and enhancers organized into domains. We analysed a subset of candidate boundary elements and show that they function as enhancer blockers in Drosophila. The functional conservation of boundary elements from mosquito in fly suggests that regulation of hox genes involving chromatin domain boundaries is an evolutionary conserved mechanism and points to an important role of such elements in key developmentally regulated loci.

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.

Evolutionary changes of the target sites of two microRNAs encoded in the Hox gene cluster of Drosophila and other insect species

MicroRNAs (miRs) are noncoding RNAs that regulate gene expression at the post-transcriptional level. In animals, the target sites of a miR are generally located in the 3′ untranslated regions (UTRs) of messenger RNAs. However, how the target sites change during evolution is largely unknown. MiR-iab-4 and miR-iab-4as are known to regulate the expression of two Hox genes, Abd-A and Ubx, in Drosophila melanogaster. We have therefore studied the evolutionary changes of these two miR genes and their target sites of the Hox genes in Drosophila, other insect species, and Daphnia. Our homology search identified a single copy of each miR gene located in the same genomic position of the Hox gene cluster in all species examined. The seed nucleotide sequence was also the same for all species. Searching for the target sites in all Hox genes, we found several target sites of miR-iab-4 and miR-iab-4as in Antp in addition to Abd-A and Ubx in most insect species examined. Our phylogenetic analysis of target sites in Abd-A, Ubx, and Antp showed that the old target sites, which existed before the divergence of the 12 Drosophila species, have been well maintained in most species under purifying selection. By contrast, new target sites, which were generated during Drosophila evolution, were often lost in some species and mostly located in unalignable regions of the 3′ UTRs. These results indicate that these regions can be a potential source of generating new target sites, which results in multiple target genes for each miR in animals.

Isolation of BAC clones containing conserved genes from libraries of three distantly related moths: a useful resource for comparative genomics of Lepidoptera

Lepidoptera, butterflies and moths, is the second largest animal order and includes numerous agricultural pests. To facilitate comparative genomics in Lepidoptera, we isolated BAC clones containing conserved and putative single-copy genes from libraries of three pests, Heliothis virescens, Ostrinia nubilalis, and Plutella xylostella, harboring the haploid chromosome number, n = 31, which are not closely related with each other or with the silkworm, Bombyx mori, (n = 28), the sequenced model lepidopteran. A total of 108–184 clones representing 101–182 conserved genes were isolated for each species. For 79 genes, clones were isolated from more than two species, which will be useful as common markers for analysis using fluorescence in situ hybridization (FISH), as well as for comparison of genome sequence among multiple species. The PCR-based clone isolation method presented here is applicable to species which lack a sequenced genome but have a significant collection of cDNA or EST sequences.

Functional analysis of Ultrabithorax in the silkworm, Bombyx mori, using RNAi

The formation of abdominal appendages in insects is suppressed by the Hox genes Ultrabithorax (Ubx) and abdominal-A (abd-A), but mechanisms of the suppression can differ among species. As the function of Ubx and abd-A has been described in only a few species, more data from various insects are necessary to elucidate the evolutionary transition of regulation on abdominal appendages. We examined the function of Ubx in the silkworm Bombyx mori (Bm-Ubx) by embryonic RNA interference (RNAi). This is the first case in which functional analysis for Ubx is performed in lepidopteran insects. Larvae treated with Bm-Ubx dsRNA displayed an additional pair of thoracic leg-like protuberances in A1, whereas the other abdominal segments had no transformation. Our results suggest that Bm-Ubx is a suppressor of leg development in A1.

Reprobing multicolor FISH preparations in lepidopteran chromosome

Multicolor fluorescence in-situ hybridization (FISH) and subsequent reprobing of chromosome preparations increase the number of chromosomes and/or anchor loci on the chromosomes simultaneously identified. Reprobing techniques have been widely applied to chromosomes of vertebrates and plants. We have developed a novel reprobing protocol that utilizes multicolor FISH and bacterial artificial chromosome (BAC) probes to examine chromosome preparations in a model lepidopteran species, the silkworm, Bombyx mori. With standard two-color BAC-FISH, routinely used to map genes on B. mori chromosomes, we could localize only two probes on one preparation, whereas our new protocol combining five-color BAC-FISH and preparation reprobing enabled us to simultaneously map 10 probes, as demonstrated with the Bombyx Z chromosome. The improved BAC-FISH technique will facilitate karyotyping and synteny analysis in Lepidoptera.

A homeotic mutation influences the wing vibration patterns during mating in males of the silkworm moth Bombyx mori

An abnormality in the wing vibration pattern in males of the E(Nc) homeotic mutant of Bombyx mori was investigated. The wild-type (+/+) males show a switching of the rhythmic wing vibrations from a sequential pattern to an intermittent pattern during mating, whereas the E(Nc) mutants show a sequential pattern both before and during mating. Wing motions in +/+ males became small during mating, but those in +/E(Nc) males did not. Ablation of the head ganglia of +/+ and +/E(Nc) males during mating caused no change in the motor patterns of wing vibrations. Ablation or cooling of the posterior abdomen in the +/+ males during mating caused sequential wing vibrations, suggesting that the change in wing vibrations is induced by signals from the posterior abdomen. The pterothoracic ganglion in the +/E(Nc) males is separated into two ganglia, in contrast to the complete ganglionic fusion in the +/+ males. The neurons in the pterothoracic ganglion stained from abdominal nerve cords are homologus in +/+ and +/E(Nc) males, but many of these in +/E(Nc) males are elongated along the anteroposterior axis. These results suggest that the wing vibration pattern is restricted by genetic factors through reconstruction of the thoracic nervous system during metamorphosis.

The bithorax complex of Drosophila an exceptional Hox cluster

In his 1978 seminal paper, Ed Lewis described a series of mutations that affect the segmental identities of the segments forming the posterior two-thirds of the Drosophila body plan. In each class of mutations, particular segments developed like copies of a more-anterior segment. Genetic mapping of the different classes of mutations led to the discovery that their arrangement along the chromosome paralleled the body segments they affect along the anteroposterior axis of the fly. As all these mutations mapped to the same cytological location, he named this chromosomal locus after its founding mutation. Thus the first homeotic gene (Hox) cluster became known as the bithorax complex (BX-C). Even before the sequencing of the BX-C, the fact that these similar mutations grouped together in a cluster, lead Ed Lewis to propose that the homeotic genes arose through a gene duplication mechanism and that these clusters would be conserved through evolution. With the identification of the homeobox in the early 1980s, Lewis' first prediction was confirmed. The two cloned Drosophila homeotic genes, Antennapedia and Ultrabithorax, were indeed related genes. Using the homeobox as an entry point, homologous genes have since been cloned in many other species. Today, Hox clusters have been discovered in almost all metazoan phyla, confirming Lewis' second prediction. Remarkably, these homologous Hox genes are also arranged in clusters with their order within each cluster reflecting the anterior boundary of their domain of expression along the anterior-posterior axis of the animal. This correlation between the genomic organization and the activity along the anteroposterior body axis is known as the principle of "colinearity." The description of the BX-C inspired decades of developmental and evolutionary biology. And although this first Hox cluster led to the identification of many important features common to all Hox gene clusters, it now turns out that the fly Hox clusters are rather exceptional when compared with the Hox clusters of other animals. In this chapter, we will review the history and salient features of bithorax molecular genetics, in part, emphasizing its unique features relative to the other Hox clusters.

Analysis of the Tribolium homeotic complex: insights into mechanisms constraining insect Hox clusters

The remarkable conservation of Hox clusters is an accepted but little understood principle of biology. Some organizational constraints have been identified for vertebrate Hox clusters, but most of these are thought to be recent innovations that may not apply to other organisms. Ironically, many model organisms have disrupted Hox clusters and may not be well-suited for studies of structural constraints. In contrast, the red flour beetle, Tribolium castaneum, which has a long history in Hox gene research, is thought to have a more ancestral-type Hox cluster organization. Here, we demonstrate that the Tribolium homeotic complex (HOMC) is indeed intact, with the individual Hox genes in the expected colinear arrangement and transcribed from the same strand. There is no evidence that the cluster has been invaded by non-Hox protein-coding genes, although expressed sequence tag and genome tiling data suggest that noncoding transcripts are prevalent. Finally, our analysis of several mutations affecting the Tribolium HOMC suggests that intermingling of enhancer elements with neighboring transcription units may constrain the structure of at least one region of the Tribolium cluster. This work lays a foundation for future studies of the Tribolium HOMC that may provide insights into the reasons for Hox cluster conservation.

Fine mapping of E(kp)-1, a locus associated with silkworm (Bombyx mori) proleg development

The silkworm homeotic mutant E(kp) has a pair of rudimentary abdominal legs, called prolegs, in its A2 segment. This phenotype is caused by a single dominant mutation at the E(kp)-1 locus, which was previously mapped to chromosome 6. To explore the possible association of Hox genes with proleg development in the silkworm, a map-based cloning strategy was used to isolate the E(kp)-1 locus. Five E(kp)-1-linked simple sequence repeat markers on chromosome 6 were used to generate a low-resolution map with a total genetic distance of 39.5 cM. Four additional cleaved amplified polymorphic sequence markers were developed based on the initial map. The closest marker to E(kp)-1 was at a genetic distance of 2.7 cM. A high-resolution genetic map was constructed using nine BC1 segregating populations consisting of 2396 individuals. Recombination suppression was observed in the vicinity of E(kp)-1. Four molecular markers were tightly linked to E(kp)-1, and three were clustered with it. These markers were used to screen a BAC library. A single bacterial artificial chromosome (BAC) clone spanning the E(kp)-1 locus was identified, and E(kp)-1 was delimited to a region less than 220 kb long that included the Hox gene abdominal-A and a non-coding locus, iab-4. These results provide essential information for the isolation of this locus, which may shed light on the mechanism of proleg development in the silkworm and possibly in Lepidoptera.

Conserved synteny of genes between chromosome 15 of Bombyx mori and a chromosome of Manduca sexta shown by five-color BAC-FISH

The successful assignment of the existing genetic linkage groups (LGs) to individual chromosomes and the second-generation linkage map obtained by mapping a large number of bacterial artificial chromosome (BAC) contigs in the silkworm, Bombyx mori , together with public nucleotide sequence databases, offer a powerful tool for the study of synteny between karyotypes of B. mori and other lepidopteran species. Conserved synteny of genes between particular chromosomes can be identified by comparatively mapping orthologous genes of the corresponding linkage groups with the help of BAC-FISH (fluorescent in situ hybridization). This technique was established in B. mori for 2 differently labeled BAC probes simultaneously hybridized to pachytene bivalents. To achieve higher-throughput comparative mapping using BAC-FISH in Lepidoptera, we developed a protocol for five-color BAC-FISH, which allowed us to map simultaneously 6 different BAC probes to chromosome 15 in B. mori. We identified orthologs of 6 B. mori LG15 genes (RpP0, RpS8, eIF3, RpL7A, RpS23, and Hsc70) for the tobacco hornworm, Manduca sexta , and selected the ortholog-containing BAC clones from an M. sexta BAC library. All 6 M. sexta BAC clones hybridized to a single M. sexta bivalent in pachytene spermatocytes. Thus, we have confirmed the conserved synteny between the B. mori chromosome 15 and the corresponding M. sexta chromosome (hence provisionally termed chromosome 15).

The rise and fall of Hox gene clusters

Although all bilaterian animals have a related set of Hox genes, the genomic organization of this gene complement comes in different flavors. In some unrelated species, Hox genes are clustered; in others, they are not. This indicates that the bilaterian ancestor had a clustered Hox gene family and that, subsequently, this genomic organization was either maintained or lost. Remarkably, the tightest organization is found in vertebrates, raising the embarrassingly finalistic possibility that vertebrates have maintained best this ancestral configuration. Alternatively, could they have co-evolved with an increased ;organization' of the Hox clusters, possibly linked to their genomic amplification, which would be at odds with our current perception of evolutionary mechanisms? When discussing the why's and how's of Hox gene clustering, we need to account for three points: the mechanisms of cluster evolution; the underlying biological constraints; and the developmental modes of the animals under consideration. By integrating these parameters, general conclusions emerge that can help solve the aforementioned dilemma.

Development of a chordate anterior–posterior axis without classical retinoic acid signaling

Developmental signaling by retinoic acid (RA) is thought to be an innovation essential for the origin of the chordate body plan. The larvacean urochordate Oikopleura dioica maintains a chordate body plan throughout life, and yet its genome appears to lack genes for RA synthesis, degradation, and reception. This suggests the hypothesis that the RA-machinery was lost during larvacean evolution, and predicts that Oikopleura development has become independent of RA-signaling. This prediction raises the problem that the anterior-posterior organization of a chordate body plan can be developed without the classical morphogenetic role of RA. To address this problem, we performed pharmacological treatments and analyses of developmental molecular markers to investigate whether RA acts in anterior-posterior axial patterning in Oikopleura embryos. Results revealed that RA does not cause homeotic posteriorization in Oikopleura as it does in vertebrates and cephalochordates, and showed that a chordate can develop the phylotypic body plan in the absence of the classical morphogenetic role of RA. A comparison of Oikopleura and ascidian evidence suggests that the lack of RA-induced homeotic posteriorization is a shared derived feature of urochordates. We discuss possible relationships of altered roles of RA in urochordate development to genomic events, such as rupture of the Hox-cluster, in the context of a new understanding of chordate phylogeny.

Patterns of conservation and change in honey bee developmental genes

The current insect genome sequencing projects provide an opportunity to extend studies of the evolution of developmental genes and pathways in insects. In this paper we examine the conservation and divergence of genes and developmental processes between Drosophila and the honey bee; two holometabolous insects whose lineages separated approximately 300 million years ago, by comparing the presence or absence of 308 Drosophila developmental genes in the honey bee. Through examination of the presence or absence of genes involved in conserved pathways (cell signaling, axis formation, segmentation and homeobox transcription factors), we find that the vast majority of genes are conserved. Some genes involved in these processes are, however, missing in the honey bee. We have also examined the orthology of Drosophila genes involved in processes that differ between the honey bee and Drosophila. Many of these genes are preserved in the honey bee despite the process in which they act in Drosophila being different or absent in the honey bee. Many of the missing genes in both situations appear to have arisen recently in the Drosophila lineage, have single known functions in Drosophila, and act early in developmental pathways, while those that are preserved have pleiotropic functions. An evolutionary interpretation of these data is that either genes with multiple functions in a common ancestor are more likely to be preserved in both insect lineages, or genes that are preserved throughout evolution are more likely to co-opt additional functions.

Are Cirripedia hopeful monsters? Cytogenetic approach and evidence for a Hox gene cluster in the cirripede crustacean Sacculina carcini

The "hopeful monster" has haunted evolutionary thinking since Richard Goldschmidt coined the phrase in 1933. The phrase is directly related to genetic mechanisms in development and evolution. Cirripedes are peculiar crustaceans in that they all lack abdomens as adults. In a previous study aimed at describing the repertoire of Hox genes of the Cirripedia, we failed to isolate the abdominal-A gene in three species representative of all three cirripede orders. To address the question of whether the cirripede ancestor could have been a "hopeful monster" arising from a rearrangement of the Hox complex, we have performed a cytogenetic analysis of the Hox complex of the cirripede Sacculina carcini. We present here molecular and cytogenetic evidence for the grouping of the Hox genes on a single chromosome. This is the first direct evidence reported for the grouping of Hox genes on the same chromosome in a non-insect arthropod species.

A second-generation integrated map of the silkworm reveals synteny and conserved gene order between lepidopteran insects

A second-generation linkage map was constructed for the silkworm, Bombyx mori, focusing on mapping Bombyx sequences appearing in public nucleotide databases and bacterial artificial chromosome (BAC) contigs. A total of 874 BAC contigs containing 5067 clones (22% of the library) were constructed by PCR-based screening with sequence-tagged sites (STSs) derived from whole-genome shotgun (WGS) sequences. A total of 523 BAC contigs, including 342 independent genes registered in public databases and 85 expressed sequence tags (ESTs), were placed onto the linkage map. We found significant synteny and conserved gene order between B. mori and a nymphalid butterfly, Heliconius melpomene, in four linkage groups (LGs), strongly suggesting that using B. mori as a reference for comparative genomics in Lepidotera is highly feasible.

Unusual gene order and organization of the sea urchin hox cluster

While the highly consistent gene order and axial colinear patterns of expression seem to be a feature of vertebrate hox gene clusters, this pattern may be less well conserved across the rest of the bilaterians. We report the first deuterostome instance of an intact hox cluster with a unique gene order where the paralog groups are not expressed in a sequential manner. The finished sequence from BAC clones from the genome of the sea urchin, Strongylocentrotus purpuratus, reveals a gene order wherein the anterior genes (Hox1, Hox2 and Hox3) lie nearest the posterior genes in the cluster such that the most 3' gene is Hox5. (The gene order is 5'-Hox1, 2, 3, 11/13c, 11/13b, 11/13a, 9/10, 8, 7, 6, 5-3'.) The finished sequence result is corroborated by restriction mapping evidence and BAC-end scaffold analyses. Comparisons with a putative ancestral deuterostome Hox gene cluster suggest that the rearrangements leading to the sea urchin gene order were many and complex.

The genesis and evolution of homeobox gene clusters

Once called the 'Rosetta stone' of developmental biology, the homeobox continues to fascinate both evolutionary and developmental biologists. The birth of the homeotic, or Hox, gene cluster, and its subsequent evolution, has been crucial in mediating the major transitions in metazoan body plan. Comparative genomics studies indicate that the more recently discovered ParaHox and NK clusters were linked to the Hox cluster early in evolution, and that together they constituted a 'megacluster' of homeobox genes that conspicuously contributed to body-plan evolution.

Evidence for a Dispersed Hox Gene Cluster in the Platyhelminth Parasite Schistosoma mansoni

In most bilaterian organisms so far studied, Hox genes are organized in genomic clusters and determine development along the anteroposterior axis. It has been suggested that this clustering, together with spatial and temporal colinearity of gene expression, represents the ancestral condition. However, in organisms with derived modes of embryogenesis and lineage-dependent mechanisms for the determination of cell fate, temporal colinearity of expression can be lost and Hox cluster organization disrupted, as is the case for the ecdysozoans Drosophila melanogaster and Caenorhabditis elegans and the urochordates Ciona intestinalis and Oikopleura dioica. We sought to determine whether a lophotrochozoan, the platyhelminth parasite Schistosoma mansoni, possesses a conserved or disrupted Hox cluster. Using a polymerase chain reaction (PCR)-based strategy, we have cloned and characterized three novel S. mansoni genes encoding orthologues of Drosophila labial (SmHox1), deformed (SmHox4), and abdominal A (SmHox8), as well as the full-length coding sequence of the previously described Smox1, which we identify as an orthologue of fushi tarazu. Quantitative reverse transcriptase-PCR showed that the four genes were expressed at all life-cycle stages but that levels of expression were differentially regulated. Phylogenetic analysis and the conservation of "parapeptide" sequences C-terminal to the homeodomains of SmHox8 and Smox1 support the grouping of platyhelminths within the lophotrochozoan clade. However, Bacterial Artificial Chromosome (BAC) library screening followed by genome walking failed to reconstitute a cluster. The BAC clones containing Hox genes were sequenced, and in no case were other Hox genes found on the same clone. Moreover, the SmHox4 and SmHox8 genes contained single very large introns (>40 kbp) further indicating that the schistosome Hox cluster is highly extended. Localization of the Hox genes to chromosomes using fluorescence in situ hybridization showed that SmHox4 and SmHox8 are on the long arm of chromosome 4, whereas SmHox1 and Smox1 are on chromosome 3. In silico screening of the available genome sequences corroborated results of Southern blotting and BAC library screening that indicate that there are no paralogues of SmHox1, SmHox4, or SmHox8. The schistosome Hox cluster is therefore not duplicated, but is both dispersed and disintegrated in the genome.

Conservation of regulatory sequences and gene expression patterns in the disintegrating Drosophila Hox gene complex

Homeotic (Hox) genes are usually clustered and arranged in the same order as they are expressed along the anteroposterior body axis of metazoans. The mechanistic explanation for this colinearity has been elusive, and it may well be that a single and universal cause does not exist. The Hox-gene complex (HOM-C) has been rearranged differently in several Drosophila species, producing a striking diversity of Hox gene organizations. We investigated the genomic and functional consequences of the two HOM-C splits present in Drosophila buzzatii. Firstly, we sequenced two regions of the D. buzzatii genome, one containing the genes labial and abdominal A, and another one including proboscipedia, and compared their organization with that of D. melanogaster and D. pseudoobscura in order to map precisely the two splits. Then, a plethora of conserved noncoding sequences, which are putative enhancers, were identified around the three Hox genes closer to the splits. The position and order of these enhancers are conserved, with minor exceptions, between the three Drosophila species. Finally, we analyzed the expression patterns of the same three genes in embryos and imaginal discs of four Drosophila species with different Hox-gene organizations. The results show that their expression patterns are conserved despite the HOM-C splits. We conclude that, in Drosophila, Hox-gene clustering is not an absolute requirement for proper function. Rather, the organization of Hox genes is modular, and their clustering seems the result of phylogenetic inertia more than functional necessity.

The Bombyx mori karyotype and the assignment of linkage groups

Lepidopteran species have a relatively high number of small holocentric chromosomes (Bombyx mori, 2n = 56). Chromosome identification has long been hampered in this group by the high number and by the absence of suitable markers like centromere position and chromosome bands. In this study, we carried out fluorescence in situ hybridization (FISH) on meiotic chromosome complements using genetically mapped B. mori bacterial artificial chromosomes (BACs) as probes. The combination of two to four either green or red fluorescence-labeled probes per chromosome allowed us to recognize unequivocally each of the 28 bivalents of the B. mori karyotype by its labeling pattern. Each chromosome was assigned one of the already established genetic linkage groups and the correct orientation in the chromosome was defined. This facilitates physical mapping of any other sequence and bears relevance for the ongoing B. mori genome projects. Two-color BAC-FISH karyotyping overcomes the problem of chromosome recognition in organisms where conventional banding techniques are not available.