Construction of a single nucleotide polymorphism linkage map for the silkworm, Bombyx mori, based on bacterial artificial chromosome end sequences

We have developed a linkage map for the silkworm Bombyx mori based on single nucleotide polymorphisms (SNPs) between strains p50T and C108T initially found on regions corresponding to the end sequences of bacterial artificial chromosome (BAC) clones. Using 190 segregants from a backcross of a p50T female x an F1 (p50T x C108T) male, we analyzed segregation patterns of 534 SNPs between p50T and C108T, detected among 3840 PCR amplicons, each associated with a p50T BAC end sequence. This enabled us to construct a linkage map composed of 534 SNP markers spanning 1305 cM in total length distributed over the expected 28 linkage groups. Of the 534 BACs whose ends harbored the SNPs used to construct the linkage map, 89 were associated with 107 different ESTs. Since each of the SNP markers is directly linked to a specific genomic BAC clone and to whole-genome sequence data, and some of them are also linked to EST data, the SNP linkage map will be a powerful tool for investigating silkworm genome properties, mutation mapping, and map-based cloning of genes of industrial and agricultural interest.

Role of the silkworm argonaute2 homolog gene in double-strand break repair of extrachromosomal DNA

The argonaute protein family provides central components for RNA interference (RNAi) and related phenomena in a wide variety of organisms. Here, we isolated, from a Bombyx mori cell, a cDNA clone named BmAGO2, which is homologous to Drosophila ARGONAUTE2, the gene encoding a repressive factor for the recombination repair of extrachromosomal double-strand breaks (DSBs). RNAi-mediated silencing of the BmAGO2 sequence markedly increased homologous recombination (HR) repair of DSBs in episomal DNA, but had no effect on that in chromosomes. Moreover, we found that RNAi for BmAGO2 enhanced the integration of linearized DNA into a silkworm chromosome via HR. These results suggested that BmAgo2 protein plays an indispensable role in the repression of extrachromosomal DSB repair.

Lipophorin receptor of Bombyx mori: cDNA cloning, genomic structure, alternative splicing, and isolation of a new isoform

The cDNA and genomic structure of a putative lipophorin receptor from the silkworm, Bombyx mori (BmLpR), indicated the presence of four isoforms, designated LpR1, LpR2, LpR3, and LpR4. The deduced amino acid sequence of each isoform showed five functional domains that are homologous to vertebrate very low density lipoprotein receptor (VLDLR). All four isoforms seem to have originated from a single gene by alternative splicing and were differentially expressed in a tissue- and stage-specific manner. BmLpR1 harbored an additional 27 amino acids in the O-linked sugar domain, resulting in an extra exon. The silkworm BmLpR gene consisted of 16 exons separated by 15 introns spanning >122 kb and was at least three times larger than the human VLDLR gene. Surprisingly, one of the isoforms, LpR4, was expressed specifically in the brain and central nervous system. Additionally, it had a unique cytoplasmic tail, leading to the proposition that it represents a new candidate LpR for possible brain-related function(s). This is the first report on the genomic characterization of an arthropod lipoprotein receptor gene and the identification of a brain-specific receptor variant from a core member of the low density lipoprotein receptor family in invertebrates.

Glycine-rich protein genes, which encode a major component of the cuticle, have different developmental profiles from other cuticle protein genes in Bombyx mori

Three types of GRP (glycine-rich proteins) cDNAs were identified in the EST database of Bombyx mori. These came from 21 ESTs in the W3-stage wing disc EST library. We named them BmGRP1, BmGRP2 and BmGRP3. BmGRP1 and BmGRP2 had 57% identity in deduced amino acid sequences. Expression of all BmGRPs was observed in the epidermis at the fourth molting stage, and in the wing at pupation and mid-pupal stage. It is suggested that BmGRPs contribute to larval, pupal and adult cuticles together with other cuticle proteins. Transcripts of BmGRP2 increased after 7 days of pupal stage. BmGRP2 is suggested to construct adult trachea in the wing. Hormonal response of BmGRPs was compared with that of another group of cuticular protein genes, BMWCPs. BmGRPs were induced by a pulse of 20E. Induction of BmGRP3 was observed in W1 wing discs in the presence of JHA which was added with 20E, whereas that of BMWCP2 was inhibited in the presence of JHA. Induction of BmGRPs was observed in the wing discs of V3 and W1 stages, while that of BMWCP2 was not observed in the V3 wing discs. These differences between BMWCPs and BmGRPs in response to hormones at different developmental stages are discussed.

