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

Notch Signaling in Insect Development: A Simple Pathway with Diverse Functions

Notch signaling is an evolutionarily conserved pathway which functions between adjacent cells to establish their distinct identities. Despite operating in a simple mechanism, Notch signaling plays remarkably diverse roles in development to regulate cell fate determination, organ growth and tissue patterning. While initially discovered and characterized in the model insect Drosophila melanogaster, recent studies across various insect species have revealed the broad involvement of Notch signaling in shaping insect tissues. This review focuses on providing a comprehensive picture regarding the roles of the Notch pathway in insect development. The roles of Notch in the formation and patterning of the insect embryo, wing, leg, ovary and several specific structures, as well as in physiological responses, are summarized. These results are discussed within the developmental context, aiming to deepen our understanding of the diversified functions of the Notch signaling pathway in different insect species.

Evolution and genomic organization of the insect sHSP gene cluster and coordinate regulation in phenotypic plasticity

Background

Conserved syntenic gene complexes are rare in Arthropods and likely only retained due to functional constraint. Numerous sHSPs have been identified in the genomes of insects, some of which are located clustered in close proximity. Previous phylogenetic analyses of these clustered sHSP have been limited to a small number of holometabolous insect species and have not determined the pattern of evolution of the clustered sHSP genes (sHSP-C) in insect or Arthropod lineages.

Results

Using eight genomes from representative insect orders and three non-insect arthropod genomes we have identified that a syntenic cluster of sHSPs (sHSP-C) is a hallmark of most Arthropod genomes. Using 11 genomes from Hymenopteran species our phylogenetic analyses have refined the evolution of the sHSP-C in Hymenoptera and found that the sHSP-C is order-specific with evidence of birth-and-death evolution in the hymenopteran lineage. Finally we have shown that the honeybee sHSP-C is co-ordinately expressed and is marked by genomic features, including H3K27me3 histone marks consistent with coordinate regulation, during honeybee ovary activation.

Conclusions

The syntenic sHSP-C is present in most insect genomes, and its conserved coordinate expression and regulation implies that it is an integral genomic component of environmental response in arthropods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01885-8.

Genome Architecture Facilitates Phenotypic Plasticity in the Honeybee (Apis mellifera)

Phenotypic plasticity, the ability of an organism to alter its phenotype in response to an environmental cue, facilitates rapid adaptation to changing environments. Plastic changes in morphology and behavior are underpinned by widespread gene expression changes. However, it is unknown if, or how, genomes are structured to ensure these robust responses. Here, we use repression of honeybee worker ovaries as a model of plasticity. We show that the honeybee genome is structured with respect to plasticity; genes that respond to an environmental trigger are colocated in the honeybee genome in a series of gene clusters, many of which have been assembled in the last 80 My during the evolution of the Apidae. These clusters are marked by histone modifications that prefigure the gene expression changes that occur as the ovary activates, suggesting that these genomic regions are poised to respond plastically. That the linear sequence of the honeybee genome is organized to coordinate widespread gene expression changes in response to environmental influences and that the chromatin organization in these regions is prefigured to respond to these influences is perhaps unexpected and has implications for other examples of plasticity in physiology, evolution, and human disease.

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.

Notch signalling mediates reproductive constraint in the adult worker honeybee

The hallmark of eusociality is the reproductive division of labour, in which one female caste reproduces, while reproduction is constrained in the subordinate caste. In adult worker honeybees (Apis mellifera) reproductive constraint is conditional: in the absence of the queen and brood, adult worker honeybees activate their ovaries and lay haploid male eggs. Here, we demonstrate that chemical inhibition of Notch signalling can overcome the repressive effect of queen pheromone and promote ovary activity in adult worker honeybees. We show that Notch signalling acts on the earliest stages of oogenesis and that the removal of the queen corresponds with a loss of Notch protein in the germarium. We conclude that the ancient and pleiotropic Notch signalling pathway has been co-opted into constraining reproduction in worker honeybees and we provide the first molecular mechanism directly linking ovary activity in adult worker bees with the presence of the queen.

