Arthropod Genomes Research Articles

Mechanisms and Dynamics of Orphan Gene Emergence in Insect Genomes

Orphan genes are defined as genes that lack detectable similarity to genes in other species and therefore no clear signals of common descent (i.e., homology) can be inferred. Orphans are an enigmatic portion of the genome because their origin and function are mostly unknown and they typically make up 10% to 30% of all genes in a genome. Several case studies demonstrated that orphans can contribute to lineage-specific adaptation. Here, we study orphan genes by comparing 30 arthropod genomes, focusing in particular on seven recently sequenced ant genomes. This setup allows analyzing a major metazoan taxon and a comparison between social Hymenoptera (ants and bees) and nonsocial Diptera (flies and mosquitoes). First, we find that recently split lineages undergo accelerated genomic reorganization, including the rapid gain of many orphan genes. Second, between the two insect orders Hymenoptera and Diptera, orphan genes are more abundant and emerge more rapidly in Hymenoptera, in particular, in leaf-cutter ants. With respect to intragenomic localization, we find that ant orphan genes show little clustering, which suggests that orphan genes in ants are scattered uniformly over the genome and between nonorphan genes. Finally, our results indicate that the genetic mechanisms creating orphan genes—such as gene duplication, frame-shift fixation, creation of overlapping genes, horizontal gene transfer, and exaptation of transposable elements—act at different rates in insects, primates, and plants. In Formicidae, the majority of orphan genes has their origin in intergenic regions, pointing to a high rate of de novo gene formation or generalized gene loss, and support a recently proposed dynamic model of frequent gene birth and death.

Variation and constraint inHoxgene evolution

Despite enormous body plan variation, genes regulating embryonic development are highly conserved. Here, we probe the mechanisms that predispose ancient regulatory genes to reutilization and diversification rather than evolutionary loss. The Hox gene fushi tarazu (ftz) arose as a homeotic gene but functions as a pair-rule segmentation gene in Drosophila. ftz shows extensive variation in expression and protein coding regions but has managed to elude loss from arthropod genomes. We asked what properties prevent this loss by testing the importance of different protein motifs and partners in the developing CNS, where ftz expression is conserved. Drosophila Ftz proteins with mutated protein motifs were expressed under the control of a neurogenic-specific ftz cis-regulatory element (CRE) in a ftz mutant background rescued for segmentation defects. Ftz CNS function did not require the variable motifs that mediate differential cofactor interactions involved in homeosis or segmentation, which vary in arthropods. Rather, CNS function did require the shared DNA-binding homeodomain, which plays less of a role in Ftz segmentation activity. The Antennapedia homeodomain substituted for Ftz homeodomain function in the Drosophila CNS, but full-length Antennapedia did not rescue CNS defects. These results suggest that a core CNS function retains ftz in arthropod genomes. Acquisition of a neurogenic CRE led to ftz expression in unique CNS cells, differentiating its role from neighboring Hox genes, rendering it nonredundant. The inherent flexibility of modular CREs and protein domains allows for stepwise acquisition of new functions, explaining broad retention of regulatory genes during animal evolution.

PhyloTreePruner: A Phylogenetic Tree-Based Approach for selection of Orthologous sequences for phylogenomics

Molecular phylogenetics relies on accurate identification of orthologous sequences among the taxa of interest. Most orthology inference programs available for use in phylogenomics rely on small sets of pre-defined orthologs from model organisms or phenetic approaches such as all-versus-all sequence comparisons followed by Markov graph-based clustering. Such approaches have high sensitivity but may erroneously include paralogous sequences. We developed PhyloTreePruner, a software utility that uses a phylogenetic approach to refine orthology inferences made using phenetic methods. PhyloTreePruner checks single-gene trees for evidence of paralogy and generates a new alignment for each group containing only sequences inferred to be orthologs. Importantly, PhyloTreePruner takes into account support values on the tree and avoids unnecessarily deleting sequences in cases where a weakly supported tree topology incorrectly indicates paralogy. A test of PhyloTreePruner on a dataset generated from 11 completely sequenced arthropod genomes identified 2,027 orthologous groups sampled for all taxa. Phylogenetic analysis of the concatenated supermatrix yielded a generally well-supported topology that was consistent with the current understanding of arthropod phylogeny. PhyloTreePruner is freely available from http://sourceforge.net/projects/phylotreepruner/.