Mechanomyographic determination of post-activation potentiation in myopathies

OBJECTIVE

There is a need to provide an index of muscle contractility in evaluating myopathies especially in the clinical setting. This study was conducted to investigate if the mechanomyogram post-activation potentiation (MMG-PAP) can be used as an index of muscle contractile force potentiation (force-PAP), if it differs between normal and myopathic muscles, and if it can reflect abnormalities in muscle fiber anatomy.

METHODS

The correlation between MMG-PAP and force-PAP was evaluated in 12 normal subjects after maximum voluntary contraction (MVC) of the biceps brachii muscle. The same method was then applied to study MMG-PAP in 16 patients with myopathies, 16 disease and 25 normal controls. Mean fiber diameters and the proportions of type 1 and 2 fibers in biopsied biceps brachii muscle were determined and compared with MMG-PAP values.

RESULTS

There was a significant positive correlation between force-PAP (197 +/- 148%) and MMG-PAP (135 +/- 68%) immediately after MVC (P < 0.05). The mean MMG-PAP in myopathies (66 +/- 53%) was significantly lower than those of the disease (128 +/- 34%; P < 0.005) and normal controls (120 +/- 56%; P < 0.005). Patients with non-dystrophic myopathies, including those with myositis, had significantly lower MMG-PAP values (38 +/- 20%; P < 0.005) than those with muscular dystrophy (148 +/- 23%). MMG-PAP did not clearly correlate with either type 2 fiber atrophy or type 2 fiber disproportion based on muscle biopsy analysis of myopathic patients.

CONCLUSIONS

This study shows that MMG-PAP can be used as an index of muscle contractility and that it is significantly lower in non-dystrophic myopathies compared to normal subjects. MMG-PAP does not seem to reflect abnormal muscle fiber anatomy.

SIGNIFICANCE

MMG-PAP may become a valuable non-invasive tool in augmenting routine clinical electrophysiologic studies especially in evaluating muscle contractility in myopathies.

Simple sequence repeat-based consensus linkage map of Bombyx mori

We established a genetic linkage map employing 518 simple sequence repeat (SSR, or microsatellite) markers for Bombyx mori (silkworm), the economically and culturally important lepidopteran insect, as part of an international genomics program. A survey of six representative silkworm strains using 2,500 (CA)n- and (CT)n-based SSR markers revealed 17-24% polymorphism, indicating a high degree of homozygosity resulting from a long history of inbreeding. Twenty-nine SSR linkage groups were established in well characterized Dazao and C108 strains based on genotyping of 189 backcross progeny derived from an F(1) male mated with a C108 female. The clustering was further focused to 28 groups by genotyping 22 backcross progeny derived from an F(1) female mated with a C108 male. This set of SSR linkage groups was further assigned to the 28 chromosomes (established linkage groups) of silkworm aided by visible mutations and cleaved amplified polymorphic sequence markers developed from previously mapped genes, cDNA sequences, and cloned random amplified polymorphic DNAs. By integrating a visible mutation p (plain, larval marking) and 29 well conserved genes of insects onto this SSR-based linkage map, a second generation consensus silkworm genetic map with a range of 7-40 markers per linkage group and a total map length of approximately 3431.9 cM was constructed and its high efficiency for genotyping and potential application for synteny studies of Lepidoptera and other insects was demonstrated.

Partial deletions of the W chromosome due to reciprocal translocation in the silkworm Bombyx mori

