Genomic sequence of a 320-kb segment of the Z chromosome of Bombyx mori containing a kettin ortholog

Sex Chromosome Dosage Compensation in Insects

Dosage compensation (DC) is of crucial importance in balancing the sex-linked gene expression between males and females. It serves to guarantee that the proteins or other enzymatic products encoded by the sex chromosome exhibit quantitative parity between the two genders. During the evolutionary process of achieving dose compensation, insects have developed a wide variety of mechanisms. There exist two primary modes of dosage compensation mechanisms, including the up-regulation of heterogametic sex chromosomes in the heterogamety and down-regulation of homogametic sex chromosomes in the homogamety. Although extensive investigations have been conducted on dosage compensation in model insects, many questions still remain unresolved. Meanwhile, research on non-model insects is attracting increasing attention. This paper systematically summarizes the current advances in the field of insect dosage compensation with respect to its types and mechanisms. The principal insects involved in this study include the Drosophila melanogaster, Tribolium castaneum, Bombyx mori, and other lepidopteran insects. This paper analyzes the controversial issues about insect dosage compensation and also provides prospects for future research.

Co-Expression Network and Time-Course Expression Analyses to Identify Silk Protein Regulatory Factors in Bombyx mori

Simple Summary

Previous studies have reported how the silk production ability of Bombyx mori can be enhanced, but the mechanism that regulates silk protein genes remains unclear. We performed co-expression network analysis using networkz, an in-house program, which led to the identification of 91 transcription factors were co-expressed with silk protein genes. Of them, 13 transcripts were identified to be novel regulatory factors by time-course expression analysis during the fifth instar larvae stage. Their expression patterns were highly relevant to those of silk protein genes. Our results suggest that the two-step expression screening was robust and highly sensitive to screen relative genes, and a complex mechanism regulates silk protein production in B. mori. The novel candidates that were identified herein can serve as key genes to develop methods to enhance the silk protein production ability of B. mori.

Abstract

Bombyx mori is an important economic insect and an animal model in pharmacomedical research. Although its physiology has been studied for many years, the mechanism via which silk protein genes are regulated remains unclear. In this study, we performed two-step expression screening, namely co-expression network and time-course expression analyses to screen silk protein regulation factors. A co-expression network analysis using RNA-seq data that were obtained from various tissues, including the silk glands of B. mori, was performed to identify novel silk protein regulatory factors. Overall, 91 transcription factors, including some known ones, were found to be co-expressed with silk protein genes. Furthermore, time-course expression analysis during the fifth instar larvae stage revealed that the expression pattern of 13 novel transcription factors was highly relevant to that of silk protein genes and their known regulatory factor genes. In particular, the expression peak of several transcription factors (TFs) was detected before the expression of silk protein genes peak. These results indicated that a larger number of genes than expected may be involved in silk protein regulation in B. mori. Functional analyses of function-unknown transcription factors should enhance our understanding of this system.

High-Quality de novo Chromosome-Level Genome Assembly of a Single Bombyx mori With BmNPV Resistance by a Combination of PacBio Long-Read Sequencing, Illumina Short-Read Sequencing, and Hi-C Sequencing

The reference genomes of Bombyx mori (B. mori), Silkworm Knowledge-based database (SilkDB) and SilkBase, have served as the gold standard for nearly two decades. Their use has fundamentally shaped model organisms and accelerated relevant studies on lepidoptera. However, the current reference genomes of B. mori do not accurately represent the full set of genes for any single strain. As new genome-wide sequencing technologies have emerged and the cost of high-throughput sequencing technology has fallen, it is now possible for standard laboratories to perform full-genome assembly for specific strains. Here we present a high-quality de novo chromosome-level genome assembly of a single B. mori with nuclear polyhedrosis virus (BmNPV) resistance through the integration of PacBio long-read sequencing, Illumina short-read sequencing, and Hi-C sequencing. In addition, regular bioinformatics analyses, such as gene family, phylogenetic, and divergence analyses, were performed. The sample was from our unique B. mori species (NB), which has strong inborn resistance to BmNPV. Our genome assembly showed good collinearity with SilkDB and SilkBase and particular regions. To the best of our knowledge, this is the first genome assembly with BmNPV resistance, which should be a more accurate insect model for resistance studies.

dmrt11E ortholog is a crucial factor for oogenesis of the domesticated silkworm, Bombyx mori

DMRT (Doublesex and Mab-3-related transcription factor) is a highly conserved transcription factor family involved in sex determination in numerous animal species. One DMRT, dmrt2/dmrt11E, has entirely different functions in invertebrate and vertebrate species, indicating unpredicted functions. Here, we performed functional analysis of the dmrt11E gene in the domesticated silkworm, Bombyx mori. This gene was preferentially expressed in ovarioles at the last larval instar stage. Its mRNA accumulated in ovarian eggs during the adult stage. CRISPR/Cas9-mediated knockout of Bombyx dmrt11E (Bmdmrt11E) caused defects in oogenesis, resulting in the production of abnormal eggs with transparent liquids. These eggs had significantly reduced fertility and lipid levels. Transcriptomic comparisons between ovaries of control and mutant insects at two developmental stages identified six genes that may be under the control of Bmdmrt11E. Finally, we provide a possible model for lipid uptake and storage in eggs of Bombyx mori.

Transgenic and knockout analyses of Masculinizer and doublesex illuminated the unique functions of doublesex in germ cell sexual development of the silkworm, Bombyx mori

Background

Masculinizer (Masc) plays a pivotal role in male sex determination in the silkworm, Bombyx mori. Masc is required for male-specific splicing of B. mori doublesex (Bmdsx) transcripts. The male isoform of Bmdsx (BmdsxM) induces male differentiation in somatic cells, while females express the female isoform of Bmdsx (BmdsxF), which promotes female differentiation in somatic cells. Our previous findings suggest that Masc could direct the differentiation of genetically female (ZW) germ cells into sperms. However, it remains unclear whether Masc directly induces spermatogenesis or if it promotes male differentiation in germ cells indirectly by inducing the expression of BmdsxM.

Results

