Cytochrome P450 CYP307A1/Spook: a regulator for ecdysone synthesis in insects

The prothoracic gland (PG) has essential roles in synthesizing and secreting a steroid hormone called ecdysone that is critical for molting and metamorphosis of insects. However, little is known about the genes controlling ecdysteroidogenesis in the PG. To identify genes functioning in the PG of the silkworm, Bombyx mori, we used differential display PCR and focused on a cytochrome P450 gene designated Cyp307a1. Its expression level positively correlates with a change in the hemolymph ecdysteroid titer. In addition, Drosophila Cyp307a1 is encoded in the spook locus, one of the Halloween mutant family members showing a low ecdysone titer in vivo, suggesting that Cyp307a1 is involved in ecdysone synthesis. While Drosophila Cyp307a1 is expressed in the early embryos and adult ovaries, the expression is not observed in the PGs of embryos or third instar larvae. These results suggest a difference in the ecdysone synthesis pathways during larval development in these insects.     

Spook and Spookier code for stage-specific components of the ecdysone biosynthetic pathway in Diptera

Ecdysteroids regulate many key developmental events in arthropods including molting and metamorphosis. Recently, members of the Drosophila Halloween group of genes, that are required for embryonic viability and cuticle deposition, have been shown to code for several cytochrome P450 enzymes that catalyze the terminal hydroxylation steps in the conversion of cholesterol to the molting hormone 20-hydroxyecdysone. These P450s are conserved in other insects and each is thought to function throughout development as the sole mediator of a particular biosynthetic step since, where analyzed, each is expressed at all stages of development and shows no closely related homolog in their respective genomes. In contrast, we show here that several dipteran genomes encode two novel, highly related, microsomal P450 enzymes, Cyp307A1 and Cyp307A2, that likely participate as stage-specific components of the ecdysone biosynthetic machinery. This hypothesis comes from the observation that Cyp307A1 is encoded by the Halloween gene spook (spo), but unlike other Halloween class genes, Dmspo is not expressed during the larval stages. In contrast, Cyp307a2, dubbed spookier (spok), is expressed primarily during larval stages within the prothoracic gland cells of the ring gland. RNAi mediated reduction in the expression of this heterochromatin localized gene leads to arrest at the first instar stage which can be rescued by feeding the larva 20E, E or ketodiol but not 7dC. In addition, spok expression is eliminated in larvae carrying mutations in molting defective (mld), a gene encoding a nuclear zinc finger protein that is required for production of ecdysone during Drosophila larval development. Intriguingly, mld is not present in the Bombyx mori genome, and we have identified only one spook homolog in both Bombyx and Manduca that is expressed in both embryos and larva. These studies suggest an evolutionary split between Diptera and Lepidoptera in how the ecdysone biosynthetic pathway is regulated during development.     

Molecular evolution of the insect Halloween family of cytochrome P450s: phylogeny, gene organization and functional conservation

The insect molting hormone, 20-hydroxyecdysone (20E), is a major modulator of the developmental processes resulting in molting and metamorphosis. During evolution selective forces have preserved the Halloween genes encoding cytochrome P450 (P450) enzymes that mediate the biosynthesis of 20E. In the present study, we examine the phylogenetic relationships of these P450 genes in holometabolous insects belonging to the orders Hymenoptera, Coleoptera, Lepidoptera and Diptera. The analyzed insect genomes each contains single orthologs of Phantom (CYP306A1), Disembodied (CYP302A1), Shadow (CYP315A1) and Shade (CYP314A1), the terminal hydroxylases. In Drosophila melanogaster, the Halloween gene spook (Cyp307a1) is required for the biosynthesis of 20E, although a function has not yet been identified. Unlike the other Halloween genes, the ancestor of this gene evolved into three paralogs, all in the CYP307 family, through gene duplication. The genomic stability of these paralogs varies among species. Intron-exon structures indicate that D. melanogaster Cyp307a1 is a mRNA-derived paralog of spookier (Cyp307a2), which is the ancestral gene and the closest ortholog of the coleopteran, lepidopteran and mosquito CYP307A subfamily genes. Evolutionary links between the insect Halloween genes and vertebrate steroidogenic P450s suggest that they originated from common ancestors, perhaps destined for steroidogenesis, before the deuterostome-arthropod split. Conservation of putative substrate recognition sites of orthologous Halloween genes indicates selective constraint on these residues to prevent functional divergence. The results suggest that duplications of ancestral P450 genes that acquired novel functions may have been an important mechanism for evolving the ecdysteroidogenic pathway.     

A cytochrome p450 conserved in insects is involved in cuticle formation