In honeybees, pheromones produced by the queen inhibit reproduction by workers and enforce a eusocial division of labour. Here, Duncan, Hyink and Dearden show that this inhibition is mediated by the Notch signalling pathway in the workers' ovaries.

Origin and evolution of the enhancer of split complex

BACKGROUND

The Enhancer of split complex is an unusual gene complex found in Arthropod genomes. Where known this complex of genes is often regulated by Notch cell signalling and is critically important for neurogenesis. The Enhancer of split complex is made up of two different classes of genes, basic helix-loop-helix-orange domain transcription factors and bearded class genes. The association of these genes has been detected in the genomes of insects and crustaceans.

RESULTS

Tracing the evolution of the Enhancer of split complex in recently sequenced Arthropod genomes indicates that enhancer of split basic helix-loop-helix orange domain genes arose before the common ancestor of insects and Crustacea, and before the formation of the complex. Throughout insect and crustacean evolution, a four-gene cluster has been present with lineage specific gene losses and duplications. The complex can be found in the vast majority of genomes, but appears to be missing from the genomes of chalcid wasps, raising questions as to how they carry out neurogenesis in the absence of these crucial genes.

CONCLUSIONS

The enhancer of split complex arose in the common ancestor of Crustacea and insects, probably through the linkage of a basic helix-loop-helix orange domain gene and a bearded class gene. The complex has been maintained, with variations, throughout insect and crustacean evolution indicating some function of the complex, such as coordinate regulation, may maintain its structure through evolutionary time.

Identification of two juvenile hormone inducible transcription factors from the silkworm, Bombyx mori

Juvenile hormone (JH) regulates many physiological processes in insects. However, the signal cascades in which JH is active have not yet been fully elucidated, particularly in comparison to another major hormone ecdysteroid. Here we identified two JH inducible transcription factors as candidate components of JH signaling pathways in the silkworm, Bombyx mori. DNA microarray analysis showed that expression of two transcription factor genes, E75 and Enhancer of split mβ (E(spl)mβ), was induced by juvenile hormone I (JH I) in NIAS-Bm-aff3 cells. Real time RT-PCR analysis confirmed that expression of four E75 isoforms (E75A, E75B, E75C and E75D) and E(spl)mβ was 3-8 times greater after JH I addition. Addition of the protein synthesis inhibitor cycloheximide did not suppress JH-induced expression of the genes, indicating that they were directly induced by JH. JH-induced expression of E75 and E(spl)mβ was also observed in four other B. mori cell lines and in larval hemocytes of final instar larvae. Notably, E75A expression was induced very strongly in larval hemocytes by topical application of the JH analog fenoxycarb; the level of induced expression was comparable to that produced by feeding larvae with 20-hydroxyecdysone. These results suggest that E75 and E(spl)mβ are general and direct target genes of JH and that the transcription factors encoded by these genes play important roles in JH signaling.

The genomes of two key bumblebee species with primitive eusocial organization

Background

The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats.

Results

We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits.

Conclusions

These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0623-3) contains supplementary material, which is available to authorized users.

Multispecies Analysis of Expression Pattern Diversification in the Recently Expanded Insect Ly6 Gene Family

Gene families often consist of members with diverse expression domains reflecting their functions in a wide variety of tissues. However, how the expression of individual members, and thus their tissue-specific functions, diversified during the course of gene family expansion is not well understood. In this study, we approached this question through the analysis of the duplication history and transcriptional evolution of a rapidly expanding subfamily of insect Ly6 genes. We analyzed different insect genomes and identified seven Ly6 genes that have originated from a single ancestor through sequential duplication within the higher Diptera. We then determined how the original embryonic expression pattern of the founding gene diversified by characterizing its tissue-specific expression in the beetle Tribolium castaneum, the butterfly Bicyclus anynana, and the mosquito Anopheles stephensi and those of its duplicates in three higher dipteran species, representing various stages of the duplication history (Megaselia abdita, Ceratitis capitata, and Drosophila melanogaster). Our results revealed that frequent neofunctionalization episodes contributed to the increased expression breadth of this subfamily and that these events occurred after duplication and speciation events at comparable frequencies. In addition, at each duplication node, we consistently found asymmetric expression divergence. One paralog inherited most of the tissue-specificities of the founder gene, whereas the other paralog evolved drastically reduced expression domains. Our approach attests to the power of combining a well-established duplication history with a comprehensive coverage of representative species in acquiring unequivocal information about the dynamics of gene expression evolution in gene families.

