Transposable element evolution in Heliconius suggests genome diversity within Lepidoptera
© Lavoie et al.; licensee BioMed Central Ltd. 2013
Received: 19 June 2013
Accepted: 27 September 2013
Published: 2 October 2013
Transposable elements (TEs) have the potential to impact genome structure, function and evolution in profound ways. In order to understand the contribution of transposable elements (TEs) to Heliconius melpomene, we queried the H. melpomene draft sequence to identify repetitive sequences.
We determined that TEs comprise ~25% of the genome. The predominant class of TEs (~12% of the genome) was the non-long terminal repeat (non-LTR) retrotransposons, including a novel SINE family. However, this was only slightly higher than content derived from DNA transposons, which are diverse, with several families having mobilized in the recent past. Compared to the only other well-studied lepidopteran genome, Bombyx mori, H. melpomene exhibits a higher DNA transposon content and a distinct repertoire of retrotransposons. We also found that H. melpomene exhibits a high rate of TE turnover with few older elements accumulating in the genome.
Our analysis represents the first complete, de novo characterization of TE content in a butterfly genome and suggests that, while TEs are able to invade and multiply, TEs have an overall deleterious effect and/or that maintaining a small genome is advantageous. Our results also hint that analysis of additional lepidopteran genomes will reveal substantial TE diversity within the group.
KeywordsHeliconius melpomene Lepidopteran Butterfly Transposable elements Genomic deletions
Transposable elements (TEs) are segments of DNA that can mobilize in a genome. They impact the structure and function of the genomes they occupy. TEs can be divided into two classes. Class I TEs are the retrotransposons, which require an RNA intermediate and use a “copy and paste” mechanism to insert themselves into a new location in the genome. Retrotransposons are further divided into two groups, the long terminal repeat elements (LTRs) and non-LTR elements. LTR retrotransposons, such as members of the Gypsy and Copia superfamilies, are similar in structure to some retroviruses. Non-LTR retrotransposons lack LTR sequences and autonomous versions (Long INterspersed Elements or LINEs) usually harbor one or two open reading frames (ORFs) that are responsible for their mobilization. Examples include the LINE1, CR1, and RTE superfamilies and can be categorized into 28 monophyletic clades . Short INterspersed Elements (SINEs) are a group of nonautonomous non-LTR retrotransposons that are mobilized via the enzymatic machinery of LINEs .
Class II elements include the DNA transposons which use a “cut and paste” mechanism to mobilize in the genomes they occupy. Typically, DNA transposons require a transposase enzyme to recognize the terminal inverted repeats (TIRs) of the transposon and then excise and reinsert the element into another location in the genome . Examples of Class II elements include members of the TcMariner, hAT, and piggyBac superfamilies. There is a second group of Class II TEs known as the rolling circle transposable elements that includes the Helitrons .
The first lepidopteran to have its whole genome sequenced, the silkworm moth Bombyx mori, has accumulated a diverse array of retrotransposons and DNA transposons . For instance, a non-LTR retrotransposon, L1Bm, is abundant in the genome with copies of the 3′ end numbering ~25,000. However, like many LINEs most copies are 5′ truncated . Multiple copies of a piggyBac-like DNA transposon that may harbor an intact transposase have also been found in B. mori and it appears to have been recently active along with other Class II elements .
Recently, the genome of Heliconius melpomene was released , providing new insights into lepidopteran genome evolution from a transposable element perspective. H. melpomene is a heliconiine butterfly that is widespread throughout Central America and South America [8, 9]. The H. melpomene genome is the third lepidopteran and second butterfly genome to be sequenced. Unfortunately, the analysis of the second genome (and the first butterfly), the monarch, Danaus plexipus, was not comprehensive . Therefore, we confine our comparisons of the H. melpomene genome to B. mori.
Our analyses indicate that H. melpomene exhibits a high rate of TE turnover, with little accumulation of older elements, especially longer, autonomous elements, suggesting that TEs have an overall deleterious effect on the genome. Furthermore, the TE landscape of H. melpomene is distinct compared to the silkworm moth, consisting of substantially higher Class II content and a distinct set of retrotransposons. This suggests that lepidopterans in general will exhibit high levels of TE diversity as additional genomes are sequenced and characterized.
Summary of overall transposable element content in Heliconius melpomene
Identification of Metulj and its subfamilies
Age analyses and relative insertion rates
Divergence values and estimated activity periods for Metulj subfamilies
Counts of intact open reading frames for full length consensus sequences of each element class
Coordinates of ORF
120 - 1208
2980 - 4896
2051 - 4969
534 - 2924
95 - 1828
55 - 1530
543 - 2975
1468 - 4407
505 - 1986
111 - 1280
1411 - 4557
119 - 4207
616 - 3636
264 - 3233
1334 - 3874
723 - 1724
323 - 1639
69 - 1130
181 - 3144
1299 - 3611
172 - 3333
590 - 3517
13 - 4542
1071 - 5060
2741 - 4402
52 - 1716
2601 - 4148
84 - 2540
3008 - 4198
49 - 1272
1694 - 3433
1260 - 3167
1525 - 3489
DNA transposons exhibit a much different pattern of succession with multiple lineages exhibiting relatively recent activity (i.e. mean activity periods estimated within the last 2 my; Additional file 3: Table S2). Only three autonomous DNA transposon families were identified in the genome but one stands out. Tc3-1_Hm exhibits an average divergence of 0.002% among 113 full length insertions. A total of 43 intact ORFs are present, suggesting that this family is a recent and active addition to the TE repertoire of H. melpomene. However, no intact transposase ORFs other than Tc3-1_Hm were evident. A second standout is the Helitron superfamily, which also appears to have undergone a relatively recent amplification and is the most prevalent Class II element, occupying ~5% of the genome. Several other element families also appear to be young and active. These include multiple nonautonomous families of the piggyBac, Mariner, hAT and Helitron superfamilies and the two autonomous piggyBac elements. For the purposes of this study, MITEs (miniature inverted repeat transposable elements) were considered a subset of non-autonomous DNA transposons.
Evidence of TE removal
Evolutionary relationships among autonomous Non-LTR retrotransposons
The other two candidates were piggyBac-1_Hm and piggyBac-2_Hm with hits matching our criteria in Manduca sexta (piggyBac-1_Hm, 99% identity over 100% of the query, E-value = 0), Bombyx mori (and piggyBac-2_Hm, 99% identity over 100% of the query, E-value = 0), and D. plexipus (and piggyBac-2_Hm, 98% identity over 85% of the query, E-value = 0). In the case of D. plexipus, the reduced coverage is due to the fact that the insertion terminates with the scaffold (AGBW01001888).
TE content in Heliconius compared to Bombyx
The genome of Heliconius melpomene is the third lepidopteran genome to be fully sequenced. Unfortunately, the authors of the monarch genome manuscript did not complete a comprehensive analysis of the TE landscape , and our comparisons were therefore limited to B. mori.
TEs make up 35% of the B. mori genome, with the largest fraction (26.6%) being non-LTR retrotransposons . Of the non-LTR content, around half is derived from SINEs, 48%. A smaller fraction, ~25%, of the H. melpomene genome is composed of TEs. 12.5% consists of non-LTR retrotransposons, and 8% of the genome is occupied by SINEs (68% of the non-LTR content). Thirty-two non-LTR families belonging to 12 clades (Jockey, RTE, CR1, CRE, R1, R2, R4, I, Vingi, Daphne, Proto2 and L2) were identified and classified from B. mori. This is two more than were identified in H. melpomene. However, despite harboring two fewer clades than B. mori, the H. melpomene genome contains more families in total and this can be attributed to higher within-family diversity in some clades. For instance, 13 families of L2 and 10 families of CR1 were identified in H. melpomene, while only one and two are present in B. mori, respectively. Of the available lepidopteran genomes (including the monarch butterfly), Metulj is restricted to Heliconius.
In H. melpomene, LTR retrotransposons make up only ~0.45% of the genome. This is within the same range as what was described for B. mori by Osanai-Futahashi et al. in 2008 , 1.7%, but substantially different from a second estimate of LTR content in B. mori by Jin-Shan et al. , 11.8%. Given that Osanai-Futahashi examined a more complete assembly of the silkworm genome, we suspect that their estimate is closer to reality. Both genomes harbor Gypsy and Copia elements. B. mori however has two additional families which include Pao and Micropia . That being said, ~2.4% of the genome consists of candidate TEs that remain unidentified by our analyses and could belong in the LTR category.
While the retrotransposon content of B. mori and H. melpomene are similar, with regard to Class II elements, the DNA transposons, the two species are strikingly different in both content and quantity. Only ~3% of the B. mori genome consists of Class II elements  while ~10% of the H. melpomene genome is derived from DNA transposons. Indeed, the butterfly genome has been the subject of considerable DNA transposon activity within the recent past. This includes massive amplification by the Helitron superfamily and very recent, if not ongoing activity, from one member of the Tc-Mariner family. At least 43 intact members of the Tc3-1_Hm autonomous element are present in the genome draft and they are 99.4% identical, indicating that these elements are likely active.
Turnover of non-LTR element families in Heliconius
The lack of intact, older LINE elements in the genome suggests that they have a high fitness cost and that they may be preferentially removed. Mechanisms to accomplish removal include ectopic recombination between similar elements and removal of individual insertions via selection. Indeed, increased rates of ectopic recombination have been suggested as a mechanism for the differences in TE accumulation in both mammals and insects . Our results suggest that this mechanism is in play in the H. melpomene genome. Figure 3 indicates that deletions of large portions of LINE elements occur at relatively high frequency.
That being said, we note that other elements families have accumulated to relatively high numbers. In particular, this is true of Metulj and many of the Helitron elements. However, those elements with high copy numbers are typically under 500 nt in length. Previous authors have noted that shorter elements are likely less prone to recombination than their longer cousins [20, 21], allowing them to remain in the genome.
Hierarchical insertion patterns (TinT) indicate short periods of activity for the longer, autonomous elements, which exhibit a clear pattern of succession (Additional file 1: Figure S1). If one ignores the wide distributions of Metulj, the only SINE, each non-LTR family occupies a relatively narrow temporal space indicating that they experience brief periods of activity before ceasing mobilization. This is similar to what has been observed in some other taxa, including the lizard Anolis carolinensis, but is distinct from mammals, which have a single lineage of LINE-1 that has accumulated high copy numbers . The same analysis was performed for B. mori, with similar results (Additional file 1: Figure S1). Like many insects, the H. melpomene genome is relatively small, ~269 Mb. These results suggest that, while TE activity occurs and novel elements can invade the genome with some success, strong selection is working against the accumulation of large TEs and that homologous recombination acts to rapidly disable elements and keep the genome compact.
Evidence of horizontal transfer
We found evidence of horizontal transfer of three DNA transposons between H. melpomene and other taxa. Multiple elements matching nhAT10_Hm were identified in three taxa, the triatomine bug, Rhodnius prolixus, a strepsipteran insect, Mengenilla moldrzyki, and the planarian, Schmidtea mediterranea. In each case, the entire nhAT10_Hm is present as part of a larger element. For example, when compared to the planarian autonomous element, hAT-11_SM, nhAT10_Hm_has the hallmarks of an internal deletion variant. The first 70 bases are essentially identical between both TEs, as are the last 420 (Figure 5). The same regions overlap with as yet unnamed repeats in R. prolixus and M. moldrzyki. The top hit for R. prolixus can be found on contig ACPB02011601.1, nt 29253–30319, and the top hits for M. moldrzyki can be found on contigs AGDA01050831.1, nt 10068–10485 and AGDA01007612.1, nt 6860–6920, respectively. In these two taxa, the overlaps are with elements that are likely nonautonomous. This suggests that a hAT-22_SM-like element has been invading multiple genomes and produced similar nonautonomous variants in each. Indeed, we subsequently used BLASTN to query the genome drafts of H. melpomene, M. moldrzyki and R. prolixus using the consensus sequence of hAT-11_SM and, while no full-length elements were obvious, we identified high scoring (E-value = 0) hits from various portions of the consensus in each. Interestingly, both S. mediterranea and R. prolixus have been implicated in horizontal transfer previously [23–25].
The other candidates are the autonomous piggyBac elements, piggyBac-1_Hm and piggyBac-2_Hm. A single instance of piggyBac-1_Hm was identified in the Manduca sexta genome draft (scaffold AIXA01012877) with 99% identity over its entire length. Two full length copies of piggyBac-2_Hm in the Dazao strain of B. mori (scaffolds AADK01008943 and AADK01013248) with the same values. The final hit, to the monarch butterfly genome, is incomplete due to the termination of the scaffold ~350 bp prior to the end of the consensus. Both moths would have diverged from the lineage leading to butterflies ~145 mya  while the monarch is thought to have diverged from Heliconius ~89.79 mya  and, given the high rate of change observed in lepidopteran genomes, it is unlikely that they would have been conserved over such an extended period. This suggests to us that horizontal transfer explains their presence in each. However, as additional genomes are characterized this interpretation could change.
In conclusion, by conducting the first full TE analysis of a butterfly we have demonstrated that TEs, specifically SINEs and Helitrons, make up a large portion of the H. melpomene genome. We identified a novel SINE family which is found only in Heliconius and demonstrated that the genome of H. melpomene has experienced recent DNA transposon activity, most notably a Tc3 element. We have also shown that older, intact LINE elements are not found within the genome and that their activity period in the genome is short due to their rapid removal. Further studies of other lepidopteran genomes will be beneficial to our understanding of TEs in lepidopterans.
The genome sequence of a male Heliconius melpomene melpomene was recently described . Briefly, the specimen was acquired from Darien, Panama and the genome was sequenced using both 454 and Illumina platforms to generate a 38X draft genome. The sequenced male was inbred for five generations of sib mating. Repeat discovery was performed as summarized elsewhere  and described briefly here. Repetitive sequences in the H. melpomene draft sequence (Genbank accession number: CAEZ01000000) were identified de novo using RepeatModeler . To infer the consensus sequences for each repeat, we used the filtered RepeatModeler output to query the entire WGS draft using BLAST v2.2.23 . Up to fifty of the top hits spanning at least 100 bases were extracted along with up to 1,000 bases of flanking sequence, and we aligned the extracted sequences with MUSCLE 4.0  to generate 50% majority rule consensus sequences. Consensus sequences were considered ‘complete’ when single copy sequence could be identified at the 5′ and 3′ ends in each component sequence. If this condition was not met, the process was repeated until single copy DNA sequence was identifiable at both ends. The resulting library was submitted to CENSOR , BLASTN and BLASTX to ascertain the identity of the consensus with regard to previously classified elements. The result was a custom library of elements, which served as our library for subsequent analyses. The library of TEs was passed through a locally implemented version of RepeatMasker  to estimate the TE content of the H. melpomene genome.
Identification of SINE subfamilies
We identified 14,196 intact insertions of Metulj between 240–294 bases in length (+/− 10% of the general consensus) and passed them to COSEG [12, 13] for subfamily identification. COSEG examines multiple instances of TE insertions and identifies significant co-segregating (2–3 bp) sites in an effort to determine subfamily structure. A perl script provided by R. Hubley was used to refine the consensus sequence for each subfamily and is available upon request. We created a custom RepeatMasker library consisting of the suggested Metulj subfamily consensus sequences and extracted the top 150 hits for each from the genome. We aligned the extracted sequences with their respective subfamily consensus sequence to confirm the presence of each in the genome.
Identification of intact ORFs
We submitted the consensus sequence of each TE to NCBI ORF finder to identify potential open reading frames (ORFs). We classified any elements with identifiable ORFs spanning 1000 bp or more as potentially full length. ORF sequences were translated and BLASTP was used to confirm identity. ORFs of BEL-1_HMM, BEL-2_HMM, Copia-1_HMM, Gypsy-10_HMM, Gypsy-1_HMM, Gypsy-2_HMM, Gypsy-3_HMM, Gypsy-4_HMM, Gypsy-5_HMM, Gypsy-6_HMM, Gypsy-7_HMM, Gypsy-8_HMM, Gypsy-9_HMM as well as RTE-1_HMe, R4-1_Hme were identified by other parties and were obtained from RepBase.
We estimated the number of intact ORFs for each family of autonomous elements by passing the ORF sequences through a local version of TBLASTN, after which, up to 50 of the top hits based on bit score were extracted with 1000 bp of buffer and aligned. Extracted sequences were trimmed so they began and ended at the same position as the ORF query sequence. We defined an intact ORF as one that is greater than or equal to 90% of the expected amino acid length, contains a single, terminal stop codon, and begins with a methionine start codon.
Age analyses and relative insertion periods
We used the TinT online server (http://www.compgen.uni-muenster.de) as a method to determine periods of relative TE activity and succession patterns . Due to low copy numbers, analysis of LTR elements could not be performed. Furthermore, DNA transposons utilize a cut-and-paste mechanism of transposition that makes a nested insertion analysis of this type less informative. Thus, we analyzed only non-LTR retrotransposons.
We also estimated activity periods based on genetic distances between individual insertions and the consensus of each subfamily as described previously [33, 34]. Briefly, we created a modified TE library consisting of the full consensus of all Metulj subfamilies and non-autonomous DNA transposons, the full ORFs of all DNA transposons and 500 bp from the 3′ end of non-LTR ORFs. This library was then used to query the genome using RepeatMasker. We estimated Kimura2-parameter  distances (including CpG sites) between each insertion and its respective consensus . A neutral mutation rate is not available for H. melpomene. We applied an estimated mutation rate of 0.01909 substitutions per site/per million years which was taken from Papilioninae, a subfamily of the butterfly family Papilionidae .
The nearly vertical succession of non-LTR retrotransposons seen in the TinT plot (Additional file 1: Figure S1) suggests a rapid turnover of longer elements. One mechanism through which elements can be removed from a genome is non-homologous recombination leading to large deletions. By taking each RepeatMasker hit from each TE subfamily and mapping its location along the consensus element, we were able to examine decay patterns among selected elements.
Evolutionary relationships among autonomous non-LTR retrotransposons
From Genbank and Repbase, we collected non-LTR retrotransposon protein sequences from diverse known clades [37, 38]. We aligned these sequences with the consensus sequences retrieved from the H. melpomene genome using Clustal W in BioEdit . The most conserved region (about 300 amino acids) from the reverse transcriptase domain was identified and used in the phylogenetic analysis. Newly identified families missing this region were excluded. We inferred a maximum-likelihood tree with 1,000 bootstrap replicates using MEGA5 .
We investigated the taxonomic distribution of all H. melpomene TEs by querying the full WGS database at NCBI with BLAST. We considered an element to be a likely candidate for HT if a BLASTN search indicated that the consensus shared >95% sequence identity over at least 80% of its length. Any hits matching these criteria were examined by extracting the highest scoring hits, alignment to the query sequence and manual examination.
Non-long terminal repeat
Long terminal repeat
Long INterspersed Elements
Short INterspersed Elements
Open reading frame
Terminal inverted repeats
Transposition in Transposition.
We thank Federico Hoffmann, Ananya Sharma, and Mike Vandewege for support and ideas. We thank Alex Suh for his help with TinT and Robert Hubley for providing help with COSEG.
This work was supported by the National Science Foundation [MCB-0841821 and DEB-1020865 (DAR), NSF DEB-1257839 (BAC)]. Additional support was provided by the College of Agriculture and Life Sciences at Mississippi State University.
- Kapitonov V, Tempel S, Jurka J: Simple and fast classification of non-LTR retrotransposons based on phylogeny of their RT domain protein sequences. Gene 2009, 448: 207-213. 10.1016/j.gene.2009.07.019PubMed CentralView ArticlePubMed
- Dewannieux M, Esnault C, Heidmann T: LINE-mediated retrotransposition of marked Alu Sequences. Nat Genet 2003, 35: 41-48. 10.1038/ng1223View ArticlePubMed
- Biemont C, Vieira C: Junk DNA as an evolutionary force. Nature 2006, 443: 521-524. 10.1038/443521aView ArticlePubMed
- Kapitonov V, Jurka J: Rolling-circle transposons in eukaryotes. Proc Natl Acad Sci U S A 2001, 98: 8714-8719. 10.1073/pnas.151269298PubMed CentralView ArticlePubMed
- 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.012View ArticlePubMed
- Ichimura S, Mita K, Sugaya K: A Major Non-LTR Retrotransposon of Bombyx mori , L1Bm. J Mol Evol 1997, 45: 253-264. 10.1007/PL00006228View ArticlePubMed
- Daimon T, Mitsuhiro M, Katsuma S, Abe H, Mita K, Shimada T: Recent transposition of yabusame, a novel piggyBac-like transposable element in the genome of the silkworm, Bombyx mori . Genome 2010, 53: 585-593. 10.1139/G10-035View ArticlePubMed
- Consortium HG: A butterfly genome reveals promiscuous exchange of mimicry adaptations among species. Nature 2012, 487: 94-98.
- Brower AVZ: Parallel race formation and the evolution of mimicry in Heliconius butterflies: A phylogenetic hypothesis from mitochondrial DNA sequences. Evolution 1996, 50: 195-221. 10.2307/2410794View Article
- Zhan S, Merlin C, Boore J, Reppert S: The Monarch butterfly genome yields insights into long-distance migration. Cell 2011, 147: 1-15.View Article
- Okada N, Shedlock A, Nikaido M: Retroposon Mapping in Molecular Systematics. In Mobile Genetic Elements: Protocols and Genomic Applications, Volume 260. Edited by: Miller W, Pierre C. NJ: Humana Press; 2004:189-226.View Article
- Price AL, Eskin E, Pevzner PA: Whole-genome analysis of Alu repeat elements reveals complex evolutionary history. Genome Res 2004, 14: 2245-2252. 10.1101/gr.2693004PubMed CentralView ArticlePubMed
- COSEG 0.2.1. http://www.repeatmasker.org
- Churakov G, Grundmann N, Kuritzin A, Brosius J, Makalowski W, Schmitz J: A novel web-based TinT application and the chronology of the primate Alu retroposon activity. BMC Evol Biol 2010, 10: 376. 10.1186/1471-2148-10-376PubMed CentralView ArticlePubMed
- Ustyugova SV, Amosova AL, Lebedev YB, Sverdlov ED: Cell line fingerprinting using retroelement insertion polymorphism. Biotechniques 2005, 38: 561-565. 10.2144/05384ST02View ArticlePubMed
- Blass E, Bell M, Boissinot S: Accumulation and rapid decay of non-LTR retrotransposons in the genome of the three-spine stickleback. Genome Biol Evol 2012, 4: 687-702. 10.1093/gbe/evs044PubMed CentralView ArticlePubMed
- Novick P, Basta H, Floumanhaft M, McClure M, Boissinot S: The evolutionary dynamics of autonomous non-LTR retrotransposons in the lizard Anolis carolinensis shows more similarity to fish than mammals. Mol Biol Evol 2009, 26: 1811-1822. 10.1093/molbev/msp090View ArticlePubMed
- Jin-Shan X, Qing-You X, Jun L, Guo-Qing P, Ze-Yang Z: Survey of long terminal repeat retrotransposons of domesticated silkworm (Bombyx mori). Insect Biochem Mol Biol 2005, 35: 921-929. 10.1016/j.ibmb.2005.03.014View ArticlePubMed
- Eickbush T, Furano A: Fruit flies and humans respond differently to retrotransposons. Genomes Evol 2002, 12: 669-674.
- Cooper J, Watanabe Y, Nurse P: Fission yeast Taz1 protein is required for meiotic telomere clustering and recombination. Nature 1998, 392: 828-831. 10.1038/33947View ArticlePubMed
- Song M, Boissinot S: Selection against LINE-1 retrotransposons results principally from their ability to mediate ectopic recombination. Gene 2007, 390: 206-213. 10.1016/j.gene.2006.09.033View ArticlePubMed
- Furano A, Duvernell D, Boissinot S: L1 (LINE-1) retrotransposon diversity differs dramatically between mammals and fish. Trends Genet 2004, 20: 9-14. 10.1016/j.tig.2003.11.006View ArticlePubMed
- Garcia-Fernandez J, Bayascas-Ramirez JR, Marfany G, Munoz-Marmol AM, Casali A, Baguna J, Salo E: High copy number of highly similar mariner-like transposons in planarian (Platyhelminthe): evidence for a trans-phyla horizontal transfer. Mol Biol Evol 1995, 12: 421-431.PubMed
- Gilbert C, Schaack S, Pace JK 2nd, Brindley PJ, Feschotte C: A role for host-parasite interactions in the horizontal transfer of transposons across phyla. Nature 2010, 464: 1347-1350. 10.1038/nature08939PubMed CentralView ArticlePubMed
- Novick P, Smith J, Ray D, Boissinot S: Independent and parallel lateral transfer of DNA transposons in tetrapod genomes. Gene 2010, 449: 85-94. 10.1016/j.gene.2009.08.017View ArticlePubMed
- Tu Z: Three novel families of miniature inverted-repeat transposable elements are associated with genes of the yellow fever mosquito, Aedes aegypti. Proc Natl Acad Sci U S A 1997, 94: 7475-7480. 10.1073/pnas.94.14.7475PubMed CentralView ArticlePubMed
- Chan D, Kim P: HIV entry and its inhibition. Cell 1998, 93: 681-684. 10.1016/S0092-8674(00)81430-0View ArticlePubMed
- RepeatModeler Open-1.0. http://www.repeatmasker.org/RMDownload.html
- Altschul S, Madden T, Schaffer A, Zhang J, Zhang Z, Miller W, Lipman D: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucl Acids Res 1997, 25: 3389-3402. 10.1093/nar/25.17.3389PubMed CentralView ArticlePubMed
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucl Acids Res 2004, 32: 1792-1797. 10.1093/nar/gkh340PubMed CentralView ArticlePubMed
- Kohany O, Gentles A, Hankus L, Jurka J: Annotation, submission and screening of repetitive elements in Repbase: RebaseSubmitter and Censor. BMC Bioinforma 2006, 7: 107. 10.1186/1471-2105-7-107View Article
- RepeatMasker Open-3.0. http://www.repeatmasker.org
- Pagan H, Smith J, Hubley R, Ray D: PiggyBac-ing on a primate genome: Novel elements, recent activity and horizontal transfer. Genome Biol Evol 2010, 2: 293-303. 10.1093/gbe/evq021PubMed CentralView ArticlePubMed
- Ray D, Feschotte C, Pagan H, Smith J, Pritham E, Arensburger P, Atkinson P, Craig N: Multiple waves of recent DNA transposon activity in the bat, Myotis lucifugus . Genome Res 2008, 18: 717-728. 10.1101/gr.071886.107PubMed CentralView ArticlePubMed
- 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/BF01731581View ArticlePubMed
- Simonsen T, Zakharov E, Djernaes M, Cotton A, Vane-Wright R, Sperling F: Phylogenetics and divergence times of Papilioninae (Lepidoptera) with special reference to the enigmatic genera Teinopalpus and Meandrusa. Cladistics 2010, 26: 1-25. 10.1111/j.1096-0031.2009.00297.xView Article
- Benson DA: Genbank. Nucleic Acids Res 2010, 38: D46-D51. 10.1093/nar/gkp1024PubMed CentralView ArticlePubMed
- 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/000084979View ArticlePubMed
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 1999, 41: 95-98.
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 2011, 28: 2731-2739. 10.1093/molbev/msr121PubMed CentralView ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.