A draft sequence for the genome of the domesticated silkworm (Bombyx mori)

Identification and expression pattern of Bmlark, a homolog of the Drosophila gene lark in Bombyx mori

Two Bombyx mori isoforms of the gene lark, which is shown to play an important role in Drosophila circadian rhythms, were identified and named Bmlark-PA and Bmlark-PB, respectively. Bmlark-PA consists of 5 exons and encodes a protein of 343 amino acid residues which contains 3 functional domains: two RRM (RNA recognization motif) domains and an RTZF (retroviral-type zinc finger) and shares 72% identity with the Drosophila gene lark at the amino acid level. Bmlark-PB lacks the sequence between 118 and 791 nt of Bmlark-PA and codes for a protein of 68 amino acid residues, which contains no distinct functional domains. Alignments of the cDNAs of Bmlark to the genomic draft sequence of B. mori showed that the gene Bmlark had a single copy in the genome, suggesting that an alternative splicing mechanism occurs in the gene Bmlark. RT-PCR analysis indicated that Bmlark-PA was expressed only in late pupae and adult but Bmlark-PB was broadly expressed in many tissues and throughout the developmental stages from embryo to adult.

SPODOBASE: an EST database for the lepidopteran crop pest Spodoptera

Background

The Lepidoptera Spodoptera frugiperda is a pest which causes widespread economic damage on a variety of crop plants. It is also well known through its famous Sf9 cell line which is used for numerous heterologous protein productions. Species of the Spodoptera genus are used as model for pesticide resistance and to study virus host interactions. A genomic approach is now a critical step for further new developments in biology and pathology of these insects, and the results of ESTs sequencing efforts need to be structured into databases providing an integrated set of tools and informations.

Description

The ESTs from five independent cDNA libraries, prepared from three different S. frugiperda tissues (hemocytes, midgut and fat body) and from the Sf9 cell line, are deposited in the database. These tissues were chosen because of their importance in biological processes such as immune response, development and plant/insect interaction. So far, the SPODOBASE contains 29,325 ESTs, which are cleaned and clustered into non-redundant sets (2294 clusters and 6103 singletons). The SPODOBASE is constructed in such a way that other ESTs from S. frugiperda or other species may be added. User can retrieve information using text searches, pre-formatted queries, query assistant or blast searches. Annotation is provided against NCBI, UNIPROT or Bombyx mori ESTs databases, and with GO-Slim vocabulary.

Conclusion

The SPODOBASE database provides integrated access to expressed sequence tags (EST) from the lepidopteran insect Spodoptera frugiperda. It is a publicly available structured database with insect pest sequences which will allow identification of a number of genes and comprehensive cloning of gene families of interest for scientific community. SPODOBASE is available from URL:

Bombyx mori cell line as a model of immune-system organs

We tested 11 Bombyx mori cell lines for induction of cecropin B gene (CecB) expression. After the immune challenge, CecB expression was induced in seven cell lines. A mixture of the cell-free supernatant from the immune-responsive cell lines and lipopolysaccharide activated a promoter of CecB in the non-immune-responsive cell line, indicating that secreted factor(s) is involved in CecB activation. The expressed sequence tags of one of the immune-responsive cell lines, NISES-BoMo-Cam1, contained genes encoding proteins similar to Relish, Cactus, clip-domain serine protease, serpin, lectin, peptidoglycan recognition protein, 6tox and gloverin, in addition to seven known B. mori immune-inducible genes. These results show that NISES-BoMo-Cam1 cells can be used as an in vitro model of the immune system organs of B. mori.

Expression of spider flagelliform silk protein in Bombyx mori cell line by a novel Bac-to-Bac/BmNPV baculovirus expression system

Bombyx mori nuclear polyhedrosis virus (BmNPV) baculovirus expression system (BES) has a lot of advantages such as high expression efficiency, convenience, and low feeding cost. In this report, we used a recently developed BmNPV bacmid, which could infect both B. mori cell lines and silkworm larvae. The results showed it takes only 7 to 10 days to generate recombinant baculovirus and permit the rapid isolation from small-scale cultures and then use it to transfect B. mori cell lines, compared to traditional homologous recombination method, which needs at least 40 days for multiple rounds of purification and amplification of viruses. Using this BES, we expressed a recombinant spider flagelliform protein in BmN cell line, which was around 37 kDa in sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and western blot analysis. The BmNPV bacmid system using silkworm would be very attractive for expression of target proteins.

A wing expressed sequence tag resource for Bicyclus anynana butterflies, an evo-devo model

BACKGROUND