Identification and classification of detoxification enzymes from Culex quinquefasciatus (Diptera: Culicidae).

Molecular characterization of the insecticide resistance has become a hot research topic ever since the first disease transmitting arthropod (Anopheles gambiae) genome sequence has unveiled in 2002. A recent publication of the Culex quinquefasciatus genome sequence has opened up new opportunities for molecular and comparative genomic analysis of multiple mosquito genomes to characterize the insecticide resistance. Here, we utilized a whole genome sequence of Cx. quinquefasciatus to identify putatively active members of the detoxification supergene families, namely cytochrome P450s (P450s), glutathione-S-transferases (GSTs), and choline/carboxylesterases (CCEs). The Culex genome analysis revealed 166 P450s, 40 GSTs, and 62 CCEs. Further, the comparative genomic analysis shows that these numbers are considerably higher than the other dipteran mosquitoes. These observed speciesspecific expansions of the detoxification super gene family members endorse the popular understanding of the involvement of these gene families in protecting the organism against multitudinous classes of toxic substances during its complex (aquatic and terrestrial) life cycle. Thus, the generated data set may provide an initial point to start with to characterize the insecticide resistance at a molecular level which could then lead the development of an easy to use molecular marker to monitor the incipient insecticide resistance in field environs.

Quantitative PCR-based genome size estimation of the astigmatid mites Sarcoptes scabiei, Psoroptes ovis and Dermatophagoides pteronyssinus

BackgroundThe lack of genomic data available for mites limits our understanding of their biology. Evolving high-throughput sequencing technologies promise to deliver rapid advances in this area, however, estimates of genome size are initially required to ensure sufficient coverage.MethodsQuantitative real-time PCR was used to estimate the genome sizes of the burrowing ectoparasitic mite Sarcoptes scabiei, the non-burrowing ectoparasitic mite Psoroptes ovis, and the free-living house dust mite Dermatophagoides pteronyssinus. Additionally, the chromosome number of S. scabiei was determined by chromosomal spreads of embryonic cells derived from single eggs.ResultsS. scabiei cells were shown to contain 17 or 18 small (< 2 μM) chromosomes, suggesting an XO sex-determination mechanism. The average estimated genome sizes of S. scabiei and P. ovis were 96 (± 7) Mb and 86 (± 2) Mb respectively, among the smallest arthropod genomes reported to date. The D. pteronyssinus genome was estimated to be larger than its parasitic counterparts, at 151 Mb in female mites and 218 Mb in male mites.ConclusionsThis data provides a starting point for understanding the genetic organisation and evolution of these astigmatid mites, informing future sequencing projects. A comparitive genomic approach including these three closely related mites is likely to reveal key insights on mite biology, parasitic adaptations and immune evasion.

A gene-based SNP resource and linkage map for the copepod Tigriopus californicus.

BackgroundAs yet, few genomic resources have been developed in crustaceans. This lack is particularly evident in Copepoda, given the extraordinary numerical abundance, and taxonomic and ecological diversity of this group. Tigriopus californicus is ideally suited to serve as a genetic model copepod and has been the subject of extensive work in environmental stress and reproductive isolation. Accordingly, we set out to develop a broadly-useful panel of genetic markers and to construct a linkage map dense enough for quantitative trait locus detection in an interval mapping framework for T. californicus--a first for copepods.ResultsOne hundred and ninety Single Nucleotide Polymorphisms (SNPs) were used to genotype our mapping population of 250 F2 larvae. We were able to construct a linkage map with an average intermarker distance of 1.8 cM, and a maximum intermarker distance of 10.3 cM. All markers were assembled into linkage groups, and the 12 linkage groups corresponded to the 12 known chromosomes of T. californicus. We estimate a total genome size of 401.0 cM, and a total coverage of 73.7%. Seventy five percent of the mapped markers were detected in 9 additional populations of T. californicus. Of available model arthropod genomes, we were able to show more colocalized pairs of homologues between T. californicus and the honeybee Apis mellifera, than expected by chance, suggesting preserved macrosynteny between Hymenoptera and Copepoda.ConclusionsOur study provides an abundance of linked markers spanning all chromosomes. Many of these markers are also found in multiple populations of T. californicus, and in two other species in the genus. The genomic resource we have developed will enable mapping throughout the geographical range of this species and in closely related species. This linkage map will facilitate genome sequencing, mapping and assembly in an ecologically and taxonomically interesting group for which genomic resources are currently under development.

