Skip to main content
  • Research article
  • Open access
  • Published:

Burst expansion, distribution and diversification of MITEs in the silkworm genome

Abstract

Background

Miniature inverted-repeat transposable elements (MITEs) are widespread in plants and animals. Although silkworm (Bombyx mori) has a large amount of and a variety of transposable elements, the genome-wide information of the silkworm MITEs is unknown.

Results

We used structure-based and homology approaches to search for MITEs in the silkworm genome. We identified 17 MITE families with a total of 5785 members, accounting for ~0.4% of the genome. 7 of 17 MITE families are completely novel based on the nucleotide composition of target site duplication (TSD) and/or terminal inverted repeats (TIR). Silkworm MITEs were widely and nonrandom distributed in the genome. One family named BmMITE-2 might experience a recent burst expansion. Network and diversity analyses for each family revealed different diversification patterns of the silkworm MITEs, reflecting the signatures of genome-shocks that silkworm experienced. Most silkworm MITEs preferentially inserted into or near genes and BmMITE-11 that encodes a germline-restricted small RNA might silence its the closest genes in silkworm ovary through a small RNA pathway.

Conclusions

Silkworm harbors 17 MITE families. The silkworm MITEs preferred to reside in or near genes and one MITE might be involved in gene silence. Our results emphasize the exceptional role of MITEs in transcriptional regulation of genes and have general implications to understand interaction between MITEs and their host genome.

Background

Since transposable elements (TEs) were first discovered in maize by McClintock [1], TEs have been a hotspot subject in genetics. As the advent of genomics era, it is now known that TEs constitute a significant component of eukaryotic genomes and that there are two major classes of TEs [2, 3]. Class I elements (retrotransposon) transpose by the mechanism of 'copy-and-paste' through a RNA intermediate, while class II elements (DNA transposons) by the mechanism of 'cut-and-paste' through a DNA intermediate.

In addition to transposition through a DNA intermediate rather than a RNA intermediate, DNA transposable element is distinguished by terminal inverted repeats (TIRs). With the action of transposase, DNA transposons usually excise from one site and reinsert elsewhere in the genome. Furthermore, DNA transposons can be also further classified into autonomous and non-autonomous. Of non-autonomous DNA transposons, a group of miniature inverted repeat transposable elements (MITEs) first discovered in maize were found to be widespread in various higher organisms [4–9]. MITEs have some common characteristics: short, terminal inverted repeats (TIR), target site duplication (TSD), high AT content, potential to form stable secondary structure, high number of copies in a genome [3, 9–11]. Importantly, recent genome-wide analyses revealed that MITEs insert preferentially into or near genes [5, 12, 13] and that several families of miRNAs in humans, Arabidopsis, rice and Solanaceae were derived from MITEs [14–16]. These results suggest that MITEs play important roles not only in genome evolution but also in transcriptional regulation of genes.

MITEs were originally discovered in plants and have two major superfamilies. They are Tourist-like and Stowaway-like on the basis of their similarity to two elements previously identified in maize and sorghum [3, 4, 10]. It is now clear that these two MITEs originated by internal deletion of the corresponding autonomous elements [4] and could transpose by other autonomous DNA elements [17, 18]. Virtually, all MITEs have been identified through computer-assisted database searches [11]. As more and more genome sequences become available, therefore, several programs have been developed for MITE identification [8, 19, 20]. Strikingly, a recently developed program, known as MUST, has been verified to be very powerful and efficient for identification of MITEs [20]. This facilitates discovery of more novel MITEs.

Silkworm, Bombyx mori, is a model insect for the order Lepidoptera and has importantly economic value for silk production and bioreactor. The draft genome sequences of silkworm were released by Mita et al. [21] and Xia et al. [22], respectively. Recently, a new assembly has been completed [23]. Analyses of the silkworm genome sequence suggested that ~40% of the genome is composed of the known TEs [24]. This number is only smaller than the TE proportion (47%) in Aedes aegypti genome that has the largest proportion of TEs in the insect genomes sequenced to date [25]. DNA transposons (Tc1-mariner, Helitron, Harbinger, hAT, P and Piggybac) are only ~3% of the genome while most of the silkworm TEs are retrotransposons [24]. Although two MITEs, Hoshidandy and Organdy, were previously discovered in silkworm [26] (also see GenBank accession no. AB455941), the genome-wide information about the silkworm MITEs remains unclear.

In this study, we scanned the new assembly of the silkworm genome sequence to identify MITEs by using a recently developed program, MUST and a strict filtering approach [20]. The results indicated that the silkworm genome harbors 17 MITE families. Estimates of insertion date and diversity for each MITE family showed that the silkworm MITE families might have experienced burst expansions at different time points of evolution and exhibited various patterns of diversification. That MITEs preferentially insert into or near genes has been also confirmed in silkworm. In addition, we found evidence that a sRNA derived from a MITE might silence its host and neighbor genes.

Results

Mining and characterization of MITEs

MUST [20], a program designed to detect MITE elements, was first used to search the silkworm genome sequence. The program identifies candidates based on common features of MITEs (short, TSDs and TIRs structure) and sequence alignment. With MUST, we mined 143333 MITE candidates in the silkworm genome, which were grouped into 1350 families. Then, we filtered out pseudo-MITEs from predicted ones by a strict approach: Those containing undetermined fragments (designated as Ns in scaffolds) as well as those solely composed of simple repeats and nested in repeats were considered as pseudo-MITEs (additional file 1). By this way, the number of the MITE families reduced to 17 (Table 1). These 17 families include 3337 intact MITEs and were designated as BmMITE-1 to BmMITE-17, respectively. TSD lengths of all MITEs range from 2 to 9 bp, TIR lengths from 8 to 59 bp, and full lengths of complete MITEs from 210 to 567 bp. Three families are flanked by TA, three families flanked by TWA, four families flanked by NNNNNNNN (the N represents A, T, C or G) and seven families flanked by WW, TDA, ATT, ATAT, ATATAT, TTCATTT, TTACTGTAT (the W represents A or T; D represents A, T or G), respectively. Based on the nucleotide composition of TSD, 3 MITE families (BmMITE-4,5,6) belong to Tourist-like family [4], 3 families (BmMITE-2,3,8) belong to Stowaway-like family [10]; based on the nucleotide composition of TSD and TIR, 4 families (BmMITE-13,14,15,16) belong to Pegasus-like family [27], and the remaining 7 families (BmMITE-1, 7, 9, 10, 11, 12, 17) were novel.

