Yuji Yasukochi

A Second-Generation Integrated Map of the Silkworm Reveals Synteny and Conserved Gene Order Between Lepidopteran Insects

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

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.

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.

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

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

Extensive Conserved Synteny of Genes Between the Karyotypes of Manduca sexta and Bombyx mori Revealed by BAC-FISH Mapping

Background Genome sequencing projects have been completed for several species representing four highly diverged holometabolous insect orders, Diptera, Hymenoptera, Coleoptera, and Lepidoptera. The striking evolutionary diversity of insects argues a need for efficient methods to apply genome information from such models to genetically uncharacterized species. Constructing conserved synteny maps plays a crucial role in this task. Here, we demonstrate the use of fluorescence in situ hybridization with bacterial artificial chromosome probes as a powerful tool for physical mapping of genes and comparative genome analysis in Lepidoptera, which have numerous and morphologically uniform holokinetic chromosomes. Methodology/Principal Findings We isolated 214 clones containing 159 orthologs of well conserved single-copy genes of a sequenced lepidopteran model, the silkworm, Bombyx mori, from a BAC library of a sphingid with an unexplored genome, the tobacco hornworm, Manduca sexta. We then constructed a BAC-FISH karyotype identifying all 28 chromosomes of M. sexta by mapping 124 loci using the corresponding BAC clones. BAC probes from three M. sexta chromosomes also generated clear signals on the corresponding chromosomes of the convolvulus hawk moth, Agrius convolvuli, which belongs to the same subfamily, Sphinginae, as M. sexta. Conclusions/Significance Comparison of the M. sexta BAC physical map with the linkage map and genome sequence of B. mori pointed to extensive conserved synteny including conserved gene order in most chromosomes. Only a few rearrangements, including three inversions, three translocations, and two fission/fusion events were estimated to have occurred after the divergence of Bombycidae and Sphingidae. These results add to accumulating evidence for the stability of lepidopteran genomes. Generating signals on A. convolvuli chromosomes using heterologous M. sexta probes demonstrated that BAC-FISH with orthologous sequences can be used for karyotyping a wide range of related and genetically uncharacterized species, significantly extending the ability to develop synteny maps for comparative and functional genomics.

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

The W chromosome of the silkworm Bombyxmori 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.

The silkworm W chromosome is a source of female-enriched piRNAs

In the silkworm, Bombyx mori, the W chromosome plays a dominant role in female determination. However, neither protein-coding genes nor transcripts have so far been isolated from the W chromosome. Instead, a large amount of functional transposable elements and their remnants are accumulated on the W chromosome. PIWI-interacting RNAs (piRNAs) are 23-30-nt-long small RNAs that potentially act as sequence-specific guides for PIWI proteins to silence transposon activity in animal gonads. In this study, by comparing ovary- and testis-derived piRNAs, we identified numerous female-enriched piRNAs. Our data indicated that female-enriched piRNAs are derived from the W chromosome. Moreover, comparative analyses on piRNA profiles from a series of W chromosome mutant strains revealed a striking enrichment of a specific set of transposon-derived piRNAs in the putative sex-determining region. Collectively, we revealed the nature of the silkworm W chromosome as a source of piRNAs.

Conserved synteny of genes between chromosome 15 of Bombyx mori and a chromosome of Manduca sexta shown by five-color BAC-FISH

The successful assignment of the existing genetic linkage groups (LGs) to individual chromosomes and the second-generation linkage map obtained by mapping a large number of bacterial artificial chromosome (BAC) contigs in the silkworm, Bombyx mori, together with public nucleotide sequence databases, offer a powerful tool for the study of synteny between karyotypes of B. mori and other lepidopteran species. Conserved synteny of genes between particular chromosomes can be identified by comparatively mapping orthologous genes of the corresponding linkage groups with the help of BAC-FISH (fluorescent in situ hybridization). This technique was established in B. mori for 2 differently labeled BAC probes simultaneously hybridized to pachytene bivalents. To achieve higher-throughput comparative mapping using BAC-FISH in Lepidoptera, we developed a protocol for five-color BAC-FISH, which allowed us to map simultaneously 6 different BAC probes to chromosome 15 in B. mori. We identified orthologs of 6 B. mori LG15 genes (RpP0, RpS8, eIF3, RpL7A, RpS23, and Hsc70) for the tobacco hornworm, Manduca sexta, and selected the ortholog-containing BAC clones from an M. sexta BAC library. All 6 M. sexta BAC clones hybridized to a single M. sexta bivalent in pachytene spermatocytes. Thus, we have confirmed the conserved synteny between the B. mori chromosome 15 and the corresponding M. sexta chromosome (hence provisionally termed chromosome 15).

Sex-linked pheromone receptor genes of the European corn borer, Ostrinia nubilalis, are in tandem arrays.

Background Tuning of the olfactory system of male moths to conspecific female sex pheromones is crucial for correct species recognition; however, little is known about the genetic changes that drive speciation in this system. Moths of the genus Ostrinia are good models to elucidate this question, since significant differences in pheromone blends are observed within and among species. Odorant receptors (ORs) play a critical role in recognition of female sex pheromones; eight types of OR genes expressed in male antennae were previously reported in Ostrinia moths. Methodology/Principal Findings We screened an O. nubilalis bacterial artificial chromosome (BAC) library by PCR, and constructed three contigs from isolated clones containing the reported OR genes. Fluorescence in situ hybridization (FISH) analysis using these clones as probes demonstrated that the largest contig, which contained eight OR genes, was located on the Z chromosome; two others harboring two and one OR genes were found on two autosomes. Sequence determination of BAC clones revealed the Z-linked OR genes were closely related and tandemly arrayed; moreover, four of them shared 181-bp direct repeats spanning exon 7 and intron 7. Conclusions/Significance This is the first report of tandemly arrayed sex pheromone receptor genes in Lepidoptera. The localization of an OR gene cluster on the Z chromosome agrees with previous findings for a Z-linked locus responsible for O. nubilalis male behavioral response to sex pheromone. The 181-bp direct repeats might enhance gene duplications by unequal crossovers. An autosomal locus responsible for male response to sex pheromone in Heliothis virescens and H. subflexa was recently reported to contain at least four OR genes. Taken together, these findings support the hypothesis that generation of additional copies of OR genes can increase the potential for male moths to acquire altered specificity for pheromone components, and accordingly, facilitate differentiation of sex pheromones.

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.