A complete full-length non-LTR retrotransposon, BMC1, on the W chromosome of the silkworm, Bombyx mori

The silkworm (Bombyx mori) neuropeptide orcokinin is involved in the regulation of pigmentation

The natural colorful cuticles of insects play important roles in many physiological processes. Pigmentation is a physiological process with a complex regulatory network whose regulatory mechanism remains unclear. Bombyx mori pigmentation mutants are ideal materials for research on pigmentation mechanisms. The purple quail-like (q-l^p) and brown quail-like (q-l^b) mutants originated from plain silkworm breeds 932VR and 0223JH respectively exhibit similar cuticle pigmentation to that of the quail mutant. The q-l^p mutant also presents a developmental abnormality. In this study, genes controlling q-l^p and q-l^b mutants were located on chromosome 8 by positional cloning. Then the neuropeptide gene orcokinin (OK) was identified to be the major gene responsible for two quail-like mutants. The B. mori orcokinin gene (BommoOK) produces two transcripts, BommoOKA and BommoOKB, by alternative splicing. The CRISPR/Cas9 system and orcokinin peptides injection were used for further functional verification. We show a novel function of BommoOKA in inhibiting pigmentation, and one mature peptide of orcokinin A, OKA_type2, is the key factor in pigmentation inhibition. These results provide a reference for studying the function of orcokinin and are of theoretical importance for studying the regulatory mechanism of pigmentation.

Positional cloning of the sex-linked giant egg (Ge) locus in the silkworm, Bombyx mori

The giant egg (Ge) locus is a Z-linked mutation that leads to the production of large eggs. Cytological observations suggest that an unusual translocation of a large fragment of the W chromosome bearing a putative egg size-determining gene, Esd, gave rise to giant egg mutants. However, there is currently no molecular evidence confirming either a W-Z translocation or the presence of Esd on the W chromosome. To elucidate the origin of giant egg mutants, we performed positional cloning. We observed that the Bombyx mori. orthologue of the human Phytanoyl-CoA dioxygenase domain containing 1 gene (PHYHD1) is disrupted in giant egg mutants. PHYHD1 is highly conserved in eukaryotes and is predicted to be a Fe(II) and 2-oxoglutarate-dependent oxygenase. Exon skipping in one of the two available Ge mutants is probably caused by the insertion of a non-long terminal repeat transposon into intron 4 in the vicinity of the 5' splice site. Segmental duplication in Ge(2) , an independent allele, was caused by unequal recombination between short interspersed elements inserted into introns 3 and 5. Our results indicate that (1) Bombyx PHYHD1 is responsible for the Ge mutants and that (2) the Ge locus is unrelated to the W-linked putative Esd. To our knowledge, this is the first report describing the phenotypic defects caused by mutations in PHYHD1 orthologues.

Antennapedia is involved in the development of thoracic legs and segmentation in the silkworm, Bombyx mori

Homeotic genes, which are associated closely with body patterning of various species, specify segment identity. The Wedge eye-spot (Wes) is a new homeotic mutant located on the sixth linkage group. Homozygous Wes/Wes embryos are lethal and display a pair of antenna-like appendages under the mouthparts as well as fused thoracic segments. These mutants also exhibit a narrower eye-spot at the larval stage compared with the wild type. By positional cloning, we identified the candidate gene of the Wes locus, Bombyx mori Antennapedia (BmAntp). Two BmAntp transcripts were identified in the homozygote of the Wes mutant, including a normal form and an abnormal form with a 1570-bp insertion. Our data showed that the insertion element was a long interspersed nuclear element (LINE)-like transposon that destroyed the original open reading frame of BmAntp. Quantitative RT-PCR analysis showed that the expression levels of normal BmAntp transcripts were increased markedly in the Wes heterozygous larvae compared with the wild type. Furthermore, we performed RNAi of BmAntp and observed fused thoracic segments and defective thoracic legs in the developing embryos. Our results indicated that BmAntp is responsible for the Wes mutant and has an important role in determining the proper development of the thoracic segments. Our identification of a homeotic mutation in the silkworm is an important contribution to our understanding of the regulation of Hox genes at different levels of expression.

Genome-wide screening and characterization of transposable elements and their distribution analysis in the silkworm, Bombyx mori