In this study, we performed genetic analyses using the transgenic line that expressed Masc, as well as various Bmdsx knockout lines. We found that Masc-expressing females with a homozygous mutation in BmdsxM showed normal development in ovaries. The formation of testis-like tissues was abolished in these females. On the other hand, Masc-expressing females carrying a homozygous mutation in BmdsxF exhibited almost complete male-specific development in gonads and germ cells. These results suggest that BmdsxM has an ability to induce male development in germ cells as well as internal genital organs, while BmdsxF inhibits BmdsxM activity and represses male differentiation. To investigate whether MASC directly controls male-specific splicing of Bmdsx and identify RNAs that form complexes with MASC in testes, we performed RNA immunoprecipitation (RIP) using an anti-MASC antibody. We found that MASC formed a complex with AS1 lncRNA, which is a testis-specific factor involved in the male-specific splicing of Bmdsx pre-mRNA.

Conclusions

Taken together, our findings suggest that Masc induces male differentiation in germ cells by enhancing the production of BmdsxM. Physical interaction between MASC and AS1 lncRNA may be important for the BmdsxM expression in the testis. Unlike in the Drosophila dsx, BmdsxM was able to induce spermatogenesis in genetically female (ZW) germ cells. To the best of our knowledge, this is the first report that the role of dsx in germ cell sexual development is different between insect species.

Transcriptome analysis of genes related to gonad differentiation and development in Muscovy ducks

BACKGROUND

Sex-related genes play a crucial role in gonadal differentiation into testes or ovaries. However, the genetic control of gonadal differentiation in Muscovy ducks remains unknown. Therefore, the objective of our study was to screen new candidate genes associated with ovarian and testicular development.

RESULTS

In this study, 24 males before gonadal differentiation (MB), 24 females before gonadal differentiation (FB), 24 males after gonadal differentiation (MA) and 24 females after gonadal differentiation (FA) were selected from Putian Muscovy ducks, forming 4 groups. RNA-Seq revealed 101.76 Gb of clean reads and 2800 differentially expressed genes (DEGs), including 46 in MB vs FB, 609 in MA vs FA, 1027 in FA vs FB, and 1118 in MA vs MB. A total of 146 signalling pathways were enriched by KEGG analysis, among which 20, 108, 108 and 116 signalling pathways were obtained in MB vs FB, MA vs MB, MA vs FA and FA vs FB, respectively. In further GO and KEGG analyses, a total of 21 candidate genes related to gonad differentiation and development in Muscovy ducks were screened. Among these, 9 genes were involved in the differentiation and development of the testes, and 12 genes were involved in the differentiation and development of the ovaries. In addition, RNA-Seq data revealed 2744 novel genes.

CONCLUSIONS

RNA-Seq data revealed 21 genes related to gonadal differentiation and development in Muscovy ducks. We further identified 12 genes, namely, WNT5B, HTRA3, RSPO3, BMP3, HNRNPK, NIPBL, CREB3L4, DKK3, UBE2R2, UBPL3KCMF1, ANXA2, and OSR1, involved in the differentiation and development of ovaries. Moreover, 9 genes, namely, TTN, ATP5A1, DMRT1, DMRT3, AMH, MAP3K1, PIK3R1, AGT and ADAMTSL1, were related to the differentiation and development of testes. Moreover, after gonadal differentiation, DMRT3, AMH, PIK3R1, ADAMTSL1, AGT and TTN were specifically highly expressed in males. WNT5B, ANXA2 and OSR1 were specifically highly expressed in females. These results provide valuable information for studies on the sex control of Muscovy ducks and reveal novel candidate genes for the differentiation and development of testes and ovaries.

Influence of Diet, Sex, and Viral Infections on the Gut Microbiota Composition of Spodoptera exigua Caterpillars

The gut microbiota plays essential roles in processes related with metabolism, physiology, and immunity in all organisms, including insects. In the present work, we performed a broad analysis of the Spodoptera exigua gut microbiota, a major agricultural pest. We analyzed the influence of multiple parameters such as diet, geographic location, sex, or viral infections on S. exigua caterpillar gut microbiota composition. Our study revealed a high variability in bacterial composition among individuals, and a major influence of environmental bacteria (including those acquired through diet) on the gut microbiota composition, supporting previous studies that claim resident microbiota are lacking in caterpillars. Previous studies with laboratory-reared insects showed that changes in caterpillar gut bacterial composition affect the insecticidal properties of entomopathogenic viruses and bacteria. Our study revealed different microbiota composition in field insects carrying a natural viral infection with Spodoptera exigua nucleopolyhedrovirus (SeMNPV) and/or Spodoptera exigua iflavirus 1 (SeIV1). Few taxa can be specifically associated with the infection, suggesting microbiota influence the infective process of these natural pathogens, and providing new strategies for insect pest management.

Estimation of long terminal repeat element content in the Helicoverpa zea genome from high-throughput sequencing of bacterial artificial chromosome pools

The lepidopteran pest insect Helicoverpa zea feeds on cultivated corn and cotton across the Americas where control remains challenging owing to the evolution of resistance to chemical and transgenic insecticidal toxins, yet genomic resources remain scarce for this species. A bacterial artificial chromosome (BAC) library having a mean genomic insert size of 145 ± 20 kbp was created from a laboratory strain of H. zea, which provides ∼12.9-fold coverage of a 362.8 ± 8.8 Mbp (0.37 ± 0.09 pg) flow cytometry estimated haploid genome size. Assembly of Illumina HiSeq 2000 reads generated from 14 pools that encompassed all BAC clones resulted in 165 485 genomic contigs (N₅₀ = 3262 bp; 324.6 Mbp total). Long terminal repeat (LTR) protein coding regions annotated from 181 contigs included 30 Ty1/copia, 78 Ty3/gypsy, and 73 BEL/Pao elements, of which 60 (33.1%) encoded all five functional polyprotein (pol) domains. Approximately 14% of LTR elements are distributed non-randomly across pools of BAC clones.

Transgenic Expression of the piRNA-Resistant Masculinizer Gene Induces Female-Specific Lethality and Partial Female-to-Male Sex Reversal in the Silkworm, Bombyx mori

In Bombyx mori (B. mori), Fem piRNA originates from the W chromosome and is responsible for femaleness. The Fem piRNA-PIWI complex targets and cleaves mRNAs transcribed from the Masc gene. Masc encodes a novel CCCH type zinc-finger protein and is required for male-specific splicing of B. mori doublesex (Bmdsx) transcripts. In the present study, several silkworm strains carrying a transgene, which encodes a Fem piRNA-resistant Masc mRNA (Masc-R), were generated. Forced expression of the Masc-R transgene caused female-specific lethality during the larval stages. One of the Masc-R strains weakly expressed Masc-R in various tissues. Females heterozygous for the transgene expressed male-specific isoform of the Bombyx homolog of insulin-like growth factor II mRNA-binding protein (ImpM) and Bmdsx. All examined females showed a lower inducibility of vitellogenin synthesis and exhibited abnormalities in the ovaries. Testis-like tissues were observed in abnormal ovaries and, notably, the tissues contained considerable numbers of sperm bundles. Homozygous expression of the transgene resulted in formation of the male-specific abdominal segment in adult females and caused partial male differentiation in female genitalia. These results strongly suggest that Masc is an important regulatory gene of maleness in B. mori.