Evidence for the plasticity of arthropod signal transduction pathways

Metazoans are known to contain a limited, yet highly conserved, set of signal transduction pathways that instruct early developmental patterning mechanisms. Genomic surveys that have compared gene conservation in signal transduction pathways between various insects and Drosophila support the conclusion that these pathways are conserved in evolution. However, the degree to which individual components of signal transduction pathways vary among more divergent arthropods is not known. Here, we report our results of a survey of the genome of the two-spotted spider mite Tetranychus urticae, using a set of 294 Drosophila orthologs of genes that function in signal transduction. We find a third of all genes surveyed absent from the spider mite genome. We also identify several novel duplications that have not been previously reported for a chelicerate. In comparison with previous insect surveys, Tetranychus contains a decrease in overall gene conservation, as well as an unusual ratio of ligands to receptors and other modifiers. These findings suggest that gene loss and duplication among components of signal transduction pathways are common among arthropods and suggest that signal transduction pathways in arthropods are more evolutionarily labile than previously hypothesized.

History of a prolific family: the Hes/Hey-related genes of the annelid Platynereis

BACKGROUND

The Hes superfamily or Hes/Hey-related genes encompass a variety of metazoan-specific bHLH genes, with somewhat fuzzy phylogenetic relationships. Hes superfamily members are involved in a variety of major developmental mechanisms in metazoans, notably in neurogenesis and segmentation processes, in which they often act as direct effector genes of the Notch signaling pathway.

RESULTS

We have investigated the molecular and functional evolution of the Hes superfamily in metazoans using the lophotrochozoan Platynereis dumerilii as model. Our phylogenetic analyses of more than 200 Metazoan Hes/Hey-related genes revealed the presence of five families, three of them (Hes, Hey and Helt) being pan-metazoan. Those families were likely composed of a unique representative in the last common metazoan ancestor. The evolution of the Hes family was shaped by many independent lineage specific tandem duplication events. The expression patterns of 13 of the 15 Hes/Hey-related genes in Platynereis indicate a broad functional diversification. Nevertheless, a majority of these genes are involved in two crucial developmental processes in annelids: neurogenesis and segmentation, resembling functions highlighted in other animal models.

CONCLUSIONS

Combining phylogenetic and expression data, our study suggests an unusual evolutionary history for the Hes superfamily. An ancestral multifunctional annelid Hes gene may have undergone multiples rounds of duplication-degeneration-complementation processes in the lineage leading to Platynereis, each gene copies ensuring their maintenance in the genome by subfunctionalisation. Similar but independent waves of duplications are at the origin of the multiplicity of Hes genes in other metazoan lineages.

The Drosophila enhancer of split gene complex: architecture and coordinate regulation by notch, cohesin, and polycomb group proteins

The cohesin protein complex functionally interacts with Polycomb group (PcG) silencing proteins to control expression of several key developmental genes, such as the Drosophila Enhancer of split gene complex [E(spl)-C]. The E(spl)-C contains 12 genes that inhibit neural development. In a cell line derived from the central nervous system, cohesin and the PRC1 PcG protein complex bind and repress E (spl)-C transcription, but the repression mechanisms are unknown. The genes in the E(spl)-C are directly activated by the Notch receptor. Here we show that depletion of cohesin or PRC1 increases binding of the Notch intracellular fragment to genes in the E(spl)-C, correlating with increased transcription. The increased transcription likely reflects both direct effects of cohesin and PRC1 on RNA polymerase activity at the E(spl)-C, and increased expression of Notch ligands. By chromosome conformation capture we find that the E(spl)-C is organized into a self-interactive architectural domain that is co-extensive with the region that binds cohesin and PcG complexes. The self-interactive architecture is formed independently of cohesin or PcG proteins. We posit that the E(spl)-C architecture dictates where cohesin and PcG complexes bind and act when they are recruited by as yet unidentified factors, thereby controlling the E(spl)-C as a coordinated domain.