Table 1 Families and characteristics of silkworm MITEs

Three families flanked respectively by TA and by TWA have no any similarities in TIR and intra-sequence; they may have independent origins but may be transposed by the same transposase. Four families flanked by NNNNNNNN have no any similarities in intra-sequence but have high similarities in TIR regions; they may derive from a common ancestor and may be transposed by the same transposase [3]. Seven families with different TSD, TIR and intra-sequences may have different origins and may be transposed by different transposases.

Next, a homology search was used to estimate number of copies for each MITE family in the silkworm genome. That is, a BLASTN (E < e-5) search [10] was used against the silkworm genome sequence. With this approach, we identified 5785 MITEs in total, which constitute ~1.86 Mb (0.4%) of the silkworm genome. Then, each MITE was classified as intact (full length) or fragmentary: individual with both complete TIRs was regarded as intact and otherwise as fragmentary. As a result, 3337 MITEs are intact whereas 2448 MITEs are fragmentary. As mentioned above, 5785 MITEs belong to 17 families. However, family size varies greatly among MITEs, which ranges from 9 to 2790 (Table 1). The information about insertion sites of the MITEs into scaffold was shown in Additional file 2. The largest family is BmMITE-2 with 2790 copies. The copies of both BmMITE-2 and BmMITE-8 occupy approximately 74% of all MITEs. The ratio of full length to all copies of each family ranges from 10% to 100%. Four MITE families (BmMITE-1, BmMITE-4, BmMITE-6 and BmMITE-8) have the ratios less than 50% and all others more than 50%. Especially, BmMITT-2, BmMITE-11, BmMITE-12, and BmMITE-15 have 85%, 94%, 86% and 100% of high ratios, respectively, suggesting that they may be very young and active recently. It should be noted that BmMITE-2 has the largest number (2371) of full length MITEs whereas BmMITE-8 has the largest number (1364) of fragmentary MITEs among the families, suggesting that the former might experience a recent burst expansion whereas the latter might undergo an old burst expansion.

The average AT content for each MITE family ranges from 44% to 73%. Only BmMITE-4 has less than 50% AT content and all others more than 50%, suggesting that the silkworm MITEs are AT rich. The average AT content of the silkworm genome is approximately 62% [23]. Eight families BmMITE-1, BmMITE-2, BmMITE-8, BmMITE-9, BmMITE-10, BmMITE-11, BmMITE-12 and BmMITE-15 have AT content more than the average of the silkworm genome.

The secondary structure and negative ΔG value of Gibbs energy were predicted by UNAFOLD [28]. The results suggested that all the silkworm MITEs have the predicted secondary structures (Figure 1 and additional file 3). Furthermore, almost all MITE families have high negative ΔG values except for BmMITE-11 (8 kcal/mol) and BmMITE-12 (14 kcal/mol), indicating that most silkworm MITEs have potential to form stable secondary structures. These results imply that they may play important roles in transcriptional regulation of genes.

Figure 1
figure 1

The secondary structure of BmMITE-2 predicted by the program UNADOLD.

The identified silkworm MITE families were annotated based on BLAST search. Using a representative member of each MITE family as query, we searched ISfinder, RepBase and NCBI nr databases, respectively. These three databases almost include all known MITEs and transposable elements. We found that BmMITE-2 perfectly matched to one silkworm MITE named Hoshidandy (No. AB455941) in NCBI nr database and BmMITE-7 perfectly matched to Organdy identified previously in silkworm [26] (Table 1). Another 15 MITE families have no match to any known MITEs or transposable elements. This suggested that the program MUST used in this study is feasible and efficient for identification of MITEs because 2 known families of the silkworm MITEs were all recovered and 15 novel families were discovered by MUST [20].

Verification of the predicted MITEs

Four predicted insertion sites of BmMITE-2 were selected for verification by PCR in 14 silkworm accessions representing four main geographic strains (Chinese, Japanese, European, and Tropical). Primers were designed based on flanking regions of each insertion site (additional file 4). The PCR results for the four insertion sites were summarized in additional file 5. Furthermore, sequencing of the corresponding PCR products confirmed the presence or absence of BmMITE-2 (additional file 6). The polymorphic distribution of indel for all four tested insertion sites suggested that they are not fixed residents of the silkworm genome and validated the prediction by the program MUST [20].

Figure 2 shows the results of amplification for insertion site 4: The BmMITE-2 is present at the lines Wu-D, DongDe-201, Wu-E, Lu-10, DaXianTuZhong, FuRongHuiLuan, Ri-9, HeiZi, ChunSi, QiongShanHaiNan, Ri-110, ShangSanHuBan and absent at strains Wu-B, PeiXianZhong. The presence/absence polymorphisms at all four insertion sites of BmMITE-2 across the 14 tested silkworm lines suggested that BmMITE-2 may be recently active.

Figure 2
figure 2