The genome of Tetranychus urticae reveals herbivorous pest adaptations

The spider mite Tetranychus urticae is a cosmopolitan agricultural pest with an extensive host plant range and an extreme record of pesticide resistance. Here we present the completely sequenced and annotated spider mite genome, representing the first complete chelicerate genome. At 90 megabases T. urticae has the smallest sequenced arthropod genome. Compared with other arthropods, the spider mite genome shows unique changes in the hormonal environment and organization of the Hox complex, and also reveals evolutionary innovation of silk production. We find strong signatures of polyphagy and detoxification in gene families associated with feeding on different hosts and in new gene families acquired by lateral gene transfer. Deep transcriptome analysis of mites feeding on different plants shows how this pest responds to a changing host environment. The T. urticae genome thus offers new insights into arthropod evolution and plant–herbivore interactions, and provides unique opportunities for developing novel plant protection strategies.

Fossil Rhabdoviral Sequences Integrated into Arthropod Genomes: Ontogeny, Evolution, and Potential Functionality

Retroelements represent a considerable fraction of many eukaryotic genomes and are considered major drives for adaptive genetic innovations. Recent discoveries showed that despite not normally using DNA intermediates like retroviruses do, Mononegaviruses (i.e., viruses with nonsegmented, negative-sense RNA genomes) can integrate gene fragments into the genomes of their hosts. This was shown for Bornaviridae and Filoviridae, the sequences of which have been found integrated into the germ line cells of many vertebrate hosts. Here, we show that Rhabdoviridae sequences, the major Mononegavirales family, have integrated only into the genomes of arthropod species. We identified 185 integrated rhabdoviral elements (IREs) coding for nucleoproteins, glycoproteins, or RNA-dependent RNA polymerases; they were mostly found in the genomes of the mosquito Aedes aegypti and the blacklegged tick Ixodes scapularis. Phylogenetic analyses showed that most IREs in A. aegypti derived from multiple independent integration events. Since RNA viruses are submitted to much higher substitution rates as compared with their hosts, IREs thus represent fossil traces of the diversity of extinct Rhabdoviruses. Furthermore, analyses of orthologous IREs in A. aegypti field mosquitoes sampled worldwide identified an integrated polymerase IRE fragment that appeared under purifying selection within several million years, which supports a functional role in the host's biology. These results show that A. aegypti was subjected to repeated Rhabdovirus infectious episodes during its evolution history, which led to the accumulation of many integrated sequences. They also suggest that like retroviruses, integrated rhabdoviral sequences may participate actively in the evolution of their hosts.

Molecular Bases of Plant Resistance to Arthropods

Arthropod-resistant crops provide significant ecological and economic benefits to global agriculture. Incompatible interactions involving resistant plants and avirulent pest arthropods are mediated by constitutively produced and arthropod-induced plant proteins and defense allelochemicals synthesized by resistance gene products. Cloning and molecular mapping have identified the Mi-1.2 and Vat arthropod resistance genes as CC-NBS-LRR (coiled coil-nucleotide binding site-leucine rich repeat) subfamily NBS-LRR resistance proteins, as well as several resistance gene analogs. Genetic linkage mapping has identified more than 100 plant resistance gene loci and linked molecular markers used in cultivar development. Rice and sorghum arthropod-resistant cultivars and, to a lesser extent, raspberry and wheat cultivars are components of integrated pest management (IPM) programs in Asia, Australia, Europe, and North America. Nevertheless, arthropod resistance in most food and fiber crops has not been integrated due primarily to the application of synthetic insecticides. Plant and arthropod genomics provide many opportunities to more efficiently develop arthropod-resistant plants, but integration of resistant cultivars into IPM programs will succeed only through interdisciplinary collaboration.

New genome size estimates of 134 species of arthropods