In the silkworm, Bombyx mori (female, ZW; male, ZZ), femaleness is determined by the presence of a single W chromosome, irrespective of the number of autosomes or Z chromosomes. The W chromosome is devoid of functional genes, except the putative female-determining gene (Fem). However, there are strains in which chromosomal fragments containing autosomal markers have been translocated on to W. In this study, we analysed the W chromosomal regions of the Zebra-W strain (T(W;3)Ze chromosome) and the Black-egg-W strain (T(W;10)+(w-2) chromosome) at the molecular level. Initially, we undertook a project to identify W-specific RAPD markers, in addition to the three already established W-specific RAPD markers (W-Kabuki, W-Samurai and W-Kamikaze). Following the screening of 3648 arbitrary 10-mer primers, we obtained nine W-specific RAPD marker sequences (W-Bonsai, W-Mikan, W-Musashi, W-Rikishi, W-Sakura, W-Sasuke, W-Yukemuri-L, W-Yukemuri-S and BMC1-Kabuki), almost all of which contained the border regions of retrotransposons, namely portions of nested retrotransposons. We confirmed the presence of eleven out of twelve W-specific RAPD markers in the normal W chromosomes of twenty-five silkworm strains maintained in Japan. These results indicate that the W chromosomes of the strains in Japan are almost identical in type. The Zebra-W strain (T(W;3)Ze chromosome) lacked the W-Samurai and W-Mikan RAPD markers and the Black-egg-W strain (T(W;10)+(w-2) chromosome) lacked the W-Mikan RAPD marker. These results strongly indicate that the regions containing the W-Samurai and W-Mikan RAPD markers or the W-Mikan RAPD marker were deleted in the T(W;3)Ze and T(W;10)+(w-2) chromosomes, respectively, due to reciprocal translocation between the W chromosome and the autosome. This deletion apparently does not affect the expression of Fem; therefore, this deleted region of the W chromosome does not contain the putative Fem gene.

Retrotransposable elements on the W chromosome of the silkworm, Bombyx mori

The sex chromosomes of the silkworm, Bombyxmori, are designated ZW(XY) for females and ZZ(XX) for males. The W chromosome of B. mori does not recombine with the Z chromosome and autosomes and no genes for morphological characters have been mapped to the W chromosome as yet. Furthermore, femaleness is determined by the presence of a single W chromosome, regardless of the number of autosomes or Z chromosomes. To understand these interesting features of the W chromosome, it is necessary to analyze the W chromosome at the molecular biology level. Initially to isolate DNA sequences specific for the W chromosome as randomly amplified polymorphic DNA (RAPD) markers, we compared the genomic DNAs between males and females by PCR with arbitrary 10-mer primers. To the present, we have identified 12 W-specific RAPD markers, and with the exception of one RAPD marker, all of the deduced amino acid sequences of these W-specific RAPD markers show similarity to previously reported amino acid sequences of retrotransposable elements from various organisms. After constructing a genomic DNA lambda phage library of B. mori we obtained two lambda phage clones, one containing the W-Kabuki RAPD sequence and one containing the W-Samurai RAPD sequence and found that these DNA sequences comprised nested structures of many retrotransposable elements. To further analyze the W chromosome, we obtained 14 W-specific bacterial artificial chromosome (BAC) clones from three BAC libraries and subjected these clones to shotgun sequencing. The resulting assembly of sequences did not produce a single contiguous sequence due to the presence of many retrotransposable elements. Therefore, we coupled PCR with shotgun sequencing. Through these analyses, we found that many long terminal repeat (LTR) and non-LTR retrotransposons, retroposons, DNA transposons and their derivatives, have accumulated on the W chromosome as strata. These results strongly indicate that retrotransposable elements are the main structural component of the W chromosome.

Identification of molting fluid carboxypeptidase A (MF-CPA) in Bombyx mori

