Genome-wide screening and characterization of transposable elements and their distribution analysis in the silkworm, Bombyx mori
Introduction
In eukaryotes, mobile elements or transposable elements (TEs) constitute a large fraction of the genome (Kazazian, 2004). TEs have been divided into two classes based on their transposition intermediate (Capy et al., 1997). Class I elements use an RNA-mediated mode of transposition and encode a reverse transcriptase (RT). They are further divided into two subclasses, the elements that are characterized by direct, long terminal repeats (LTR retrotransposons), and the elements that lack terminal repeats (non-LTR retrotransposons). Autonomous non-LTR retrotransposons are also called LINEs (long interspersed nuclear elements) and are thought to be responsible for mobilization of non-autonomous short interspersed nuclear elements (SINEs). On the other hand, class II elements (transposons) use a DNA-mediated mode of ‘cut and paste’ transposition. In the silkworm (Bombyx mori), many TEs have been identified [e.g. Pao (Xiong et al., 1993), Kamikaze and Yamato (Abe et al., 2001), BMC1 (Ogura et al., 1994), L1Bm (Ichimura et al., 1997), R1Bm (Xiong and Eickbush, 1988), R2Bm (Burke et al., 1987), SART1 (Takahashi et al., 1997), TRAS1 (Okazaki et al., 1995), Bm1 (Adams et al., 1986) and mariner (Robertson and Asplund, 1996, Tomita et al., 1997)].
With a recent availability of genomic sequences in many species, it has been suggested that TEs have a major impact upon genome evolution in various scales e.g. genome size, genome rearrangement and deletion/insertion of a single nucleotide (Kazazian, 2004, Kidwell, 2002, Petrov, 2001). Although the elucidation of interaction of TEs with their host genome is undoubtedly an important work, it had been impracticable in the silkworm genome due to the lack of a sufficient genomic resource. However, the currently available silkworm genomic resource (ISGSC, in preparation), e.g. the new genome assembly by integrating two independent whole genome shotgun sequences data (Mita et al., 2004, Xia et al., 2004) and comprehensive catalogue of silkworm repetitive sequences, provides powerful new opportunities for rigorous analysis of TEs distribution in the genome. The preliminary analysis revealed that TEs make up about 35% of the silkworm genome, and that SINEs and non-LTR retrotransposons constitute about 75% of the TEs-derived sequences (ISGSC, in preparation). This tendency contrasts the genomic structure of B. mori with that of insects such as Drosophila melanogaster (Bergman et al., 2006), Anopheles gambiae (Holt et al., 2002) and Apis mellifera (HGSC, 2006), but similar to Aedes aegypti in which TEs share about 47% of the genome (Nene et al., 2007), and also to the human genome, in which L1 non-LTR retrotransposon and an SINE named Alu make up more than 27% of the genome (IHGSC, 2001).
In this paper, we aimed at detailed characterization of TEs, especially non-LTR retrotransposons, and their distribution in the silkworm genome in order to clarify their contribution to the genome structure and evolution.
Section snippets
Bombyx genomic sequence and transposable elements
We used silkworm scaffold build2 (ISGSC, in preparation) as genomic sequences data. The total length of all scaffolds with and without gaps amounted to 480.78 and 431.76 Mb, respectively. Total length of scaffolds mapped onto chromosomes with and without gaps (Ns) amounted to 419.98 and 377.18 Mb, respectively. The GC% of scaffolds in total was 37.7. The library of the silkworm TEs was same as used in the related paper (ISGSC, in preparation). To state briefly, to analyze the share of
Fraction and classes of TE-derived sequences in Bombyx genome
Table 1 summarizes the number of copies and fraction of genome for each TE classes. As mentioned above, the TEs constitute about 35% of the silkworm genome. Compared to other insects in which genomic sequences have been released, the silkworm genome was found to be second repetitive. The most TE-rich insect is A. aegypti, in which TEs share 47% of the genome (Nene et al., 2007). The TE content in D. melanogaster euchromatin is about 5.5% (Bergman et al., 2006), and about 77% including other
Conclusion
The genome-wide screening and characterization of the silkworm TEs revealed that a large amount of the silkworm genome is due to TEs, especially non-LTRs and SINEs. From the distribution analysis, it was shown that the silkworm TEs have positive correlations with regional GC content and recombination rate, indicating that TEs have greatly marked the genomic structure of the silkworm genome. These data are valuable in considering how TEs contribute to genome evolution.
Acknowledgements
This work was supported in part by grants in aid for scientific research from the Ministry of Education, Science, and Culture of Japan (to HF), the Insect Technology Program of MAFF (to KM), PROBRAIN, the Basic Research Program (to KM and HF), and Research Fellowship of Japan Society for the Promotion of Science for Young Scientists (to MO). We thank Kazutoshi Yoshitake for helping with the data analysis.
References (62)
- et al.
A highly reiterated family of transcribed oligo(A)-terminated, interspersed DNA elements in the genome of Bombyx mori
J. Mol. Biol.
(1986) - et al.
Tandem UAA repeats at the 3′-end of the transcript are essential for the precise initiation of reverse transcription of the I factor in Drosophila melanogaster
J. Biol. Chem.
(2002) Transposable elements: how non-LTR retrotransposons do it
Curr. Biol.
(1997)Evolution of genome size: new approaches to an old problem
Trends Genet.
(2001)- et al.
Bmmar1: a basal lineage of the mariner family of transposable elements in the silkworm moth, Bombyx mori
Insect Biochem. Mol. Biol.
(1996) - et al.
A complete full-length non-LTR retrotransposon, BMC1, on the W chromosome of the silkworm, Bombyx mori
Genes Genet. Syst.
(1998) - et al.
Two novel Pao-like retrotransposons (Kamikaze and Yamato) from the silkworm species Bombyx mori and B. mandarina: common structural features of Pao-like elements
Mol. Genet. Genomics
(2001) - et al.
Nested retrotransposons on the W chromosome of the wild silkworm Bombyx mandarina
Insect Mol. Biol.
(2002) - et al.
Functional roles of 3′-terminal structures of template RNA during in vivo retrotransposition of non-LTR retrotransposon, R1Bm
Nucleic Acids Res.
(2005) - et al.
A distal enhancer and an ultraconserved exon are derived from a novel retroposon
Nature
(2006)
Recurrent insertion and duplication generate networks of transposable element sequences in the Drosophila melanogaster genome
Genome Biol.
Non-LTR retrotransposons in the African malaria mosquito, Anopheles gambiae: unprecedented diversity and evidence of recent activity
Mol. Biol. Evol.
The site-specific ribosomal insertion element type II of Bombyx mori (R2Bm) contains the coding sequence for a reverse transcriptase-like enzyme
Mol. Cell. Biol.
The domain structure and retrotransposition mechanism of R2 elements are conserved throughout arthropods
Mol. Biol. Evol.
Dynamics and Evolution of Transposable Elements
Retrotransposition of the I factor, a non-long terminal repeat retrotransposon of Drosophila, generates tandem repeats at the 3′ end
Nucleic Acids Res.
Multiple sequence alignment with the Clustal series of programs
Nucleic Acids Res.
Transposons but not retrotransposons are located preferentially in regions of high recombination rate in Caenorhabditis elegans
Genetics
Recombination rates between adjacent genic and retrotransposon regions in maize vary by 2 orders of magnitude
Proc. Natl. Acad. Sci. U.S.A.
Telomere-specific non-LTR retrotransposons and telomere maintenance in the silkworm, Bombyx mori
Chromosome Res.
A rapidly rearranging retrotransposon within the miniexon gene locus of Crithidia fasciculata
Mol. Cell. Biol.
The genetics and genomics of the silkworm, Bombyx mori
Annu. Rev. Entomol.
L1-like non-LTR retrotransposons in the yeast Candida albicans
Curr. Genet.
The modulation of DNA content: proximate causes and ultimate consequences
Genome Res.
The genome sequence of the malaria mosquito Anopheles gambiae
Science
Insights into social insects from the genome of the honeybee Apis mellifera
Nature
A major non-LTR retrotransposon of Bombyx mori, L1Bm
J. Mol. Evol.
Initial sequencing and analysis of the human genome
Nature
Comparative recombination rates in the rat, mouse, and human genomes
Genome Res.
Cited by (92)
Masculinizer is not post-transcriptionally regulated by female-specific piRNAs during sex determination in the Asian corn borer, Ostrinia furnacalis
2023, Insect Biochemistry and Molecular BiologyCRISPR/Cas9-mediated mutagenesis of Ago2 and Siwi in silkworm cultured cells
2021, GeneCitation Excerpt :Reads shorter than 20 nt or longer than 32 nt were excluded, thereby obtaining reads of 20–32 nt. Small RNAs were mapped to the Bombyx genome (Kawamoto et al., 2019), transposons (Osanai-Futahashi et al., 2008), and BmLV genome (Iwanaga et al., 2012) using bowtie (Langmead et al., 2009). Reads that could be aligned to the genome with up to two mismatches were used to calculate the mapping rate of each library.
High-quality genome assembly of the silkworm, Bombyx mori
2019, Insect Biochemistry and Molecular BiologyCitation Excerpt :This suggests that the basic results are not entirely dependent on the source of TE libraries. The percentage of repetitive elements in the old genome assembly was previously estimated to be 43.6% (International Silkworm Genome Consortium, 2008; Osanai-Futahashi et al., 2008). The differences between the current and previous data may be due to the version of RepeatMasker and the search algorism used.
Evolutionary dynamics of genome size and transposable elements in crickets (Ensifera: Gryllidea)
2024, Systematic Entomology
- 1
These authors contributed equally to this work.