To elucidate the contribution of transposable elements (TEs) to the silkworm genome structure and evolution, we have conducted genome-wide analysis of TEs using the newly released genome assembly. The TEs made up 35% of the genome and contributed greatly to the genome size. Non-long terminal repeat retrotransposons (non-LTRs) and short interspersed nuclear elements (SINEs) were the predominant TE classes. From characterization of the TE distribution in the genome, it was revealed that non-LTRs, especially R1 clade elements, are frequently inserted into GC-rich regions. The GC content of non-LTRs themselves was over 40%, which indicate their contribution to the GC content of the insertion region. TEs accumulated in regions with low gene density, and there were relatively strong positive correlations between TE density and chromosomal recombination rate. We also characterized the clade distribution of the non-LTRs. The silkworm non-LTRs represented 10 of the 16 previously defined clades, which had the most variety than that reported for other genomes. Two partial CRE clade elements were found, which is one of the most ancient lineages of non-LTRs, and have been only found in Trypanosoma and fungi before. This analysis suggests that Bombyx genome is influenced by numerous amounts and variety of TEs.

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.

Probing the W chromosome of the codling moth, Cydia pomonella, with sequences from microdissected sex chromatin

The W chromosome of the codling moth, Cydia pomonella, like that of most Lepidoptera species, is heterochromatic and forms a female-specific sex chromatin body in somatic cells. We collected chromatin samples by laser microdissection from euchromatin and W-chromatin bodies. DNA from the samples was amplified by degenerate oligonucleotide-primed polymerase chain reaction (DOP-PCR) and used to prepare painting probes and start an analysis of the W-chromosome sequence composition. With fluorescence in situ hybridization (FISH), the euchromatin probe labelled all chromosomes, whereas the W-chromatin DNA proved to be a highly specific W-chromosome painting probe. For sequence analysis, DOP-PCR-generated DNA fragments were cloned, sequenced, and tested by Southern hybridization. We recovered single-copy and low-copy W-specific sequences, a sequence that was located only in the W and the Z chromosome, multi-copy sequences that were enriched in the W chromosome but occurred also elsewhere, and ubiquitous multi-copy sequences. Three of the multi-copy sequences were recognized as derived from hitherto unknown retrotransposons. The results show that our approach is feasible and that the W-chromosome composition of C. pomonella is not principally different from that of Bombyx mori or from that of Y chromosomes of several species with an XY sex-determining mechanism. The W chromosome has attracted repetitive sequences during evolution but also contains unique sequences.

APE-type non-LTR retrotransposons: determinants involved in target site recognition

Non-long terminal repeat (Non-LTR) retrotransposons represent a diverse and widely distributed group of transposable elements and an almost ubiquitous component of eukaryotic genomes that has a major impact on evolution. Their copy number can range from a few to several million and they often make up a significant fraction of the genomes. The members of the dominating subtype of non-LTR retrotransposons code for an endonuclease with homology to apurinic/apyrimidinic endonucleases (APE), and are thus termed APE-type non-LTR retrotransposons. In the last decade both the number of identified non-LTR retrotransposons and our knowledge of biology and evolution of APE-type non-LTR retrotransposons has increased tremendously.

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.

Molecular Analysis of Divergence in Tachinid Uzi (Exorista Sorbillans) Populations in India

