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

The Draft Genome of Yellow Stem Borer, an Agriculturally Important Pest, Provides Molecular Insights into Its Biology, Development and Specificity Towards Rice for Infestation

Simple Summary

Yellow stem borer (YSB), is the most destructive and widely occurring pest that attacks rice throughout the growing season. Rice (Oryza sativa L.) is a major staple cereal worldwide, providing essential caloric requirements for more than half of the world’s population. Annual losses to rice borers are approximately 5–10%, but losses in individual fields may reach up to 50–60%. The use of traditional pest management strategies in controlling YSB is somewhat challenging due to its unique internal feeding habit. Genome sequence information of economically important crop pests is important for designing or developing pest-resistant rice varieties. In an approach to achieve this, we present our first-ever study on the draft genome sequence of YSB. The information provided from our current study might be useful in developing genome-based approaches for the management of pest species.

Abstract

Yellow stem borer (YSB), Scirpophaga incertulas (Walker) (Lepidoptera: Crambidae), a major monophagous insect pest of rice, causes significant yield losses. The rice–YSB interaction is very dynamic, making it difficult for management. The development of resistant lines has been unsuccessful as there are no effective resistant sources in the germplasm. Genome information is necessary for a better understanding of interaction with rice in terms of its recognition, response, and infestation mechanism. The draft genome of YSB is predicted to have 46,057 genes with an estimated size of 308 Mb, being correlated with the flow cytometry analysis. The existence of complex metabolic mechanisms and genes related to specific behavior was identified, being conditioned by a higher level of regulation. We deciphered the possible visual, olfactory, and gustatory mechanisms responsible for its evolution as a monophagous pest. Comparative genomic analysis revealed that YSB is unique in the way it has evolved. The obvious presence of high-immunity-related genes, well-developed RNAi machinery, and diverse effectors provides a means for developing genomic tools for its management. The identified 21,696 SSR markers can be utilized for diversity analysis of populations across the rice-growing regions. We present the first draft genome of YSB. The information emanated paves a way for biologists to design novel pest management strategies as well as for the industry to design new classes of safer and specific insecticide molecules.

LTR-retrotransposons in R. exoculata and other crustaceans: the outstanding success of GalEa-like copia elements

Transposable elements are major constituents of eukaryote genomes and have a great impact on genome structure and stability. They can contribute to the genetic diversity and evolution of organisms. Knowledge of their distribution among several genomes is an essential condition to study their dynamics and to better understand their role in species evolution. LTR-retrotransposons have been reported in many diverse eukaryote species, describing a ubiquitous distribution. Given their abundance, diversity and their extended ranges in C-values, environment and life styles, crustaceans are a great taxon to investigate the genomic component of adaptation and its possible relationships with TEs. However, crustaceans have been greatly underrepresented in transposable element studies. Using both degenerate PCR and in silico approaches, we have identified 35 Copia and 46 Gypsy families in 15 and 18 crustacean species, respectively. In particular, we characterized several full-length elements from the shrimp Rimicaris exoculata that is listed as a model organism from hydrothermal vents. Phylogenic analyses show that Copia and Gypsy retrotransposons likely present two opposite dynamics within crustaceans. The Gypsy elements appear relatively frequent and diverse whereas Copia are much more homogeneous, as 29 of them belong to the single GalEa clade, and species- or lineage-dependent. Our results also support the hypothesis of the Copia retrotransposon scarcity in metazoans compared to Gypsy elements. In such a context, the GalEa-like elements present an outstanding wide distribution among eukaryotes, from fishes to red algae, and can be even highly predominant within a large taxon, such as Malacostraca. Their distribution among crustaceans suggests a dynamics that follows a "domino days spreading" branching process in which successive amplifications may interact positively.

Establishment of a novel in vivo sex-specific splicing assay system to identify a trans-acting factor that negatively regulates splicing of Bombyx mori dsx female exons

The Bombyx mori homolog of doublesex, Bmdsx, plays an essential role in silkworm sexual development. Exons 3 and 4 of Bmdsx pre-mRNA are specifically excluded in males. To explore how this occurs, we developed a novel in vivo sex-specific splicing assay system using sexually differentiated cultured cells. A series of mutation analyses using a Bmdsx minigene with this in vivo splicing assay system identified three distinct sequences (CE1, CE2, and CE3) positioned in exon 4 as exonic splicing silencers responsible for male-specific splicing. Gel shift analysis showed that CE1 binds to a nuclear protein from male cells but not that from female cells. Mutation of UAA repeats within CE1 inhibited the binding of the nuclear protein to the RNA and caused female-specific splicing in male cells. We have identified BmPSI, a Bombyx homolog of P-element somatic inhibitor (PSI), as the nuclear factor that specifically binds CE1. Down-regulation of endogenous BmPSI by RNA interference significantly increased female-specific splicing in male cells. This is the first report of a PSI homolog implicated in the regulated sex-specific splicing of dsx pre-mRNA.

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.

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

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

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

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

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.

Role of the male BmDSX protein in the sexual differentiation of Bombyx mori