Butterfly wing color patterns are a key model for integrating evolutionary developmental biology and the study of adaptive morphological evolution. Yet, despite the biological, economical and educational value of butterflies they are still relatively under-represented in terms of available genomic resources. Here, we describe an Expression Sequence Tag (EST) project for Bicyclus anynana that has identified the largest available collection to date of expressed genes for any butterfly.

RESULTS

By targeting cDNAs from developing wings at the stages when pattern is specified, we biased gene discovery towards genes potentially involved in pattern formation. Assembly of 9,903 ESTs from a subtracted library allowed us to identify 4,251 genes of which 2,461 were annotated based on BLAST analyses against relevant gene collections. Gene prediction software identified 2,202 peptides, of which 215 longer than 100 amino acids had no homology to any known proteins and, thus, potentially represent novel or highly diverged butterfly genes. We combined gene and Single Nucleotide Polymorphism (SNP) identification by constructing cDNA libraries from pools of outbred individuals, and by sequencing clones from the 3' end to maximize alignment depth. Alignments of multi-member contigs allowed us to identify over 14,000 putative SNPs, with 316 genes having at least one high confidence double-hit SNP. We furthermore identified 320 microsatellites in transcribed genes that can potentially be used as genetic markers.

CONCLUSION

Our project was designed to combine gene and sequence polymorphism discovery and has generated the largest gene collection available for any butterfly and many potential markers in expressed genes. These resources will be invaluable for exploring the potential of B. anynana in particular, and butterflies in general, as models in ecological, evolutionary, and developmental genetics.

Identification and characterization of the telomerase reverse transcriptase of Bombyx mori (silkworm) and Tribolium castaneum (flour beetle)

Chromosomal ends of most eukaryotes are composed of simple telomeric repeats. Arthropod telomeres are generally constituted by TTAGG pentanucleotide repeats; however, some insect species including Drosophila melanogaster do not have telomeric repeats. In contrast, the domestic silkworm Bombyx mori contains TTAGG-type telomeric repeats, but the telomerase activity has not been detected in all investigated tissues. To search for a cause of unusual telomere structure in insects, we here identified telomerase reverse transcriptase (TERT) subunit from the domestic silkworm B. mori and the flour beetle Tribolium castaneum. This is the first report of telomerase genes from arthropods. The domestic silkworm TERT gene (BmoTERT) and the flour beetle TERT gene (TcasTERT) both did not have the N-terminal GQ motif. Comparison between cDNA and genomic DNA of BmoTERT revealed that it includes no introns. BmoTERT contains five ATG codons in its 5'UTR, which could reduce the translation of BmoTERT proteins. Also, Northern hybridization indicated that BmoTERT is transcribed at a very low level. These unique features of BmoTERT possibly explain the undetectable Bombyx telomerase activity.

Identification and characterization of piggyBac-like elements in the genome of domesticated silkworm, Bombyx mori

piggyBac is a short inverted terminal repeat (ITR) transposable element originally discovered in Trichoplusia ni. It is currently the preferred vector of choice for enhancer trapping, gene discovery and identifying gene function in insects and mammals. Many piggyBac-like sequences have been found in the genomes of phylogenetically species from fungi to mammals. We have identified 98 piggyBac-like sequences (BmPBLE1-98) from the genome data of domesticated silkworm (Bombyx mori) and 17 fragments from expressed sequence tags (ESTs). Most of the BmPBLE1-98 probably exist as fossils. A total of 21 BmPBLEs are flanked by ITRs and TTAA host dinucleotides, of which 5 contain a single ORF, implying that they may still be active. Interestingly, 16 BmPBLEs have CAC/GTG not CCC/GGG as the characteristic residues of ITRs, which is a surprising phenomenon first observed in the piggyBac families. Phylogenetic analysis indicates that many BmPBLEs have a close relation to mammals, especially to Homo sapiens, only a few being grouped with the T. ni piggyBac element. In addition, horizontal transfer was probably involved in the evolution of the piggyBac-like elements between B. mori and Daphnia pulicaria. The analysis of the BmPBLEs will contribute to our understanding of the characteristic of the piggyBac family and application of piggyBac in a wide range of insect species.

Gene expression profiling between embryonic and larval stages of the silkworm, Bombyx mori