PCR verification in 14 silkworm strains for BmMITE-2. Lanes from 1 to M are: Wu-B, Wu-D, DongDe-201, Wu-E, Lu-10, DaXianTuZhong, FuRongHuiLuan, PeiXianZhong, Ri-9, HeiZi, ChunSi, QiongShanHaiNan, Ri-110, ShangSanHuBan, DNA marker, respectively. The black arrow points to the BmMITE-2 occupied at this genomic location. The white arrow represents the BmMITE-2 lack at this genomic location.

Estimates of insertion date and diversity

We estimated age of each full length MITE by the method used for maize MITE families [13]. Briefly, we first determined divergence between each MITE and the family consensus sequence, and then estimated insertion time based on the divergence. Figure 3 showed that insertion date varies greatly among members of each family as well as among families, which ranges from 0 to 4 million years ago (mya). Strikingly, BmMITE-2 might be dramatically expanded during a period from 0 mya to 1 mya (Figure 3B) and accumulated up to 2173 copies during this short period while the rest 16 families might experience major expansion events within 2 mya (Figure 3).

Figure 3
figure 3

Insertion dates of each MITE family in the silkworm genome. The level of nucleotide substitutions (k) between each MITE element and the family consensus sequence was estimated using Kimura 2-parameter distance. Then, the insertion time was calculated by the method in [13].

To examine intrafamily diversity pattern for each family, we performed network analyses on the basis of the alignment of full length sequences. Topology of a network reflects corresponding demographic events that each MITE family experienced [13, 29]. For instance, a network topology showing numerous nodes distributed around its centre and separated by long branches implies that this family might experience an old population expansion. In contrast, a network topology characterized by a central node surrounded by many, almost identical, and short branches, indicates a very recent expansion from an ancestral element [13]. 8 (BmMITE-4, BmMITE-6, BmMITE-7, BmMITE-11, BmMITE-14, BmMITE-15, BmMITE-16 and BmMITE-17) of 17 MITE families presented the topologies that these MITE families might experience old population expansions. The topologies of the rest 9 families indicated that these families might undergo recent expansions (Figure 4). These observations are basically consistent with the insertion time estimates.

Figure 4
figure 4

Median-joining networks of 17 silkworm MITE families. 1 to 17 represent BmMITE-1 to BmMITE-17 family, respectively. The circle, circle area and branch length represent MITE sequence, proportional to the number of identical copies and proportional to the number of nucleotide changes, respectively.

Pattern of regional variation for each family provides details of divergence between copies. We performed diversity (pi) sliding window analysis for the alignment of full length sequence of each MITE family. Window size is 20 nucleotides and step size is 20 nucleotides [13]. The regions with diversity equal or higher than average sequence diversity (pi) +2SD was defined as variable whereas those with diversity less than average sequence diversity (pi)-2SD were considered as conserved regions (additional file 7). It was found that regional variation pattern varies greatly among these 17 silkworm MITE families. Almost all families have highly conserved TIRs, however, 6 families (BmMITE-4, BmMITE-5, BmMITE-9, BmMITE-12, BmMITE-14, and BmMITE-15) have only one conserved TIR. Because TIRs must be first recognized by corresponding transposase for transposition, it is not surprising to see more conserved in TIR than intra-sequence.

Distributions of MITEs on chromosomes and Estimation of MITE richness in genes

We recorded information about position of each MITE on a chromosome to look at distribution of MITEs in the genome. The results suggested that all MITE families are widely distributed on all 28 silkworm chromosomes. Then, we examined whether MITEs were randomly distributed among 28 silkworm chromosomes using χ2 test. The null hypothesis was rejected (P < 0.01), suggesting that the distribution of MITEs is nonrandom in the silkworm genome (additional file 8).

Next, we examined whether insertion site of each MITE preferentially is in or close to genes. If a MITE inserts into within the 5 kb flanking regions of a genes, this MITE is regarded as close to a gene [13, 30–33]. The results indicated that a larger number of MITEs inserted into gene regions (exon, intron) and flanking regions of genes. 3794 (66%) of the 5785 predicted MITEs inserted into gene regions (Table 2). Of 3794 MITEs inserted into gene regions, 962 (25%) were located in 5` flanking regions of the closest genes, 60 (2%) in exons, 1427 (38%) in introns and 1343 (35%) in 3`flanking regions of the closest genes, respectively. It appears that the silkworm MITEs preferentially inserted into introns and 3'-flanking regions rather than 5'-flanking regions and exons. To determine whether insertions of MITEs into gene regions are due to chance, we performed a computer simulation as a negative control (see Methods for details). It was found that the silkworm MITEs have significantly higher insertion frequencies into gene regions than those in control (66% vs 39%: χ2 test, P = 0.0, Table 2), implying that the silkworm MITEs preferentially do insert into gene regions. When a MITE insertion into within respective 500, 500-3000, and 3000-5000 bp flanking regions of a gene was assumed to be close to a gene, the patterns of MITE insertion in the genome are also similar to that a MITE within 5 kb flanking regions of a genes (additional file 9).

Table 2 Characteristics of insertion sites of the silkworm MITEs.

To investigate whether silkworm MITEs preferentially insert in a special gene family, we annotated all the closest genes using WEGO, an online gene ontology website, and compared the closest genes with all silkworm genes (additional file 10). The results indicated that functionalities of the closest genes were randomly selected.

A MITE containing small RNA and expressions of its closest genes

We scanned nr database at NCBI using a representative member of each silkworm MITE family as query and found a small RNA known as RNA-36850 (No. AB423040) perfectly matched to BmMITE-11 (Figure 5A). This small RNA is 28 bp long and was first discovered in silkworm ovary by Kawaoka et al. [30]. Previous studies showed that this small RNA (RNA-36850) belonged to a new class of germline-restricted small RNAs (26-33 bp) called piRNAs, which is thought to defend the host genome against transposons [31].

Figure 5
figure 5