The sex determination pathway is different between Drosophila melanogaster and Bombyx mori in the initial signal. Here we show evidence that the sex determination pathway in B. mori is similar to that of D. melanogaster at the level of the terminal regulator, doublesex (dsx), which is essential for the proper differentiation of the sexually dimorphic somatic features of D. melanogaster. In B. mori, a homolog of dsx (Bmdsx) is expressed in various tissues, and its primary transcript is alternatively spliced in males and females to yield sex-specific mRNAs that encode male-specific (BmDSXM) and female-specific (BmDSXF) polypeptides. In the studies reported here, transgenic silkworms carrying a construct with a Bmdsx male cDNA placed under the control of either an hsp70 promoter or a Bombyx actin3 promoter were generated by piggyBac-mediated germline transformation. Ectopic expression of the male cDNA in females resulted in abnormal differentiation of certain female-specific genital organs and caused partial male differentiation in female genitalia. Transgenic analysis also revealed that the expression of BmDSXM in females caused repression of the female-specifically expressed gene, the vitellogenin gene, and also resulted in activation of the pheromone-binding protein gene that is dominantly expressed in males. These results provide evidence that the role of BmDSXM includes the activation of some aspects of male differentiation as well as the repression of female differentiation. Taken together with our previous data on the function of BmDSXF, we can conclude that Bmdsx is a double-switch gene at the final step in the sex-determination cascade of B. mori.

The genetics and genomics of the silkworm, Bombyx mori

We review progress in applying molecular genetic and genomic technologies to studies in the domesticated silkworm, Bombyx mori, highlighting its use as a model for Lepidoptera, and in sericulture and biotechnology. Dense molecular linkage maps are being integrated with classical linkage maps for positional cloning and marker-assisted selection. Classical mutations have been identified by a candidate gene approach. Cytogenetic and sequence analyses show that the W chromosome is composed largely of nested full-length long terminal repeat retrotransposons. Z-chromosome-linked sequences show a lack of dosage compensation. The downstream sex differentiation mechanism has been studied via the silkworm homolog of doublesex. Expressed sequence tagged databases have been used to discover Lepidoptera-specific genes, provide evidence for horizontal gene transfer, and construct microarrays. Physical maps using large-fragment bacterial artificial chromosome libraries have been constructed, and whole-genome shotgun sequencing is underway. Germline transformation and transient expression systems are well established and available for functional studies, high-level protein expression, and gene silencing via RNA interference.

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.

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.

Molecular cloning of two prophenoloxidase genes from the mosquito Aedes aegypti

The biosynthesis of melanotic materials is an important process in the life of a mosquito. Melanin production is critical for many diverse processes such as egg chorion tanning, cuticular sclerotization, and melanotic encapsulation of metazoan parasites. Prophenoloxidase plays a critical role in this biochemical cascade. Two cDNAs, one full length and one partial clone, and two genomic clones encoding prophenoloxidase (pro-PO) were isolated from the yellow fever mosquito, Aedes aegypti. The full-length cDNA, pAaProPO1, is 2286 bp long with a 2055 bp open reading frame encoding a 685 amino acid protein that shares 89% identity with Armigeres subalbatus pro-PO. It contains two putative copper binding domains (amino acids 197-243 and 346-423) that are homologous to other insect pro-POs. AaProPO1 messenger RNA (mRNA) was detected by reverse transcription polymerase chain reaction (RT-PCR) only from third-stage larvae and not in adult mosquitoes after blood feeding, during the melanotic encapsulation of Dirofilaria immitis microfilariae or following exposure to bacteria. A 750 bp fragment of the second cDNA (pAaProPO2) was cloned using RT-PCR from mRNA obtained from 14-day postovipostional eggs. AaProPO2 mRNA was not found in any other life stages, and may be in low abundance or transiently expressed. AaProPO2 and AaProPO1 each contain three introns that are 60, 68 and 58 bp and 61, 69 and 59 bp long, respectively, and the intron sequences of these two genes are not similar.

A homologue of the Drosophila doublesex gene is transcribed into sex-specific mRNA isoforms in the silkworm, Bombyx mori

The doublesex (dsx) gene is known as the final gene of the sex-determining cascade in Drosophila melanogaster. We have isolated a homologue of dsx in the silkworm, Bombyx mori, which has an epistatic feminizing gene located on the W chromosome. RT-PCR analysis indicated that B. mori dsx (Bmdsx) was transcribed in all the examined tissues, and the size of the amplified products was different between males and females. In Northern blot hybridization of poly(A)(+) RNA, the Bmdsx probe also detected a band with a sex-specific size difference. The male-specific cDNA lacked the sequence between 713 and 961nt of the female-specific cDNA. An RNase protection assay indicated that this sequence was male-specifically removed from the Bmdsx pre-mRNA. Southern blot analysis showed that Bmdsx is present at a single copy in the genome. These results suggested that the primary Bmdsx transcript is alternatively spliced to yield male- and female-specific mRNA isoforms. These sex-specific isoforms encode polypeptides with a common amino-terminal sequence but sex-specific carboxyl termini. DNA binding domain (DM domain) of BmDSX has 80% identity with D. melanogaster DSX proteins. These results suggest the Bmdsx would also regulate sexual differentiation, as does the Drosophila dsx gene.

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

In the silkworm, Bombyx mori, a non-long terminal repeat (non-LTR) retrotransposon, BMC1, is considered to be a LINE (long interspersed nuclear element)-like element. So far, a BMC1 containing two intact open reading frames (ORFs) has not been found. However, we discovered a complete full-length BMC1 on the W chromosome. This BMC1 is 5091 bp and contains a 5' untranslated region (5'-UTR), two intact ORFs, and 3'-UTR which terminates in a poly(A) tail. ORF1 encodes a putative nucleic acid-binding protein, while ORF2 encodes a protein containing an endonuclease domain and a reverse transcriptase domain.