To elucidate the molecular mechanisms associated with metamorphic phenomenon relating to Bombyx mori, an important organism in the sericulture industry, we identified genes that are expressed in the different developmental stages, specifically the embryonic (ES) and larval (LS) stages of B. mori. Of 8230 high-quality ESTs from two full-length enriched cDNA libraries, 3442 of the ES ESTs were coalesced into 1325 clusters, while 4788 were coalesced into 927 clusters. The functional classification of these ESTs based on Gene Ontology showed that the types of genes that are associated with oxidoreductase activity, enzyme inhibition, and larval development were highly observed in LS, whereas the types of genes that are involved in nucleotide binding, enzyme activity, and protein transport activity were highly observed in ES. In addition, when the gene expression profile between ES and LS was examined by counting the EST frequencies in each library, 69 genes were identified as being either up- or down-regulated in the larval stage compared to the embryonic stage (P>0.99) and this was confirmed by semi-quantitative RT-PCR. The results show that genes involved in proteolysis and peptidolysis, and lipid and carbohydrate metabolism were dramatically up-regulated in LS, while those related to protein metabolism, DNA/RNA, and coenzymes were highly down-expressed. In particular, a GO analysis of these genes revealed that genes that are involved in hydrolase activity were observed to be highly expressed in amount as well as diversity in LS, while those involved in nucleic acid binding were highly expressed in ES. These data may contribute to elucidating genetic events that distinguish the developmental stage and to our understanding of the metamorphosis of B. mori.

Activation of epidermal growth factor receptor mediates receptor axon sorting and extension in the developing olfactory system of the moth Manduca sexta

During development of the adult olfactory system of the moth Manduca sexta, olfactory receptor neurons extend axons from the olfactory epithelium in the antenna into the brain. As they arrive at the brain, interactions with centrally derived glial cells cause axons to sort and fasciculate with other axons destined to innervate the same glomeruli. Here we report studies indicating that activation of the epidermal growth factor receptor (EGFR) is involved in axon ingrowth and targeting. Blocking the EGFR kinase domain pharmacologically leads to stalling of many axons in the sorting zone and nerve layer as well as abnormal axonal fasciculation in the sorting zone. We also find that neuroglian, an IgCAM known to activate the EGFR through homophilic interactions in other systems, is transiently present on olfactory receptor neuron axons and on glia during the critical stages of the sorting process. The neuroglian is resistant to extraction with Triton X-100 in the sorting zone and nerve layer, possibly indicating its stabilization by homophilic binding in these regions. Our results suggest a mechanism whereby neuroglian molecules on axons and possibly sorting zone glia bind homophilically, leading to activation of EGFRs, with subsequent effects on axon sorting, pathfinding, and extension, and glomerulus development.

Proteome analysis of silk gland proteins from the silkworm,Bombyx mori

The silk gland of Bombyx mori is an organ specialized for the synthesis and secretion of silk proteins. We report here the resolution of silk gland proteins by 2-DE and the identification of many of those proteins. This was accomplished by dissecting the glands into several sections, with each exhibiting more than 400 protein spots by 2-DE, of which 100 spots were excised and characterized by in-gel digestion followed by PMF. Ninety-three proteins were tentatively identified. These were then categorized into groups involved in silk protein secretion, transport, lipid metabolism, defense, etc. Western blotting of a 2-DE gel using an antibody of the carotenoid binding protein confirmed the presence of this protein in the silk gland. Proteins including fibroin L-chain and P25 were found as multiple isoforms, some of which contained differential amounts of phosphate residues as analyzed by on-probe dephosphorylation. The current analysis contributes to our understanding of proteins expressed by the silk gland not only of the model lepidopteran B. mori, but also to proteins from other silk-producing insects such as Philosamia cynthia ricini.

A novel droplet-tap sample-loading method for two-dimensional gel electrophoresis

A novel procedure, droplet-tap mode, has been devised for sample application for two-dimensional gel electrophoresis (2DE) expression profiles. The sample was loaded by evenly distributed tapping of droplets of the sample on to the rehydration buffer (RB) and then lowering the strip on to the solution surface. At normal loading concentrations, the number of spots obtained was increased by approximately one-third by this new approach compared with the rehydration loading procedure. The method also resulted in significantly improved resolution compared with cup loading when high concentrations of proteins were present, indicating its potential usefulness in micropreparative separation. In addition, recovery of the proteins confirmed that protein uptake was enhanced by use of this method. By enabling improved performances in 2DE, the proposed procedure has much potential for sample loading to meet the requirements of global proteome analysis.

A second-generation integrated map of the silkworm reveals synteny and conserved gene order between lepidopteran insects