BmMITE-11 containing the small RNA and expression profiles of its the closest genes. (A) The alignment of BmMITE-11 and small RNA (RNA-36850, Genbank accession no. AB423040). (B) Expression profiles of the closest genes to which BmMITE-11 is located. Red indicates upregulated genes, green down-regulated genes, black no expression change. a, b, c represent the closest genes with BmMITE-11 inserted into their introns, 3` flanking and 5`flanking regions, respectively.

To determine whether BmMITE-11 can regulate expressions of the closest genes through the small RNA, we examined expression profiles of the closest genes using the available microarray data [32]. First, we found the six closest genes of BmMITE-11: BGIBMGA009213, BGIBMGA001511, BGIBMGA013723, BGIBMGA013305, and BGIBMGA014223 in SilkDB http://silkworm.swu.edu.cn/silkdb/. With information of these six genes, corresponding expression profiles were extracted from the microarray data. Because BGIBMGA002049 has no information in microarray data, so it was not included in the study. Figure 5B showed that the five closest genes have different expression patterns. However, they have a common feature that all the five genes were not expressed in the ovary and testis. There are 10393 active genes in the silkworm microarray data. 8431 of 10393 genes were expressed in ovary or testis [32]. Thus, the probability of 5 random selected genes that were not expressed in ovary or testis is 0.000241 (P = (1- 8431/10393)5 = 0.000241). Chance cannot explain the observed expression pattern. Therefore, BmMITE-11 might self-silence germline-specifically by small RNA. As a result, BmMITE-11 might silence these five closest genes in silkworm ovary through small RNA pathway.

Discussion

Discovery and characterization of the silkworm MITEs

In this study, we carried out a systematic and genome-wide analysis to search for MITEs in silkworm using a novel computational approach, MUST [20]. First, 143333 MITEs were predicted and were grouped into 1350 families. Having further examined these predicted MITEs, however, we found that a large number of MITEs were pseudo-MITEs. After removing pseudo-MITEs by a strict filtering process, we identified only 17 silkworm MITE families. Of 17 families, Hoshidandy and Organdy have been reported previously [26]. However, 7 of 17 families are novel according to the classification of nucleotide composition of TSD. This suggests that the procedure for MITE identification used in this study is reliable.

Previous studies identified various MITEs in other insects. For example, eight MITE families were found in mosquito [8]. Our results indicated that silkworm harbors 17 MITE families and 7 of which are completely novel. Furthermore, the silkworm MITEs show huge diversification in both TSD composition and full length sequence (Table 1). However, all families have MITE-specific features: high number of copies, high AT content and potential to form a stable secondary structure. This also suggested validity of MUST in predicting the silkworm MITE. In addition, the verification of predicted MITEs by PCR and sequencing of the PCR products showed existence of MITEs.

It should be pointed out that a large number of candidate MITEs were predicted by MUST at first. However, it was found that most of these candidates were false positive after filtering out by a strict process. Probably, a high rate of false positive in MITE identification by MUST is due to incompleteness of the silkworm genome sequence, that is, there are many gaps to be sequenced in the genome. Again, that silkworm has a large proportion of repetitive sequences may affect correct identification of MITEs. In addition, a high sensitivity of MUST in identifying MITEs may be a reason. It can be known from considerations above that 17 MITE families contained in silkworm should be a conserved number. More silkworm MITEs are to be identified in future.

Expansion and diversity pattern

The ages of 17 MITE families vary greatly (Figure 3), ranging from 0 to 4 mya ago. Strikingly, BmMITE-2 might be dramatically amplified during a period from 0 mya to 1 mya (Figure 3B). A similar result has been reported that mPing underwent a dramatic amplification during rice domestication [33]. It is interesting to determine whether silkworm MITEs experienced recent expansions is due to bottleneck of silkworm domestication in future.

It should also be noted that the method used to estimate time of insertion of MITE copies in this study has an implicit assumption that all full sequences evolve neutrally and rate-equally [13]. However, this may not be true: Not all sequences are under neutral evolution; there may be differential rate of substitution among nucleotides as some regions may be important for transposition and there may be different rate of change in different genomic regions. All these will affect the accuracy of estimate time of MITE insertion. To solve these problems, a sophisticated estimation method needs to be established in future. Nevertheless, the results of MITE insertion time presented in this study should provide some useful information about the evolutionary dynamics of the silkworm MITEs.

Our diversity analyses suggested that members of each MITE family have a high level of similarity and similar sequence lengths. This shows that they can be considered as a population that experienced several successive steps of amplifications from a handful of master copies [13]. Right after amplification, copies are almost identical but will diverge in sequence and length over time due to random mutation. Therefore, a large number of identical or very similar copies within a MITE family imply a recent burst. Our network analyses revealed that 9 MITE families (BmMITE-1-3, BmMITE-5, BmMITE-8-10, BmMITE-12, and BmMITE-13) might experience several successive recent expansions (Figure 4). If these recent expansions were associated with the demographic events that silkworm experienced such as domestication bottleneck, then network topologies of above 9 MITE families can be explained by the genome-shock theory proposed by McClintock [34]. Thus, it is interesting to investigate dynamics of these silkworm MITE families during domestication in future.

Distribution of MITEs in the genome and their contribution to gene regulation

Our observations suggested that the silkworm MITEs are widely distributed in the genome. However, they are not randomly distributed on chromosomes (additional file 8). Since the silkworm MITEs preferentially insert into gene regions like MITEs in other higher organisms, a nonrandom distribution of the silkworm MITEs may be due to different gene densities on chromosomes. However, we did not find the significant correlation between the densities of genes and MITEs among the silkworm chromosomes (R2 = 0.055, P > 0.05, df = 26). Thus, the causes of nonrandom distribution of the silkworm MITEs on chromosomes are to be explored.