Insect genome size diversity remains poorly sampled, with sparse and sporadic sampling of a few select orders and with many orders unrepresented or underrepresented in the genome size database. Here, we present 134 genome size estimates for 18 orders, including the first ever genome size estimates for eight orders, 38 families, 102 genera, and 131 species. Also reported here are three insect species genome size estimates that are corrections for unpublished genome size values that made it into the literature, including the smallest arthropod genome of the two spot spider mite (1C = 91 Mb). These estimates range from 91 to 7,752 Mb and provide a broader picture of genome size variation within Insecta and among all Arthropods. Proposed developmental constraints for holometabolous insect genome sizes are supported, with the majority of the species examined falling well under the hypothesized 1,978 Mb (2 pg) limit. The only exceptions occur in the highly diverse beetles (Coleoptera) (154 < 1C < 2,578) and the Mecoptera (1,890 < 1C < 2,134). Hemimetabolous orders include species with larger genomes. Of these, the Orthopteran genomes remain the largest (675 < 1C < 7,752). Hemipteran genomes are smaller (407 < 1C < 1,230), but with a very notable exception, the Cicadidae family, with bloated genomes (5,238 < 1C < 7,000 Mb) of a size previously known to exist only within Orthoptera. The genome sizes in all other Hemimetabolous orders (521 Mb < 1C < 3,860 Mb) fall within the range exhibited by the Hemiptera. New preparations and procedures for insect genome estimates with frozen material are described and discussed. To further improve and extend the number and quality of arthropod genome size estimates, a new flow cytometry insect standard, Blattaria: Periplaneta americana♂ (1C = 3,338 Mb) is established for genome size determination of insect genomes in excess of 2,000 Mb.

The complexity of Rhipicephalus (Boophilus) microplus genome characterised through detailed analysis of two BAC clones

BackgroundRhipicephalus (Boophilus) microplus (Rmi) a major cattle ectoparasite and tick borne disease vector, impacts on animal welfare and industry productivity. In arthropod research there is an absence of a complete Chelicerate genome, which includes ticks, mites, spiders, scorpions and crustaceans. Model arthropod genomes such as Drosophila and Anopheles are too taxonomically distant for a reference in tick genomic sequence analysis. This study focuses on the de-novo assembly of two R. microplus BAC sequences from the understudied R microplus genome. Based on available R. microplus sequenced resources and comparative analysis, tick genomic structure and functional predictions identify complex gene structures and genomic targets expressed during tick-cattle interaction.ResultsIn our BAC analyses we have assembled, using the correct positioning of BAC end sequences and transcript sequences, two challenging genomic regions. Cot DNA fractions compared to the BAC sequences confirmed a highly repetitive BAC sequence BM-012-E08 and a low repetitive BAC sequence BM-005-G14 which was gene rich and contained short interspersed elements (SINEs). Based directly on the BAC and Cot data comparisons, the genome wide frequency of the SINE Ruka element was estimated. Using a conservative approach to the assembly of the highly repetitive BM-012-E08, the sequence was de-convoluted into three repeat units, each unit containing an 18S, 5.8S and 28S ribosomal RNA (rRNA) encoding gene sequence (rDNA), related internal transcribed spacer and complex intergenic region.In the low repetitive BM-005-G14, a novel gene complex was found between to 2 genes on the same strand. Nested in the second intron of a large 9 Kb papilin gene was a helicase gene. This helicase overlapped in two exonic regions with the papilin. Both these genes were shown expressed in different tick life stage important in ectoparasite interaction with the host. Tick specific sequence differences were also determined for the papilin gene and the protein binding sites of the 18S subunit in a comparison to Bos taurus.ConclusionIn the absence of a sequenced reference genome we have assembled two complex BAC sequences, characterised novel gene structure that was confirmed by gene expression and sequencing analyses. This is the first report to provide evidence for 2 eukaryotic genes with exon regions that overlap on the same strand, the first to describe Rhipicephalinae papilin, and the first to report the complete ribosomal DNA repeated unit sequence structure for ticks. The Cot data estimation of genome wide sequence frequency means this research will underpin future efforts for genome sequencing and assembly of the R. microplus genome.

Cloning and characterization ofmariner-like elements in the soybean aphid,Aphis glycinesMatsumura

Soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae), is currently the most important insect pest of soybean (Glycine max (L.) Merr.) in the United States and causes significant economic damage worldwide, but little is known about the aphid at the molecular level. Mariner-like transposable elements (MLEs) are ubiquitous within the genomes of arthropods and various other invertebrates. In this study, we report the cloning of MLEs from the soybean aphid genome using degenerate PCR primers designed to amplify conserved regions of mariner transposases. Two of the ten sequenced clones (designated as Agmar1 and Agmar2) contained partial but continuous open reading frames, which shared high levels of homology at the protein level with other mariner transposases from insects and other taxa. Phylogenetic analysis revealed Agmar1 to group within the irritans subfamily of MLEs and Agmar2 within the mellifera subfamily. Southern blot analysis and quantitative PCR analysis indicated a low copy number for Agmar1-like elements within the soybean aphid genome. These results suggest the presence of at least two different putative mariner-like transposases encoded by the soybean aphid genome. Both Agmar1 and Agmar2 could play influential roles in the architecture of the soybean aphid genome. Transposable elements are also thought to potentially mediate resistance in insects through changes in gene amplification and mutations in coding sequences. Finally, Agmar1 and Agmar2 may represent useful genetic tools and provide insights on A. glycines adaptation.

Tracking the Ancestry of a Deeply Conserved Eumetazoan SINE Domain

Transposable elements (TEs), such as short interspersed elements (SINEs), evolve rapidly and are generally restricted to specific lineages. Here, we demonstrate that a central core of the previously described Deu-domain located within DeuSINEs (Nishihara et al. 2006) is widely distributed throughout the Metazoa. We characterize five new SINEs with this core sequence from the genomes of cnidarians, molluscs, annelids, and arthropods. Because this domain can be traced back to the cnidarian-bilaterian split >600 Ma, we propose naming it the "Nin" domain (the meaning of the Japanese character "Nin" is to endure and hide). Given that conserved noncoding DNA, such as that derived from the activity of SINEs, can be functionally relevant for the host genome (Sasaki et al. 2008), our findings highlight the need to understand these functions and the roles they may have played in supporting the evolution of multicellular genomes.

Not just another genome

Sequence analysis of the Daphnia pulex genome holds some surprises that could not have been anticipated from what was learned so far from other arthropod genomes. It establishes Daphnia as an eco-genetical model organism par excellence.

Insulin/FOXO Signaling Regulates Ovarian Prostaglandins Critical for Reproduction

Abnormalities in insulin/IGF-1 signaling are associated with infertility, but the molecular mechanisms are not well understood. Here we use liquid chromatography with electrospray ionization tandem mass spectrometry to show that the C. elegans insulin/FOXO pathway regulates the metabolism of locally acting lipid hormones called prostaglandins. C. elegans prostaglandins are synthesized without prostaglandin G/H synthase homologs, the targets of nonsteroidal anti-inflammatory drugs. Our results support the model that insulin signaling promotes the conversion of oocyte polyunsaturated fatty acids (PUFAs) into F-series prostaglandins that guide sperm to the fertilization site. Reduction in insulin signaling activates DAF-16/FOXO, which represses the transcription of germline and intestinal genes required to deliver PUFAs to oocytes in lipoprotein complexes. Nutritional and neuroendocrine cues target this mechanism to control prostaglandin metabolism and reproductive output. Prostaglandins may be conserved sperm guidance factors and widespread downstream effectors of insulin actions that influence both reproductive and nonreproductive processes.

Genomic survey of the ectoparasitic mite Varroa destructor, a major pest of the honey bee Apis mellifera.

BackgroundThe ectoparasitic mite Varroa destructor has emerged as the primary pest of domestic honey bees (Apis mellifera). Here we present an initial survey of the V. destructor genome carried out to advance our understanding of Varroa biology and to identify new avenues for mite control. This sequence survey provides immediate resources for molecular and population-genetic analyses of Varroa-Apis interactions and defines the challenges ahead for a comprehensive Varroa genome project.ResultsThe genome size was estimated by flow cytometry to be 565 Mbp, larger than most sequenced insects but modest relative to some other Acari. Genomic DNA pooled from ~1,000 mites was sequenced to 4.3× coverage with 454 pyrosequencing. The 2.4 Gbp of sequencing reads were assembled into 184,094 contigs with an N50 of 2,262 bp, totaling 294 Mbp of sequence after filtering. Genic sequences with homology to other eukaryotic genomes were identified on 13,031 of these contigs, totaling 31.3 Mbp. Alignment of protein sequence blocks conserved among V. destructor and four other arthropod genomes indicated a higher level of sequence divergence within this mite lineage relative to the tick Ixodes scapularis. A number of microbes potentially associated with V. destructor were identified in the sequence survey, including ~300 Kbp of sequence deriving from one or more bacterial species of the Actinomycetales. The presence of this bacterium was confirmed in individual mites by PCR assay, but varied significantly by age and sex of mites. Fragments of a novel virus related to the Baculoviridae were also identified in the survey. The rate of single nucleotide polymorphisms (SNPs) in the pooled mites was estimated to be 6.2 × 10-5per bp, a low rate consistent with the historical demography and life history of the species.ConclusionsThis survey has provided general tools for the research community and novel directions for investigating the biology and control of Varroa mites. Ongoing development of Varroa genomic resources will be a boon for comparative genomics of under-represented arthropods, and will further enhance the honey bee and its associated pathogens as a model system for studying host-pathogen interactions.