Using microarray analyses, we identified carboxypeptidase A (MF-CPA), which was induced during pupal ecdysis in the wing discs of Bombyx mori. Here, we report the functional characterization of MF-CPA. MF-CPA has amino acid sequence similarities with the proteins in the carboxypeptidase A/B subfamily, from human to nematode. The MF-CPA gene is expressed during the molting periods in the epithelial tissues. MF-CPA is detected in the molting fluid, which fills the space between the old and new cuticle during molting. By Western blot analysis, we show that MF-CPA is secreted as a zymogen and processed in the molting fluid. Recombinant MF-CPA expressed in the insect cells has carboxypeptidase A activity. We propose that MF-CPA degrades the proteins from the old cuticle during the molting periods and contributes to recycling of the amino acids.

Novel macula-like virus identified in Bombyx mori cultured cells

We identified a novel, 6,513-bp-long RNA, termed Bombyx mori macula-like latent virus (BmMLV) RNA, which abundantly expressed in B. mori cultured BmN cells. BmMLV RNA potentially encodes two proteins, putative RNA replicase and coat protein, which share structural features and sequence similarities with those of a plant RNA virus, the genus Maculavirus. Northern blot analysis showed that two transcripts were expressed in BmN cells: a 6.5-kb-long RNA, which contains both putative RNA replicase and coat protein genes, and a 1.2-kb-long RNA, which contains only a coat protein gene. Southern blot analysis showed that BmMLV RNA is not carried by the B. mori genome. RT-PCR analysis also revealed the presence of BmMLV RNA in several B. mori cell lines other than BmN cells, suggesting that BmMLV RNA latently exists in B. mori cultured cells. Infection studies showed that BmMLV virions were able to infect BmMLV-negative Spodoptera frugiperda Sf-9 cells and B. mori larvae. Electron microscopy and Northern blot analysis of a purified BmMLV revealed that isometric virions appear to be 28 to 30 nm in diameter and contain a 6.5-kb genomic RNA. These results showed that BmMLV is a novel macula-like virus infectious to and replicable in B. mori-derived cells.

Characteristics of two genes encoding proteins with an ADAM-type metalloprotease domain, which are induced during the molting periods in Bombyx mori

By microarray analyses, we identified two genes (BmADAMTS-1 and BmADAMTS-like) encoding a protein, which are induced during the pupal ecdysis in the wing discs of Bombyx mori; these genes are homologous to ADAMTS family members (a disintegrin and metalloproteinase domain, with thrombospondin type-1 repeats). A complete metal-binding motif of the ADAM-type metalloprotease domain (HEXXHXXGXXHD) was contained in both amino acid sequences. However, thrombospondin type 1 (TSP-1) repeats were observed only in BmADAMTS-1. The BmADAMTS-1 gene was expressed in the hemocyte and midgut of the larvae at day 2 of wandering stage (W2), and strongly induced during the pupal ecdysis in the hemolymph. The BmADAMTS-like gene was expressed in the epithelial tissues of the larvae at W2, and had expression peaks slightly later than the BmADAMTS-1 gene. Our results indicate that BmADAMTS-1 and BmADAMTS-like may cleave the extracellular matrix (ECM) in the degenerating and remodeling tissues during the molting periods.

Hormonal control of a metamorphosis-specific transcriptional factor Broad-Complex in silkworm

Insect metamorphosis is induced by the steroid 20-hydroxyecdysone (20E) in the absence of sesquiterpenoid juvenile hormone (JH). In Drosophila melanogaster, the Broad-Complex (BR-C) transcriptional factor plays critical roles during metamorphosis. We isolated cDNAs encoding BR-C in the silkworm Bombyx mori and examined their mRNA expression. cDNAs for three BR-C isoforms with zinc finger pairs (Z1, Z2 and Z4) and four isoforms lacking them were cloned. Their mRNAs were expressed in multiple tissues at the larval-pupal metamorphosis. In the anterior silk gland, BR-C mRNAs were expressed at the end of the last larval instar but not expressed during the penultimate instar. 20E administration induced BR-C mRNA expression and JH suppressed 20E-induced BR-C expression in this tissue both in vivo and in vitro. Thus, BR-C mRNAs are inducible by 20E only in the absence of JH, a finding that explains their metamorphosis-specific expression.