We observed that the silkworm MITEs in introns, 3'flank and 5'flank are in total much more than in exons (Table 2). It should be noted that the silkworm MITEs inserted into introns and 3' flanking regions significantly more than 5' flanking regions (χ2 test: P < 0.01, respectively). These observations could be explained by two reasons: MITE insertions in exons were rapidly purged out from a population because they are deleterious [12]; because many MITEs have been found to contain poly (A) signal [10], these MITEs are likely to be maintained in 3' flanking regions to act in regulation. Since introns have been suggested to harbor regulatory elements [35], much more silkworm MITE insertions into introns imply that the MITEs may play important roles in gene expression by changing regulatory motifs. In addition, that MITEs preferentially insert into gene regions provides a material basis for establishment of a TE-derived genetic regulatory network [36] and this is also an interesting topic to be studied in future.

A recent study tried to relate MITEs to biogenesis of their siRNAs in Solanaceae [7]. However, evidence for functional implications of MITEs in gene regulation through small RNA pathways is still lacking. Most importantly, we found a BmMITE-11 from which a silkworm ovarian small RNA (RNA-36850: No. AB423040) derived (Figure 5A). The expression data further suggested that the closest genes of BmMITE-11 were germline-specifically silenced (Figure 5B). Because this small RNA (RNA-36850) belongs to a new class of germline-restricted small RNAs (26-33 bp) called piRNAs [30, 31], our result may be the first evidence for effects of a MITE on its neighbor genes in transcriptional regulation through the sRNA pathway. The validation of this mechanism is underway. Given the high copy number of MITEs, many small RNAs and miRNAs derived from MITEs [14–16] and their preferential insertion into gene regions, it will be important to systematically account for the different mechanisms of MITEs and their potential functional roles in transcriptional regulation of genes.

Conclusions

Although MITEs in various higher organisms including mosquito and Drosophila as well as beetle have been extensively investigated, little is known about the genome-wide information of MITEs in the silkworm genome. We identified the 17 silkworm MITE families by using a recently developed algorism to scan the genome sequence. Silkworm has 17 MITE families, and furthermore, 7 of 17 families are completely novel based on the nucleotide composition of TSD. These results added new knowledge for understanding the evolution of MITEs. Importantly, we not only corroborated the preference of MITEs inserted into or near genes seen in the other genomes but also found that BmMITE-11 might silence its closest genes through sRNA pathway. These results emphasize the exceptional role of MITEs in the transcriptional regulation of genes and have general implications to understand the interaction between MITEs and their host genome.

Methods

Mining and characterization of MITEs

The new assembly of the silkworm genome sequence was downloaded from SilkDB http://silkworm.swu.edu.cn/silkdb[23]. First, a structure-based approach implemented in the MUST program was used to search the silkworm genome sequence for candidate MITEs that have characteristics [20]: TIR with 8-50 bp of length; TSD with 2-30 bp of length; 100-600 bp of the sequence length. The parameters were selected based on the common features of known MITEs. Candidate MITEs were grouped by all-blast-all. The MITEs for which any pair has an identity >0.8 were defined as a family. Families that have <3 members were excluded to reduce false positive rate for identification of MITEs.

Next, a homology BLAST search was used to scan the silkworm genome sequence with a representative member for estimating copy number of each MITE family. At this step, a MITE family was defined by a sequence similarity with E < e-5 for BLASTN. This standard has been recently used by Kuang et al. [7]. The BLAST results were filtered out using a Perl script according to criteria: a minimum nucleotide identity rate >90% and query coverage >80%. Those including both complete TIRs were defined as intact MITEs and others as fragmentary MITEs.

UNAFOLD http://www.bioinfo.rpi.edu/applications/hybrid/twostate-fold.php was used to predict secondary structures of MITEs [28]. Then, all MITE families were used against the ISfinder database http://www-is.biotoul.f/is.html[37], RepBase database (version 14.11) [38] and NCBI nr database respectively to find known families.

Verification of predicted MITEs

Fourteen silkworm accessions (Wu-B, Wu-D, DongDe-201, Wu-E, Lu-10, DaXianTuZhong, FuRongHuiLuan, PeiXianZhong, Ri-9, HeiZi, ChunSi, QiongShanHaiNan, Ri-110, ShangSanHuBan) that represent four main geographic strains (Chinese, Japanese, European, and Tropical) were used in insertion validation of a predicted MITE at four sites. A MITE used in insertion validation was randomly selected and the primers were designed based on flanking regions of each of insertion sites (additional file 4).

Estimates of insertion time and diversity of MITE families

To estimate the age of MITEs, DNA sequences of each MITE family were aligned using Clustal W [39] and the family consensus sequences were constructed using the program DAMBE [40], then the level of nucleotide substitutions (k) between each MITE element and the family consensus sequence was estimated using Kimura 2-parameter distance [41]. The each MITE age was estimated using the formula T = k/2r, assumed r = 1.56 × 10-8, it's the fruitfly neutral rate of substitutions per year and has been used in silkworm [22].

All full length sequences of each MITE family were aligned using Clustal W [39]. Then, using the alignment of all full length sequences of each family and the program Network4.5 http://www.fluxus-engineering.com/sharenet.htm[42], we constructed median-joining (MJ) networks to estimate BmMITE intrafamily diversity. This method has been used in maize [13]. Additionally, a sliding window analysis with a window size 20 bp and a step size 20 bp was used to look at pattern of regional nucleotide diversity (pi) for each family and completed by DnaSP (version5.10) [43].

Distributions of MITEs on chromosomes and Estimation of MITE richness in genes

All identified members of each MITE family were mapped to chromosomes by SilkMap http://silkworm.swu.edu.cn/silksoft/silkmap.html and copy numbers of each family on chromosomes were counted.