Construction, complete sequence, and annotation of a BAC contig covering the silkworm chorion locus

The silkmoth chorion was studied extensively by F.C. Kafatos’ group for almost 40 years. However, the complete structure of the chorion locus was not obtained in the genome sequence of Bombyx mori published in 2008 due to repetitive sequences, resulting in gaps and an incomplete view of the locus. To obtain the complete sequence of the chorion locus, expressed sequence tags (ESTs) derived from follicular epithelium cells were used as probes to screen a bacterial artificial chromosome (BAC) library. Seven BACs were selected to construct a contig which covered the whole chorion locus. By Sanger sequencing, we successfully obtained complete sequences of the chorion locus spanning 871,711 base pairs on chromosome 2, where we annotated 127 chorion genes. The dataset reported here will recruit more researchers to revisit one of the oldest model systems which has been used to study developmentally regulated gene expression. It also provides insights into egg development and fertilization mechanisms and is relevant to applications related to improvements in breeding procedures and transgenesis.

A comprehensive analysis of the chorion locus in silkmoth

Despite more than 40 years of intense study, essential features of the silkmoth chorion (eggshell) are still not fully understood. To determine the precise structure of the chorion locus, we performed extensive EST analysis, constructed a bacterial artificial chromosome (BAC) contig, and obtained a continuous genomic sequence of 871,711 base pairs. We annotated 127 chorion genes in two segments interrupted by a 164 kb region with 5 non-chorion genes, orthologs of which were on chorion bearing scaffolds in 4 ditrysian families. Detailed transcriptome analysis revealed expression throughout choriogenesis of most chorion genes originally categorized as “middle”, and evidence for diverse regulatory mechanisms including cis-elements, alternative splicing and promoter utilization, and antisense RNA. Phylogenetic analysis revealed multigene family associations and faster evolution of early chorion genes and transcriptionally active pseudogenes. Proteomics analysis identified 99 chorion proteins in the eggshell and micropyle localization of 1 early and 6 Hc chorion proteins.

Sexually biased transcripts at early embryonic stages of the silkworm depend on the sex chromosome constitution

In the silkworm, Bombyx mori, females are heterogametic (WZ) whereas males have two Z chromosomes. Femaleness of B. mori is determined by the presence of the W chromosome, suggesting that there is a dominant feminizing gene on this chromosome. Recently, by transcriptome analysis of B. mori embryos, we discovered that a single W-chromosome-derived PIWI-interacting RNA (piRNA) is the long-sought primary determinant of femaleness in B. mori. However, sexual bias in the transcriptome of B. mori early embryos has not yet been well characterized. Using deep sequencing data from molecularly sexed RNA of B. mori embryos, we identified and characterized 157 transcripts that are statistically differentially expressed between male and female early embryos. Most of the female-biased transcripts were transposons or repeat sequences that are produced presumably from the W chromosome. Bioinformatic analysis revealed that these repetitive sequences are piRNA precursors. In contrast, male-biased genes were frequently transcribed from the Z chromosome, suggesting that dosage compensation in Z-linked genes does not occur or is incomplete at early embryonic stages. Our analysis has drawn a picture of a global landscape of sexually biased transcriptome during early B. mori embyogenesis and has suggested for the first time that most sexually biased embryonic transcripts depend on sex chromosomes.

Effects of RNAi-mediated knockdown of histone methyltransferases on the sex-specific mRNA expression of Imp in the silkworm Bombyx mori

Sexual differentiation in Bombyx mori is controlled by sex-specific splicing of Bmdsx, which results in the omission of exons 3 and 4 in a male-specific manner. In B. mori, insulin-like growth factor II mRNA-binding protein (Imp) is a male-specific factor involved in male-specific splicing of Bmdsx. Male-specific Imp mRNA results from the male-specific inclusion of exon 8. To verify the link between histone methylation and alternative RNA processing in Imp, we examined the effects of RNAi-mediated knockdown of several histone methyltransferases on the sex-specific mRNA expression of Imp. As a result, male-specific expression of Imp mRNA was completely abolished when expression of the H3K79 methyltransferase DOT1L was repressed to <10% of that in control males. Chromatin immunoprecipitation-quantitative PCR analysis revealed a higher distribution of H3K79me2 in normal males than in normal females across Imp. RNA polymerase II (RNAP II) processivity assays indicated that RNAi knockdown of DOT1L in males caused a twofold decrease in RNAP II processivity compared to that in control males, with almost equivalent levels to those observed in normal females. Inhibition of RNAP II-mediated elongation in male cells repressed the male-specific splicing of Imp. Our data suggest the possibility that H3K79me2 accumulation along Imp is associated with the male-specific alternative processing of Imp mRNA that results from increased RNAP II processivity.

A single sex-linked dominant gene does not fully explain the codling moth's resistance to granulovirus

BACKGROUND

In 2004, resistance to a commercial formulation of the Cydia pomonella granulovirus (CpGV) was identified in a field population of Cydia pomonella from an organic orchard in southern France. The genetic inheritance of this resistance was analysed in the resistant laboratory strain RGV. This strain was obtained using successive crosses between the resistant field population and a susceptible laboratory strain, SV, with selection for CpGV resistance at each generation.

RESULTS

After eight generations of introgression of the resistant trait into SV, the RGV-8 strain exhibited 7000-fold higher resistance than SV. Mass-crossing experiments showed that resistance to CpGV is strongly dominant, sex dependent and under the control of a single major gene. However, the contribution of other genes is required to explain all of the data obtained in this study. These additional genes do not follow the laws of classical Mendelian transmission.

CONCLUSION

Transmission of granulovirus resistance in the RGV-8 strain of C. pomonella cannot be fully explained by the effect of a locus located on the Z chromosome. The action of other factors needs to be considered.

Molecular analysis of the muscle protein projectin in Lepidoptera