Exorista sorbillans is a tachinid endoparasitoid of silkworm, Bombyx mori, and is globally known as uzi. It causes economic injury to the cocoon crop in silkworm cultivating areas of India, except those above 400 m above mean sea level (AMSL) in the foothills of the Himalayas (Darjeeling). It is reported that the sericulture tract of south India became infected with this pest only since 1980 through an accidental transportation of cocoons from West Bengal. To ascertain whether the genome of this parasitoid is differentiating into discrete gene pools in contrasting geo-climatic conditions, molecular profiling of four populations (Es (Annatapur), Es(Ramanagaram), Es (Channapatna) and Es(Kodathi) from south India and Es(Murshidabad) from Murshidabad, West Bengal was undertaken with 13 ISSR, 3 RAPD and six non-random primers designed from various repeat sequences of B. mori . MANOVA indicated significance for the Roy's largest root estimate (55.4; F =18.47; p = 0.002) for the variability contributed by the replication. Further, hierarchical clustering done on the basis of Euclidean distance matrix and Nei's unbiased Phylip clustering put Es(Murshidabad) at the maximum distance from those of south India and 29 markers could also be identified which significantly differentiateEs(Murshidabad) from others. However, Nei's statistics for gene diversity in sub-populations reveal considerably high gene-flow (3.44 and 2.51) among the populations around Bangalore. The gene-flow between Es(Murshidabad) and other population is lowest but cannot be ignored. The comparison of endosymbiont specific 16SrRNA and fts Z gene (partial) sequences through clustalW (gcgMSF) revealed a closer relationship of Es(Murshidabad) with Es(Annatapur) and Es (Ramanagaram) and is not congruent with the relationships discussed above. The significance of this maiden study with a tachinid fly-pest is discussed in the context of understanding the diversification of Uzi fly-pest and also establishing this pest as a relevant biological material for studying microevolution in future.

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.

Cytogenetic and molecular characterization of the MBSAT1 satellite DNA in holokinetic chromosomes of the cabbage moth, Mamestra brassicae (Lepidoptera)

Digestion of Mamestra brassicae DNA with DraI produced a prominent fragment of approximately 200 bp and a ladder of electrophoretic bands with molecular weights which are a multiple of 200 bp. Southern blotting revealed that this ladder is composed of DNA fragments that are multimers of the 200-bp DraI band suggesting that DraI isolated a satellite that has been called Mamestra brassicae satellite DNA1 (MBSAT1). MBSAT1 is the first satellite DNA isolated in Lepidoptera. In-situ DraI digestion of chromosome spreads, together with fluorescent in-situ hybridization, showed that MBSAT1 sequences are clustered in heterochromatin of the sex chromosomes, Z and W. MBSAT1 was 234 bp long with an AT content of 60.7%. The curvature-propensity plot suggested a curvature in the MBSAT1 structure.

W-derived BAC probes as a new tool for identification of the W chromosome and its aberrations in Bombyx mori

We isolated four W chromosome-derived bacterial artificial chromosome (W-BAC) clones from Bombyx mori BAC libraries by the polymerase chain reaction and used them as probes for fluorescence in situ hybridization (FISH) on chromosome preparations from B. mori females. All four W-BAC probes surprisingly highlighted the whole wild-type W sex chromosome and also identified the entire original W-chromosomal region in W chromosome-autosome translocation mutants. This is the first successful identification of a single chromosome by means of BAC-FISH in species with holokinetic chromosomes. Genomic in situ hybridization (GISH) by using female-derived genomic probes highlighted the W chromosome in a similar chromosome-painting manner. Besides the W, hybridization signals of W-BAC probes also occurred in telomeric and/or subtelomeric regions of the autosomes. These signals coincided well with those of female genomic probes except one additional GISH signal that was observed in a large heterochromatin block of one autosome pair. Our results support the opinion that the B. mori W chromosome accumulated transposable elements and other repetitive sequences that also occur, but scattered, elsewhere in the respective genome.

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

The sex chromosome constitution of the silkworm, Bombyx mori, is ZW in the female and ZZ in the male. Very little molecular information is available about the Z chromosome in Lepidoptera, although the topic is interesting because of the absence of gene dosage compensation in this chromosome. We constructed a 320-kb BAC contig around the Bmkettin gene on the Z chromosome in Bombyx and determined its nucleotide sequence by the shotgun method. We found 13 novel protein-coding sequences in addition to Bmkettin. All the transposable elements detected in the region were truncated, and no LTR retrotransposons were found, in stark contrast to the situation on the W chromosome. In this 320-kb region, four genes for muscle proteins (Bmkettin, Bmtitin1, Bmtitin2, and Bmprojectin) are clustered, together with another gene (Bmmiple) on the Z chromosome in B. mori; their orthologs are also closely linked on chromosome 3 in Drosophila, suggesting a partial synteny. Real-time RT-PCR experiments demonstrated that transcripts of 13 genes of the 14 Z-linked genes found accumulated in larger amounts in males than in female moths, indicating the absence of gene dosage compensation. The implications of these findings for the evolution and function of the Z chromosome in Lepidoptera are discussed.

Nested retrotransposons on the W chromosome of the wild silkworm Bombyx mandarina

The W chromosome of the silkworms Bombyx mori or B. mandarina is recombinationally isolated from the Z chromosome and the autosomes. We previously characterized a female-specific randomly amplified polymorphic DNA (RAPD), designated W-Yamato, derived from the W chromosome of the wild silkworm Bombyx mandarina. To further analyse the W chromosome of B. mandarina, we obtained a lambda phage clone that contains the W-Yamato RAPD sequence and sequenced the 16.7 kb DNA insert. We found that this DNA comprises a nested structure of at least seven elements: six retrotransposons and one transposable element-like sequence. The transposable element-like sequence is inserted into a micropia-like retrotransposon (Karate). The Karate and the non-long terminal repeat (non-LTR) retrotransposon BMC1 are inserted into a 412-like retrotransposon (Judo). Furthermore, this Judo, and two non-LTR retrotransposons (Kurosawa and Kendo) are inserted into a Pao-like retrotransposon (Yamato). These results indicate that the retrotransposons inserted into the W chromosome are not efficiently removed but accumulate gradually as strata without recombination.

Meiotic pairing of sex chromosome fragments and its relation to atypical transmission of a sex-linked marker in Ephestia kuehniella (Insecta: Lepidoptera)

The physical basis of non-Mendelian segregation of a sex-linked marker was studied in sex- chromosome mutant females of eight ASF ('abnormal segregating females') lines in the flour moth, Ephestia kuehniella. Electron microscopical analysis of microspread synaptonemal complexes revealed that in one line, the Z chromosome segment that contained the dz+ allele was translocated onto an autosome. The resulting quadrivalent visible in early female meiosis was 'corrected' into two bivalents in later stages. This explains autosomal inheritance of the sex chromosome marker in this strain. In the other seven ASF lines, the type of meiotic pairing of an additional fragment (Zdz+) of the Z chromosome was responsible for abnormal segregation of the marker gene. In several of these lines, Zdz+ contained a piece of the W chromosome in addition to the Z segment, as was confirmed by comparative genomic hybridization (CGH). Zdz+ formed three alternative pairing configurations with the original sex chromosomes: (i) a WZZdz+ trivalent, (ii) a WZ bivalent and a Zdz+ univalent or (iii) a ZZdz+ bivalent and a W univalent. In the most frequent WZZdz+ configuration, Zdz+ synapsed with Z and, consequently, segregated with W, simulating W linkage. This explains the predominant occurrence of the parental phenotypes in the progeny. Zdz+ univalents or W univalents, on the other hand, segregated randomly, resulting in both parental and nonparental phenotypes. In two of these lines, the Zdz+ was transmitted only to females. The results suggest that the W chromosome segment in Zdz+ of these lines contains a male-killing factor which makes it incompatible with male development. Our data provide direct evidence for the regular transmission of radiation-induced fragments from lepidopteran chromosomes through more than 50 generations. This is facilitated by the holokinetic nature of lepidopteran chromosomes. We conclude that Zdz+ fragments may persist as long as they possess active kinetochore elements.

A novel hybrid open reading frame formed by multiple cellular gene transductions by a plant long terminal repeat retroelement

The discovery that vertebrate retroviruses could transduce cellular sequences was central to cancer etiology and research. Although not well documented, transduction of cellular sequences by retroelements has been suggested to modify cellular functions. The maize Bs1 transposon was the first non-vertebrate retroelement reported to have transduced a portion of a cellular gene (c-pma). We show that Bs1 has, in addition, transduced portions of at least two more maize cellular genes, namely for 1,3-beta-glucanase (c-bg) and 1,4-beta-xylan endohydrolase (c-xe). We also show that Bs1 has maintained a truncated gag domain with similarity to the magellan gypsy-like long terminal repeat retrotransposon and a region that may correspond to an env-like domain. Our findings suggest that, like oncogenic retroviruses, the three transduced gene fragments and the Bs1 gag domain encode a fusion protein that has the potential to be expressed. We suggest that transduction by retroelements may facilitate the formation of novel hybrid genes in plants.

Molecular structure of the copia-like retrotransposable element Yokozuna on the W chromosome of the silkworm, Bombyx mori

We discovered a novel retrotransposable element, designated Yokozuna, on the W chromosome of Bombyx mori. The size of this element is 4738 bp, including a 208-bp long terminal repeat (LTR) on one side and a 183-bp LTR on the other. This retrotransposable element is flanked by a 5-bp target site duplication, TAATT. Yokozuna contains a single long open reading frame (ORF) and the whole deduced amino acid sequence of ORF reveals strong homology with copia of Drosophila. Moreover, an alignment analysis of the reverse transcriptase (RT) sequences suggested that the Yokozuna element is the first Bombyx retrotransposable element belonging to the Ty1-copia group. The number of Yokozuna per haploid genome is approximately four copies. Although Yokozuna was discovered on the W chromosome, it is not specific for the W chromosome.