Two files for positions of predicted genes in scaffolds and for lengths of scaffolds were downloaded from SilkDB http://silkworm.swu.edu.cn/silkdb. Then, Perl script was written to scan the files for extracting the information of MITEs close to or in predicted genes. To determine whether insertions of MITEs close to gene regions are due to chance, a computer simulation, which is similar to that used by Naito et al. [33], was performed. Briefly, the fragments of up to 10 kb were randomly sampled from the silkworm genome sequence, the middle of each 10 kb sequence was presumed as the insertion site and the information about the insertion site close to or in predicted genes (i.e., in intron and exon) was recorded accordingly. The genes that have less than 5 kb distances to MITEs or contain MITEs were defined as the closest genes.

The closest genes were annotated by WEGO http://silkworm.swu.edu.cn/cgi-bin/wego/index.pl, and gene expression profiles were examined by using the available microarray data as described in Xia et al. [32].

References

  1. McClintock B: Chromosome organization and genic expression. Cold Spring Harbour Sym Quant Biol. 1951, 16: 13-47.

    Article  CAS  Google Scholar 

  2. Finnegan DJ: Eukaryotic transposable elements and genome evolution. Trends Genet. 1989, 5: 103-107. 10.1016/0168-9525(89)90039-5.

    Article  CAS  PubMed  Google Scholar 

  3. Feschotte C, Jiang N, Wessler SR: Plant transposable elements: where genetics meets genomics. Nat Rev Genet. 2002, 3: 329-341. 10.1038/nrg793.

    Article  CAS  PubMed  Google Scholar 

  4. Bureau TE, Wessler SR: Tourist: A large family of small inverted repeat elements frequently associated with maize genes. Plant Cell. 1992, 4: 1283-1294. 10.1105/tpc.4.10.1283.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  5. Bureau TE, Ronald PC, Wessler SR: A computer-based systematic survey reveals the predominance of small inverted repeat elements in wild-type rice genes. Proc Natl Acad Sci USA. 1996, 93: 8524-8529. 10.1073/pnas.93.16.8524.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  6. Santiago N, Herraiz C, Goni JR, Messeguer X, Casacuberta JM: Genome-wide analysis of the Emigrant family of MITEs of Arabidopsis thaliana. Mol Biol Evol. 2002, 19: 2285-2293.

    Article  CAS  PubMed  Google Scholar 

  7. Kuang H, Padmanabhan C, Li F, Kamei A, Bhaskar PB, Ouyang S, Jiang J, Buell CR, Baker B: Identification of miniature inverted-repeat transposable elements (MITEs) and biogenesis of their siRNAs in the Solanaceae: New functional implications for MITEs. Genome Res. 2009, 19: 42-56. 10.1101/gr.078196.108.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  8. Tu Z: Eight novel families of miniature inverted repeat transposable elements in the African malaria mosquito, Anopheles gambiae. Proc Natl Acad Sci USA. 2001, 98: 1699-1704. 10.1073/pnas.041593198.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  9. Smit AF, Riggs AD: Tiggers and DNA transposon fossils in the human genome. Proc Natl Acad Sci USA. 1996, 93: 1443-1448. 10.1073/pnas.93.4.1443.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. Bureau TE, Wessler SR: Stowaway: A new family of inverted repeat elements associated with the genes of both monocotyledonous and dicotyledonous plants. Plant Cell. 1994, 6: 907-916. 10.1105/tpc.6.6.907.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Jiang N, Bao Z, Zhang X, Hirochika H, Eddy SR, McCouch SR, Wessler SR: An active DNA transposon family in rice. Nature. 2003, 421: 163-167. 10.1038/nature01214.

    Article  CAS  PubMed  Google Scholar 

  12. Oki N, Yano K, Okumoto Y, Tsukiyama T, Teraishi M, Tanisaka T: A genome-wide view of miniature inverted-repeat transposable elements (MITEs) in rice, Oryza sativa ssp. Japonica. Genes Genet Syst. 2008, 83: 321-329. 10.1266/ggs.83.321.

    Article  CAS  PubMed  Google Scholar 

  13. Zerjal T, Joets J, Alix K, Grandbastien MA, Tenaillon MI: Contrasting evolutionary patterns and target specificities among three Tourist-like MITE families in the maize genome. Plant Mol Biol. 2009, 71: 99-114. 10.1007/s11103-009-9511-0.

    Article  CAS  PubMed  Google Scholar 

  14. Piriyapongsa J, Marino-Ramirez L, Jordan IK: Origin and evolution of human microRNAs from transposable elements. Genetics. 2007, 176: 1323-1337. 10.1534/genetics.107.072553.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  15. Piriyapongsa J, Jordan IK: A family of human microRNA genes from miniature inverted-repeat transposable elements. PLoS One. 2007, 2: e203-10.1371/journal.pone.0000203.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Piriyapongsa J, Jordan IK: Dual coding of siRNAs and miRNAs by plant transposable elements. RNA. 2008, 14: 814-821. 10.1261/rna.916708.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Feschotte C, Osterlund MT, Peeler R, Wessler SR: DNA-binding specificity of rice mariner-like transposases and interactions with Stowaway MITEs. Nucleic Acids Res. 2005, 33: 2153-2165. 10.1093/nar/gki509.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  18. Yang G, Nagel DH, Feschotte C, Hancock CN, Wessler SR: Tuned for transposition: Molecular determinants underlying the hyperactivity of a Stowaway MITE. Science. 2009, 325: 1391-1394. 10.1126/science.1175688.

    Article  CAS  PubMed  Google Scholar 

  19. Bao Z, Eddy SR: Automated de novo identification of repeat sequence families in sequenced genomes. Genome Res. 2002, 12: 1152-1155. 10.1101/gr.453102.

    Article  Google Scholar 

  20. Chen Y, Zhou F, Li G, Xu Y: MUST: A system for identification of miniature inverted-repeat transposable elements and applications to Anabaena variabilis and Haloquadratum walsbyi. Gene. 2009, 436: 1-7. 10.1016/j.gene.2009.01.019.

    Article  CAS  PubMed  Google Scholar 

  21. Mita K, Kasahara M, Sasaki S, Nagayasu Y, Yamada T, Kanamori H, Namiki N, Kitagawa M, Yamashita H, Yasukochi Y, Kadono-Okuda K, Yamamoto K, Ajimura M, Ravikumar G, Shimomura M, Nagamura Y, Shin-I T, Abe H, Shimada T, Morishita S, Sasaki T: The genome sequence of silkworm, Bombyx mori. DNA Res. 2004, 11: 27-35. 10.1093/dnares/11.1.27.

    Article  CAS  PubMed  Google Scholar 

  22. Xia Q, Zhou Z, Lu C, Cheng D, Dai F, Li B, Zhao P, Zha X, Cheng T, Chai C, Pan G, Xu J, Liu C, Lin Y, Qian J, Hou Y, Wu Z, Li G, Pan M, Li C, Shen Y, Lan X, Yuan L, Li T, Xu H, Yang G, Wan Y, Zhu Y, Yu M, Shen W, Wu D, Xiang Z, Yu J, Wang J, Li R, Shi J, Li H, Li G, Su J, Wang X, Li G, Zhang Z, Wu Q, Li J, Zhang Q, Wei N, Xu J, Sun H, Dong L, Liu D, Zhao S, Zhao X, Meng Q, Lan F, Huang X, Li Y, Fang L, Li C, Li D, Sun Y, Zhang Z, Yang Z, Huang Y, Xi Y, Qi Q, He D, Huang H, Zhang X, Wang Z, Li W, Cao Y, Yu Y, Yu H, Li J, Ye J, Chen H, Zhou Y, Liu B, Wang J, Ye J, Ji H, Li S, Ni P, Zhang J, Zhang Y, Zheng H, Mao B, Wang W, Ye C, Li S, Wang J, Wong GK, Yang H: A draft sequence for the genome of the domesticated silkworm (Bombyx mori). Science. 2004, 306: 1937-1940. 10.1126/science.1102210.

    Article  PubMed  Google Scholar 

  23. The International Silkworm Genome Consortium: The genome of a lepidopteran model insect, the silkworm Bombyx mori. Insect Biochem Mol Biol. 2008, 38: 1036-1045. 10.1016/j.ibmb.2008.11.004.

    Article  CAS  Google Scholar 

  24. Osanai-Futahashi M, Suetsugu Y, Mita K, Fujiwara H: Genome-wide screening and characterization of transposable elements and their distribution analysis in the silkworm, Bombyx mori. Insect Biochem Mol Biol. 2008, 38: 1046-1057. 10.1016/j.ibmb.2008.05.012.

    Article  CAS  PubMed  Google Scholar 

  25. Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Tu ZJ, Loftus B, Xi Z, Megy K, Grabherr M, Ren Q, Zdobnov EM, Lobo NF, Campbell KS, Brown SE, Bonaldo MF, Zhu J, Sinkins SP, Hogenkamp DG, Amedeo P, Arensburger P, Atkinson PW, Bidwell S, Biedler J, Birney E, Bruggner RV, Costas J, Coy MR, Crabtree J, Crawford M, Debruyn B, Decaprio D, Eiglmeier K, Eisenstadt E, El-Dorry H, Gelbart WM, Gomes SL, Hammond M, Hannick LI, Hogan JR, Holmes MH, Jaffe D, Johnston JS, Kennedy RC, Koo H, Kravitz S, Kriventseva EV, Kulp D, Labutti K, Lee E, Li S, Lovin DD, Mao C, Mauceli E, Menck CF, Miller JR, Montgomery P, Mori A, Nascimento AL, Naveira HF, Nusbaum C, O'leary S, Orvis J, Pertea M, Quesneville H, Reidenbach KR, Rogers YH, Roth CW, Schneider JR, Schatz M, Shumway M, Stanke M, Stinson EO, Tubio JM, Vanzee JP, Verjovski-Almeida S, Werner D, White O, Wyder S, Zeng Q, Zhao Q, Zhao Y, Hill CA, Raikhel AS, Soares MB, Knudson DL, Lee NH, Galagan J, Salzberg SL, Paulsen IT, Dimopoulos G, Collins FH, Birren B, Fraser-Liggett CM, Severson DW: Genome sequence of Aedes aegypti, a major arbovirus vector. Science. 2007, 316: 1718-1723. 10.1126/science.1138878.

    Article  CAS  PubMed  Google Scholar 

  26. Komoto N, Sezutsu H, Yukuhiro K, Banno Y, Fujii Y: Mutations of the silkworm molybdenum cofactor sulfurase gene, og, cause translucent larval skin. Insect Biochem Mol Biol. 2003, 33: 417-427. 10.1016/S0965-1748(03)00006-7.

    Article  CAS  PubMed  Google Scholar 

  27. Besansky NJ, Mukabayire O, Bedell JA, Lusz H: Pegasus, a small terminal inverted repeat transposable element found in the white gene of Anopheles gambiae. Genetica. 1996, 98: 119-129. 10.1007/BF00121360.

    Article  CAS  PubMed  Google Scholar 

  28. Markham NR, Zuker M: UNAFold: Software for nucleic acid folding and hybridization. Methods Mol Biol. 2008, 453: 3-31. full_text.

    Article  CAS  PubMed  Google Scholar 

  29. Jobling MA, Hurles ME, Tyler-Smith C: Human evolutionary genetics. Garland Publishing 2004, Abingdon and New York

  30. Kawaoka S, Hayashi N, Katsuma S, Kishino H, Kohara Y, Mita K, Shimada T: Bombyx small RNAs: Genomic defense system against transposons in the silkworm, Bombyx mori. Insect Biochem Mol Biol. 2008, 38: 1058-1065. 10.1016/j.ibmb.2008.03.007.

    Article  CAS  PubMed  Google Scholar 

  31. Saito K, Nishida KM, Mori T, Kawamura Y, Miyoshi K, Nagami T, Siomi H, Siomi MC: Specific association of Piwi with rasiRNAs derived from retrotransposon and heterochromatic regions in the Drosophila genome. Genes Dev. 2006, 20: 2214-2222. 10.1101/gad.1454806.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Xia Q, Cheng D, Duan J, Wang G, Cheng T, Zha X, Liu C, Zhao P, Dai F, Zhang Z, He N, Zhang L, Xiang Z: Microarray-based gene expression profiles in multiple tissues of the domesticated silkworm, Bombyx mori. Genome Biol. 2007, 8: R162-10.1186/gb-2007-8-8-r162.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Naito K, Cho E, Yang G, Campbell MA, Yano K, Okumoto Y, Tanisaka T, Wessler SR: Dramatic amplification of a rice transposable element during recent domestication. Proc Natl Acad Sci USA. 2006, 103: 17620-17625. 10.1073/pnas.0605421103.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  34. McClintock B: The significance of responses of the genome to challenge. Science. 1984, 226: 792-801. 10.1126/science.15739260.

    Article  CAS  PubMed  Google Scholar 

  35. Le Hir H, Nott A, Moore MJ: How introns influence and enhance eukaryotic gene expression. Trends Biochem Sci. 2003, 28: 215-220. 10.1016/S0968-0004(03)00052-5.

    Article  CAS  PubMed  Google Scholar 

  36. Jordan IK: Evolutionary tinkering with transposable elements. Proc Natl Acad Sci USA. 2006, 103: 7941-7942. 10.1073/pnas.0602656103.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M: ISfinder: The reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006, 34: D32-D36. 10.1093/nar/gkj014.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110: 462-467. 10.1159/000084979.

    Article  CAS  PubMed  Google Scholar 

  39. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001, 92: 371-373. 10.1093/jhered/92.4.371.

    Article  CAS  PubMed  Google Scholar 

  41. Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.

    Article  CAS  PubMed  Google Scholar 

  42. Bandelt HJ, Forster P, Rohl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.

    Article  CAS  PubMed  Google Scholar 

  43. Librado P, Rozas J: DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Dr. Tian Li for his help in drafting Perl scripts used in the study. We have also benefited from the discussion with Dr. Jun Duan when the study was launched. We also thank three reviewers for their constructive comments that improved the manuscript. This work was supported by the Hi-Tech Research and Development (863) Program of China (2006AA10A117), a grant from National Science Foundation of China (No. 30970409), and the Programme of Introducing Talents of Discipline to Universities (B07045).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Li-Yong Chen or Ze Zhang.