Striated muscles of both vertebrates and insects contain a third filament composed of the giant proteins, namely kettin and projectin (insects) and titin (vertebrates). All three proteins have been shown to contain several domains implicated in conferring elasticity, in particular a PEVK segment. In this study, the characterization of the projectin protein in the silkmoth, Bombyx mori L. (Lepidoptera: Bombycidae), and the monarch butterfly, Danaus plexippus L. (Lepidoptera: Nymphalidae), as well as a partial characterization in the Carolina sphinx, Manduca sexta L. (Lepidoptera: Sphingidae), are presented. This study showed that, similar to other insects, projectin's overall modular organization was conserved, but in contrast, the PEVK region had a highly divergent sequence. The analysis of alternative splicing in the PEVK region revealed a small number of possible isoforms and the lack of a flight-muscle specific variant, both characteristics being in sharp contrast with findings from other insects. The possible correlation with difference in flight muscle stiffness and physiology between Lepidoptera and other insect orders is discussed.

CD36 homolog divergence is responsible for the selectivity of carotenoid species migration to the silk gland of the silkworm Bombyx mori

Dietary carotenoids are absorbed in the intestine and delivered to various tissues by circulating lipoproteins; however, the mechanism underlying selective delivery of different carotenoid species to individual tissues remains elusive. The products of the Yellow cocoon (C) gene and the Flesh (F) gene of the silkworm Bombyx mori determine the selectivity for transport of lutein and β-carotene, respectively, to the silk gland. We previously showed that the C gene encodes Cameo2, a CD36 family member, which is thought to function as a transmembrane lipoprotein receptor. Here, we elucidated the molecular identity of the F gene product by positional cloning, as SCRB15, a paralog of Cameo2 with 26% amino acid identity. In the F mutant, SCRB15 mRNA structure was severely disrupted, due to a 1.4 kb genomic insertion in a coding exon. Transgenic expression of SCRB15 in the middle silk gland using the binary GAL4-UAS expression system enhanced selective β-carotene uptake by the middle silk gland, while transgenic expression of Cameo2 enhanced selective lutein uptake under the same GAL4 driver. Our findings indicate that divergence of genes in the CD36 family determines the selectivity of carotenoid species uptake by silk gland tissue and that CD36-homologous proteins can discriminate among carotenoid species.

Retrogenes moved out of the z chromosome in the silkworm

Previous studies on organisms with well-differentiated X and Y chromosomes, such as Drosophila and mammals, consistently detected an excess of genes moving out of the X chromosome and gaining testis-biased expression. Several selective evolutionary mechanisms were shown to be associated with this nonrandom gene traffic, which contributed to the evolution of the X chromosome and autosomes. If selection drives gene traffic, such traffic should also exist in species with Z and W chromosomes, where the females are the heterogametic sex. However, no previous studies on gene traffic in species with female heterogamety have found any nonrandom chromosomal gene movement. Here, we report an excess of retrogenes moving out of the Z chromosome in an organism with the ZW sex determination system, Bombyx mori. In addition, we showed that those "out of Z" retrogenes tended to have ovary-biased expression, which is consistent with the pattern of non-retrogene traffic recently reported in birds and symmetrical to the retrogene movement in mammals and fruit flies out of the X chromosome evolving testis functions. These properties of gene traffic in the ZW system suggest a general role for the heterogamety of sex chromosomes in determining the chromosomal locations and the evolution of sex-biased genes.

Identification and expression analysis of nervous wreck, which is preferentially expressed in the brain of the male silkworm moth, Bombyx mori

Sexually dimorphic neural circuits are essential for reproductive behaviour. The molecular basis of sexual dimorphism in the silkworm moth (Bombyx mori) brain, however, is unclear. We conducted cDNA subtraction screening and identified nervous wreck (Bmnwk), a synaptic growth regulatory gene, whose expression is higher in the male brain than in the female brain of the silkworm. Bmnwk was preferentially expressed in the brain at the late pupae and adult stages. In situ hybridization revealed that Bmnwk is highly expressed in the optic lobe of the male moth brain. These findings suggest that Bmnwk has a role in the development and/or maintenance of the optic lobe in the male silkworm brain.

Amphioxus connectin exhibits merged structure as invertebrate connectin in I-band region and vertebrate connectin in A-band region

Connectin is an elastic protein found in vertebrate striated muscle and in some invertebrates as connectin-like proteins. In this study, we determined the structure of the amphioxus connectin gene and analyzed its sequence based on its genomic information. Amphioxus is not a vertebrate but, phylogenetically, the lowest chordate. Analysis of gene structure revealed that the amphioxus gene is approximately 430 kb in length and consists of regions with exons of repeatedly aligned immunoglobulin (Ig) domains and regions with exons of fibronectin type 3 and Ig domain repeats. With regard to this sequence, although the region corresponding to the I-band is homologous to that of invertebrate connectin-like proteins and has an Ig-PEVK region similar to that of the Neanthes sp. 4000K protein, the region corresponding to the A-band has a super-repeat structure of Ig and fibronectin type 3 domains and a kinase domain near the C-terminus, which is similar to the structure of vertebrate connectin. These findings revealed that amphioxus connectin has the domain structure of invertebrate connectin-like proteins at its N-terminus and that of vertebrate connectin at its C-terminus. Thus, amphioxus connectin has a novel structure among known connectin-like proteins. This finding suggests that the formation and maintenance of the sarcomeric structure of amphioxus striated muscle are similar to those of vertebrates; however, its elasticity is different from that of vertebrates, being more similar to that of invertebrates.

Molecular analysis of sex chromosome-linked mutants in the silkworm Bombyx mori

In Bombyx mori, the W chromosome determines the female sex. A few W chromosome-linked mutations that cause masculinization of the female genitalia have been found. In female antennae of a recently isolated mutant, both female-type and male-type Bmdsx mRNAs were expressed, and BmOr1 (bombykol receptor) and BmOr3 (bombykal receptor), which are predominantly expressed in the antennae of male moths, were expressed about 50 times more abundantly in the antennae of mutant females than in those of normal females. These mutants are valuable resources for the molecular analysis of the sexdetermination system. Besides the Fem gene, the quantitative egg size-determining gene Esd is thought to be present on the W chromosome, based on the observation that ZWW triploid moths produce larger eggs than ZZW triploids. The most recently updated B. mori genome assembly comprises 20.5 Mb of Z chromosome sequence. Using these sequence data, responsible genes or candidate genes for four Z-linked mutants have been reported. The od (distinct oily) and spli (soft and pliable) are caused by mutation in BmBLOS2 and Bmacj6, respectively. Bmap is a candidate gene for Vg (vestigial). Similarly, Bmprm is a candidate gene for Md (muscle dystrophy), causing abnormal development of indirect flight muscle.

