The Enhancer of split and Achaete-Scute complexes of Drosophilids derived from simple ur-complexes preserved in mosquito and honeybee
© Schlatter and Maier; licensee BioMed Central Ltd. 2005
Received: 15 July 2005
Accepted: 17 November 2005
Published: 17 November 2005
In Drosophila melanogaster the Enhancer of split-Complex [E(spl)-C] consists of seven highly related genes encoding basic helix-loop-helix (bHLH) repressors and intermingled, four genes that belong to the Bearded (Brd) family. Both gene classes are targets of the Notch signalling pathway. The Achaete-Scute-Complex [AS-C] comprises four genes encoding bHLH activators. The question arose how these complexes evolved with regard to gene number in the evolution of insects concentrating on Diptera and the Hymenoptera Apis mellifera.
In Drosophilids both gene complexes are highly conserved, spanning roughly 40 million years of evolution. However, in species more diverged like Anopheles or Apis we find dramatic differences. Here, the E(spl)-C consists of one bHLH (mβ) and one Brd family member (mα) in a head to head arrangement. Interestingly in Apis but not in Anopheles, there are two more E(spl) bHLH like genes within 250 kb, which may reflect duplication events in the honeybee that occurred independently of that in Diptera. The AS-C may have arisen from a single sc/l'sc like gene which is well conserved in Apis and Anopheles and a second ase like gene that is highly diverged, however, located within 50 kb.
E(spl)-C and AS-C presumably evolved by gene duplication to the nowadays complex composition in Drosophilids in order to govern the accurate expression patterns typical for these highly evolved insects. The ancestral ur-complexes, however, consisted most likely of just two genes: E(spl)-C contains one bHLH member of mβ type and one Brd family member of mα type and AS-C contains one sc/l'sc and a highly diverged ase like gene.
The Notch pathway is one of the best studied cell to cell communication systems in the animal kingdom. It is highly conserved and used from worm to man. This pathway is needed whenever cell decisions are influenced by cell-cell communication, and also during proliferation or pathway crosstalk . Using Drosophila melanogaster as a model system, the Notch pathway was intensely studied over many years. The best defined process governed by Notch is called "lateral inhibition": cells of a given fate are singled out from an equivalence group of the same fate, whereas the differentiation of the other cells is suppressed by the Notch signal. This happens for example during neurogenesis, where neuroblasts are selected from proneural clusters; they keep neural fate, whereas the surrounding cells eventually differentiate as epidermoblasts. The name giving transmembrane Notch-receptor interacts physically with the extracellular domain of the transmembrane ligands Delta or Serrate of the signalling cell. After this activation, the intracellular domain of Notch is cleaved and travels into the nucleus where it transcriptionally activates together with Suppressor of Hairless [Su(H)] genes of the Enhancer of split complex [E(spl)-C]. E(spl) gene products in turn repress the activity of proneural genes encoded for example by the Achaete-Scute-Complex [AS-C]. As consequence these cells stay undifferentiated to become epidermoblasts later on, whereas the signalling cell enters into the programmed neural cell fate. In the focus of our studies are these two complexes, E(spl)-C and AS-C, since in D. melanogaster they are composed of several genes with complex expression patterns and specific yet partly redundant functions.
The AS-C in D. melanogaster comprises four genes, achaete (ac), scute (sc), lethal of scute (l'sc) and asense (ase) that all encode transcriptional activators of the bHLH class. They determine proneural fate and are thus required for the development of the central and peripheral nervous system [13, 26–29]. These genes have been also studied in vertebrates. In the mouse there are three achaete-scute family members abbreviated ASH for Ac-Sc-Homolog : MASH-1 and XASH-3/CASH-4 are two members that are involved in the development of the nervous system . Both complexes are therefore good candidates to look for the changes that occurred during insect evolution.
In this work, we studied the evolution of the E(spl)-C and AS-C by making use of the recent advances in the genome projects of the Diptera D. pseudoobscura, D. virilis, Anopheles gambiae and the Hymenoptera Apis mellifera. The estimated distances are nearly 30 million years (Myr) separating D. melanogaster and D. pseudoobscura which both belong to the Sophophora subgenus and around 40 Myr separating the Sophophora from the Drosophila subgenus where D. virilis belongs to . The distance between these modern dipterans and the more ancient ones like the Culicidae A. gambiae is estimated 200–250 Myr and that between Diptera and Hymenoptera like the honeybee A. mellifera 250–300 Myr . Our study shows that both gene complexes are highly conserved in Drosophilids with regard to overall size, gene number and structural similarity of the encoded proteins. In contrast, more ancient dipterans like the mosquito and similarly also the honeybee have much simpler gene complex structures: the E(spl)-C consists of just one mβ like bHLH gene and one mα like Brd-type gene. However, in Apis two more bHLH/WRPW coding genes are found within about 200 kb and may reflect an enlargement of the E(spl)-C in this species. The AS-C consists of only one sc/l'sc like member and one further ase-like gene that is, however, widely diverged. These data suggest that the evolution of modern Drosophilids included an enlargement of these complexes, notably a multiplication of the genes in the E(spl)-C that seem to have subsequently specified their roles in Notch signalling pathway. Their strict conservation in Drosophilids argues for a diversification presumably driven by their highly specified expression patterns and regulatory activities.
The Enhancer of split complex in Drosophilids
The Enhancer of split complex [E(spl)-C] consists of 13 transcription units (Fig. 1): seven genes (m3, m5, m7, m8, mβ, mγ, mδ) encode highly related basic helix-loop-helix proteins, four have been grouped to the Brd-family (m2, m4, m6, mα); m1 which encodes a serine protease inhibitor and l(3) groucho (gro) which encodes a co-repressor of the E(spl) bHLH protein family. Although the E(spl) bHLH proteins are partly redundant and, therefore, a loss or addition of genes could be without consequences, the complex is highly conserved in all studied Drosophilids with respect to gene order and number, transcription orientation and overall size. As expected, the evolution rate of the orthologs is different. The best conservation is found between the Gro orthologs (more than 96% identity, Fig. 1), whereas M1 displays the highest evolutionary rate (less than 61% identity between D. melanogaster and D. virilis; Fig. 1).
The bHLH proteins of the E(spl)-C
In the Drosophilids, all E(spl) bHLH genes are without intron, and the proteins contain the typical bHLH and orange domains and end with the WRPW motif. The best conserved bHLH ortholog is Mγ followed by M3, Mβ, M5, Mδ and M7/8 (Fig. 1). However, the evolutionary rate varies quite strongly. Comparing D. melanogaster with D. virilis, the highest identity score is found for Mγ with ~92% (Fig. 1A,B) and the lowest for M8 with ~81%. An even more striking difference in the identity scores is observed when comparing the D. melanogaster and D. pseudoobscura orthologs M3 and M7 (93% versus 73%; Fig. 1A,C). The low degree of conservation of the M7 proteins is rather surprising. It is based on one hand on peculiar size variations: 206 residues in D. pseudoobscura, 186 in D. melanogaster and 197 in D. virilis (Fig. 1C). On the other hand, the amino acid composition of the bHLH and orange domains is much better conserved between D. melanogaster and D. virilis than between D. melanogaster and D. pseudoobscura. This is different from all the other bHLH orthologs: the bHLH domains of M3 and Mβ are identical in all three species and also the other ones are extremely similar. Only one conservative change is detected in the Mγ bHLH domain of D. melanogaster compared with D. virilis, and just two in M5. The bHLH domains of the M8 and Mδ proteins contain also single non-conservative replacements, apart from a few conservative changes. However, D. pseudoobscura M7 shows an unusual high number of changes – six replacements and five conservative changes – within the bHLH domain if compared to the melanogaster ortholog (Fig. 1C).
The Bearded-protein family in the E(spl)-C
In general, the Brd-type proteins evolve faster than the bHLH proteins: M4 and Mα are the best conserved members with ~82% identity between virilis and melanogaster and are thus within the range of the fastest evolving m7/8 bHLH coding genes (Fig. 1). The overall structure of M4 and Mα orthologs is quite similar in all studied Drosophila species. The so-called 'bearded'-domain is completely identical in the Mα orthologs, M4 has only a few gaps. However, the predicted D. virilis m4 gene has an extended open reading frame of novel 132 residues at the 5' end, whereas the remaining 156 residues are conserved. The melanogaster m4 5' region reveals similarity at the DNA level, however, has no open reading frame. Therefore it remains questionable whether the larger open reading frame in D. virilis is indeed translated. The other two Brd-family members, M2 and M6, are much less conserved: Only approximately 68% identity is found between the respective orthologs of D. melanogaster and D. virilis. Despite this little conservation, the typical Brd protein domains can still be recognized. The most prominent is the predicted basic amphipathic α-helix domain in the N-terminal protein region .
The m1 gene in the E(spl)-C
The m1 gene encodes a protein that has the signatures of serine protease inhibitors . Despite a low degree of conservation which ranges between 60 and 70% identity between the M1 orthologs (Fig. 1A), the structurally important cysteines residues are conserved in number and spacing . Notably, the orthologous genes of D. pseudoobscura and D. virilis have a significantly longer open reading frame at the 5'end that extends the proteins for approximately 50 residues to 203 in the case of D. virilis. The extended protein parts share ~68% similarities within the first 30 residues between the two orthologs. Furthermore, the first nine residues have only one conservative exchange arguing for its translation in vivo. In melanogaster all three reading frames at the 5' end are blocked by several stop codons, excluding a likewise 5' extension. However, at the DNA level there are identities of 69% to the virilis ortholog and 79% to the pseudoobscura ortholog which could be also interpreted as conserved regulatory sequence.
The E(spl)-C in Anopheles gambiae
Identity matrix of E(spl) bHLH proteins in D. melanogaster (% similarity/identity)
Approximately 100 kb from the 3' end of A.g.mβ, we detected another sequence that might encode an E(spl) bHLH type protein (Fig. 2A; ENSANGG00000017548). However, the presumptive gene product is more highly related to D. melanogaster Deadpan (Dpn; 57% identity) than to E(spl) bHLH Mβ (52% identity). The conservation extends beyond the amino acid sequence: we find the same intron/exon structure in this Anopheles gene as in the deadpan and hairy genes from Drosophila (see also Fig. 8a). We believe that this Anopheles protein corresponds to Dpn (A.g.Dpn) rather than to Hairy, since it shares little more than 47% identity with D. melanogaster Hairy. The best hit with the D. melanogaster Hairy protein is found on the second chromosome in Anopheles (A.g.h; 72% identity at protein level). Other genes of the Drosophila E(spl)-C were not detected nearby: maybe they are not conserved enough to be discovered like e.g. m6 or they are located at totally different positions in the genome like gro.
The E(spl)-C in Apis mellifera
For example, the A.m.mα homolog is predicted to consist of five exons that code for a protein with 402 residues. This is considerably larger than the 138 residues of D.m.Mα. Again, we propose a different gene structure based on the analysis of the genomic DNA, where the translation extends into the first postulated intron and ends shortly afterwards. The encoded protein than consists of only 168 amino acids and terminates with residues that are very similar to the Mα Drosophila homolog (M Q V A) (Fig. 3A). Moreover, it shows the typical amphipathic α-helix domain in the N-terminal part (Fig. 3B). The other motifs are also conserved, the extreme C-terminus is, however, only similar (Fig. 3A). This protein shows also high similarity to Drosophila Twin of m4 (Tom) that might slightly exceed that to D.m.Mα depending on the parameters used. As in mosquito, no other Brd-like protein was found in the honeybee database such that a Brd-complex seems non-existent, unlike in Drosophilids .
The Apis database annotates a single intron for the A.m.mβ gene that conforms to the GT AG rule. This is in contrast to all the Drosophila E(spl) bHLH genes that are intronless. Moreover, there are several possible start sites and it remains unclear which one is used. In D. melanogaster, there are other bHLH-WRPW encoding genes that contain introns, like deadpan (dpn) or Side. However, the respective protein sequences do not align well with A.m.Mβ, and the Apis genome contains predicted orthologs to these genes at other locations (see below). The A.m.Mβ protein is also highly similar to melanogaster Mγ with 72.3 versus 77.7% similarity and 65 versus 67.5% identity compared to D.m.Mβ (Fig. 4B). These differences are extremely small. In fact, the Apis database proposes this gene as mγ based on the similarity within the bHLH domain. However, we used for our comparison the entire protein sequences and calculated identity overall. Moreover, based on the E(spl)-C structure in Anopheles, we favour the hypothesis that mβ is the ancestral gene.
Within the ~250 kb contig (GroupUn. 159), there are two further stretches that might encode E(spl) type bHLH proteins. One is located about 150 kb apart. This gene contains no introns and the encoded protein shows homology to both Mβ and Mγ. In this case, the similarity seems slightly higher to Mγ than to Mβ (72.2 vs. 69.3% similarity and 66 vs. 62.7% identity). Therefore, we call the gene A.m.mγ (Fig. 4A,B). Both FlyBase and BlastN give a higher score to D. melanogaster Mβ than to Mγ. Presumably, both use similar paradigms based on an alignment of only the best conserved sequences, whereas we used less stringent parameters (see Methods) to allow an alignment of the complete sequences. As to be expected, the two E(spl) honeybee proteins A.m.Mβ and A.m.Mγ are highly related to each other with a similarity of 73% and an identity of 67%. Interestingly, this numbers are very similar to those from a likewise comparison of the D. melanogaster Mβ and Mγ homologs (Table 1).
Another 50 kb further up at position 50 (Fig. 4A), we found a short alignment to the E(spl) bHLH domain. However, there was no predicted gene, nor a start codon, nor a WRPW motif. Nearby at position ~46 kb there is an open reading frame of 133 residues split by one intron that belongs to a predicted database gene consisting of five exons (GENSCAN00000025907; black in Fig. 4C). However, lacking any similarity to known genes of Apis or other species, this gene remained without functional prediction in the Apis database. In agreement, our searches in the FlyBase did not detect any similar sequences. However, we predict an E(spl) bHLH-type protein encoded by this gene region: extension of the open reading frame into the adjacent intron ends in WRPW (Fig. 4C, blue exons). The bHLH encoding sequences (purple in Fig. 4C) are located within the second predicted intron of the putative gene shown in black. There are respective exon/intron boundary consensus sequences to allow for a single transcript that contains the bHLH domain, an orange domain and ends with the WRPW motif (Fig. 4C). This third E(spl) bHLH gene would comprise five exons. In fact, in the update May 2005 of the Ensembl honeybee database, the Ensembl automatic analysis pipeline predicts a very similar protein, however, with a different N-terminus and slightly smaller third and fourth exons. This gene than consists of four exons which would be similar to D. melanogaster hairy and dpn that contain two introns. However, the encoded protein is most similar to the E(spl) protein Mβ, so we call the gene A.m.mβ'. Since Apis hairy and dpn are found elsewhere in the genome (see below), we propose that the E(spl)-C in Apis mellifera consists of the ur-complex plus two further E(spl) bHLH genes most closely related to mγ and mβ. No other genes of the Drosophila E(spl)-C are present in that of the honeybee. We find a highly conserved Groucho ortholog, however, at a completely different position in the genome.
Conservation of other Hairy/E(spl)-like proteins known from Drosophila
Gene annotation used by the respective databases
(contig 4374 Contig 4847)
D.m. mδ (CG8328)
D.p. mδ (178 652-178 089)
D.m. mγ (CG8333)
D.p. mγ (175 594-174 953)
A.m. mγ (ENSAPMG0000004887)
D.m. mβ (CG14548)
D.p. mβ (169 335–169 934)
A.m. mβ (ENSAPMG0000004881)
A.g. mβ (ENSANGG00000017601)
D.m. mα (CG8337)
D.p. mα (163 786-163 358)
A.m. mα (GENSCAN0000001764)
A.g. mα (SNA00000011401)
D.m. m1 (CG8342)
D.p. m1 (158 207–158 734)
D.m. m2 (CG6104)
D.p. m2 (157 012–157 407)
D.m. m3 (CG8346)
D.p. m3 (151 792-151 136)
D.m. m4 (CG6099)
D.p. m4 (149 049–149 516)
D.m. m5 (CG6096)
D.p. m5 (142 321–142 881)
D.m. m6 (CG8354)
D.p. m6 (137 883-137 656)
D.m. m7 (CG8361)
D.p. m7 (134 008-133 391)
D.m. m8 (CG8365)
D.p. m8 (130 185-129 628)
D.m. gro (CG8384)
D.p.gro (148 473-117 489)
A.m. mβ' (ENSAPMG00000016895)
D.m. stich1 (CG17100)
A.m. stich1 (ENSAPMG00000005857)
A.g. stich1 (ENSANGG00000016365)
D.m. side (CG10446)
A.m. side (ENSAPMG0000000088)
A.g. side (ENSANGG00000014329)
D.m. dpn (CG8704)
A.m. dpn (ENSAPMG00000004551)
A.g. dpn (ENSANGG00000017548)
D.m. Hey (CG11194)
A.m. Hey (ENSAPMG0000000726)
A.g. Hey (ENSANGG00000021744)
D.m. Her (CG5927)
D.m. h (CG6494)
A.m. h (ENSAPMG00000004545)
A.g. h (ENSANGG00000018369)
D.m. ac (CG3796)
D.m. sc (CG3827)
D.m. l'sc (CG3839)
A.m. ash (ENSAPMG00000003261)
A.g. ash (ENSANGG00000010650(Q95VY6)
D.m. ase (CG3258)
A.m. ase (ENSAPMG00000003265)
A.g. ase (ENSANGG00000015341)
D.m. da (CG5102)
A.m. da (ENSAPMP00000005673)
A.g. da (ENSANGEST00000361691/SNAP000000012539)
D.m. Ocho (CG5138)
D.m. Tom (CG5185)
D.m. Brd (CG3096)
The Achaete-Scute complex in Drosophilids
Of the four AS-C gene members in D. melanogaster, ase stands out because it is much larger than the other three. In D. virilis, the size increase is even more striking: D.v.Ase is predicted to comprise 619 residues, whereas D.m.Ase is only 486 residues in length (Fig. 5B). This extension of more than 20% additional residues is caused by multiple insertions of repetitive sequences that code for poly-glutamine (Q), poly-alanine (A) and poly-asparagine (N) stretches (Fig. 5B). Like in D. melanogaster the unrelated gene pepsinogen-like (pcl) is located between l'sc and ase (Fig. 5A).
The AS-C in Anopheles gambiae
The second presumptive gene has two predicted introns. The derived amino acid sequence shares between 54% and 57% identity with all four D. melanogaster AS-C proteins. By searching the D. melanogaster genome with the predicted protein sequence, the best hit was found to Ase protein (Fig. 6B). Accordingly, the Anopheles database defined it as ase homolog and so we named it A.g.ase. We note, however, that a precise appointment is difficult based on the lack of a significant similarity at the C-terminus which normally allows the distinction between the AS-C members.
The AS-C in honeybee
Interpretation of the second gene is much more difficult. The annotation by Ensembl automatic pipeline using GeneWise model based on either protein or aligned EST's resulted in a coding sequence within a single exon that lacked both start and stop codons, however, aligned well with AS-C bHLH domains. We thus extended our studies into the surrounding genomic DNA, where we detected one further open reading frame. The deduced amino acid sequence matched well with the N-terminal part of D. melanogaster Ac and Ase proteins. However, the C-terminus did not align convincingly (Fig. 7B/GW). Another gene prediction using Chris Burge's Genscan program  gave a transcript of 8 exons spanning over 8 kb of genomic DNA. The translation gave a larger protein again without start methionine that contained the single exon predicted by the GenWise model (Fig. 7B/GS). Again there was very little similarity in the C-terminal part of this and Drosophila AS-C proteins (Fig 7B). Because Genscan could not predict the protein start, we propose a combination of both models with the N-terminus as shown in the A.m.Ase/GW sequence and the C-terminus as in A.m.Ase/GS that might, however, end shorter than shown in Fig. 7B.
The high divergence from the Drosophila AS-C proteins renders precise predictions very difficult. In fact, under standard conditions like FlyBase BlastN only parts of the bHLH domain can be identified. Comparison with the D. melanogaster AS-C proteins gave minimally different scores, with the highest score found with D.m.Ac followed by D.m.Ase, dependent on the parameters. For example, the Apis database finds best scores with melanogaster Sc and L'sc. We named this gene A.m.ase by its similarity notably in the N-terminal part and based on the arrangement of the mosquito AS-C.
Conservation of predicted regulatory elements in E(spl)-C and AS-C
Neurogenesis in Drosophila is subjected to various levels of regulation. As described in the introduction, E(spl) bHLH proteins repress proneural gene activity. Negative regulation is brought about by repression of transcription as well as at the protein level [36, 39–42]. A further mode of regulation involves RNA:RNA duplexes . These are formed by small sequence stretches located in the 3' untranslated region (3' UTR) of the mRNAs of proneural AS-C genes and different members of the E(spl)-C. For example, l'sc mRNA contains so-called proneural boxes (PB-box) (GGAAGAC) which bind to the GY-boxes (GUCUUCC) of E(spl) m4 RNA . We searched the genomic sequences adjacent to the coding sequences for respective regulatory elements and found them in D. virilis: there are two PB boxes in D.v.l'sc and GY-boxes in the 3'UTR of the predicted virilis E(spl) genes m3, m4, m5 and mγ just as in D. melanogaster and meanwhile published by Lai et al. . Moreover, it has been shown that E(spl) genes are direct targets of the Notch signal involving the DNA binding protein Su(H) [11, 13, 15, 19]. Accordingly, there are potential Su(H) binding sites (C/TGTGA/GGA) in all D. melanogaster E(spl) genes including m2 and m6 [14, 15], which are also present in the respective virilis orthologs . However, the predicted binding sites for proneural bHLH activators (E box: GCAGGTG)  are less well conserved during evolution. Whereas the D. virilis m2 and m4 orthologs contain such a regulatory element, we found no sequence fitting the E-box consensus in either the D.v.m6 or D.v.mα sequence.
In mosquito and honeybee, regulatory elements of RNA:RNA duplex-type were not detected. None of the two AS-C genes of either Apis or Anopheles contained PB-box like sequences in the 3' UTR, albeit the highly diverged sequence and gene structure of A.g.ase does not allow to definitely exclude their presence. Since there is no predicted m4 ortholog in Anopheles and Apis, we looked at the 3' end of the mα gene as the single Brd-family member. However, in none of these gene sequences did we find GY type boxes like in Drosophila. We note that some of the predicted gene structures are still incomplete and a search for small sequence stretches is notoriously difficult if one allows for variations. Therefore, there is the formal possibility that we have missed these sites.
The Enhancer of split complex
Extensive genome analyses in the recent years revealed that there are not many examples of large gene complexes that are widely conserved. Prominent examples are the HOX (h omeobox) complexes, which contain homeotic genes in Drosophila. HOX complexes are well conserved in metazoans despite some variations in gene number. HOX-genes encode regulatory proteins with specific individual functions and mutations affect different aspects of the body plan . Not surprisingly, it is almost only the homeodomain, which serves as sequence-specific DNA binding motif that is conserved amongst different species . In contrast, similarity amongst bHLH proteins encoded by the E(spl)-C extends over the entire length, even within the same species indicating rather recent duplication events. The D. melanogaster proteins M8/M5 and Mβ/M3 are most similar with over 70% identity, whereas Mδ is the most diverged. However, Mδ still shares at least 50% identity with other E(spl) bHLH protein members (see pair wise comparison in the identity matrix of Mδ with M8, M5 and M3; Table 1). More interesting is the analysis of the overall similarity amongst these proteins. Here, any one of the proteins is compared with the other six and the result is averaged. Clearly, Mβ (73/64%, similarity/identity) closely followed by Mγ (72/63%) is most similar to all others, whereas Mδ (66/55%) shows the lowest values (Table 1). One interpretation might be that the different bHLH genes evolved by duplication out of mβ or mγ. Remarkably, these two bHLH proteins besides M3 are the best conserved in the three Drosophila species (Fig. 1A). We would like to postulate that these are the most ancient proteins with the most general function and, therefore, the highest selection pressure. This hypothesis is supported by the finding that mβ has the most general expression pattern from which the others can be derived by a decrease of gene activity . The conspicuous conservation of M3 might hint to an important function during egg development as this gene is expressed also maternally [17, 22]. The high degree of conservation of all E(spl) bHLH orthologous proteins in Drosophilids, which is clearly higher than the similarity within this protein family in D. melanogaster, indicates specific and non-redundant roles during development (see also ). Some of these functions have been identified in the past [19, 39]. It is conceivable that regulatory sequences were not duplicated or evolved more rapidly so that we now find highly dynamic expression patterns of these genes.
The ancestral E(spl)-C is composed of mβ and mα
As outlined above, mβ appears to be the ancestral bHLH gene of the E(spl)-C in Drosophilids based on its great similarity with all the other bHLH proteins. This assumption is strongly supported by the sequence conservation of the E(spl) bHLH proteins in A. gambiae and A. mellifera. The single E(spl) bHLH protein encoded by the mosquito genome has the highest identity to Mβ. The genome of honeybee contains three prospective genes that encode proteins most highly related to E(spl) D.m.Mβ and D.m.Mγ. All three are clustered within a single sequence contig, albeit they span a large segment of about 250 kb, whereas the whole E(spl)-C in D. melanogaster comprises roughly 50 kb. Despite the fact that two of these genes possess introns just within the bHLH domain and at positions close to the ones found in the D. melanogaster genes dpn, hairy or Her (Figs. 4, 8), the amino acid sequence similarity classifies them clearly as E(spl) bHLH proteins. A comparison of Anopheles and Apis proteins reveals, that the presumptive Mβ homologs have highest similarity (83%) and identity (76%), whereas the protein that we classified as A.m.Mγ is just 70% similar and 66% identical to A.g.Mβ.
In Drosophilids, mα is located close to mβ and is transcribed in the opposite direction (head to head; Fig. 1A). This arrangement is likewise found in Anopheles and Apis (Figs. 2A, 4A). Notably, A.m.mα is next to A.m.mβ, whereas the two Apis A.m.mγ and A.m.mβ' genes are much further apart (Fig. 4). We find this arrangement to be very ancient. In the beetle Tribolium, which on the tree of evolution is found even more deeply rooted (~300 Myr to Dipterans ), two similar genes coding for Mβ-like proteins (~65% and ~67% identity to D.m.Mβ) are found and one is within ~18 kb to a gene coding for an Mα-like protein (~52% identity to D.m.Mα) (unpublished data derived from the Tribolium database). We postulate that the ur-complex consisted of these two ancestral genes, mα and mβ. It is intriguing that they belong to the two different classes of Notch-responsive genes in the E(spl)-C, the bHLH and the Brd-class. In the fly, mα and bHLH genes are similarly expressed [14, 15, 22]. It is not unlikely that they share common regulatory elements that could explain their co-segregation in the process of evolution.
Comparison between Apis mellifera, D. melanogaster and Anopheles gambia (% similarity/identity)
The Achaete-Scute complex
Genes related to achaete or scute have been identified in a large number of species, from hydra  to mouse , and so we expect these also in the different insects. The AS-C was most intensely studied in various species of Schizophora flies, apart from Drosophila [28, 37, 48–53]. The number of genes varies between one and four, however, is not strictly correlated with the position in the phylogenetic tree. For example, AS-C of Calliphora vicina contains three genes, whereas other dipteran flies like Drosophila contain four. Two genes are found in the branchiopod crustacean Triops longicaudatu s  and only one in hydra . In Dipteran flies the expression patterns of the proneural genes are largely varied. This is regulated by positional information through the Iroquois Complex and pannier and in addition by a transcriptional feed-back loop involving AS-C proteins. Eventually, neural precursors are selected by the repressive activity of E(spl) bHLH proteins [55, 56]. Thereby, location and number of the large bristles on the notum is precisely controlled. The mosquito is covered with rows of large sensory bristle, where number and position varies between individuals . This is in accordance with the fact that there is only one scute-like gene, A.g.ash that is expressed all over the presumptive notum in a modular pattern . Recently it was shown that the Anopheles A.g.ash gene can mimic the endogenous Drosophila genes and that overexpression leads to many ectopic bristles .
Albeit the bristle pattern on the notum of different Drosophilids varies slightly, bristle number and position is highly stereotyped . Therefore, it is not surprising to find the AS-C highly conserved within Drosophilids. Yet, the rate of change came unexpected and is quite remarkable outside of the bHLH domain. Compared to E(spl) bHLH proteins, those encoded by AS-C have a rather low degree of similarity, most notably Ac. In fact, the big flesh fly Calliphora vicina, which like Drosophila belongs to the Schizophora, is totally lacking the ac gene and is covered with bristles . In agreement, we were unable to find ac in Anopheles or Apis, arguing for rapid evolution. The best conservation rate is found in Sc and L'sc suggesting high evolutionary pressure and maybe common ancestry. Not only the bHLH domain, but also two small stretches outside (aa 203; SPTPS in D. melanogaster L'sc) and also the C-terminus are of high similarity, the latter found identical in Calliphora . Presumably these protein domains are of functional importance. Indeed, the C-terminus acts as transcriptional activation domain and is also used to recruit E(spl) bHLH proteins . Although the alignments of the respective genes of honeybee and mosquito to sc and l'sc are very similar, the tendency goes to a closer relationship to l'sc. However, we propose that this gene pair arose by duplication in the course of Drosophilid evolution, such that we may be looking at a common ancestor in the other two species.
The rate of conservation is very limited for the Ase homologs. Decent conservation is found within the bHLH domain, and moreover, a further well-conserved box is present (NGxQYxRIPGTNTxQxL; x are differences between A. gambiae and D. melanogaster). This sequence is likewise detected in the Ase protein of C. vicina, which however shares many more similarities with D.m.Ase . In Apis, there is no such conservation outside of the bHLH domain, which itself is highly diverged. The overall degree of conservation is so poor that further statements about the relationship are difficult. We argue that this gene represents A.m.ase by its close proximity to A.m.ash, although other interpretations are similarly possible. An analysis of its expression pattern in honeybee may help to solve these questions.
We aimed towards an understanding of the evolution of E(spl)-C and AS-C complexes which in D. melanogaster comprise genes of apparent redundant functions. Our analysis covered insect species that belong to the orders Hymenoptera (honeybee) and Diptera and there to the suborders Nematocera (mosquito) and Brachycera (three species of the genus Drosophila) and thus spans an about 300 Myr window of evolution. We find that both E(spl)-C and AS-C expanded rather recently as they are only present in their nowadays complex structures in Drosophilids. In Apis and in Anopheles, we find very similar arrangements indicative of an ancient ur-complex. The E(spl)-C seems to have evolved from two genes, one HES-like and one Brd-like gene that are arranged in a head to head orientation. Both types of genes are responsive to Notch signalling in Drosophila. Our data suggest that the most ancient genes are E(spl) bHLH mβ and E(spl) m α from which the other E(spl)-C genes derived by duplication and subsequent change. Moreover, an E(spl) ur-complex is likewise detected in Tribolium castaneum that belongs to the order Coleoptera. In Drosophila the complex also gained unrelated genes like m1 and gro. The latter is highly conserved, however, located at different genomic positions. Whereas in Anopheles the ur-complex seems to exist in its original form, two additional mβ-like bHLH genes are found in the Apis genome that possess introns. These introns are at similar positions as the introns of two other HES-like genes, dpn and h which themselves are highly conserved in the three insect species, arguing for a common evolutionary history. Presumably, the introns are evolutionarily ancient as they are also found in the C. elegans E(spl)/h like gene lin-22. The AS-C seems to originate from a single sc/l'sc like bHLH gene and a second largely diverged bHLH gene that shares similarity with Drosophila ase. The high degree of variation in the latter makes it difficult to conclusively decide on the original arrangement of this gene complex.
Drosophila melanogaster gene and protein sequences were accessed in FlyBase . The D. pseudoobscura database is found at the Human Genome Sequencing Center . The genome of D. virilis is sequenced by Agencourt Bioscience Corporation and can be downloaded  or searched . The Apis mellifera and Anopheles gambiae databases can be accessed with the Ensembl genome browser . The honeybee genomic sequence is also available at the Human Genome Sequencing Center  and the EST sequences were obtained from the honeybee brain EST project . Table 2 lists accession numbers of genes, contig numbers and positions of genes therein as well as accession numbers to gene predictions. Due to updates of the databases, gene annotation and names may be different now than reported here.
Sequences were downloaded and further studied applying the HUSAR programs of the Deutsche Krebsforschungszentrum [64, 65]. Genomic DNA was translated with MAP. BESTFIT and GAP programs were used for alignments and calculation of similarity and identity scores. DOMAIN SWEEP was applied to define bHLH and orange domains, respectively. The amphipathic α-helices were predicted and drawn using WHEEL.
D. pseudoobscura sequence acquisition
The annotation of the D. pseudoobscura database is well advanced. TBlastN searches starting from FlyBase  or using D. melanogaster genes allow easy access to both E(spl)-C and AS-C gene complexes. The database also gives information on transcript length, structure and orientation.
D. virilis sequence acquisition
Searches were done with BlastN. However, the D. virilis database allows this scan only with DNA against DNA. We therefore searched first for a characteristic part of a D. melanogaster gene (i.e. the bHLH domain) using lowest percent identity option (75% identity) to find the respective virilis ortholog. We then used the FIND program to identify the respective contig containing these sequences within the downloaded genomic sequence. Subsequently, the entire contig was translated in all six possible reading frames and manually analysed for E(spl)-C and AS-C genes and proteins. The complete E(spl)-C was located in one contig. To identify the whole sequence of the AS-C, three overlapping contigs had to be investigated. However, the contigs covering the AS-C still contained many large unsequenced or uncertain stretches. Therefore, the exact size of the complex could not be defined. The sequence gaps do not affect the coding sequences of the studied genes.
A. gambiae and A. mellifera sequence acquisition
The other databases were screened with D. melanogaster protein sequences for members of E(spl)-C and AS-C against genomic DNA. After detecting convincing similarities, the surrounding genomic DNA was downloaded for further studies with the HUSAR programs. Most of the genes that we describe here have also been annotated as transcribed regions. However, the majority of the database predictions were either incomplete or inconsistent. For our predictions, we carefully analysed genomic DNA and possible coding sequences. All six reading frames were searched for bHLH domain sequences and in the case of the E(spl)-C bHLH proteins also for WRPW motifs. Predicted introns were scanned for open reading frames and their borders reinvestigated. Exon/intron boundaries were defined by obeying to the GT AG rule. Afterwards it was analysed whether the newly predicted introns affected the open reading frame.
Classification of gene homology
Gene homology and orthology was classified based on sequence identity with the respective D. melanogaster protein. In the studied Drosophila species all analysed genes have been identified in a 1:1 ratio and could therefore be classified as orthologs. The degree of similarity and identity between two related proteins was determined with the BestFit program . Comparison of proteins from the three Drosophila species was done under the pre-configured standard conditions of the program (Gap weight: 8, length weight: 2). A TBlastN search in the Apis and Anopheles databases always returned several sequences, and we analysed up to 15 of the best hits. If the hit was in a region without predicted gene, we translated the respective sequence and analysed the open reading frames manually. They were analysed for expected protein motifs, like bHLH, orange or WRPW domain. The predicted protein sequences were then used for a BlastP search in FlyBase. In case of a 1:1 relationship, the genes were classified as orthologs. Since there are several E(spl)-C bHLH and AS-C bHLH genes in Drosophila, a 1:1 allocation was not possible, therefore, we classified them as homolog (see results). Protein sequences derived from the Apis or Anopheles genome projects are so diverged, that standard conditions only align the best conserved domains. We therefore changed the conditions to Gap weight: 1, length weight: 1, max. penalized length: 30. These relaxed conditions had little influence on either alignment of well-conserved sequences, the similarity or identity values. For example, the D.m.Mβ and D.v.Mβ orthologs share 89% identity under stringent conditions versus 93% identity under relaxed conditions. The reduced stringency however, allowed an alignment of the entire protein sequence also of the diverged proteins with the consequence that the identity values increased considerably compared to standard conditions. For example an alignment of A.g.Ash (371 residues) with D.m.L'sc (257 residues) under standard conditions aligns only the residues A.g.Ash 50–233 with D.m.L'sc 40–218 with an identity of 46%. Under relaxed conditions the whole protein sequences align with 64% identity.
We are indebted to the different Genome Projects for free and easy access to genomic data resources, the FlyBase, the Human Genome Sequencing Center, the Agencourt Bioscience Corporation, the Ensembl Genome Browser and the Honeybee Brain EST Project at the University of Illinois. We thank A. Preiss for constant support and AP, A.C. Nagel and an anonymous reviewer for critically reading the manuscript.
- Artavanis-Tsakonas S, Rand MR, Lake RJ: Notch signaling: cell fate control and signal integration in development. Science. 1999, 284: 770-776. 10.1126/science.284.5415.770.View ArticlePubMedGoogle Scholar
- Welshons WJ: Dosage sensitive interactions with split mutations in the presence of an Enhancer of split. Drosophila Inf Service. 1956, 30: 157-158.Google Scholar
- Dawson SR, Turner DL, Weintraub H, Parkhurst SM: Specificity for the hairy/Enhancer of split basic helix-loop-helix (bHLH) proteins maps outside the bHLH domain and suggests two separable modes of transcriptional repression. Mol Cell Biol. 1995, 15: 6923-6931.PubMed CentralView ArticlePubMedGoogle Scholar
- Delidakis C, Artavanis-Tsakonas S: The Enhancer of split locus of Drosophila encodes seven independent helix-loop-helix proteins. Proc Natl Acad Sci USA. 1992, 89: 8731-8735.PubMed CentralView ArticlePubMedGoogle Scholar
- Hartley D, Preiss A, Artavanis-Tsakonas S: A deduced gene product from the Drosophila neurogenic locus, Enhancer of split, shows homology to mammalian G-protein β subunit. Cell. 1988, 55: 785-795. 10.1016/0092-8674(88)90134-1.View ArticlePubMedGoogle Scholar
- Klämbt C, Knust E, Tietze K, Campos-Ortega JA: Closely related transcripts encoded by the neurogenic gene complex Enhancer of split of Drosophila melanogaster. EMBO J. 1989, 8: 203-210.PubMed CentralPubMedGoogle Scholar
- Knust E, Tietze K, Campos-Ortega JA: Molecular analysis of the neurogenic locus Enhancer of split of Drosophila melanogaster. EMBO J. 1987, 6: 4113-4123.PubMed CentralPubMedGoogle Scholar
- Knust E, Schrons H, Grawe F, Campos-Ortega JA: Seven genes of the Enhancer of split complex of Drosophila melanogaster encode Helix-loop-helix proteins. Genetics. 1992, 132: 505-518.PubMed CentralPubMedGoogle Scholar
- Paroush Z, Finley RL, Kidd T, Wainwright SM, Ingham PW, Brent R, Ish-Horowitz D: groucho is required for Drosophila neurogenesis, segmentation, and sex determination and interacts directly with hairy-related bHLH proteins. Cell. 1994, 79: 805-815. 10.1016/0092-8674(94)90070-1.View ArticlePubMedGoogle Scholar
- Preiss A, Hartley DA, Artavanis-Tsakonas S: The molecular genetics of Enhancer of split, a gene required for embryonic neural development in Drosophila. EMBO J. 1988, 7: 3917-3927.PubMed CentralPubMedGoogle Scholar
- Eastman DS, Slee R, Skoufos E, Bangalore L, Bray S, Delidakis C: Synergy between Suppressor of Hairless and Notch in regulation of Enhancer of split mγ and mδ expression. Mol Cell Biol. 1997, 17: 5620-5628.PubMed CentralView ArticlePubMedGoogle Scholar
- Jennings B, Preiss A, Delidakis C, Bray S: The Notch signaling pathway is required for Enhancer of split bHLH protein expression during neurogenesis in the Drosophila embryo. Development. 1994, 120: 3537-3548.PubMedGoogle Scholar
- Jennings B, de Celis J, Delidakis C, Preiss A, Bray S: Role of Notch and achaete-scute complex in the expression of Enhancer of split bHLH proteins. Development. 1995, 121: 3745-3752.Google Scholar
- Lai EC, Bodner R, Posakony JW: The Enhancer of split Complex of Drosophila includes four Notch-regulated members of the Bearded family. Development. 2000, 127: 3441-3455.PubMedGoogle Scholar
- Wurmbach E, Wech I, Preiss A: The Enhancer of split complex of Drosophila melanogaster harbors three classes of Notch responsive genes. Mech Dev. 1999, 80: 171-180. 10.1016/S0925-4773(98)00212-3.View ArticlePubMedGoogle Scholar
- Delidakis C, Preiss A, Hartley DA, Artavanis-Tsakonas S: Two genetically and molecularly distinct functions involved in early neurogenesis reside within the Enhancer of split locus of Drosophila melanogaster. Genetics. 1991, 129: 803-823.PubMed CentralPubMedGoogle Scholar
- Nagel AC, Maier D, Krauss S, Mezger M, Preiss A: Neurogenic phenotypes induced by RNA interference with bHLH genes of the Enhancer of split complex of Drosophila melanogaster. Genesis. 2004, 39: 105-114. 10.1002/gene.20033.View ArticlePubMedGoogle Scholar
- Schrons H, Knust E, Campos-Ortega JA: The Enhancer of split complex and adjacent genes in the 96F region of Drosophila melanogaster are required for segregation of neural and epidermal progenitor cells. Genetics. 1992, 132: 481-503.PubMed CentralPubMedGoogle Scholar
- Cooper MT, Tyler DM, Furriols M, Chalkiadaki A, Delidakis C, Bray S: Spatially restricted factors cooperate with Notch in the regulation of Enhancer of split genes. Dev Biol. 2000, 221: 390-403. 10.1006/dbio.2000.9691.View ArticlePubMedGoogle Scholar
- de Celis JF, de Celis J, Ligoxygakis P, Preiss A, Delidakis C, Bray S: Functional relationships between Notch, Su(H) and the bHLH genes of the E(spl) complex: the E(spl) genes mediate only a subset of Notch activities during imaginal development. Development. 1996, 122: 2719-2728.PubMedGoogle Scholar
- Tata F, Hartley DA: Inhibition of cell fate in Drosophila by Enhancer of split genes. Mech Dev. 1995, 51: 305-315. 10.1016/0925-4773(95)00377-0.View ArticlePubMedGoogle Scholar
- Wech I, Bray S, Delidakis C, Preiss A: Distinct expression patterns of different Enhancer of split bHLH genes during Drosophila embryogenesis. Dev Genes Evol. 1999, 209: 370-375. 10.1007/s004270050266.View ArticlePubMedGoogle Scholar
- Maier D, Marte BM, Schäfer W, Yu Y, Preiss A: Drosophila evolution challenges postulated redundancy in the E(spl) gene complex. Proc Natl Acad Sci USA. 1993, 90: 5464-5468.PubMed CentralView ArticlePubMedGoogle Scholar
- Fisher AL, Caudy M: Groucho proteins: transcriptional corepressors for specific subsets of DNA-binding transcription factors in vertebrates and invertebrates. Genes Dev. 1998, 12: 1931-1940.View ArticlePubMedGoogle Scholar
- Fisher A, Caudy M: The function of hairy-related bHLH repressor proteins in cell fate decisions. Bioassays. 1998, 20: 298-306. 10.1002/(SICI)1521-1878(199804)20:4<298::AID-BIES6>3.0.CO;2-M.View ArticleGoogle Scholar
- Artavanis-Tsakonas S, Simpson P: Choosing a cell fate: a view from the Notch locus. Trends Genet. 1991, 7: 403-8. 10.1016/0168-9525(91)90220-K.View ArticlePubMedGoogle Scholar
- Campuzano S, Modolell J: Patterning of the Drosophila nervous system: the achaete-scute gene complex. Trends Genet. 1992, 8: 202-208.View ArticlePubMedGoogle Scholar
- Cubas P, de Celis J, Campuzano S, Modolell J: Proneural clusters of achaete-scute expression and the generation of sensory organs in the Drosophila imaginal wing disc. Genes Dev. 1991, 5: 996-1008.View ArticlePubMedGoogle Scholar
- Heitzler P, Bourois M, Ruel L, Cateret C, Simpson P: Genes of the Enhancer of split and achaete-scute complexes are required for a regulatory loop between Notch and Delta during lateral signalling in Drosophila. Development. 1996, 122: 161-171.PubMedGoogle Scholar
- Bertrand N, Castro D, Guillemot F: Proneural genes and the specification of neural cell fates. Nature Rev Neurosci. 2002, 3: 517-530. 10.1038/nrn874.View ArticleGoogle Scholar
- Flybase, A Database of the Drosophila Genom. [http://flybase.bio.indiana.edu/]
- Gerhart J, Kirschner M: Arthropod phylogeny. Cells, embryos, and evolution. 1997, Blackwell Science. Malden, 305-Google Scholar
- Lai EC, Tam B, Rubin GM: Pervasive regulation of Drosophila Notch target genes by GY-box, Brd-box, and K-box-class microRNAs. Genes Dev. 2005, 19: 1067-1080. 10.1101/gad.1291905.PubMed CentralView ArticlePubMedGoogle Scholar
- Moore AD, Barbel S, Jan LY, Jan YN: A genomewide survey of basic helix-loop-helix factors in Drosophila. Proc Natl Acad Sci USA. 2000, 97: 10436-10441. 10.1073/pnas.170301897.PubMed CentralView ArticlePubMedGoogle Scholar
- Prokopenko SN, He Y, Lu Y, Bellen HJ: Mutations affecting the development of the peripheral nervous system in Drosophila: A molecular screen for novel proteins. Genetics. 2000, 156: 1691-1715.PubMed CentralPubMedGoogle Scholar
- Giagtzoglou N, Koumbanakis KA, Fullard J, Zarifi I, Delidakis C: Role of the Sc C-terminus in transcription activation and E(spl) repressor recruitment. J Biol Chem. 2005, 280: 1299-1305. 10.1074/jbc.M408949200.View ArticlePubMedGoogle Scholar
- Wülbeck C, Simpson P: The expression of pannier and achaete-scute homologues in a mosquito suggests an ancient role of pannier as a selector gene in the regulation of the dorsal body plan. Development. 2002, 129: 3861-3871.PubMedGoogle Scholar
- Burge C, Karlin S: Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997, 268: 78-94. 10.1006/jmbi.1997.0951.View ArticlePubMedGoogle Scholar
- Giagtzoglou N, Alifragis P, Koumbanakis KA, Delidakis C: Two modes of recruitment of E(spl) repressors onto target genes. Development. 2003, 130: 259-270. 10.1242/dev.00206.View ArticlePubMedGoogle Scholar
- Gibert J-M, Simpson P: Evolution of cis-regulation of the proneural genes. Int J Dev Biol. 2003, 47: 643-651.PubMedGoogle Scholar
- Massari ME, Murre C: Helix-Loop-Helix proteins: regulators of transcription in eucaryotic organisms. Mol Cell Biol. 2000, 20: 429-440. 10.1128/MCB.20.2.429-440.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Oellers N, Dehio M, Knust E: bHLH proteins encoded by the Enhancer of split complex of Drosophila negatively interfere with transcriptional activation mediated by proneural genes. Mol Gen Genet. 1994, 244: 465-473. 10.1007/BF00583897.View ArticlePubMedGoogle Scholar
- Lai EC, Posakony JW: Regulation of Drosophila neurogenesis by RNA:RNA duplexes?. Cell. 1998, 93: 1103-1104. 10.1016/S0092-8674(00)81454-3.View ArticlePubMedGoogle Scholar
- Castro B, Barolo S, Bailey AM, Posakony JM: Lateral inhibition in proneural clusters: cis-regulatory logic and default repression by Suppressor of Hairless. Development. 2005, 132: 3333-3344. 10.1242/dev.01920.View ArticlePubMedGoogle Scholar
- Carroll SB, Grenier JK, Weatherbee SD: From DNA to diversity: molecular genetics and the evolution of animal design. 2001, Blackwell Science. Malden, 46: 148-153.Google Scholar
- McGinnis W, Garber RL, Wirz J, Kuroiwa A, Gehring WJ: A homologous protein-coding sequence in Drosophila homeotic genes and its conservation in other metazoans. Cell. 1984, 37: 403-408. 10.1016/0092-8674(84)90370-2.View ArticlePubMedGoogle Scholar
- Grens A, Mason E, Marsh JL, Bode HR: Evolutionary conservation of a cell fate specification gene: the Hydra achaete scute homolog has proneural activity in Drosophila. Development. 1995, 121: 4027-4035.PubMedGoogle Scholar
- Alonso MC, Cabrera CV: The achaete-scute gene complex of Drosophila melanogaster comprises four homologous genes. EMBO J. 1988, 7: 2585-2591.PubMed CentralPubMedGoogle Scholar
- Ghysen A, Dambly-Chaudiere C: From DNA to form: the achaete-scute complex. Genes Dev. 1988, 2: 495-501.View ArticlePubMedGoogle Scholar
- Gonzalez F, Romani S, Cubas P, Modolell J, Campuzano S: Molecular analysis of the asense gene, a member of the achaete-scute complex of Drosophila melanogaster, and its novel role in optic lobe development. EMBO J. 1989, 8: 3553-3562.PubMed CentralPubMedGoogle Scholar
- Pistillo D, Skaer N, Simpson P: scute expression in Calliphora vicina reveals an ancestral pattern of longitudinal stripes on the thorax of higher Diptera. Development. 2002, 129: 563-572.PubMedGoogle Scholar
- Skeath JB, Carroll SB: Regulation of the achaete-scute gene expression and sensory organ pattern formation in the Drosophila wing. Genes Dev. 1991, 5: 984-995.View ArticlePubMedGoogle Scholar
- Villares R, Cabrera CV: The achaete-scute gene complex of D. melanogaster: conserved domains in a subset of genes required for neurogenesis and their homology to myc. Cell. 1987, 50: 415-424. 10.1016/0092-8674(87)90495-8.View ArticlePubMedGoogle Scholar
- Wheeler SR, Skeath JB: The identification and expression of achaete-scute genes in the branchiopod crustacean Triops longicaudatus. Gene Expr Patterns. 2005, 5: 695-700. 10.1016/j.modgep.2005.02.005.View ArticlePubMedGoogle Scholar
- Garcia-Garcia MJ, Ramain P, Simpson P, Modolell J: Different contributions of pannier and wingless to the patterning of the dorsal mesothorax of Drosophila. Development. 1999, 126: 3523-3532.PubMedGoogle Scholar
- Gomez-Skarmeta JL, Diez del Corral R, de la Calle-Mustienes E, Ferre-Marco D, Modolell J: Araucan and caupolican, two members of the novel iroquois complex, encode homeoproteins that control proneural and vein-forming genes. Cell. 1996, 85: 95-105. 10.1016/S0092-8674(00)81085-5.View ArticlePubMedGoogle Scholar
- Simpson P, Woehl R, Usui K: The development and evolution of bristle patterns in Diptera. Development. 1999, 126: 1349-1364.PubMedGoogle Scholar
- Bächli G, Bural H: Taxonomische Merkmale von Drosophila-Arten. Diptera Drosophilidae. 1985, Schweizerische Entomologische Gesellschaft, Zürich: ETH-Zentrum, 20-30.Google Scholar
- Human Genome Sequencing Center. [http://www.hgsc.bcm.tmc.edu/]
- D. virilis (Drosophila virilis) Genome Browser Gateway. [http://genome.ucsc.edu/cgi-bin/hgGateway?db=droVir1]
- NCBI, Trace Archive database Mega BLAST search. [http://www.ncbi.nlm.nih.gov./blast/mmtrace.shtml]
- Ensembl genome browser. [http://www.ensembl.org/]
- Honeybee brain EST project. [http://titan.biotec.uiuc.edu/bee/honeybee_project.htm]
- HUSAR – Sequence Analysis Package. [http://genius.embnet.dkfz-heidelberg.de/menu/w2h/w2hdkfz/]
- Senger M, Flores T, Glatting K-H, Ernst P, Hotz-Wagenblatt A, Suhai S: W2H: WWW interface to the GCG sequence analysis package. Bioinformatics. 1998, 14: 452-457. 10.1093/bioinformatics/14.5.452.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.