Additional information

Authors' contributions

ZZ and MJH designed the study. MJH did the data analyses and drafted the manuscript. YHS and YHG performed the experiments. LYC revised the manuscript. ZHX supervised the study. ZZ drafted and revised the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2010_3114_MOESM1_ESM.PDF

Additional file 1: The examples for Pseudo-MITEs: (A) undetermined fragments (designated as Ns in scaffolds), (B) solely composed of simple repeats, and (C) nested in repeats, arrow represents TSD and underline represents TIR. (PDF 1006 KB)

Additional file 2: The information about insertion sites of the MITEs into scaffolds. (PDF 185 KB)

12864_2010_3114_MOESM3_ESM.PDF

Additional file 3: The secondary structures of silkworm MITEs predicted by the program of UNADOLD: 1: BmMITE-1,2-16: BmMITE-3-17, respectively. (PDF 4 MB)

Additional file 4: Primers for PCR verification of BmMITE-2; (DOC 28 KB)

Additional file 5: The results for PCR verification of predicted BmMITE-2. (DOC 36 KB)

12864_2010_3114_MOESM6_ESM.PDF

Additional file 6: Sequencing result of PCR products for the presence or absence of BmMITE-2 at insertion site 3. "DXTZ", "SSHB" represent DaXianTuZhong and ShangSanHuBan strains, respectively. The box represents the sequence of BmMITE-2. (PDF 343 KB)