Classification and evolutionary analysis of the basic helix-loop-helix gene family in the green anole lizard, Anolis carolinensis

Helix-loop-helix (bHLH) proteins play essential regulatory roles in a variety of biological processes. These highly conserved proteins form a large transcription factor superfamily, and are commonly identified in large numbers within animal, plant, and fungal genomes. The bHLH domain has been well studied in many animal species, but has not yet been characterized in non-avian reptiles. In this study, we identified 102 putative bHLH genes in the genome of the green anole lizard, Anolis carolinensis. Based on phylogenetic analysis, these genes were classified into 43 families, with 43, 24, 16, 3, 10, and 3 members assigned into groups A, B, C, D, E, and F, respectively, and 3 members categorized as "orphans". Within-group evolutionary relationships inferred from the phylogenetic analysis were consistent with highly conserved patterns observed for introns and additional domains. Results from phylogenetic analysis of the H/E(spl) family suggest that genome and tandem gene duplications have contributed to this family's expansion. Our classification and evolutionary analysis has provided insights into the evolutionary diversification of animal bHLH genes, and should aid future studies on bHLH protein regulation of key growth and developmental processes.

The pea aphid (Acyrthosiphon pisum) genome encodes two divergent early developmental programs

The pea aphid (Acyrthosiphon pisum) can reproduce either sexually or asexually (parthenogenetically), giving rise, in each case, to almost identical adults. These two modes of reproduction are accompanied by differences in ovarian morphology and the developmental environment of the offspring, with sexual forms producing eggs that are laid, whereas asexual development occurs within the mother. Here we examine the effect each mode of reproduction has on the expression of key maternal and axis patterning genes; orthodenticle (otd), hunchback (hb), caudal (cad) and nanos (nos). We show that three of these genes (Ap-hb, Ap-otd and Ap-cad) are expressed differently between the sexually and asexually produced oocytes and embryos of the pea aphid. We also show, using immunohistochemistry and cytoskeletal inhibitors, that Ap-hb RNA is localized differently between sexually and asexually produced oocytes, and that this is likely due to differences in the 3' untranslated regions of the RNA. Furthermore, Ap-hb and Ap-otd have extensive expression domains in early sexually produced embryos, but are not expressed at equivalent stages in asexually produced embryos. These differences in expression likely correspond with substantial changes in the gene regulatory networks controlling early development in the pea aphid. These data imply that in the evolution of parthenogenesis a new program has evolved to control the development of asexually produced embryos, whilst retaining the existing, sexual, developmental program. The patterns of modification of these developmental processes mirror the changes that we see in developmental processes between species, in that early acting pathways in development are less constrained, and evolve faster, than later ones. We suggest that the evolution of the novel asexual development pathway in aphids is not a simple modification of an ancestral system, but the evolution of two very different developmental mechanisms occurring within a single species.

Canonical terminal patterning is an evolutionary novelty

Patterning of the terminal regions of the Drosophila embryo is achieved by an exquisitely regulated signal that passes between the follicle cells of the ovary, and the developing embryo. This pathway, however, is missing or modified in other insects. Here we trace the evolution of this pathway by examining the origins and expression of its components. The three core components of this pathway: trunk, torso and torso-like have different evolutionary histories and have been assembled step-wise to form the canonical terminal patterning pathway of Drosophila and Tribolium. Trunk, torso and a gene unrelated to terminal patterning, prothoraciotrophic hormone (PTTH), show an intimately linked evolutionary history, with every holometabolous insect, except the honeybee, possessing both PTTH and torso genes. Trunk is more restricted in its phylogenetic distribution, present only in the Diptera and Tribolium and, surprisingly, in the chelicerate Ixodes scapularis, raising the possibility that trunk and torso evolved earlier than previously thought. In Drosophila torso-like restricts the activation of the terminal patterning pathway to the poles of the embryo. Torso-like evolved in the pan-crustacean lineage, but based on expression of components of the canonical terminal patterning system in the hemimetabolous insect Acyrthosiphon pisum and the holometabolous insect Apis mellifera, we find that the canonical terminal-patterning system is not active in these insects. We therefore propose that the ancestral function of torso-like is unrelated to terminal patterning and that torso-like has become co-opted into terminal patterning in the lineage leading to Coleoptera and Diptera. We also show that this co-option has not resulted in changes to the molecular function of this protein. Torso-like from the pea aphid, honeybee and Drosophila, despite being expressed in different patterns, are functionally equivalent. We propose that co-option of torso-like into restricting the activity of trunk and torso facilitated the final step in the evolution of this pathway; the capture of transcriptional control of target genes such as tailless and huckebein by this complex and novel patterning pathway.