A conserved gene cluster as a putative functional unit in insect innate immunity

The Nimrod gene superfamily is an important component of the innate immune response. The majority of its member genes are located in close proximity within the Drosophila melanogaster genome and they lie in a larger conserved cluster (“Nimrod cluster”), made up of non-related groups (families, superfamilies) of genes. This cluster has been a part of the Arthropod genomes for about 300–350 million years. The available data suggest that the Nimrod cluster is a functional module of the insect innate immune response.

The Distribution of GYR- and YLP-Like Motifs in Drosophila Suggests a General Role in Cuticle Assembly and Other Protein-Protein Interactions

BackgroundArthropod cuticle is composed predominantly of a self-assembling matrix of chitin and protein. Genes encoding structural cuticular proteins are remarkably abundant in arthropod genomes, yet there has been no systematic survey of conserved motifs across cuticular protein families.Methodology/Principal FindingsTwo short sequence motifs with conserved tyrosines were identified in Drosophila cuticular proteins that were similar to the GYR and YLP Interpro domains. These motifs were found in members of the CPR, Tweedle, CPF/CPFL, and (in Anopheles gambiae) CPLCG cuticular protein families, and the Dusky/Miniature family of cuticle-associated proteins. Tweedle proteins have a characteristic motif architecture that is shared with the Drosophila protein GCR1 and its orthologs in other species, suggesting that GCR1 is also cuticular. A resilin repeat, which has been shown to confer elasticity, matched one of the motifs; a number of other Drosophila proteins of unknown function exhibit a motif architecture similar to that of resilin. The motifs were also present in some proteins of the peritrophic matrix and the eggshell, suggesting molecular convergence among distinct extracellular matrices. More surprisingly, gene regulation, development, and proteolysis were statistically over-represented ontology terms for all non-cuticular matches in Drosophila. Searches against other arthropod genomes indicate that the motifs are taxonomically widespread.ConclusionsThis survey suggests a more general definition for GYR and YLP motifs and reveals their contribution to several types of extracellular matrix. They may define sites of protein interaction with DNA or other proteins, based on ontology analysis. These results can help guide experimental studies on the biochemistry of cuticle assembly.

Widespread Horizontal Gene Transfer from Double-Stranded RNA Viruses to Eukaryotic Nuclear Genomes

Horizontal gene transfer commonly occurs from cells to viruses but rarely occurs from viruses to their host cells, with the exception of retroviruses and some DNA viruses. However, extensive sequence similarity searches in public genome databases for various organisms showed that the capsid protein and RNA-dependent RNA polymerase genes from totiviruses and partitiviruses have widespread homologs in the nuclear genomes of eukaryotic organisms, including plants, arthropods, fungi, nematodes, and protozoa. PCR amplification and sequencing as well as comparative evidence of junction coverage between virus and host sequences support the conclusion that these viral homologs are real and occur in eukaryotic genomes. Sequence comparison and phylogenetic analysis suggest that these genes were likely transferred horizontally from viruses to eukaryotic genomes. Furthermore, we present evidence showing that some of the transferred genes are conserved and expressed in eukaryotic organisms and suggesting that these viral genes are also functional in the recipient genomes. Our findings imply that horizontal transfer of double-stranded RNA viral genes is widespread among eukaryotes and may give rise to functionally important new genes, thus entailing that RNA viruses may play significant roles in the evolution of eukaryotes.

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.