A second-generation linkage map was constructed for the silkworm, Bombyx mori, focusing on mapping Bombyx sequences appearing in public nucleotide databases and bacterial artificial chromosome (BAC) contigs. A total of 874 BAC contigs containing 5067 clones (22% of the library) were constructed by PCR-based screening with sequence-tagged sites (STSs) derived from whole-genome shotgun (WGS) sequences. A total of 523 BAC contigs, including 342 independent genes registered in public databases and 85 expressed sequence tags (ESTs), were placed onto the linkage map. We found significant synteny and conserved gene order between B. mori and a nymphalid butterfly, Heliconius melpomene, in four linkage groups (LGs), strongly suggesting that using B. mori as a reference for comparative genomics in Lepidotera is highly feasible.

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.

Structures, regulatory regions, and inductive expression patterns of antimicrobial peptide genes in the silkworm Bombyx mori

Antimicrobial peptides (AMPs) are a group of immune proteins that protect the host from infection. In Drosophila, seven groups of inducible AMPs have been identified, with activities against fungi and gram-positive and gram-negative bacteria. On the basis of the silkworm genome sequence and expressed sequence tags, we identified 35 AMP genes, mostly belonging to the cecropin, moricin, and gloverin gene families. We predicted the core promoters required for gene transcription and the cis-regulatory elements for NF-kappaB/Rel and GATA transcription factors. The expression profiles of these genes after an immune challenge with lipopolysaccharide were examined by reverse transcription PCR. Members of the cecropin B and gloverin A subfamilies were intensely expressed in the fat body after induction. In contrast, those of the moricin B subfamily were not expressed under the same conditions. Such results suggest that these regulatory elements and their positions in the upstream regions play an important role in regulating the transcription of these defense genes.

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.

Synergy between sequence and size in large-scale genomics

Until recently the study of individual DNA sequences and of total DNA content (the C-value) sat at opposite ends of the spectrum in genome biology. For gene sequencers, the vast stretches of non-coding DNA found in eukaryotic genomes were largely considered to be an annoyance, whereas genome-size researchers attributed little relevance to specific nucleotide sequences. However, the dawn of comprehensive genome sequencing has allowed a new synergy between these fields, with sequence data providing novel insights into genome-size evolution, and with genome-size data being of both practical and theoretical significance for large-scale sequence analysis. In combination, these formerly disconnected disciplines are poised to deliver a greatly improved understanding of genome structure and evolution.

Small RNA Viruses of Insects: Expression in Plants and RNA Silencing

Interest in insect small RNA viruses (SRVs) has grown slowly but steadily. A number of new viruses have been analyzed at the sequence level, adding to our knowledge of their diversity at the level of both individual virus species and families. In particular, a number of possible new virus families have emerged. This research has largely been driven by interest in their potential for pest control, as well as in their importance as the causal agents of disease in beneficial arthropods. At the same time, research into known viruses has made valuable contributions to our understanding of an emerging new field of central importance to molecular biology-the existence of RNA-based gene silencing, developmental control, and adaptive immune systems in eukaryotes. Subject to RNA-based adaptive immune responses in their hosts, viruses have evolved a variety of genes encoding proteins capable of suppressing the immune response. Such genes were first identified in plant viruses, but the first examples known from animal viruses were identified in insect RNA viruses. This chapter will address the diversity of insect SRVs, and attempts to harness their simplicity in the engineering of transgenic plants expressing viruses for resistance to insect pests. We also describe RNA interference and antiviral pathways identified in plants and animals, how they have led viruses to evolve genes capable of suppressing such adaptive immunity, and the problems presented by these pathways for the strategy of expressing viruses in transgenic plants. Approaches for countering these problems are also discussed.

EST-based profiling and comparison of gene expression in the silkworm fat body during metamorphosis

The fat body plays key roles in metabolism and sustenance of growth throughout the life of a silkworm, and thus represents a model tissue for studying development and metamorphosis. Analysis of 18,480 ESTs derived from larval and pupal fat body cDNA libraries allowed characterization of gene expression patterns in the fat body during metamorphosis. By Phrap assembling, 6,814 estimated transcription units (TUs) were generated, 2,673 of which contained 60% of the total ESTs and represented 1,366 distinctive known genes. These genes were classified into 14 categories based on their known or putative functions, and most were found to be involved in metabolism and cell defense. Further comparative analysis of the expression profiles revealed that the gene expression patterns in the larval fat body strikingly differed from those in the pupal fat body. From this, we identified a subset of genes potentially associated with metamorphic events such as the formation or destruction of specific tissues, and simultaneously confirmed the stage-specific expression patterns of several identified genes using RT-PCR examination. This work will provide a valuable resource for studying regulatory mechanisms associated with fat body developmental changes and silkworm metamorphosis.