The sequencing of numerous insect genomes has revealed dynamic changes in the number and identity of cytochrome P450 genes in different insects. In the evolutionary sense, the rapid birth and death of many P450 genes is observed, with only a small number of P450 genes showing orthology between insects with sequenced genomes. It is likely that these conserved P450s function in conserved pathways. In this study, we demonstrate the P450 gene, Cyp301a1, present in all insect genomes sequenced to date, affects the formation of the adult cuticle in Drosophila melanogaster. A Cyp301a1 piggyBac insertion mutant and RNAi of Cyp301a1 both show a similar cuticle malformation phenotype, which can be reduced by 20-hydroxyecdysone, suggesting that Cyp301a1 is an important gene involved in the formation of the adult cuticle and may be involved in ecdysone regulation in this tissue.     

Genome-wide analysis of cytochrome P450 monooxygenase genes in the silkworm, Bombyx mori

Based on the advances in the silkworm genome project, a new genome-wide analysis of cytochrome P450 genes was performed. A total of 84 CYP-related sequences were identified and could be classified into 26 families and 47 subfamilies according to standard nomenclature. Seventy eight of the eighty four genes appear to be functional and six are probable pseudogenes. The distribution of Bombyx mori P450s in the genome shows that most of them are tandem arranged on chromosomes, only 34 genes are present as singletons, with 8 clusters including 3 or more than 3 genes. Sequence alignments were used to reconstruct phylogenetic trees and to analyze the intron-exon organizations of the functional genes. The conserved intron positioning agrees perfectly with their common grouping on the tree. The presence of three extremely ancient introns which are conserved across different clans indicates that a few introns are still highly conserved after they have undergone extensive evolutionary changes of B. mori P450 duplication and divergence. Comparison of the P450s from B. mori to the P450s from Drosophila melanogaster shows that the expansion is not uniform across the gene families. Remarkably, two mitochondrial families, the B. mori CYP333 and D. melanogaster Cyp12, formed two orthologous groups in the phylogenetic tree. All CYP333s can be proposed to be related to xenobiotic metabolism in accordance with the D. melanogaster Cyp12s. The characterization and evolutionary analysis of P450s from B. mori in the current study provide useful information for understanding the characteristics and diversity of P450s from B. mori and the baseline for functional analyses of individual P450s in this model Lepidopteran insect.     

CYP306A1, a cytochrome P450 enzyme, is essential for ecdysteroid biosynthesis in the prothoracic glands of Bombyx and Drosophila

Ecdysteroids mediate a wide variety of developmental and physiological events in insects. In the postembryonic development of insects, ecdysone is synthesized in the prothoracic gland (PG). Although many studies have revealed the biochemical and physiological properties of the enzymes for ecdysteroid biosynthesis, most of the molecular identities of these enzymes have not been elucidated. Here we describe an uncharacterized cytochrome P450 gene, designated Cyp306a1, that is essential for ecdysteroid biosynthesis in the PGs of the silkworm Bombyx mori and fruit fly Drosophila melanogaster. Using the microarray technique for analyzing gene expression profiles in PG cells during Bombyx development, we identified two PG-specific P450 genes whose temporal expression patterns are correlated with changes in ecdysteroid titer during development. Amino acid sequence analysis showed that one of the Bombyx P450 genes belongs to the CYP306A1 subfamily. The temporal and spatial expression pattern of the Drosophila Cyp306a1 homolog is essentially the same as that of Bombyx Cyp306a1. We also found that Drosophila Cyp306a1 is disrupted in the phantom (phm) mutant, known also as the Halloween mutant. The morphological defects and decreased expression of ecdysone-inducible genes in phm suggest that this mutant cannot produce a high titer of ecdysone. Finally we demonstrate that S2 cells transfected with Cyp306a1 convert ketodiol to ketotriol via carbon 25 hydroxylation. These results strongly suggest that CYP306A1 functions as a carbon 25 hydroxylase and has an essential role in ecdysteroid biosynthesis during insect development.     

Drosophila genes that affect meiosis duration are among the meiosis related genes that are more often found duplicated

Using a phylogenetic approach, the examination of 33 meiosis/meiosis-related genes in 12 Drosophila species, revealed nine independent gene duplications, involving the genes cav, mre11, meiS332, polo and mtrm. Evidence is provided that at least eight out of the nine gene duplicates are functional. Therefore, the rate at which Drosophila meiosis/meiosis-related genes are duplicated and retained is estimated to be 0.0012 per gene per million years, a value that is similar to the average for all Drosophila genes. It should be noted that by using a phylogenetic approach the confounding effect of concerted evolution, that is known to lead to overestimation of the duplication and retention rate, is avoided. This is an important issue, since even in our moderate size sample, evidence for long-term concerted evolution (lasting for more than 30 million years) was found for the meiS332 gene pair in species of the Drosophila subgenus. Most striking, in contrast to theoretical expectations, is the finding that genes that encode proteins that must follow a close stoichiometric balance, such as polo, mtrm and meiS332 have been found duplicated. The duplicated genes may be examples of gene neofunctionalization. It is speculated that meiosis duration may be a trait that is under selection in Drosophila and that it has different optimal values in different species.     