A unique red fluorescent protein of silkworm bearing two photochromic moieties

A silkworm excretory red fluorescent protein (SE-RFP) having light-dependent activity against Bombyx mori nucleopolyhedrovirus (BmNPV) was purified. Light was observed to be essential also for the SE-RFP synthesis as it was produced only when silkworms were reared in light. SE-RFP has exhibited a high fluorescence quantum yield of 0.86. The apparent mass of native SE-RFP was about 1100 kDa as analysed by gel filtration chromatography. Two photochromic moieties associated with the SE-RFP, namely tetrapyrrole-I (TP-I) and tetrapyrrole-II (TP-II), were isolated by employing TLC and HPTLC techniques. The purified tetrapyrroles were characterized by UV-absorption, fluorescence, atomic absorption and FT-IR spectral analyses. The molecular masses of TP-I and TP-II were 535 and 870 Da, respectively, as determined by ESI-MS and MALDI-TOF-MS. The molar ratio of TP-I to TP-II was 1.14 : 1.00, and a total of 7.251 micromol tetrapyrroles (TP-I + TP-II) were found to be present per mg of SE-RFP. TP-I and TP-II were identified as chlorophyll derivatives, namely, pyropheophorbide a and pheophytin a, respectively. Hence, the SE-RFP was concluded to be a unique insect red fluorescent protein having two photochromic moieties and potent photobiological activity.

Dosage analysis of Z chromosome genes using microarray in silkworm, Bombyx mori

In many organisms, dosage compensation is needed to equalize sex-chromosome gene expression in males and females. Several genes on silkworm Z chromosome were previously detected to show a higher expression level in males and lacked dosage compensation. Whether silkworm lacks global dosage compensation still remains poorly known. Here, we analyzed male:female (M:F) ratios of expression of chromosome-wide Z-linked genes in the silkworm using microarray data. The expression levels of genes on Z chromosome in each tissue were significantly higher in males compared to females, which indicates no global dosage compensation in silkworm. Interestingly, we also found some genes with no bias (M:F ratio: 0.8-1.2) on the Z chromosome. Comparison of male-biased (M:F ratio more than 1.5) and unbiased genes indicated that the two sets of the genes have functional differences. Analysis of gene expression by sex showed that M:F ratios were, to some extent, associated with their expression levels. These results provide useful clues to further understanding roles of dosage of Z chromosome and some Z-linked sexual differences in silkworms.

The silkworm Z chromosome is enriched in testis-specific genes

The role of sex chromosomes in sex determination has been well studied in diverse groups of organisms. However, the role of the genes on the sex chromosomes in conferring sexual dimorphism is still being experimentally evaluated. An unequal complement of sex chromosomes between two sexes makes them amenable to sex-specific evolutionary forces. Sex-linked genes preferentially expressed in one sex over the other offer a potential means of addressing the role of sex chromosomes in sexual dimorphism. We examined the testis transcriptome of the silkworm, Bombyx mori, which has a ZW chromosome constitution in the female and ZZ in the male, and show that the Z chromosome harbors a significantly higher number of genes expressed preferentially in testis compared to the autosomes. We hypothesize that sexual antagonism and absence of dosage compensation have possibly led to the accumulation of many male-specific genes on the Z chromosome. Further, our analysis of testis-specific paralogous genes suggests that the accumulation on the Z chromosome of genes advantageous to males has occurred primarily by translocation or tandem duplication.

RECK negatively regulates matrix metalloproteinase-9 transcription

RECK, a glycosylphosphatidylinositol-anchored glycoprotein, inhibits the enzymatic activities of some matrix metalloproteinases (MMP), thereby suppressing tumor cell metastasis; however, the detailed mechanism is still obscure. In this study, we compared the gene expression profiles between mock- and RECK-transfected HT1080 cells and showed that RECK decreases MMP-9 mRNA levels but not other MMP mRNA levels. Moreover, treatment with RECK-specific siRNA increased MMP-9 mRNA in RECK-expressing cells. The promoter assay showed that MMP-9 promoter activity was suppressed by RECK and that RECK-mediated suppression of MMP-9 promoter activity requires 12-O-tetradecanoylphorbol-13-acetate-responsive element (TRE) and kappaB sites. Moreover, the binding ability of Fra-1 and c-Jun to TRE within the MMP-9 promoter region was suppressed by RECK. Thus, these results show that RECK is a negative regulator of MMP-9 transcription.

A bird's-eye view of sex chromosome dosage compensation

Intensive study of a few genetically tractable species with XX/XY sex chromosomes has produced generalizations about the process of sex chromosome dosage compensation that do not fare well when applied to ZZ/ZW sex chromosome systems, such as those in birds. The inherent sexual imbalance in dose of sex chromosome genes has led to the evolution of sex-chromosome-wide mechanisms for balancing gene dosage between the sexes and relative to autosomal genes. Recent advances in our knowledge of avian genomes have led to a reexamination of sex-specific dosage compensation (SSDC) in birds, which is less effective than in known XX/XY systems. Insights about the mechanisms of SSDC in birds also suggest similarities to and differences from those in XX/XY species. Birds are thus offering new opportunities for studying dosage compensation in a ZZ/ZW system, which should shed light on the evolution of SSDC more broadly.

Diversity of aminopeptidases, derived from four lepidopteran gene duplications, and polycalins expressed in the midgut of Helicoverpa armigera: identification of proteins binding the delta-endotoxin, Cry1Ac of Bacillus thuringiensis

Helicoverpa armigera midgut proteins that bind the Bacillus thuringiensis (Bt) delta-endotoxin Cry1Ac were purified by affinity chromatography. SDS-PAGE showed that several proteins were eluted with N-acetylgalactosamine and no further proteins were detected after elution with urea. Tandem mass spectral data for tryptic peptides initially indicated that the proteins resembled aminopeptidases (APNs) from other lepidopterans and cDNA sequences for seven APNs were isolated from H. armigera through a combination of cloning with primers derived from predicted peptide sequences and established EST libraries. Phylogenetic analysis showed lepidopteran APN genes in nine clades of which five were part of a lepidopteran-specific radiation. The Cry1Ac-binding proteins were then identified with four of the seven HaAPN genes. Three of those four APNs are likely orthologs of APNs characterised as Cry1Ac-binding proteins in other lepidopterans. The fourth Cry1Ac-binding APN has orthologs not previously identified as Cry1Ac-binding partners. The HaAPN genes were expressed predominantly in the midgut through larval development. Each showed consistent expression along the length of the midgut but five of the genes were expressed at levels about two orders of magnitude greater than the remaining two. The remaining mass spectral data identified sequences encoding polycalin proteins with multiple lipocalin-like domains. A polycalin has only been previously reported in another lepidopteran, Bombyx mori, but polycalins in both species are now linked with binding of Bt Cry toxins. This is the first report of hybrid, lipocalin-like domains in shorter polycalin sequences that are not present in the longest sequence. We propose that these hybrid domains are generated by alternative splicing of the mRNA.