BmStart1, a novel carotenoid-binding protein isoform from Bombyx mori, is orthologous to MLN64, a mammalian cholesterol transporter

Carotenoid-binding protein (CBP) from the silkworm Bombyx mori is an essential molecule for carotenoid dependent cocoon pigmentation. We identified a novel isoform of CBP, Start1 of B. mori (BmStart1). BmStart1 contains a membrane-spanning MENTAL domain in its N-terminus and a lipid-binding START domain in its C-terminus. This domain architecture is identical to the mammalian MLN64 and Start1 of Drosophila melanogaster (DmStart1), both of which have been implicated to function in cholesterol transport and regulation of steroidogenesis. BmStart1 is expressed in both white and yellow cocoon strains of B. mori, while CBP is only detected in the yellow cocoon strain. BmStart1 mRNA abundance in the prothoracic gland, the main ecdysteroidogenic tissue, positively correlates with changes in the hemolymph ecdysteroid level. Genomic sequence analysis revealed that BmStart1 and CBP are generated from the same gene locus by alternative splicing. Splice site comparison and homology search indicate that BmStart1 is orthologous to both MLN64 and DmStart1. This study implies that alternative splicing of the BmStart1/CBP gene generates unique protein isoforms whose endogenous ligands, sterol or carotenoid, are structurally different.

Evolutionary change - patterns and processes

The present review considered: (a) the factors that conditioned the early transition from non-life to life; (b) genome structure and complexity in prokaryotes, eukaryotes, and organelles; (c) comparative human chromosome genomics; and (d) the Brazilian contribution to some of these studies. Understanding the dialectical conflict between freedom and organization is fundamental to give meaning to the patterns and processes of organic evolution.

An expansion of the dual clip-domain serine proteinase family in Manduca sexta: gene organization, expression, and evolution of prophenoloxidase-activating proteinase-2, hemolymph proteinase 12, and other related proteinases

Prophenoloxidase-activating proteinases (PAPs) take part in insect defense responses including melanotic encapsulation and wound healing. To understand their gene structure and regulation, we screened a genomic library and isolated overlapping lambda clones for Manduca sexta PAP-2, hemolymph proteinase 12 (HP12), and HP24. Complete nucleotide sequence analysis indicated that all three genes encode polypeptides with two regulatory clip domains at the amino terminus, a linker region, and a catalytic serine proteinase domain at the carboxyl terminus. Each gene contains eight exons, with introns located at equivalent positions. Similar sequences are present in introns as well as exons, indicating that these genes arose from recent gene duplication and sequence divergence. We analyzed their 5' flanking sequences and identified putative immune and hormone responsive elements. Reverse transcription-polymerase chain reactions confirmed that PAP-2 and HP12 mRNA levels in the larval fat body and hemocytes increased after a bacterial challenge. However, HP24 expression was barely detected. PAP-2 transcripts in cultured fat body became less abundant after 20-hydroxyecdysone treatment. Thus, PAP-2, HP12, and HP24 mRNA levels are differentially regulated by immune and developmental signals. Comparison with HP15, HP23, and PAP-3 sequences suggested an evolutionary pathway of the dual clip-domain serine proteinases in M. sexta.

Radiation-induced changes in gene expression in the silkworm revealed by serial analysis of gene expression (SAGE)

Serial analysis of gene expression (SAGE) was used to examine the profile of expressed genes during embryonic development in the domesticated silkworm, Bombyx mori, after irradiation with Cobalt-60. A comparison of the SAGE sequence tags derived from irradiated embryos with those from normal embryos revealed 673 differentially expressed genes (P < 0.01 and at least three folds change). Of these, 292 genes were highly expressed in normal embryos and 381 genes were highly expressed in irradiated embryos. These results provide valuable information for understanding the mechanisms of radiation-induced changes in gene expression. In addition, it was noted that the generation of longer cDNA fragments from SAGE tags is an efficient way to identify genes, thereby facilitating the analysis of large numbers of unknown genes.

Mutually exclusive splicing of the insect Dscam pre-mRNA directed by competing intronic RNA secondary structures

Drosophila Dscam encodes 38,016 distinct axon guidance receptors through the mutually exclusive alternative splicing of 95 variable exons. Importantly, known mechanisms that ensure the mutually exclusive splicing of pairs of exons cannot explain this phenomenon in Dscam. I have identified two classes of conserved elements in the Dscam exon 6 cluster, which contains 48 alternative exons--the docking site, located in the intron downstream of constitutive exon 5, and the selector sequences, which are located upstream of each exon 6 variant. Strikingly, each selector sequence is complementary to a portion of the docking site, and this pairing juxtaposes one, and only one, alternative exon to the upstream constitutive exon. The mutually exclusive nature of the docking site:selector sequence interactions suggests that the formation of these competing RNA structures is a central component of the mechanism guaranteeing that only one exon 6 variant is included in each Dscam mRNA.