Interspecific Comparison of the Transformer Gene of Drosophila Reveals an Unusually High Degree of Evolutionary Divergence

The transformer (tra) gene of Drosophila melanogaster occupies an intermediate position in the regulatory pathway controlling all aspects of somatic sexual differentiation. The female-specific expression of this gene's function is regulated by the Sex lethal (Sxl) gene, through a mechanism involving sex-specific alternative splicing of tra pre-mRNA. The tra gene encodes a protein that is thought to act in conjunction with the transformer-2 (tra-2) gene product to control the sex-specific processing of doublesex (dsx) pre-mRNA. The bifunctional dsx gene carries out opposite functions in the two sexes, repressing female differentiation in males and repressing male differentiation in females. Here we report the results from an evolutionary approach to investigate tra regulation and function, by isolating the tra-homologous genes from selected Drosophila species, and then using the interspecific DNA sequence comparisons to help identify regions of functional significance. The tra-homologous genes from two Sophophoran subgenus species, Drosophila simulans and Drosophila erecta, and two Drosophila subgenus species, Drosophila hydei and Drosophila virilis, were cloned, sequenced and compared to the D. melanogaster tra gene. This comparison reveals an unusually high degree of evolutionary divergence among the tra coding sequences. These studies also highlight a highly conserved sequence within intron one that probably defines a cis-acting regulator of the sex-specific alternative splicing event.     

Evolution of nucleotide substitutions and gene regulation in the amylase multigenes in Drosophila kikkawai and its sibling species

In order to determine evolutionary changes in gene regulation and the nucleotide substitution pattern in a multigene family, the amylase multigenes were characterized in Drosophila kikkawai and its sibling species. The nucleotide substitution pattern was investigated. Drosophila kikkawai has four amylase genes. The Amy1 and Amy2 genes are a head-to-head duplication in the middle of the B arm of the second chromosome, while the Amy3 and Amy4 genes are a tail-to-tail duplication near the centromere of the same chromosome. In the sibling species of D. kikkawai (Drosophila bocki, Drosophila leontia, and Drosophila lini), sequencing of the Amy1, Amy2, Amy3, and Amy4 genes revealed that the Amy1 and Amy2 gene group diverged from Amy3 and Amy4 after duplication. In the Amy1 and Amy2 genes, the divergent evolution occurred in the flanking regions; in contrast, the coding regions have evolved in concerted fashion. The electrophoretic pattern of AMY isozymes was also examined. In D. kikkawai and its siblings, two or three electrophoretically different isozymes are encoded by the Amy1 and Amy2 genes (S isozyme) and by the Amy3 and Amy4 genes (F (M) isozymes). The S and F (M) isozymes show different patterns of band intensity when larvae and flies were fed in different media. Amy1 and Amy2, which encode the S isozyme, are more strikingly regulated than Amy3 and Amy4, which encode the F (M) isozyme. The GC content and codon usage bias were higher for the Amy1 and Amy2 genes than for the Amy3 and Amy4 genes. Although the ratio of synonymous and replacement substitutions within the Amy1 and Amy2 gene group was not significantly different from that within the Amy3 and Amy4 gene group, the synonymous substitution rate in the lineage of Amy1 and Amy2 was lower than that of Amy3 and Amy4. In conclusion, after the first duplication but before speciation of four species, the synonymous substitution rate between the two lineages and the electrophoretic pattern of the isozymes encoded by them changed, although we do not know whether there was any evolutionary relationship between the two.     

Evolutionary Changes in Gene Expression,Coding Sequence and Copy-Number at the Cyp6g1 Locus Contribute to Resistance to Multiple Insecticides in Drosophila

Widespread use of insecticides has led to insecticide resistance in many populations of insects. In some populations, resistance has evolved to multiple pesticides. In Drosophila melanogaster, resistance to multiple classes of insecticide is due to the overexpression of a single cytochrome P450 gene, Cyp6g1. Overexpression of Cyp6g1 appears to have evolved in parallel in Drosophila simulans, a sibling species of D. melanogaster, where it is also associated with insecticide resistance. However, it is not known whether the ability of the CYP6G1 enzyme to provide resistance to multiple insecticides evolved recently in D. melanogaster or if this function is present in all Drosophila species. Here we show that duplication of the Cyp6g1 gene occurred at least four times during the evolution of different Drosophila species, and the ability of CYP6G1 to confer resistance to multiple insecticides exists in D. melanogaster and D. simulans but not in Drosophila willistoni or Drosophila virilis. In D. virilis, which has multiple copies of Cyp6g1, one copy confers resistance to DDT and another to nitenpyram, suggesting that the divergence of protein sequence between copies subsequent to the duplication affected the activity of the enzyme. All orthologs tested conferred resistance to one or more insecticides, suggesting that CYP6G1 had the capacity to provide resistance to anthropogenic chemicals before they existed. Finally, we show that expression of Cyp6g1 in the Malpighian tubules, which contributes to DDT resistance in D. melanogaster, is specific to the D. melanogaster–D. simulans lineage. Our results suggest that a combination of gene duplication, regulatory changes and protein coding changes has taken place at the Cyp6g1 locus during evolution and this locus may play a role in providing resistance to different environmental toxins in different Drosophila species.     