SilkSatDb: a microsatellite database of the silkworm, Bombyx mori

The SilkSatDb (silkmoth microsatellite database) (http://www.cdfd.org.in/silksatdb) is a relational database of microsatellites extracted from the available expressed sequence tags and whole genome shotgun sequences of the silkmoth, Bombyx mori. The database has been rendered with a simple and robust web-based search facility, developed using PHP. The SilkSatDb also stores information on primers developed and validated in the laboratory. Users can retrieve information on the microsatellite and the protocols used, along with informative figures and polymorphism status of those microsatellites. In addition, the interface is coupled with Autoprimer, a primer-designing program, using which users can design primers for the loci of interest.

Change in the expressed gene patterns of the wing disc during the metamorphosis of Bombyx mori

The present study was conducted to clarify what occurs during the metamorphosis of the imaginal disc in insects. To understand the metamorphosis on a molecular level, the changes in expression profiles in the imaginal disc during metamorphosis were investigated. For this purpose, we constructed cDNA libraries from four different stages of wing discs of Bombyx mori, sequenced about 1000 cDNAs randomly collected from each library, and constructed a database of expressed sequence tags (EST). The morphological changes and expression profiles from EST were compared during those four stages. Microarray analysis was applied to quantify the expression of each gene in each stage in order to confirm whether the expression of the genes identified from EST was induced by 20-hydroxyecdysone (20E) in a stage-specific manner. Wing discs showed dynamic morphogenesis in 4-5 days during the preparatory stage of metamorphosis, which was under the control of an ecdysteroid. Different expressed profiles were observed in each of the four different stages by comparison of each EST clone. These profiles reflected the morphological changes of the Bombyx wing disc during metamorphosis. The results of expression profiles from the four stages suggested that the V4 stage was cell proliferating; W0, proliferating and the beginning of differentiation; W2, morphologically changing; W3, cuticle secreting. Microarray analysis showed the effectiveness of its application on 20E induction of genes in wing discs. The wing disc of B. mori is an exceptionally suitable system for understanding the relationship between morphological changes and the distribution of mRNA.

A baculovirus-encoded protein tyrosine phosphatase gene induces enhanced locomotory activity in a lepidopteran host

Enhanced locomotory activity (ELA), such as wandering, is a normal behavior that occurs at the end of the larval stage in lepidopteran (butterflies and moths) insects. Baculovirus infection can also induce ELA in lepidopteran larvae. The belief is that the virus induces this behavior to increase its transmission [Goulson, D. (1997) Oecologia 109, 219-228]. Here we show that a baculovirus-encoded protein tyrosine phosphatase (PTP) gene (ptp) induces ELA that is activated by light. ELA was induced in silkworm Bombyx mori infected with the baculovirus B. mori nucleopolyhedrovirus (BmNPV) beginning at approximately 3.75 days postinfection (p.i.) and continued until 4.75 days p.i. The intensity of the ELA was dramatically reduced immediately before death at 5.25 days p.i. Light activated the intensity of the ELA by approximately 3-fold, and larvae with ELA showed positive phototropism. ELA was not induced in larvae of B. mori infected with a BmNPV ptp knockout mutant (BmPTPD). However, when a silkworm-derived ptp gene (Bmptp-h) was inserted into BmPTPD, ELA was partially recovered. Bmptp-h was identified from silkworms at 2 days after the start of the natural wandering stage. The deduced amino acid sequence of Bmptp-h showed 48.2% identify (80.7% similarity) to the deduced amino acid sequence of BmNPV ptp. On the basis of the high homology and larval stage at which Bmptp-h was isolated, we postulate that the modern baculovirus may have acquired its ptp gene from an ancestral host and that this gene was selectively maintained because it increases virus transmission.

Interferon-beta induction/interferon-alpha2b plus ribavirin therapy in patients with chronic hepatitis C

Treatment of chronic hepatitis C virus (HCV) infection with interferon (IFN) and ribavirin improves the rate of eradication of the virus by less than 20% in patients with genotype 1b and a high viral load. In this study we assessed whether IFN-beta induction/IFN-alpha2b plus ribavirin enhances the efficacy of the therapy in patients with chronic hepatitis C. The efficacy of IFN-beta induction/IFN-alpha2b plus ribavirin therapy (group A, n=7) was compared with that of IFN-alpha2b plus ribavirin (group B, n=7) in 14 patients with high levels of HCV-RNA (> 100 K/U/ml). No significant differences were observed in the clearance of HCV-RNA between the two groups (A and B, respectively) 2 weeks after the start of the treatment (0% and 14.3%), at the end of the treatment (71.4% and 100%) and 6 months after the end of the treatment (28.6% and 14.3%). Recovery was complete in 28.6% and 14.3%, transient in 42.9% and 85.7% and absent in 28.6% and 0% in groups A and B, respectively. Early log changes in the viral load from the baseline after 2 weeks of treatment were 2.41 +/- 0.91 and 2.77 +/- 0.20 in groups A and B, respectively, with no significant difference between the two groups. In the present study, we were not able to demonstrate that IFN-beta induction/IFN-alpha2b plus ribavirin therapy was superior to IFN-alpha2b plus ribavirin therapy in patients with genotype 1b and high viral loads.

Electromyographic and mechanomyographic estimation of motor unit activation strategy in voluntary force production

Electromyographic and mechanomyographic estimation of motor unit activation strategy in voluntary force production. In order to determine whether electromyogram (EMG) and mechanomyogram (MMG) are suitable for the noninvasive estimation of the motor unit (MU) activation strategy, the EMG/force and MMG/force relationships were examined simultaneously during isometric ramp contractions in biceps brachii muscle. The highest mean power frequency (MPF) of the EMG, which reflects the full MU recruitment, was determined at 51% MVC. Two obvious inflection points were identified on the MMG-amplitude/force relationship that showed an initial slow increase followed by a rapid increase and a progressive decrease at higher force levels. Our results suggest that the MMG amplitude allows the estimation of the beginning of recruitment of MUs that innervate the first-twitch fibers in addition to identification of the full MU recruitment. The rate coding strategy was qualitatively reflected by the MMG-MPF/force relationship. We conclude that the MU activation strategy is estimated in more detail by the MMG than by the EMG.

The genome sequence of silkworm, Bombyx mori

We performed threefold shotgun sequencing of the silkworm (Bombyx mori) genome to obtain a draft sequence and establish a basic resource for comprehensive genome analysis. By using the newly developed RAMEN assembler, the sequence data derived from whole-genome shotgun (WGS) sequencing were assembled into 49,345 scaffolds that span a total length of 514 Mb including gaps and 387 Mb without gaps. Because the genome size of the silkworm is estimated to be 530 Mb, almost 97% of the genome has been organized in scaffolds, of which 75% has been sequenced. By carrying out a BLAST search for 50 characteristic Bombyx genes and 11,202 non-redundant expressed sequence tags (ESTs) in a Bombyx EST database against the WGS sequence data, we evaluated the validity of the sequence for elucidating the majority of silkworm genes. Analysis of the WGS data revealed that the silkworm genome contains many repetitive sequences with an average length of <500 bp. These repetitive sequences appear to have been derived from truncated transposons, which are interspersed at 2.5- to 3-kb intervals throughout the genome. This pattern suggests that silkworm may have an active mechanism that promotes removal of transposons from the genome. We also found evidence for insertions of mitochondrial DNA fragments at 9 sites. A search for Bombyx orthologs to Drosophila genes controlling sex determination in the WGS data revealed 11 Bombyx genes and suggested that the sex-determining systems differ profoundly between the two species.

Organization of the Hox gene cluster of the silkworm, Bombyx mori: a split of the Hox cluster in a non-Drosophila insect

A bacterial artificial chromosome (BAC) contig was constructed by chromosome walking, starting from the Hox genes of the silkworm, Bombyx mori. Bombyx orthologues of the labial (lab) and zerknult (zen) genes were newly identified. The size of the BAC contig containing the Hox gene cluster-except the lab and Hox 2 genes-was estimated to be more than 2 Mb. The Bombyx Hox cluster was mapped to linkage group (LG) 6. The lab gene was mapped on the same LG, but far apart from the cluster. Fluorescence in situ hybridization analysis confirmed that the major Hox gene cluster and lab were at different locations on the same chromosome in B. mori.

Survey and analysis of microsatellites in the silkworm, Bombyx mori: frequency, distribution, mutations, marker potential and their conservation in heterologous species

We studied microsatellite frequency and distribution in 21.76-Mb random genomic sequences, 0.67-Mb BAC sequences from the Z chromosome, and 6.3-Mb EST sequences of Bombyx mori. We mined microsatellites of >/=15 bases of mononucleotide repeats and >/=5 repeat units of other classes of repeats. We estimated that microsatellites account for 0.31% of the genome of B. mori. Microsatellite tracts of A, AT, and ATT were the most abundant whereas their number drastically decreased as the length of the repeat motif increased. In general, tri- and hexanucleotide repeats were overrepresented in the transcribed sequences except TAA, GTA, and TGA, which were in excess in genomic sequences. The Z chromosome sequences contained shorter repeat types than the rest of the chromosomes in addition to a higher abundance of AT-rich repeats. Our results showed that base composition of the flanking sequence has an influence on the origin and evolution of microsatellites. Transitions/transversions were high in microsatellites of ESTs, whereas the genomic sequence had an equal number of substitutions and indels. The average heterozygosity value for 23 polymorphic microsatellite loci surveyed in 13 diverse silkmoth strains having 2-14 alleles was 0.54. Only 36 (18.2%) of 198 microsatellite loci were polymorphic between the two divergent silkworm populations and 10 (5%) loci revealed null alleles. The microsatellite map generated using these polymorphic markers resulted in 8 linkage groups. B. mori microsatellite loci were the most conserved in its immediate ancestor, B. mandarina, followed by the wild saturniid silkmoth, Antheraea assama.

Annotation pattern of ESTs from Spodoptera frugiperda Sf9 cells and analysis of the ribosomal protein genes reveal insect-specific features and unexpectedly low codon usage bias

MOTIVATION

A whole set of Expressed Sequence Tags (ESTs) from the Sf9 cell line of Spodoptera frugiperda is presented here for the first time. By this way we want to identify both conserved and specific genes of this pest species. We also expect from this analysis to find a class of protein sequences providing a tool to explore genomic features and phylogeny of Lepidoptera.

RESULTS

The ESTs display both housekeeping as well as developmentally regulated genes, and a high percentage of sequences with unknown function. Among the identified ORFs, almost all ribosomal proteins (RPs) were found with high EST redundancy and hence sequence accuracy. The codon usage found among RP genes is in average surprisingly much less biased in Lepidoptera than in other organisms. Other Spodoptera genes also displayed a low bias, suggesting a general genome expression feature in this Lepidoptera. We also found that the L35A and L36 RP sequences, respectively, display 40 and 10 amino-acid insertions, both being present only in insects. Sequence analysis suggests that they are probably not subjected to a strong selective pressure and may be good phylogenetic markers for Lepidoptera. Most interestingly, the Lepidoptera sequences of 9 RP genes displayed a specific signature different from the canonical one. We conclude that the RP family allows valuable comparative genomics and phylogeny of Lepidoptera.

AVAILABILITY

All EST sequence data are available from the private 'Spodo-Base' upon request.

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.