Telomere-specific non-LTR retrotransposons and telomere maintenance in the silkworm, Bombyx mori

Most insects have telomeres that consist of pentanucleotide (TTAGG) telomeric repeats, which are synthesized by telomerase. However, all species in Diptera so far examined and several species in other orders of insect have lost the (TTAGG)n repeats, suggesting that some of them recruit telomerase-independent telomere maintenance. The silkworm, Bombyx mori, retains the TTAGG motifs in the chromosomal ends but expresses quite a low level of telomerase activity in all stages of various tissues. Just proximal to a 6-8-kb stretch of the TTAGG repeats in B. mori, more than 1000 copies of non-LTR retrotransposons, designated TRAS and SART families, occur among the telomeric repeats and accumulate. TRAS and SART are abundantly transcribed and actively retrotransposed into TTAGG telomeric repeats in a highly sequence-specific manner. They have three possible mechanisms to ensure specific integration into the telomeric repeats. This article focuses on the telomere structure and telomere-specific non-LTR retrotransposons in B. mori and discusses the mechanisms for telomere maintenance in this insect.

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.

Juvenile hormone diol kinase, a calcium-binding protein with kinase activity, from the silkworm, Bombyx mori

Juvenile hormone (JH) diol kinase (JHDK) is an important enzyme involved in the JH degradation pathway. Bombyx mori (Bommo)-JHDK cDNA (637bp) contains an open reading frame encoding a 183-amino acid protein, which reveals a high degree of identity to the two previously reported JHDKs. JHDK is similar to GTP-binding proteins with three conserved sequence elements involved in purine nucleotide binding, contains eight alpha-helices and three EF-hand motifs, and resembles the three-dimensional model of 2SCP and some other calcium-binding proteins. The Bommo-JHDK gene has only a single copy in the silkworm haploid genome, contains only one exon, and its 5'-upstream sequence does not have a JH response element. Although Bommo-JHDK is highly expressed in the gut of the silkworm, its mRNA expression remains at a constant level during larval development suggesting this enzyme is constitutive and not regulated by JH, at least at the transcriptional level. Recombinant Bommo-JHDK catalyzed the conversion of 10S-JH diol into JH diol phosphate, confirming its enzymatic function. Recombinant enzyme formed a dimer and had biochemical characteristics similar to other JHDKs. Bommo-JHDK, a calcium-binding protein with kinase activity, provides unique insights on how JH levels are regulated in the silkworm.

Kelp fly virus: a novel group of insect picorna-like viruses as defined by genome sequence analysis and a distinctive virion structure

The complete genomic sequence of kelp fly virus (KFV), originally isolated from the kelp fly, Chaetocoelopa sydneyensis, has been determined. Analyses of its genomic and structural organization and phylogeny show that it belongs to a hitherto undescribed group within the picorna-like virus superfamily. The single-stranded genomic RNA of KFV is 11,035 nucleotides in length and contains a single large open reading frame encoding a polypeptide of 3,436 amino acids with 5' and 3' untranslated regions of 384 and 343 nucleotides, respectively. The predicted amino acid sequence of the polypeptide shows that it has three regions. The N-terminal region contains sequences homologous to the baculoviral inhibitor of apoptosis repeat domain, an inhibitor of apoptosis commonly found in animals and in viruses with double-stranded DNA genomes. The second region contains at least two capsid proteins. The third region has three sequence motifs characteristic of replicase proteins of many plant and animal viruses, including a helicase, a 3C chymotrypsin-like protease, and an RNA-dependent RNA polymerase. Phylogenetic analysis of the replicase motifs shows that KFV forms a distinct and distant taxon within the picorna-like virus superfamily. Cryoelectron microscopy and image reconstruction of KFV to a resolution of 15 A reveals an icosahedral structure, with each of its 12 fivefold vertices forming a turret from the otherwise smooth surface of the 20-A-thick capsid. The architecture of the KFV capsid is unique among the members of the picornavirus superfamily for which structures have previously been determined.

Chromosome evolution in eukaryotes: a multi-kingdom perspective