The molecular clockwork of the fire ant Solenopsis invicta

The circadian clock is a core molecular mechanism that allows organisms to anticipate daily environmental changes and adapt the timing of behaviors to maximize efficiency. In social insects, the ability to maintain the appropriate temporal order is thought to improve colony efficiency and fitness. We used the newly sequenced fire ant (Solenopsis invicta) genome to characterize the first ant circadian clock. Our results reveal that the fire ant clock is similar to the clock of the honeybee, a social insect with an independent evolutionary origin of sociality. Gene trees for the eight core clock genes, period, cycle, clock, cryptochrome-m, timeout, vrille, par domain protein 1 & clockwork orange, show ant species grouping closely with honeybees and Nasonia wasps as an outgroup to the social Hymenoptera. Expression patterns for these genes suggest that the ant clock functions similar to the honeybee clock, with period and cry-m mRNA levels increasing during the night and cycle and clockwork orange mRNAs cycling approximately anti-phase to period. Gene models for five of these genes also parallel honeybee models. In particular, the single ant cryptochrome is an ortholog of the mammalian-type (cry-m), rather than Drosophila-like protein (cry-d). Additionally, we find a conserved VPIFAL C-tail region in clockwork orange shared by insects but absent in vertebrates. Overall, our characterization of the ant clock demonstrates that two social insect lineages, ants and bees, share a similar, mammalian-like circadian clock. This study represents the first characterization of clock genes in an ant and is a key step towards understanding socially-regulated plasticity in circadian rhythms by facilitating comparative studies on the organization of circadian clockwork.

Essential roles of Da transactivation domains in neurogenesis and in E(spl)-mediated repression

E proteins are a special class of basic helix-loop-helix (bHLH) proteins that heterodimerize with many bHLH activators to regulate developmental decisions, such as myogenesis and neurogenesis. Daughterless (Da) is the sole E protein in Drosophila and is ubiquitously expressed. We have characterized two transcription activation domains (TADs) in Da, called activation domain 1 (AD1) and loop-helix (LH), and have evaluated their roles in promoting peripheral neurogenesis. In this context, Da heterodimerizes with proneural proteins, such as Scute (Sc), which is dynamically expressed and also contributes a TAD. We found that either one of the Da TADs in the Da/Sc complex is sufficient to promote neurogenesis, whereas the Sc TAD is incapable of doing so. Besides its transcriptional activation role, the Da AD1 domain serves as an interaction platform for E(spl) proteins, bHLH-Orange family repressors which antagonize Da/Sc function. We show that the E(spl) Orange domain is needed for this interaction and strongly contributes to the antiproneural activity of E(spl) proteins. We present a mechanistic model on the interplay of these bHLH factors in the context of neural fate assignment.

Comparative and evolutionary analysis of the HES/HEY gene family reveal exon/intron loss and teleost specific duplication events

Background

HES/HEY genes encode a family of basic helix-loop-helix (bHLH) transcription factors with both bHLH and Orange domain. HES/HEY proteins are direct targets of the Notch signaling pathway and play an essential role in developmental decisions, such as the developments of nervous system, somitogenesis, blood vessel and heart. Despite their important functions, the origin and evolution of this HES/HEY gene family has yet to be elucidated.

Methods and Findings