12864_2010_3114_MOESM7_ESM.PDF

Additional file 7: Nucleotide variation (π) along MITE sequence for each family. Nucleotide diversity π and its standard deviation were calculated by the program of DnaSP version 5.10. Then, the same program was used to define conserved and variable of the MITE sequences by sliding window analysis with both window size and step size 20 nucleotides. Windows with diversity equal or higher than average sequence diversity (π) + 2SD were defined as variable. Those with diversity less than average sequence diversity (π) - 2SD were considered as conserved. (PDF 26 KB)

12864_2010_3114_MOESM8_ESM.PDF

Additional file 8: Distribution of the silkworm MITEs on the 28 chromosomes. The observed distribution is significantly different from the expected one based on the total length of 28 chromosomes (Chi-square = 297, df = 27, P < 0.01). *The chromosomes that show the observed copies more than expected. (PDF 65 KB)

Additional file 9: MITE distances to the nearest genes. (DOC 56 KB)

Additional file 10: The annotation of the closest genes using WEGO. (PDF 559 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Han, MJ., Shen, YH., Gao, YH. et al. Burst expansion, distribution and diversification of MITEs in the silkworm genome. BMC Genomics 11, 520 (2010). https://doi.org/10.1186/1471-2164-11-520

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-11-520

Keywords