Positional cloning of a Bombyx wingless locus flugellos (fl) reveals a crucial role for fringe that is specific for wing morphogenesis

Mutations at the flügellos (fl) locus in Bombyx mori produce wingless pupae and moths because of the repressed response of wing discs to ecdysteroid. Four recessive fl alleles occurred spontaneously and were mapped at 13.0 of the silkworm genetic linkage group 10. By positional cloning, we confirmed that the gene responsible for fl is fringe (fng) encoding Fng glycosyltransferase, which is involved in regulating the Notch signaling pathway. In four different fl alleles, we detected a large deletion of the fng gene in fl(k) and nonsense mutations in fl, fl(o), and fl(n). In the wild-type (WT) silkworm, fng is expressed actively in the wing discs, brain, and reproductive organs from the fourth to final instars but barely in the other tissues tested. In situ hybridization showed that fng mRNA is expressed in the dorsal layer of the WT wing discs. The wingless (wg) mRNA, a downstream marker of Fng-mediated Notch signaling, is localized at the dorsoventral boundary in the WT wing discs but repressed markedly in the fl wing discs. Although null mutants of Drosophila fng result in postembryonic lethality, loss of fng function in Bombyx affects only wing morphogenesis, suggesting different essential roles for fng in tissue differentiation among insects.

Deletion of a gene encoding an amino acid transporter in the midgut membrane causes resistance to a Bombyx parvo-like virus

Bombyx mori densovirus type 2 (BmDNV-2), a parvo-like virus, replicates only in midgut columnar cells and causes fatal disease. The resistance expressed in some silkworm strains against the virus is determined by a single gene, nsd-2, which is characterized as nonsusceptibility irrespective of the viral dose. However, the responsible gene has been unknown. We isolated the nsd-2 gene by positional cloning. The virus resistance is caused by a 6-kb deletion in the ORF of a gene encoding a 12-pass transmembrane protein, a member of an amino acid transporter family, and expressed only in midgut. Germ-line transformation with a wild-type transgene expressed in the midgut restores susceptibility, showing that the defective membrane protein is responsible for resistance. Cumulatively, our data show that the membrane protein is a functional receptor for BmDNV-2. This is a previously undescribed report of positional cloning of a mutant gene in Bombyx and isolation of an absolute virus resistance gene in insects.

Hybrid zone origins, species boundaries, and the evolution of wing-pattern diversity in a polytypic species complex of North American admiral butterflies (Nymphalidae: Limenitis)

Hybrid zones present opportunities to study the effects of gene flow, selection, and recombination in natural populations and, thus, provide insights into the genetic and phenotypic changes that occur early in speciation. Here we investigate a hybrid zone between mimetic (Limenitis arthemis astyanax) and nonmimetic (Limenitis arthemis arthemis) populations of admiral butterflies using DNA sequence variation from mtDNA and seven nuclear gene loci. We find three distinct mitochondrial clades within this complex, and observe a strong overall concordance between wing-pattern phenotypes and mitochondrial variation. Nuclear gene genealogies, in contrast, revealed no evidence of exclusivity for either wing-pattern phenotype, suggesting incomplete barriers to gene exchange and/or insufficient time for lineage sorting. Coalescent simulations indicate that gene flow between these two subspecies is highly asymmetric, with the majority of migration occurring from mimetic into nonmimetic populations. Selective sweeps of alleles responsible for mimetic phenotypes may have occurred more than once when mimetic and nonmimetic Limenitis occurred together in the presence of the model (Battus philenor).

A BAC-based integrated linkage map of the silkworm Bombyx mori

An integrated map of the Bombyx mori genome has been constructed using 361.1 Mb of BAC contigs and singletons together with a genetic map containing 1689 independent genes and synteny among Apis, Tribolium, and Bombyx was examined.

Background

In 2004, draft sequences of the model lepidopteran Bombyx mori were reported using whole-genome shotgun sequencing. Because of relatively shallow genome coverage, the silkworm genome remains fragmented, hampering annotation and comparative genome studies. For a more complete genome analysis, we developed extended scaffolds combining physical maps with improved genetic maps.

Results

We mapped 1,755 single nucleotide polymorphism (SNP) markers from bacterial artificial chromosome (BAC) end sequences onto 28 linkage groups using a recombining male backcross population, yielding an average inter-SNP distance of 0.81 cM (about 270 kilobases). We constructed 6,221 contigs by fingerprinting clones from three BAC libraries digested with different restriction enzymes, and assigned a total of 724 single copy genes to them by BLAST (basic local alignment search tool) search of the BAC end sequences and high-density BAC filter hybridization using expressed sequence tags as probes. We assigned 964 additional expressed sequence tags to linkage groups by restriction fragment length polymorphism analysis of a nonrecombining female backcross population. Altogether, 361.1 megabases of BAC contigs and singletons were integrated with a map containing 1,688 independent genes. A test of synteny using Oxford grid analysis with more than 500 silkworm genes revealed six versus 20 silkworm linkage groups containing eight or more orthologs of Apis versus Tribolium, respectively.

Conclusion

The integrated map contains approximately 10% of predicted silkworm genes and has an estimated 76% genome coverage by BACs. This provides a new resource for improved assembly of whole-genome shotgun data, gene annotation and positional cloning, and will serve as a platform for comparative genomics and gene discovery in Lepidoptera and other insects.

Molecular divergence of the W chromosomes in pyralid moths (Lepidoptera)