In eukaryotes, chromosomal rearrangements, such as inversions, translocations and duplications, are common and range from part of a gene to hundreds of genes. Lineage-specific patterns are also seen: translocations are rare in dipteran flies, and angiosperm genomes seem prone to polyploidization. In most eukaryotes, there is a strong association between rearrangement breakpoints and repeat sequences. Current data suggest that some repeats promoted rearrangements via non-allelic homologous recombination, for others the association might not be causal but reflects the instability of particular genomic regions. Rearrangement polymorphisms in eukaryotes are correlated with phenotypic differences, so are thought to confer varying fitness in different habitats. Some seem to be under positive selection because they either trap favorable allele combinations together or alter the expression of nearby genes. There is little evidence that chromosomal rearrangements cause speciation, but they probably intensify reproductive isolation between species that have formed by another route.

Genomic tools and cDNA derived markers for butterflies

The Lepidoptera have long been used as examples in the study of evolution, but some questions remain difficult to resolve due to a lack of molecular genetic data. However, as technology improves, genomic tools are becoming increasingly available to tackle unanswered evolutionary questions. Here we have used expressed sequence tags (ESTs) to develop genetic markers for two Müllerian mimic species, Heliconius melpomene and Heliconius erato. In total 1363 ESTs were generated, representing 330 gene objects in H. melpomene and 431 in H. erato. User-friendly bioinformatic tools were used to construct a nonredundant database of these putative genes (available at http://www.heliconius.org), and annotate them with blast similarity searches, InterPro matches and Gene Ontology terms. This database will be continually updated with EST sequences for the Papilionideae as they become publicly available, providing a tool for gene finding in the butterflies. Alignments of the Heliconius sequences with putative homologues derived from Bombyx mori or other public data sets were used to identify conserved PCR priming sites, and develop 55 markers that can be amplified from genomic DNA in both H. erato and H. melpomene. These markers will be used for comparative linkage mapping in Heliconius and will have applications in other phylogenetic and genomic studies in the Lepidoptera.

Survey of long terminal repeat retrotransposons of domesticated silkworm (Bombyx mori)

Long terminal retrotransposons are major components of eukaryotic transposable elements. We have surveyed the long terminal repeats (LTR) retrotransposons of domesticated silkworm (Bombyx mori) by mining the data produced by Bombyx mori Genome Sequencing Project. At least 29 separate families of LTR retrotransposons are identified in this survey, comprising of 11.8% of the complete sequence. Families of domesticated silkworm LTR retrotransposons can be mainly classified into three groups: gypsy-like, copia-like, Pao-Bel. Fourteen families identified consist of gypsy-like elements, four families consist of copia-like elements and seven families consist of Pao-Bel elements. In addition to the three groups of LTR retrotransposons, two families of unusual non-coding elements are identified in the genome of this species. Further phylogenetic analysis of RT domain indicates that the elements of B.mori show high diversity and can form different clades in each group. An analysis of sequence variation from different families reveals distinct patterns of variation for the elements belonging to three groups. The analysis of the domesticated silkworm LTR retrotransposons should assist in our understanding of the roles of retroelement in lepidopteron insect genome evolution.

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.

Biosynthesis and cocoon-export of a recombinant globular protein in transgenic silkworms

A gene construct was made by fusing the coding sequence of the red fluorescent protein (DsRed) to the exon 2 of the fibrohexamerin gene (fhx), that encodes a subunit of fibroin, the major silk protein of the silkworm Bombyx mori. The fusion gene was inserted into a piggyBac vector to establish a series of transgenic lines. The expression of the transgene was monitored during the course of larval life and was found restricted to the posterior silk gland cells as the endogenous fhx gene, in all the selected transgenic lines. The exogenous polypeptide was secreted into the lumen of the posterior silk gland together with fibroin, and further exported with the silk proteins as a foreign constituent of the cocoon fiber. The capacity of DsRed to emit fluorescence in the air-dried silk thread led to show that the recombinant protein was distributed over the whole length of the fiber. A remarkable property of the system lies in the localization of the globular protein at the periphery of the silk thread, allowing its rapid and easy recovery in aqueous solutions, without dissolving fibroin. The procedure represents a novel and promising strategy for the production of massive recombinant proteins of biomedical and pharmaceutical interest, with reduced cost.

Generating phenotypic variation: prospects from "evo-devo" research on Bicyclus anynana wing patterns

Understanding the generation of phenotypic variation is an important challenge for modern evolutionary biology, and butterfly wing patterns are an exciting system that can shed some light on this issue. Here, we report on recent advances in the genetics of Bicyclus anynana butterflies. This system provides the potential for a fully integrated study of the evolutionary and developmental processes underlying diversity in morphology.

Conserved genomic organisation of Group B Sox genes in insects

Background

Sox domain containing genes are important metazoan transcriptional regulators implicated in a wide rage of developmental processes. The vertebrate B subgroup contains the Sox1, Sox2 and Sox3 genes that have early functions in neural development. Previous studies show that Drosophila Group B genes have been functionally conserved since they play essential roles in early neural specification and mutations in the Drosophila Dichaete and SoxN genes can be rescued with mammalian Sox genes. Despite their importance, the extent and organisation of the Group B family in Drosophila has not been fully characterised, an important step in using Drosophila to examine conserved aspects of Group B Sox gene function.

Results

We have used the directed cDNA sequencing along with the output from the publicly-available genome sequencing projects to examine the structure of Group B Sox domain genes in Drosophila melanogaster, Drosophila pseudoobscura, Anopheles gambiae and Apis mellifora. All of the insect genomes contain four genes encoding Group B proteins, two of which are intronless, as is the case with vertebrate group B genes. As has been previously reported and unusually for Group B genes, two of the insect group B genes, Sox21a and Sox21b, contain introns within their DNA-binding domains. We find that the highly unusual multi-exon structure of the Sox21b gene is common to the insects. In addition, we find that three of the group B Sox genes are organised in a linked cluster in the insect genomes. By in situ hybridisation we show that the pattern of expression of each of the four group B genes during embryogenesis is conserved between D. melanogaster and D. pseudoobscura.

Conclusion

The DNA-binding domain sequences and genomic organisation of the group B genes have been conserved over 300 My of evolution since the last common ancestor of the Hymenoptera and the Diptera. Our analysis suggests insects have two Group B1 genes, SoxN and Dichaete, and two Group B2 genes. The genomic organisation of Dichaete and another two Group B genes in a cluster, suggests they may be under concerted regulatory control. Our analysis suggests a simple model for the evolution of group B Sox genes in insects that differs from the proposed evolution of vertebrate Group B genes.

Connecting the Navigational Clock to Sun Compass Input in Monarch Butterfly Brain

Migratory monarch butterflies (Danaus plexippus) use a time-compensated sun compass to navigate to their overwintering grounds in Mexico. Although polarized light is one of the celestial cues used for orientation, the spectral content (color) of that light has not been fully explored. We cloned the cDNAs of three visual pigment-encoding opsins (ultraviolet [UV], blue, and long wavelength) and found that all three are expressed uniformly in main retina. The photoreceptors of the polarization-specialized dorsal rim area, on the other hand, are monochromatic for the UV opsin. Behavioral studies support the importance of polarized UV light for flight orientation. Next, we used clock protein expression patterns to identify the location of a circadian clock in the dorsolateral protocerebrum of butterfly brain. To provide a link between the clock and the sun compass, we identified a CRYPTOCHROME-staining neural pathway that likely connects the circadian clock to polarized light input entering brain.

Detection and analysis of alternative splicing in the silkworm by aligning expressed sequence tags with the genomic sequence

We identified 277 alternative splice forms in silkworm genes based on aligning expressed sequence tags with genomic sequences, using a transcript assembly program. A large fraction (74%) of these alternative splices are located in protein-coding regions and alter protein products, whereas only 26% are in untranslated regions. From the alternative splices located in protein-coding regions, some (43%) affect protein domains that bind various biological molecules. The vast majority of the detected alternative forms in this study appear to be novel, and potentially affect biologically meaningful control of function in silkworm genes. Our results indicate that alternative splicing in silkworm largely produces protein diversity and functional diversity, and is a widely used mechanism for regulating gene expression.

The Bombyx mori karyotype and the assignment of linkage groups

Lepidopteran species have a relatively high number of small holocentric chromosomes (Bombyx mori, 2n = 56). Chromosome identification has long been hampered in this group by the high number and by the absence of suitable markers like centromere position and chromosome bands. In this study, we carried out fluorescence in situ hybridization (FISH) on meiotic chromosome complements using genetically mapped B. mori bacterial artificial chromosomes (BACs) as probes. The combination of two to four either green or red fluorescence-labeled probes per chromosome allowed us to recognize unequivocally each of the 28 bivalents of the B. mori karyotype by its labeling pattern. Each chromosome was assigned one of the already established genetic linkage groups and the correct orientation in the chromosome was defined. This facilitates physical mapping of any other sequence and bears relevance for the ongoing B. mori genome projects. Two-color BAC-FISH karyotyping overcomes the problem of chromosome recognition in organisms where conventional banding techniques are not available.

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.