Daphnia Halloween genes that encode cytochrome P450s mediating the synthesis of the arthropod molting hormone: Evolutionary implications
© Rewitz and Gilbert; licensee BioMed Central Ltd. 2008
Received: 15 November 2007
Accepted: 25 February 2008
Published: 25 February 2008
In crustaceans and insects, development and reproduction are controlled by the steroid hormone, 20-hydroxyecdysone (20E). Like other steroids, 20E, is synthesized from cholesterol through reactions involving cytochrome P450s (CYPs). In insects, the CYP enzymes mediating 20E biosynthesis have been identified, but evidence of their probable presence in crustaceans is indirect, relying solely on the ability of crustaceans to synthesize 20E.
To investigate the presence of these genes in crustaceans, the genome of Daphnia pulex was examined for orthologs of these genes, the Halloween genes, encoding those biosynthetic CYP enzymes. Single homologs of spook-CYP307A1, phantom-CYP306A1, disembodied-CYP302A1, shadow-CYP315A1 and shade-CYP314A1 were identified in the Daphnia data base. Phylogenetic analysis indicates an orthologous relationship between the insect and Daphnia genes. Conserved intron/exon structures and microsynteny further support the conclusion that these steroidogenic CYPs have been conserved in insects and crustaceans through some 400 million years of evolution.
Although these arthropod steroidogenic CYPs are related to steroidogenic CYPs in Caenorhabditis elegans and vertebrates, the data suggest that the arthropod steroidogenic CYPs became functionally specialized in a common ancestor of arthropods and are unique to these animals.
Steroid hormones, regulate essential processes during development and reproduction, and are synthesized from cholesterol under the control of steroidogenic enzymes in the cytochrome P450 (CYP) family . In Caenorhabditis elegans, insects and vertebrates, different steroids are produced to control developmental processes, suggesting that steroidogenic CYPs evolved and became functionally specialized in different lineages during evolution. In insects, a specific biosynthetic pathway yielding 20-hydroxyecdysone (20E), the arthropod molting hormone, evolved, whereas in the line leading to vertebrates, biosynthetic CYPs that produce the vertebrate-type steroids evolved . Since there is some evidence of the presence of vertebrate-type sex steroids in invertebrates such as echinoderms and mollusks, although no unequivocal evidence that they can synthesize these steroids , the possibility remains that CYPs with the capacity to produce vertebrate-type sex steroids were present in the common ancestor even before the protostome-deuterostome split. Thus, the evolution of steroidogenic CYPs is still an open question.
Although one would expect that orthologs of the insect Halloween genes are present in crustaceans, there is no molecular evidence for the existence of these genes in crustaceans. We have tried for several years to probe hexapod crustaceans for Halloween gene orthologs under various hormonal regimens using degenerate primers based on the Drosophila and Bombyx genes, but we have been unsuccessful (K.F. Rewitz, J.T. Warren, E. Chang and L.I. Gilbert). The development of the genome data base of the more primitive crustacean, Daphnia pulex, allowed us to survey this genome and conduct phylogenetic analyses that suggest strongly that orthologs of spo, phm, dib, sad and shd do exist in Daphnia and thus, in a crustacean i.e. the genes appeared in arthropods before the radiation of insects.
Results and Discussion
By searching the Daphnia data bases [18, 19] we obtained candidate sequences for orthologs of the insect Halloween genes in Daphnia. Single orthologs of phm, dib, sad and shd were retrieved and only one sequence exhibited significant similarity to the spo-like genes in the CYP307 family. We also searched the genomes of non-arthropod invertebrates including the cnidarian Nematostella vectensis, the nematodes C. elegans and Brugia malayi, the annelid Capitella capitata, the mollusk Lottia gigantea, the echinoderm Strongylocentrotus purpuratus [19, 20] for Halloween gene orthologs using TBLASTN. Analyses of these invertebrate genomes did not result in any significant hits indicating the absence of Halloween orthologs and neither have we been able to identify orthologs of these genes in any of several vertebrate species.
The genes obtained from Daphnia encode approximately 500 amino acid open reading frames (ORFs) which is typically for proteins belonging to the CYP family . Alignment of Daphnia sequences with orthologs of the insect Halloween genes shows that the genes, in addition to being conserved in areas that comprise canonical structural CYP motifs, exhibit considerable conservation in regions that are believed to determine substrate specificity (Additional file 1). This indicates that the genes are functionally conserved in Daphnia as they are in insects . The overall amino acid identity between deduced orthologous proteins from eight insects belonging to four different orders and the Daphnia orthologs ranges from an average of 55.2% ± 6.8 SD (standard deviation) for Spo-like proteins to somewhat lower values for Sad proteins (38.7% ± 8.8 SD). Thus, spo-like genes are the most highly conserved of these genes and this appears to be true for Daphnia like it is for insects. Conservation of these genes from insect to crustaceans, which separated ~400 millions years ago , shows that selection has preserved the genes because of their function. The reason that spo-like genes are more conserved than the other arthropod steroidogenic CYP enzymes is not known, although it may be related to the possibility that Spo acts in the rate-limiting Black Box reaction(s) [5, 22]. If Spo is involved in control of the flux through the pathway it may have been a particular target for selection because mutations altering its enzyme activity would have had increased consequences compared to mutations altering the activity of enzymes that are less rate-limiting. A phylogenetic analysis was performed using the Daphnia sequences.
Phylogenetic analysis of Daphnia candidate orthologs
A phylogenetic tree was constructed with Daphnia sequences retrieved from the BLAST searches and orthologs of the insect Halloween genes from Drosophila, red flour beetle Tribolium castaneum and honey bee Apis mellifera, which represent three different orders (Diptera, Coleoptera and Hymenoptera, respectively) of holometabolous insects. The Daphnia gene products separate with the insect orthologs of Spo/Spok, Phm, Dib, Sad and Shd in this phylogenetic tree with high bootstrap support (Fig. 1). In examining the genome of Daphnia for orthologs of these genes, we noticed no closely related paralogs to the obtained candidate orthologs. This supports the orthology of the genes and indicates that the genes became functionally specialized before the split between crustaceans and insects and have been under heavy selection pressure ever since. Only one spo-like gene was obtained and this gene is phylogenetically a CYP307A subfamily gene i.e. most closely related to spo/spok. No ortholog of spookiest (spot: CYP307B1) was found indicating that this gene does not exist in Daphnia, which implies that it is insect-specific i.e. originating from a duplication occurring after insects arose from crustaceans  or that it was lost in Daphnia as it has been in some insects e.g. lepidopterans and Drosophila species .
In the phylogenetic analysis, sequences of steroidogenic CYPs from vertebrates and C. elegans were included to probe ancestral relationships and the origin of CYPs involved in steroid biogenesis (Fig. 1). Figure 1 includes major groups of metazoans CYPs i.e. mitochondrial, CYP2-related, CYP3-related, CYP4-related, CYP19s and CYP51s. The steroidogenic CYPs of insects and their Daphnia orthologs are evolutionarily related to vertebrate and C. elegans steroidogenic CYPs since they cluster in two major groups, those related to CYP2 enzymes (Spo, Phm, CYP17, CYP21 and DAF-9) and those that are mitochondrial (Dib, Sad, Shd and CYP11A1). Therefore, it is likely that different steroidogenic CYP enzymes are derived from common ancestors and were recruited for steroid biosynthesis prior to the protostome-deuterostome split. Evolution of these ancestral CYPs in C. elegans, arthropods and vertebrates, by gene duplication and divergence, likely shaped the biosynthetic pathways yielding the different groups of steroids in these animals.
Conservation of gene structure
Although several introns are conserved among the mitochondrial dib, sad and shd, introns that are unique to each group of orthologs (e.g. the first phase 1 intron in dib and shd) indicate an orthologous relationship to the respective groups. The conserved gene structure of these genes suggests that they evolved from a common ancestor of mitochondrial CYPs by duplications and diverged into the steroidogenic enzymes that are functionally specialized to carry out the last three steps in which C22, C2 and C20 hydroxylations form 20E. This implies, at least in part, that the biosynthetic pathway evolved step-wise by duplication and functional divergence of the biosynthetic enzymes. In vertebrates, functional divergence and specialization of steroidogenic CYPs are known from studies of the CYP11 paralogs. CYP11A1 is the enzyme that mediates the removal of the side chain of cholesterol . CYP11B1 and CYP11B2 are paralogs, from a recent duplication, that became functionally specialized as the steroid 11β-hydroxylase and the aldosterone synthase, respectively . dib, sad and shd also exhibit conservation of gene structure when compared to vertebrate mitochondrial CYPs involved in steroidogenesis .
Microsynteny supports orthology
The conserved arrangement of genes on the chromosome of different species (microsynteny) can be used to infer phylogenetic relationships that support the orthology of genes [2, 17]. Analysis of local genome structures surrounding the Halloween genes shows little conservation of microsynteny between insects and crustaceans. This is likely due to the relatively great evolutionary time separating insects from crustaceans. However, remnants showing microsynteny can be found in at least two cases. CYP18A1 is a paralog of phm which is found in Drosophila and most other insects . In Drosophila and Apis, phm and CYP18A1 are arranged tail-to-tail and adjacent to these genes are CG6696 and fused (Fig. 2). In Daphnia, the microsyntenic relationship of phm and CYP18A1 is conserved, although CG6696 and fused appear to have been rearranged to different chromosomal locations.
Comparison of gene organization in a region of approximately 100 kb surrounding the spo(k) locus reveals that spo is adjacent to neverland (nvd)  in Daphnia and the mosquito Anopheles gambiae. Analysis of the genome region of spo-like genes in species of Drosophila (available at FlyBase ) shows that in Drosophila pseudoobscura, Drosophila willistoni and Drosophila mojavensis nvd is located within 60 kb of spok. This is in agreement with the view that spok is the ancestral gene in Drosophilidae. Intriguingly, Nvd is a conserved Rieske-like oxygenase believed to be involved in the conversion of cholesterol to 7-dehydrocholesterol in insects and C. elegans [26, 28]. In insects, this is presumably the first step in the biosynthetic pathway of 20E from cholesterol via 7-dehydrocholesterol, the latter being the substrate that enters the Black Box, in which Spo/Spok likely participate (although at different developmental stages) to produce the substrate for Phm, the ketodiol (Fig. 1). Although unequivocal evidence for the exact function of Nvd and Spo/Spok is lacking, experiments utilizing ecdysteroid precursors to rescue mutant phenotypes (RNAi-induced for nvd and spok) suggest that in Drosophila nvd and spo/spok genes function in these two succeeding steps in the biochemical pathway [16, 26]. However, it must be emphasized that the Black Box may hold more than one reaction and more than one enzyme may be participating in this conversion .
Currently, it is clear that CYPs are involved in steroid biosynthesis in vertebrates and invertebrates. In C. elegans, insects and vertebrates, the majority of steroidogenic CYPs are related to one of two groups, the CYP2 enzymes and the mitochondrial CYPs. It is therefore likely that these steroidogenic CYPs evolved from common ancestors which had the capacity to modify cholesterol before the split of these animals. This implies that the distinct biosynthetic pathways in these metazoans, which convert cholesterol to different types of steroids, likely result from lineage-specific evolution of ancestral steroidogenic CYPs. The results suggest strongly that steroidogenic CYP enzymes mediating the biosynthesis of 20E are present in crustaceans and are relegated to the Arthropoda. Therefore, these enzymes probably became functionally specialized in the line leading to arthropods. Although functional characterization of the gene products, such as that done for the insect Halloween genes [1, 29], is required to obtain unequivocal evidence of their role in 20E production in Daphnia, the present data provide the first evidence for steroidogenic CYP genes in crustaceans as well as the groundwork for future functional genomic analyses in the field of crustacean endocrinology.
Data base search
The Daphnia pulex genome sequence (v1.0) made available by the Daphnia Genomics Consortium and the DOE Joint Genome Institute  was mined for orthologs of the insect Halloween genes using TBLASTN searches. Apis and Tribolium sequences, previously described in Rewitz et al. , were used as queries to identify homologous sequences in the Daphnia data base. Sequences of insect Halloween genes, previously described in Rewitz et al. , were gleaned from FlyBase  and Rene Feyereisen's website for insect CYPs . C. elegans and vertebrate CYP protein sequences were acquired from NCBI .
Phylogenetic tree construction, intron positions and microsynteny analysis
Predicted Daphnia orthologs of the insect Halloween genes were analyzed using a phylogenetic tree with insect Halloween genes from insects belonging to three different orders, Apis (Hymenoptera), Tribolium (Coleoptera) and Drosophila (Diptera). To infer phylogenetic relationships of these steroidogenic CYPs with other members of this multigene family, steroidogenic and non-steroidogenic CYPs representing some of the major groups of metazoan CYPs were included in this analysis i.e. mitochondrial, CYP2-related, CYP3-related and CYP4s, CYP19 and CYP51. Deduced Daphnia protein sequences used in the alignment can be found in Additional file 2. A multiple alignment of the protein sequences was constructed with ClustalX (1.83)  and manually edited using SeaView . A phylogenetic tree was constructed from this alignment using the Maximum Likelihood method under the Jones-Taylor-Thornton (JTT) substitution model using PHYML (v2.4.5) . Support values were obtained by bootstrapping 100 replicates. Branches with bootstrap support below 50 were collapsed to form polytomies.
Gene structure, that is, intron position and phase (the nucleotide position of the intron within a codon: phase 0 between codons, phase 1 after the first base and phase 2 after the second base) were predicted from manual annotation of the Daphnia genes with support from expressed sequence tags (ESTs) data and by homology to insect Halloween genes for which gene structures are known. Intron position and phase were mapped onto a multiple alignment to investigate conservation of intron/exon structure in relation to phylogeny.
Additional support for orthology of Daphnia and insect genes was sought by analyses of microsynteny (conserved order of genes). Genome regions surrounding the candidate Halloween orthologs of Daphnia and several insects were inspected for ORFs encoding conserved genes. Putative gene orthologs exhibiting conserved arrangement in relation to the Halloween orthologs were tested for orthology by reciprocal best hit Blast searches in which orthologs of each species were BLAST searched against the genome of the other species to look for best hits (putative orthologs).
(open reading frame)
We offer special thanks to John Colbourne and Joseph Shaw for their help in gaining access to the Daphnia data base and Michael B. O'Connor and James T. Warren for their continued collaboration on the Halloween gene project. The Daphnia sequence data were produced by the US Department of Energy Joint Genome Institute  in collaboration with the Daphnia Genomics Consortium . Our work benefits from, and contributes to, the Daphnia Genomics Consortium. The present work was supported by grant 0516623 from the National Science Foundation.
- Gilbert LI, Warren JT: A molecular genetic approach to the biosynthesis of the insect steroid molting hormone. Vitam Horm. 2005, 73: 31-57.View ArticlePubMedGoogle Scholar
- Rewitz KF, O'Connor MB, Gilbert LI: Molecular evolution of the insect Halloween family of cytochrome P450s: phylogeny, gene organization and functional conservation. Insect Biochem Mol Biol. 2007, 37: 741-753. 10.1016/j.ibmb.2007.02.012.View ArticlePubMedGoogle Scholar
- Lafont R, Mathieu M: Steroids in aquatic invertebrates. Ecotoxicology. 2007, 16: 109-130. 10.1007/s10646-006-0113-1.View ArticlePubMedGoogle Scholar
- Glenner H, Thomsen PF, Hebsgaard MB, Sørensen MV, Willerslev E: The origin of insects. Science. 2006, 314: 1883-1884. 10.1126/science.1129844.View ArticlePubMedGoogle Scholar
- Gilbert LI, Rybczynski R, Warren JT: Control and biochemical nature of the ecdysteroidogenic pathway. Ann Rev Entomol. 2002, 47: 883-916. 10.1146/annurev.ento.47.091201.145302.View ArticleGoogle Scholar
- Lachaise F, Le Roux A, Hubert M, Lafont R: The molting gland of crustaceans: localization, activity, and endocrine control. J Crustacean Biol. 1993, 13: 198-234. 10.2307/1548970.View ArticleGoogle Scholar
- Chavez VM, Marques G, Delbecque JP, Kobayashi K, Hollingsworth M, Burr J, Natzle JE, O'Connor MB: The Drosophila disembodied gene controls late embryonic morphogenesis and codes for a cytochrome P450 that regulates embryonic ecdysone levels. Development. 2000, 127 (19): 4115-4126.PubMedGoogle Scholar
- Petryk A, Warren JT, Marques G, Jarcho MP, Gilbert LI, Parvy J-P, Dauphin-Villemant C, O'Connor MB: Shade is the Drosophila P450 enzyme that mediates the hydroxylation of ecdysone to the steroid insect molting hormone 20-hydroxyecdysone. Proc Natl Acad Sci USA. 2003, 100: 13773-13778. 10.1073/pnas.2336088100.PubMed CentralView ArticlePubMedGoogle Scholar
- Niwa R, Matsuda T, Yoshiyama T, Namiki T, Mita K, Fujimoto Y, Kataoka H: CYP306A1, a cytochrome P450 enzyme is essential for ecdysteroid biosynthesis in the prothoracic glands of Bombyx and Drosophila. J Biol Chem. 2004, 279: 35942-35949. 10.1074/jbc.M404514200.View ArticlePubMedGoogle Scholar
- Niwa R, Sakudoh T, Namiki T, Saida K, Fujimoto Y, Kataoka H: The ecdysteroidogenic CYP CYP302A1/disembodied from the silkworm, Bombyx mori, is transcriptionally regulated by prothoracicotropic hormone. Insect Mol Biol. 2005, 14: 563-571. 10.1111/j.1365-2583.2005.00587.x.View ArticlePubMedGoogle Scholar
- Rewitz KF, Rybczynski R, Warren JT, Gilbert LI: Identification, characterization and developmental expression of Halloween genes encoding P450 enzymes mediating ecdysone biosynthesis in the tobacco hornworm, Manduca sexta. Insect Biochem Mol Biol. 2006, 36: 188-199. 10.1016/j.ibmb.2005.12.002.View ArticlePubMedGoogle Scholar
- Rewitz KF, Rybczynski R, Warren JT, Gilbert LI: Developmental expression of Manduca shade, the P450 mediating the final step in molting hormone synthesis. Mol Cell Endocrinol. 2006, 247: 166-174. 10.1016/j.mce.2005.12.053.View ArticlePubMedGoogle Scholar
- Warren JT, Petryk A, Marques G, Jarcho M, Parvy J-P, Dauphin-Villemaont C, O'Connor MB, Gilbert LI: Molecular and biochemical characterization of two P450 enzymes in the ecdysteroidogenic pathway of Drosophila melanogaster. Proc Natl Acad Sci USA. 2002, 99: 11043-11048. 10.1073/pnas.162375799.PubMed CentralView ArticlePubMedGoogle Scholar
- Warren JT, Petryk A, Marques G, Jarcho M, Parvy JP, Shinoda T, Itoyama K, Kobayashi J, Jarcho M, Li Y, O'Connor MB, Dauphin-Villemant C, Gilbert LI: Phantom encodes the 25-hydroxylase of Drosophila melanogaster and Bombyx mori: a P450 enzyme critical in ecdysone biosynthesis. Insect Biochem Mol Biol. 2004, 34: 991-1010. 10.1016/j.ibmb.2004.06.009.View ArticlePubMedGoogle Scholar
- Namiki T, Niwa R, Sakudoh T, Shirai K, Takeuchi H, Kataoka H: Cytochrome P450 CYP307A1/Spook: a regulator for ecdysone synthesis in insects. Biochem Biophys Res Comm. 2005, 337: 367-374. 10.1016/j.bbrc.2005.09.043.View ArticlePubMedGoogle Scholar
- Ono H, Rewitz KF, Shinoda T, Itoyama K, Petryk A, Rybczynski R, Jarcho M, Warren JT, Marques G, Shimell MJ, Gilbert LI, O'Connor MB: Spook and spookier code for stage-specific components of the biosynthetic pathway in Diptera. Dev Bio. 2006, 298 (2): 555-5780. 10.1016/j.ydbio.2006.07.023. Epub 2006 Jul 29.View ArticleGoogle Scholar
- Sztal T, Chung H, Gramzow L, Daborn PJ, Batterham P, Robin C: Two independent duplications forming the Cyp307a genes in Drosophila. Insect Biochem Mol Biol. 2007, 37: 1044-1053. 10.1016/j.ibmb.2007.05.017.View ArticlePubMedGoogle Scholar
- Daphnia genome database. [http://wFleaBase.org]
- DOE Joint Genome Institute. [http://jgi.doe.gov/]
- NCBI. [http://ncbi.nlm.nih.gov/]
- Feyereisen R: Insect cytochrome P450. Comprehensive Molecular Insect Science. Edited by: Gilbert LI, Iatrou K. 2005, Gill S: Elsevier; Oxford, 4: 1-77.View ArticleGoogle Scholar
- Lafont R, Dauphin-Villemant C, Warren JT, Rees H: Ecdysteroid chemistry and biochemistry. Comprehensive Molecular Insect Science. Edited by: Gilbert LI, Iatrou K. 2005, Gill S: Elsevier; Oxford, 3: 125-195.View ArticleGoogle Scholar
- Lisurek M, Bernhardt R: Modulation of aldosterone and cortisol synthesis on the molecular level. Mol Cell Endocrinol. 2004, 215: 149-159. 10.1016/j.mce.2003.11.008.View ArticlePubMedGoogle Scholar
- Bülow HE, Bernhardt R: Analyses of the CYP11B gene family in the guinea pig suggest the existence of a primordial CYP11B gene with aldosterone synthase activity. Eur J Biochem. 2002, 269: 3838-3846. 10.1046/j.1432-1033.2002.03076.x.View ArticlePubMedGoogle Scholar
- Claudianos C, Ranson H, Johnson RM, Biswas S, Schuler MA, Berenbaum MR, Feyereisen R, Oakeshott JG: A deficit of detoxification enzymes: pesticide sensitivity and environmental response in the honeybee. Insect Mol Biol. 2006, 15: 615-636. 10.1111/j.1365-2583.2006.00672.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoshiyama T, Namiki T, Mita K, Kataoka H, Niwa R: Neverland is an evolutionarily conserved Rieske-domain protein that is essential for ecdysone synthesis and insect growth. Development. 2006, 133: 565-2574. 10.1242/dev.02428.View ArticleGoogle Scholar
- FlyBase. [http://flybase.bio.indiana.edu/]
- Rottiers V, Motola DL, Gerisch B, Cummins CL, Nishiwaki K, Mangelsdorf DJ, Antebi1 A: Hormonal control of C. elegans Dauer formation and life span by a Rieske-like oxygenase. Dev Cell. 2006, 10: 473-482.View ArticlePubMedGoogle Scholar
- Rewitz KF, Rybczynski R, Warren JT, Gilbert LI: The Halloween genes code for cytochrome P450 enzymes mediating synthesis of the insect moulting hormone. Biochem Soc Trans. 2006, 34: 1256-1260. 10.1042/BST0341256.View ArticlePubMedGoogle Scholar
- Rene Feyereisen's website for insect CYPs. [http://p450.antibes.inra.fr/]
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins D: The CLUSTALX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMedGoogle Scholar
- Galtier N, Gouy M, Gautier C: SEAVIEW and PHYLO_WIN: two graphic tools for sequence alignment and molecular phylogeny. Comput Appl Biosci. 1996, 12: 543-548.PubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Daphnia Genomics Consortium. [http://daphnia.cgb.indiana.edu/]
- Motola DL, Cummins CL, Rottiers V, Sharma KK, Li T, Li Y, Suino-Powell K, Xu HE, Auchus RJ, Antebi A, Mangelsdorf1 DJ: Identification of ligands for DAF-12 that govern Dauer formation and reproduction in C. elegans. Cell. 2006, 124: 1209-1223. 10.1016/j.cell.2006.01.037.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.