Most Lepidoptera have a WZ/ZZ sex chromosome system. We compared structure of W chromosomes in four representatives of the family Pyralidae--Ephestia kuehniella, Cadra cautella, Plodia interpunctella, and Galleria mellonella--tracing pachytene bivalents which provide much higher resolution than metaphase chromosomes. In each species, we prepared a W-chromosome painting probe from laser-microdissected W-chromatin of female polyploid nuclei. The Ephestia W-probe was cross-hybridized to chromosomes of the other pyralids to detect common parts of their W chromosomes, while the species-specific W-probes identified the respective W chromosome. This so-called Zoo-FISH revealed a partial homology of W-chromosome regions between E. kuehniella and two other pyralids, C. cautella and P. interpunctella, but almost no homology with G. mellonella. The results were consistent with phylogenetic relationships between the species. We also performed comparative genomic hybridization, which indicated that the W chromosome of C. cautella is composed mainly of repetitive DNA common to both sexes but accumulated in the W chromosome, whereas E. kuehniella, P. interpunctella, and G. mellonella W chromosomes also possess a large amount of female specific DNA sequences, but differently organized. Our results support the hypothesis of the accelerated molecular divergence of the lepidopteran W chromosomes in the absence of meiotic recombination.

Identification of the female-determining region of the W chromosome in Bombyx mori

The W chromosome of the silkworm Bombyx mori is devoid of functional genes, except for the putative female-determining gene (Fem). To localize Fem, we investigated the presence of W-specific DNA markers on strains in which an autosomal fragment containing dominant marker genes was attached to the W chromosome. We produced new W-chromosomal fragments from the existing Zebra-W strain (T(W;3)Ze chromosome) by X-irradiation, and then carried out deletion mapping of these and sex-limited yellow cocoon strains (T(W;2)Y-Chu, -Abe and -Ban types) from different Japanese stock centers. Of 12 RAPD markers identified in the normal W chromosomes of most silkworm strains in Japan, the newly irradiated W(B-YL-YS)Ze chromosome contained three, the T(W;2)Y-Chu chromosome contained six, and the T(W;2)Y-Abe and -Ban chromosomes contained only one (W-Rikishi). To investigate the ability of the reduced W-chromosome translocation fragments to form heterochromatin bodies, which are found in nuclei of normal adult female sucking stomachs, we examined cells of the normal type p50 strain and the T(W;2)Y-Chu and -Abe strains. A single sex heterochromatin body was found in nuclei of p50 females, whereas we detected only small sex heterochromatin bodies in the T(W;2)Y-Chu strain and no sex heterochromatin body in the T(W;2)Y-Abe strain. Since adult females of all strains were normal and fertile, we conclude that only extremely limited region, containing the W-Rikishi RAPD sequence of the W chromosome, is required to determine femaleness. Based on a comparison of the normal W-chromosome and 7 translocation and W-deletion strains we present a map of Fem relative to the 12 W-specific RAPD markers.

End-sequencing and characterization of silkworm (Bombyx mori) bacterial artificial chromosome libraries

Background

We performed large-scale bacterial artificial chromosome (BAC) end-sequencing of two BAC libraries (an EcoRI- and a BamHI-digested library) and conducted an in silico analysis to characterize the obtained sequence data, to make them a useful resource for genomic research on the silkworm (Bombyx mori).

Results

More than 94000 BAC end sequences (BESs), comprising more than 55 Mbp and covering about 10.4% of the silkworm genome, were sequenced. Repeat-sequence analysis with known repeat sequences indicated that the long interspersed nuclear elements (LINEs) were abundant in BamHI BESs, whereas DNA-type elements were abundant in EcoRI BESs. Repeat-sequence analysis revealed that the abundance of LINEs might be due to a GC bias of the restriction sites and that the GC content of silkworm LINEs was higher than that of mammalian LINEs. In a BLAST-based sequence analysis of the BESs against two available whole-genome shotgun sequence data sets, more than 70% of the BESs had a BLAST hit with an identity of ≥ 99%. About 14% of EcoRI BESs and about 8% of BamHI BESs were paired-end clones with unique sequences at both ends. Cluster analysis of the BESs clarified the proportion of BESs containing protein-coding regions.

Conclusion

As a result of this characterization, the identified BESs will be a valuable resource for genomic research on Bombyx mori, for example, as a base for construction of a BAC-based physical map. The use of multiple complementary BAC libraries constructed with different restriction enzymes also makes the BESs a more valuable genomic resource. The GenBank accession numbers of the obtained end sequences are DE283657–DE378560.

Sex determination in the silkworm, Bombyx mori: a female determinant on the W chromosome and the sex-determining gene cascade

In insects, the sex is determined completely by genetic mechanisms, which at least in somatic tissues, are cell autonomous. The sex of the silkworm, Bombyx mori, is strongly controlled by the presence of the W chromosome. Genetic studies using translocations and deletions of W suggested that a presumptive feminizing gene (Fem) is located in a limited region of the W chromosome. Recent genomic studies revealed a small number of potential candidates for the Fem gene in this region. In addition, a Bombyx homologue of the Drosophila sex determining gene doublesex has been identified on an autosome and analyzed. Whereas the Drosophila doublesex gene is regulated by activation of splicing in females, the Bombyx doublesex gene (Bmdsx) encodes female- and male-specific mRNAs regulated via male-specific repression of splicing. The vitellogenin gene (Vg) is a target of the BmDSX protein, which directly binds to the Vg promoter. Furthermore, as ectopic expression of the male-type Bmdsx induces male-like transformation of the sexual organs, BmDSX may control sex-specific morphological characteristics in Bombyx. This suggests that although upstream events in Drosophila and Bombyx sex determination differ, similarities between the two species do exist in downstream genetic control of sex determination.

SAGE analysis of early oogenesis in the silkworm, Bombyx mori

To identify genes involved in the differentiation of Bombyx cystoblast, we constructed two 3' long serial analysis of gene expression (Long SAGE) libraries from stage 1-3 or stage 2-3 egg chambers and compared their gene expression profiles. In both libraries, the most frequent tags were derived from the same novel transcript. The transcript does not have any open reading frame capable of encoding a protein with over 100 amino acids in length. RNA blot analysis revealed that this transcript is specifically and abundantly expressed in the Bombyx ovary, mainly the germ line cells in the ovarioles. These results suggest that Bombyx oogenesis may be regulated by a previously unidentified non-coding RNA. Comparison of the gene expression profiles between the stage 1-3 and stage 2-3 egg chamber libraries revealed that 272 tags were significantly more abundant in stage 1-3 egg chambers (p<0.05 and at least two-fold change) than in library 2. Among the differentially expressed transcripts were the sequences that correspond to ATP synthase subunit d (3.1-fold enriched) and ATP synthase coupling factor 6 (9.1-fold enriched), suggesting that they are involved in regulation of cell cycle of cystocytes.

The female-killing chromosome of the silkworm, Bombyx mori, was generated by translocation between the Z and W chromosomes