In this study, we identified genes of the HES/HEY family in representative species and performed evolutionary analysis to elucidate their origin and evolutionary process. Our results showed that the HES/HEY genes only existed in metazoans and may originate from the common ancestor of metazoans. We identified HES/HEY genes in more than 10 species representing the main lineages. Combining the bHLH and Orange domain sequences, we constructed the phylogenetic trees by different methods (Bayesian, ML, NJ and ME) and classified the HES/HEY gene family into four groups. Our results indicated that this gene family had undergone three expansions, which were along with the origins of Eumetazoa, vertebrate, and teleost. Gene structure analysis revealed that the HES/HEY genes were involved in exon and/or intron loss in different species lineages. Genes of this family were duplicated in bony fishes and doubled than other vertebrates. Furthermore, we studied the teleost-specific duplications in zebrafish and investigated the expression pattern of duplicated genes in different tissues by RT-PCR. Finally, we proposed a model to show the evolution of this gene family with processes of expansion, exon/intron loss, and motif loss.

Conclusions

Our study revealed the evolution of HES/HEY gene family, the expression and function divergence of duplicated genes, which also provide clues for the research of Notch function in development. This study shows a model of gene family analysis with gene structure evolution and duplication.

Ancestral and conserved cis-regulatory architectures in developmental control genes

Among developmental control genes, transcription factor-target gene "linkages"--the direct connections between target genes and the factors that control their patterns of expression--can show remarkable evolutionary stability. However, the specific binding sites that mediate and define these regulatory connections are themselves often subject to rapid turnover. Here we describe several instances in which particular transcription factor binding motif combinations have evidently been conserved upstream of orthologous target genes for extraordinarily long evolutionary periods. This occurs against a backdrop in which other binding sites for the same factors are coming and going rapidly. Our examples include a particular Dpp Silencer Element upstream of insect brinker genes, in combination with a novel motif we refer to as the Downstream Element; combinations of a Suppressor of Hairless Paired Site (SPS) and a specific proneural protein binding site associated with arthropod Notch pathway target genes; and a three-motif combination, also including an SPS, upstream of deuterostome Hes repressor genes, which are also Notch targets. We propose that these stable motif architectures have been conserved intact from a deep ancestor, in part because they mediate a special mode of regulation that cannot be supplied by the other, unstable motif instances.

The Enhancer of split complex arose prior to the diversification of schizophoran flies and is strongly conserved between Drosophila and stalk-eyed flies (Diopsidae)

Background

In Drosophila, the Enhancer of split complex (E(spl)-C) comprises 11 bHLH and Bearded genes that function during Notch signaling to repress proneural identity in the developing peripheral nervous system. Comparison with other insects indicates that the basal state for Diptera is a single bHLH and Bearded homolog and that the expansion of the gene complex occurred in the lineage leading to Drosophila. However, comparative genomic data from other fly species that would elucidate the origin and sequence of gene duplication for the complex is lacking. Therefore, in order to examine the evolutionary history of the complex within Diptera, we reconstructed, using several fosmid clones, the entire E(spl)-complex in the stalk-eyed fly, Teleopsis dalmanni and collected additional homologs of E(spl)-C genes from searches of dipteran EST databases and the Glossina morsitans genome assembly.

Results

Comparison of the Teleopsis E(spl)-C gene organization with Drosophila indicates complete conservation in gene number and orientation between the species except that T. dalmanni contains a duplicated copy of E(spl)m5 that is not present in Drosophila. Phylogenetic analysis of E(spl)-complex bHLH and Bearded genes for several dipteran species clearly demonstrates that all members of the complex were present prior to the diversification of schizophoran flies. Comparison of upstream regulatory elements and 3' UTR domains between the species also reveals strong conservation for many of the genes and identifies several novel characteristics of E(spl)-C regulatory evolution including the discovery of a previously unidentified, highly conserved SPS+A domain between E(spl)mγ and E(spl)mβ.

Conclusion

Identifying the phylogenetic origin of E(spl)-C genes and their associated regulatory DNA is essential to understanding the functional significance of this well-studied gene complex. Results from this study provide numerous insights into the evolutionary history of the complex and will help refine the focus of studies examining the adaptive consequences of this gene expansion.