Duplicative and conservative transpositions of larval serum protein 1 genes in the genus Drosophila

Interspecific comparative molecular analyses of transposed genes and their flanking regions can help to elucidate the time, direction, and mechanism of gene transposition. In the Drosophila melanogaster genome, three Larval serum protein 1 (Lsp1) genes (alpha, beta and gamma) are present and each of them is located on a different chromosome, suggesting multiple transposition events. We have characterized the molecular organization of Lsp1 genes in D. buzzatii, a species of the Drosophila subgenus and in D. pseudoobscura, a species of the Sophophora subgenus. Our results show that only two Lsp1 genes (beta and gamma) exist in these two species. The same chromosomal localization and genomic organization, different from that of D. melanogaster, is found in both species for the Lsp1beta and Lsp1gamma genes. Overall, at least two duplicative and two conservative transpositions are necessary to explain the present chromosomal distribution of Lsp1 genes in the three Drosophila species. Clear evidence for implication of snRNA genes in the transposition of Lsp1beta in Drosophila has been found. We suggest that an ectopic exchange between highly similar snRNA sequences was responsible for the transposition of this gene. We have also identified the putative cis-acting regulatory regions of these genes, which seemingly transposed along with the coding sequences.     

Evolutionary profiling of the U49 snoRNA gene

Small nucleolar RNAs (snoRNAs) are involved in precursor ribosomal RNA (pre-rRNA) processing and rRNA base modifications (2'-O-ribose methylation and pseudouridylation). Their genomic organization show great flexibility: some are individually or polycistronically transcribed, while others are encoded within introns of other genes. Here, we present an evolutionary analysis of the U49 gene in seven species. In all species analyzed, U49 contains the typical hallmarks of C and D box motifs, and a conserved 12-15 nt sequence complementary to rRNA that define them as homologs. In mouse, human, and Drosophila U49 is found encoded within introns of different genes, and in plants it is transcribed polycistronically from four different locations. In addition, U49 has two copies in two different introns of the RpL14 gene in Drosophila. The results indicate a substantial degree of duplication and translocation of the U49 gene in evolution. In light of its variable organization we discuss which of the two proposed mechanisms of rearrangement has acted upon the U49 snoRNA gene: chromosomal duplication or transposition through an RNA intermediate.     

Molecular phylogeny and divergence times of drosophilid species

The phylogenetic relationships and divergence times of 39 drosophilid species were studied by using the coding region of the Adh gene. Four genera--Scaptodrosophila, Zaprionus, Drosophila, and Scaptomyza (from Hawaii)--and three Drosophila subgenera--Drosophila, Engiscaptomyza, and Sophophora--were included. After conducting statistical analyses of the nucleotide sequences of the Adh, Adhr (Adh-related gene), and nuclear rRNA genes and a 905-bp segment of mitochondrial DNA, we used Scaptodrosophila as the outgroup. The phylogenetic tree obtained showed that the first major division of drosophilid species occurs between subgenus Sophophora (genus Drosophila) and the group including subgenera Drosophila and Engiscaptomyza plus the genera Zaprionus and Scaptomyza. Subgenus Sophophora is then divided into D. willistoni and the clade of D. obscura and D. melanogaster species groups. In the other major drosophilid group, Zaprionus first separates from the other species, and then D. immigrans leaves the remaining group of species. This remaining group then splits into the D. repleta group and the Hawaiian drosophilid cluster (Hawaiian Drosophila, Engiscaptomyza, and Scaptomyza). Engiscaptomyza and Scaptomyza are tightly clustered. Each of the D. repleta, D. obscura, and D. melanogaster groups is monophyletic. The splitting of subgenera Drosophila and Sophophora apparently occurred about 40 Mya, whereas the D. repleta group and the Hawaiian drosophilid cluster separated about 32 Mya. By contrast, the splitting of Engiscaptomyza and Scaptomyza occurred only about 11 Mya, suggesting that Scaptomyza experienced a rapid morphological evolution. The D. obscura and D. melanogaster groups apparently diverged about 25 Mya. Many of the D. repleta group species studied here have two functional Adh genes (Adh-1 and Adh-2), and these duplicated genes can be explained by two duplication events.