Bombyx mori is a female-heterogametic organism (female, ZW; male, ZZ) that appears to have a putative feminizing gene (Fem) on the W chromosome. The paternally transmitted mutant W chromosome, Df(p ( Sa ) + ( p )W + ( od ))Fem, derived from the translocation-carrying W chromosome (p ( Sa ) + ( p )W + ( od )), is inert as femaleness determinant. Moreover, this Df(p ( Sa ) + ( p )W + ( od ))Fem chromosome has been thought to have a female-killing factor because no female larvae having the Df(p ( Sa ) + ( p )W + ( od ))Fem chromosome are produced. Initially, to investigate whether the Df(p ( Sa ) + ( p )W + ( od ))Fem chromosome contains any region of the W chromosome or not, we analyzed the presence or absence of 12 W-specific RAPD markers. The Df(p ( Sa ) + ( p )W + ( od ))Fem chromosome contained 3 of 12 W-specific RAPD markers. These results strongly indicate that the Df(p ( Sa ) + ( p )W + ( od ))Fem chromosome contains the region of the W chromosome. Moreover, by using phenotypic and molecular markers, we confirmed that the Df(p ( Sa ) + ( p )W + ( od ))Fem chromosome is connected with a partially deleted Z chromosome and that this fused chromosome behaves as a Z chromosome during male meiosis. Furthermore, we demonstrated that the ZZW-type triploid female having the Df(p ( Sa ) + ( p )W + ( od ))Fem chromosome is viable. Therefore, we concluded that the Df(p ( Sa ) + ( p )W + ( od ))Fem chromosome does not have a female-killing factor but that partial deletion of the Z chromosome causes the death of the ZW-type diploid female having the Df(p ( Sa ) + ( p )W + ( od ))Fem chromosome. Additionally, our results of detailed genetic analyses strongly indicate that the female-killing chromosome composed of the Df(p ( Sa ) + ( p )W + ( od ))Fem chromosome and deleted Z chromosome was generated by translocation between the Z chromosome and the translocation-carrying W chromosome, p ( Sa ) + ( p )W + ( od ).

Isolation and characterization of sex chromosome rearrangements generating male muscle dystrophy and female abnormal oogenesis in the silkworm, Bombyx mori

In deletion-mapping of W-specific RAPD (W-RAPD) markers and putative female determinant gene (Fem), we used X-ray irradiation to break the translocation-carrying W chromosome (W( Ze )). We succeeded in obtaining a fragment of the W( Ze ) chromosome designated as Ze (W), having 3 of 12 W-RAPD markers (W-Bonsai, W-Yukemuri-S, W-Yukemuri-L). Inheritance of the Ze (W) fragment by males indicates that it does not include the Fem gene. On the basis of these results, we determined the relative positions of W-Yukemuri-S and W-Yukemuri-L, and we narrowed down the region where Fem gene is located. In addition to the Ze (W) fragment, the Z chromosome was also broken into a large fragment (Z(1)) having the +( sch ) (1-21.5) and a small fragment (Z(2)) having the +( od ) (1-49.6). Moreover, a new chromosomal fragment (Ze (W)Z(2)) was generated by a fusion event between the Ze (W) and the Z(2) fragments. We analyzed the genetic behavior of the Z(1) fragment and the Ze (W)Z(2) fragment during male (Z/Z(1) Ze (W)Z(2)) and female (Z(1) Ze (W)Z(2)/W) meiosis using phenotypic markers. It was observed that the Z(1) fragment and the Z or the W chromosomes separate without fail. On the other hand, non-disjunction between the Ze (W)Z(2) fragment and the Z chromosome and also between the Ze (W)Z(2) fragment and the W chromosome occurred. Furthermore, the females (2A: Z/Ze (W)Z(2)/W) and males (2A: Z/Z(1)) resulting from non-disjunction between the Ze (W)Z(2) fragment and the W chromosome had phenotypic defects: namely, females exhibited abnormal oogenesis and males were flapless due to abnormal indirect flight muscle structure. These results suggest that Z(2) region of the Z chromosome contains dose-sensitive gene(s), which are involved in oogenesis and indirect flight muscle development.

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.

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.

Genetic mapping of Z chromosome and identification of W chromosome-specific markers in the silkworm, Bombyx mori

In the silkworm, Bombyx mori, the female is the heterogametic (ZW) sex and the male is homogametic (ZZ). The female heterogamety is a typical situation in the insect order Lepidoptera. Although the W chromosome in silkworm is strongly female determining, no W-linked gene for a morphological character has been found on it. The Z chromosome carries important traits of economic value as well as genes for various phenotypic traits, but only 2% of molecular information based on its relative size is known. Studies conducted so far indicate that the Z-linked genes are not dosage compensated. In the present study, we constructed a genetic map of randomly amplified polymorphic DNA fragments (RAPD), simple sequence repeats (SSR), and fluorescent intersimple sequence repeat PCR (FISSR) markers for the Z chromosome using a backcross mapping population. A total of 16 Z-linked markers were identified, characterized, and mapped using od, a recessive trait for translucent skin as an anchor marker yielding a total recombination map of 334.5 cM. The linkage distances obtained suggested that the markers were distributed throughout the Z chromosome. Four RAPD and four SSR markers that were linked to W chromosome were also identified. The proposed mapping approach should be useful to identify and map sex-linked traits in the silkworm. The economic and evolutionary significance of Z- and W-linked genes in silkworm, in particular, and lepidopterans, in general, is discussed.

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.

Invertebrate Muscles: Muscle Specific Genes and Proteins

This is the first of a projected series of canonic reviews covering all invertebrate muscle literature prior to 2005 and covers muscle genes and proteins except those involved in excitation-contraction coupling (e.g., the ryanodine receptor) and those forming ligand- and voltage-dependent channels. Two themes are of primary importance. The first is the evolutionary antiquity of muscle proteins. Actin, myosin, and tropomyosin (at least, the presence of other muscle proteins in these organisms has not been examined) exist in muscle-like cells in Radiata, and almost all muscle proteins are present across Bilateria, implying that the first Bilaterian had a complete, or near-complete, complement of present-day muscle proteins. The second is the extraordinary diversity of protein isoforms and genetic mechanisms for producing them. This rich diversity suggests that studying invertebrate muscle proteins and genes can be usefully applied to resolve phylogenetic relationships and to understand protein assembly coevolution. Fully achieving these goals, however, will require examination of a much broader range of species than has been heretofore performed.