Functional characterization of the active Mutator-like transposable element, Muta1 from the mosquito Aedes aegypti
© The Author(s). 2017
Received: 21 October 2016
Accepted: 19 December 2016
Published: 11 January 2017
Mutator-like transposable elements (MULEs) are widespread with members in fungi, plants, and animals. Most of the research on the MULE superfamily has focused on plant MULEs where they were discovered and where some are extremely active and have significant impact on genome structure. The maize MuDR element has been widely used as a tool for both forward and reverse genetic studies because of its high transposition rate and preference for targeting genic regions. However, despite being widespread, only a few active MULEs have been identified, and only one, the rice Os3378, has demonstrated activity in a non-host organism.
Here we report the identification of potentially active MULEs in the mosquito Aedes aegypti. We demonstrate that one of these, Muta1, is capable of excision and reinsertion in a yeast transposition assay. Element reinsertion generated either 8 bp or 9 bp target site duplications (TSDs) with no apparent sequence preference. Mutagenesis analysis of donor site TSDs in the yeast assay indicates that their presence is important for precise excision and enhanced transposition. Site directed mutagenesis of the putative DDE catalytic motif and other conserved residues in the transposase protein abolished transposition activity.
Collectively, our data indicates that the Muta1 transposase of Ae. aegypti can efficiently catalyze both excision and reinsertion reactions in yeast. Mutagenesis analysis reveals that several conserved amino acids, including the DDE triad, play important roles in transposase function. In addition, donor site TSD also impacts the transposition of Muta1.
KeywordsTransposable elements Mutator-like elements (MULE) Aedes aegypti Yeast assay Target site duplication (TSD) Transposase
Transposable elements (TEs) are mobile fragments of DNA that can move from one locus to another in the host genome, often replicating in the process. TEs usually make up the largest fraction of eukaryotic genomes; accounting for almost half of the human genome, and more than 70% of the genomes of some grass species . Based on their transposition intermediate, eukaryotic TEs are divided into two classes. Class 1 elements utilize an RNA intermediate in the transposition reaction while the intermediate for most class 2 elements is the element itself that is mobilized by a 'cut and paste' mechanism . A TE family is composed of one or more transposase-encoding autonomous elements and up to several thousand nonautonomous elements that do not encode functional transposase. Family members share the same terminal inverted repeats (TIRs) and target site duplications (TSDs) of the same length allowing them to move in trans by utilizing the transposase encoded by the autonomous family member(s) .
Prior studies have classified class 2 TEs into upwards of 19 superfamilies on the basis of the relatedness of the element-encoded transposase . The original Mutator element, now called MuDR, was first isolated from a maize strain as the agent responsible for its high forward mutation rate [5, 6]. Subsequently, members of this superfamily, called collectively Mutator-like transposable elements (MULEs), were found in other plants, and in fungi, protozoans, and in a variety of animals, (from insects to fish to other metazoans) [3, 7–9]. Typical features of MULEs include long terminal inverted repeats (TIRs) (>100 bp) and an 8-10 bp TSD . Nonautonomous family members often contain a variety of sequences between the TIRs, including fragments from host genes; such elements are called Pack-MULEs . To date, only a few active MULEs have been identified, including MuDR, Hop1, Jittery, TED, and Os3378 [6, 12–15]. Importantly, only the rice Os3378 element has been shown to transpose in a heterologous host .
Most MULE transposases contain a N-terminal zinc finger DNA binding motif and a conserved C-terminal DDE domain, which has been shown to be the catalytic core for transposition reactions in other superfamilies but not to date for MULEs [16, 17]. Phylogenetic studies indicate that the DDE domain of MULE transposases is closely related to the IS256 family, which is present in diverse prokaryotes . Each residue of the predicted DDE motif of IS256 has been shown to be necessary for transposition . Most MULE transposases also harbor a CH motif and a W residue between the second D and E of the DDE domain, which are also conserved in the hAT superfamily. Mutagenesis analyses in hAT elements demonstrated the importance of these residues for transposition [20, 21]. Whether these residues are also important for the MULE transposase has not as yet been determined.
TE-based genetic tools have facilitated our deep understanding of the biology of both plants and animals, however, such tools are not currently available in the mosquito species, including Aedes aegypti, Aedes albopictus, Culex quinquefasciatus and Anopheles gambiae, which can spread dengue fever, yellow fever, chikungunya, zika, and many other diseases that are responsible for over one million deaths annually [22–25]. Four exogenous TEs that transpose at high frequencies in Drosophila melanogaster (Hermes, Mos1, Minos, piggyBac) were found to move rarely or not at all in the germ line of mosquitoes as very few germinal mutations were detected [26, 27]. To explain this result, it was hypothesized that mosquitoes have a strong genome defense system that could effectively recognize and silence foreign TEs . Therefore endogenous active TEs could be effective mutagens from generation to generation as they might be able to evade genome surveillance. Availability of the genome sequence enabled the search for potentially active TEs in the mosquito genomes [29–32]. In this study we show that Muta1, identified by computer-assisted analysis of the Ae. aegypti genome, is an active TE capable of both insertion and excision in a yeast transposition assay. Site-directed mutagenesis analysis revealed that transposition activity in yeast was influenced by disruption of several conserved residues and by the presence of TSDs at the donor site. With characteristics such as high transposition activity, precise excision, and no target sequence preference, Muta1 could be crafted into an effective tool for forward mutation analyses in mosquitoes.
Distribution of the MULE superfamily in mosquito genomes
To identify potentially active MULEs, we first assessed the abundance and distribution of MULEs in the genomes of Ae. aegtpti, Ae albopictus, C quinquefasciatus, and An. gambiae [29–32]. The conserved MULE DDE domain sequences were used as queries in TBlastN searches against the Ae. aegypti genome through the TARGeT program, which is designed for TE discovery [20, 33]. After removal of duplicate hits, sequences with significant similarity (e- value <10-15) to the MULE DDE domain were identified and used to build a phylogenetic tree. As a result, no significant hits were identified in the An. gambiae genome, whereas 141, 105, and 10 putative MULE DDE domain sequences were identified from the Ae. aegtpti, Ae albopictus, and C quinquefasciatus genomes, respectively (Additional file 1: Figure S1).
Characterization of the MULE superfamily in Ae. aegypti
All of the 14 MULE families identified from Ae. aegypti have 8-9 bp TSD, 12 of 14 have TIRs >100 bp and two have TIRs <60 bp. Ten of the 12 long TIR families contain short subterminal tandem repeats (9-15 bp). All families except Muta6 also contain derivative nonautonomous elements (Additional file 2: Table S1), which share high sequence similarity in their TIR and subterminal regions, but carry heterogeneous internal sequences including fragments of host genes.
The phylogenetic tree also reveals the evolutionary relationships of the 14 families. Among the 5 major groups resolved on the tree, group B includes half of the 14 families, while group D and E contains one family each (Fig. 1). Furthermore, long branches indicate extensive sequence differences between families. For example, although Muta5 and Muta8 belong to group A and locate to adjacent clades, no significant nucleotide similarity was detected between these elements. In another example, nucleotide similarity between Muta3 and Muta11 is restricted to the DDE region (71% identity, ~390 bp) and the TIRs (76% identity, ~ 110 bp).
After removal of other TE insertions and manual correction of frameshifts caused by small insertions and deletions, the 31 full-length elements were predicted to encode transposase proteins ranging in size from 416 to 554 residues (Additional file 2: Table S1). Comparison of conceptually translated transposases identified conserved regions, other than the DDE domain, that could have functional significance. The transposases of 13 families were predicted to harbor a FLYWCH type zinc finger DNA binding domain in the N-terminus, while Muta6 was predicted to harbor a SWIM type zinc finger DNA binding domain in the C-terminus. Each family contains at least one putative full-length copy in which the coding region is not interrupted by frameshifts or stop codons. Of particular interest to this study, 3 families (Muta1, Muta3 and Muta5) include at least 2 members that are identical or nearly identical (>99% identity). These features indicate recent and possible ongoing activity of multiple MULE families in Ae. aegypti.
Identification of the active Muta1 family
Muta1 can transpose in yeast
A yeast transposition assay was employed to determine whether Muta1 transposase is able to catalyze the movement of natural and/or artificial nonautonomous elements. In prior studies members of five superfamilies (MULE, Tc1/mariner, hAT, PIF/Harbinger and piggyBac) were shown to transpose in yeast [15, 35–38]. Our yeast assay consisted of two plasmid constructs: an expression vector containing the Muta1 transposase coding sequence downstream of the galactose inducible GAL1 promoter, and a reporter vector containing a nonautonomous element inserted in, and blocking expression of, the ADE2 gene (Additional file 4: Figure S3A). We first tested the artificial Muta1AR element which contains 250 bp from each end of Muta1 (Fig. 3). Transposase mediated TE excision restored ADE2 expression and resulted in reversion of adenine auxotrophy to prototrophy (Additional file 4: Figure S3B, Additional file 5: Figure S4A). Excision events were validated by PCR amplification of the ADE2 empty site from revertants (Additional file 6: Figure S5A).
In subsequent experiments, the reporter plasmid was modified by inserting several natural nonautonomous elements (Muta1NA1-5, Fig. 3, Additional file 5: Figure S4B-F) into the ADE2 coding region. Despite differences in TIR length and the length and sequence of internal regions, Muta1 transposase was able to mobilize all of these elements, but none with a frequency as high as the artificial element Muta1AR.
Muta1 integration in yeast
2.44 × 10−5 Ade+ cells/total cells
His+ 5-FOAR excisant
2.71 × 10−6 His+ 5FOAR cells/total cells
ratio of reintegration = 11%
His+ 5-FOAR excisant
ratio of reintegration = 14%
To determine the precise insertion sites, polymorphic bands on transposon display gels were recovered , sequenced and mapped to yeast chromosomes (Additional file 6: Figure S5B). Sixty of 62 sites had significant matches with yeast genomic sequence, while two matched plasmid sequences. Fourteen of the 60 genomic sites were in gene bodies while 27 were within 1kb of genes. Thus 41/60 insertions (68%) were in gene-rich regions . Amplification and sequencing of each site revealed the presence of TSDs with 8 bp for 21 events and 9 bp for 39 events (Additional file 7: Table S2). Consensus sequences generated from the yeast insertion sites and the 171 sites for Muta1 derivative elements in the Ae. aegypti genome indicated little or no sequence preference for insertion (Additional file 8: Figure S6).
The Muta5 element does not transpose in yeast
The successful transposition of Muta1 in yeast prompted us to perform a similar analysis of the Muta5 family which contains 3 identical copies, a 4th with >99% sequence identity, and 3 copies with large deletions or insertions (Fig. 1). Muta5 is 3496 bp, flanked by TIRs of 151 bp and TSDs of 8 bp or 9 bp and is predicted to contain 3 exons that encode a 554-amino acid transposase. The TIR of Muta5 contains 9 copies of a 15 bp subterminal tandem repeat and an 8bp terminal motif (Additional file 9: Figure S7A). There are over 200 Muta5 derivative elements (e- value <10-15) in the Ae. aegypti genome, with 93 flanked by 8 bp or 9 bp TSD. Taken together the features of the Muta5 family strongly suggested that it was an active element. However, in the yeast transposition assay the Muta5 transposase was unable to catalyze transposition of Muta5AR (Additional file 9: Figure S7A), which contains 250 bp from the ends of Muta5 (Additional file 9: Figure S7B).
TSD at donor site affects Muta1 transposition in yeast
Impact of TSD on transposition
frequency (X 10 -5)
11.03 ± 1.13
13.9 ± 1.8
63.34 ± 9.12
20.88 ± 1.53
76.33 ± 11.53
18.22 ± 2.90
68.53 ± 12.16
19.15 ± 2.39
31.52 ± 5.27
27.83 ± 3.01
36.43 ± 6.12
24.35 ± 1.98
28.38 ± 3.11
29.01 ± 2.76
We next addressed the question of whether the presence or absence of TSDs at the donor site impacts the quality of excision by analyzing so-called transposon footprints. Specifically, class 2 TEs often leaves a footprint upon excision consisting of a few nucleotides or small rearrangements at the site of excision site . Formation of footprints involves DNA repair of sequences flanking the excised element. To assess the impact of TSDs on the repair of excision sites, the donor element construct was modified so that all excision footprints (not only those that maintain the reading frame) could be analyzed. First, nonautonomous elements were inserted in the 5’ UTR of ADE2 (Additional file 4: Figure S3C). Second, because insertion of Muta1AR in the 5’UTR resulted in leaky ADE2 expression, we substituted the longer Muta1NA1 (Fig. 3), which, in the absence of transposase, blocked ADE2 expression. When either 8 bp or 9 bp TSDs were added to Muta1NA1, about 90% of revertants were precise (Table 2), meaning that the element was removed as well as a single copy of the TSD. The quality of excision appeared to be independent of TSD sequence (Additional file 10: Figure S8B-G). In contrast, absence of donor site TSDs reduced perfect excision to only 10% of the ADE2 revertant colonies sequenced (Table 2). Most excision sites reflected loss of a few nucleotides from either side of the flanking DNA. Occasionally, part of the TIR (up to 13 bp) was left after excision and repair (Additional file 10: Figure. S8A).
Mutagenesis of Muta1 transposase
The functional significance of two additional highly conserved residues in Muta1, H307 and W401, were tested. In a prior study these residues were shown to be essential for transposition of the hAT superfamily member Hermes [21, 43]. Mutation of the corresponding Hermes H and W residues (H268 and W319) to alanine completely abolished transposition activity. Similarly for the Muta1 transposase, mutation of H307 or W401 to alanine abolished activity (Fig. 4b). Analysis of Hermes transposase also showed that the W319 residue was likely necessary for the correct positioning of flanking DNA during the excision reaction, and that other aromatic residues can partially substitute for this function . When W401 of Muta1 was mutated to phenylalanine, transposition activity was reduced by 79% (Fig. 4b, Additional file 5: Figure S3L) and the frequency of precise excision dropped from 90% to approximately 48% (8 bp donor site TSD, 14/29 events) and 50% (9 bp donor site TSD, 15/30 events) (Additional file 10: Figure S8H&I), suggesting that the Muta1 W401F mutation also led to inaccurate positioning. These results confirmed the importance of the putative DDE motif, the conserved H307, W401 and identified the nonconserved E373 as a potentially important residue for transposase function.
The MULE superfamily is widespread in eukaryotic genomes and is closely related to prokaryotic IS256 elements. However, although it is also found in the genomes of many insects, no active elements have been reported in insect species. In this study, we performed a thorough search for potentially active MULEs in the Ae. aegypti genome and demonstrate that Muta1 encodes a transposase that catalyzes the excision and reinsertion of nonautonomous derivative elements in yeast. With the availability of this heterologous transposition assay, the function of the conserved MULE DDE domain and the role of TSD in transposition were tested.
The DDE/D domain is proposed to be the catalytic core involved in transposition of class 2 TEs and has been identified in all superfamilies . Prior to this study, the functional significance of this domain had been experimentally validated only for members of the piggyBac, Tc1/Mariner, and hAT superfamilies [38, 42, 43]. Results of this study provide the first experimental evidence for the importance of the DDE motif in the transposition reaction of Muta1, a member of the MULE superfamily. Specifically, transposition was completely abolished when any of the three residues were mutated to alanine (Fig. 4b, Additional file 5: Figure S3I-K).
In addition to the DDE triad, other residues critical for transposition were identified. W401 of Muta1 is a conserved residue that is also found in the hAT superfamily [20, 43]. Crystallographic analysis and in vitro biochemical assays showed that the corresponding W318 residue in Hermes functions in the positioning of flanking DNA, which ensures that the double strand break occurs at the correct position when an element excises from flanking DNA . Other aromatic residues partially substitute for its function, however the mutant transposase generated additional species of intermediates in double strand break repair than the wild-type transposase, suggesting that inaccuracy in the position of the cleavage site may be the cause by these mutations . For Muta1, the W401A mutation completely abolished transposition activity while the W401F mutation resulted in a 79% reduction of transposition frequency (Fig. 4b, Additional file 5: Figure S3L) and caused inaccurate excision as the frequency of precise excision dropped from 90% to ~50% (Additional file 10: Figure S8H&I). Taken together these data suggest that this conserved tryptophan residue is likely playing a similar role in hAT and MULE transposases, which is to correctly position flanking DNA for the excision reaction. In addition to the W residue, a CxxH motif is also shared between MULE and hAT elements, the Hermes H268 was found to be located close to the DDE active center and involved in the interaction with TIRs . Mutation of the corresponding H307 to alanine resulted in a 99.5% reduction of Muta1 transposition activity (Fig. 4b), suggesting the importance of this residue for Muta1 transposase function. Taken together, our study provides experimental evidence to support the close evolutional relationship reported previously between the MULE and hAT superfamilies .
Comparison of Muta1 and Os3378 transposition in yeast
Excision frequency no TSD (X 10-7)a
Excision frequency 8bp TSD (X 10-7)a
Excision frequency 9bp TSD (X 10-7)a
Reintegration frequency No TSDb
Reintegration frequency 8bp TSDb
Reintegration frequency 9bp TSDb
Percentage of reinsertion in gene rich regions
The presence of donor site TSDs impacts the quality of Muta1-mediated excision events as various footprints were generated without donor TSD (Additional file 10: Figure S8A). The predominance of small deletions (1-4 bp) suggests that the Muta1 transposase cuts outside the TIR. In contrast, with the presence of either 8 bp or 9 bp TSD at the donor site, most excision events were precise and the actual TSD sequence did not seem to matter (Additional file 10: Figure S8 B-G). Similar behavior was also observed for IS256, the prokaryotic TE family related to MULEs, and for the one other MULE tested, Os3378. Reduction of TSD from 8bp to 6bp eliminated precise excision of IS256 and reduced the Os3378 precise excision frequency from 97.44% to 82.05% [15, 19]. For IS256 it was hypothesized that precise excision is achieved through a transposase-independent replication slippage mechanism that requires a short stretch of homologous DNA with a minimum length of 8 bp . In our assay, the absence of donor TSDs resulted in a 90% reduction in precise excision (Additional file 10: Figure S8A-G), which suggests a role for TSDs in promoting precise excision.
Thirty-one full-length MULEs that group into 14 families were identified in the Ae. aegypti genome. Several families have identical or nearly identical full-length copies including Muta1 (7 identical), Muta5 (3 identical), and Muta3 (2 with only 2 noncoding SNPs). Although the existence of identical genomic copies is a feature of active TEs, Muta5 was unable to catalyze the movement of nonautonomous derivative elements in yeast (Additional file 9: Figure S7). One explanation for our success with Muta1 but failure with Muta5 is that the latter has 3 predicted exons while the former has 2. More predicted exons would increase the chances of incorrectly assembling the actual/functional Muta5 transposase.
Accumulation of seven identical copies of Muta1 in Ae. aegypti suggests that this element may still be active or that it has some success evading the genome surveillance system shown previously to effectively silence exogenous TEs [26–28]. In this regard, it may be possible to engineer Muta1 to make it an effective endogenous mutagen. Like the MuDR system in maize, where the genome has numerous copies of nonautonomous Mu elements, there are over 300 Muta1 derivative elements in the Ae. aegypti genome (Additional file 3: Figure S2A). Although mobility in the yeast assay does not guarantee mobility in the host, if Muta1 was able to mobilize even a subset of these elements as it does in yeast, it could be an effective tool for high frequency insertional mutagenesis, especially when coupled with its preference for genic insertions and a lack of target sequence preference (Additional file 3: Figure S2B, Additional file 8: Figure S6).
This is the first report of the transposition of a non-plant MULE, Muta1, in a heterologous system and provides the first experimental evidence for the functional significance of the DDE domain in the transposition reaction in the MULE superfamily. High frequency transposition in a yeast assay facilitated the determination of Muta1 transposition features including precise excision, genic targeting with no sequence preference and the impact of TIR and TSD for insertion and excision. Taken together, Muta1 may be a valuable tool for forward genetics in mosquitoes.
Identification of MULEs in Ae. aegypti
The mosquito genomes used in this study: AaegL3 build for Ae. aegypti (https://www.vectorbase.org/organisms/aedes-aegypti/liverpool/aaegl3); AaloF1 build for Ae. albopictus (https://www.vectorbase.org/organisms/aedes-albopictus/foshan/aalof1); CpipJ2 build for C. quinquefasciatus (https://www.vectorbase.org/organisms/culex-quinquefasciatus/johannesburg/cpipj2) and AgamP4 build for An. gambiae (https://www.vectorbase.org/organisms/anopheles-gambiae/pest/agamp4). The conserved MULE DDE domain from all eukaryotes  was used as query to search the mosquito genomes by TBLASTN, as implemented in the TARGeT pipeline  with an E-value cutoff of 0.001. Flanking DNA sequences with 10 kb upstream and downstream of the matched region were retrieved. The ends of a putative element were determined by aligning two closely related elements with their 20 kb flanking sequences, TIR boundaries and TSDs were manually identified. Coding capacity of each element was predicted by the GENSCAN program (http://genes.mit.edu/GENSCANinfo.html).
To identify Muta1 derivative nonautonomous elements, 50 bp from each end of Muta1 was used in a BLASTN search with TARGeT  using default parameters. One hundred bp of flanking DNA sequences were retrieved for manual verification of the TIR and TSD of each derivative element. Fifty bp flanking each element were used for BLASTN searches against the Ae. aegypti genome (AaegL3 build) to determine the genomic location and compared to the genome annotation (release AaegL3.3, https://www.vectorbase.org/organisms/aedes-aegypti/liverpool/aaegl3) to determine the adjacent genes. For the control data set, 171 genome coordinates across the 4,757 scaffolds were randomly generated, and compared to the genome annotation (release AaegL3.3) to determine the surrounding sequences and genes. The random insertion sites generation used 1,000 replicates to estimate the expected number of insertions (and standard deviations) in each category.
Yeast construct construction
Genomic DNA of individual Ae. aegypti mosquito (Liverpool strain, obtained from Dr. Atkinson, UC Riverside) was extracted using the DNeasy Blood & Tissue Kit (Qiagen). The two exons of Muta1 predicted by GENSCAN program were cloned from genomic DNA and fused through overlap PCR (all primer sequences are available in Additional file 11: Table S3). The complete transposase coding sequence was then cloned into the Gateway cloning vector pENTR and transferred to destination vector pGAL415-ccdb  with LR Clonase (Invitrogen) to generate the pGAL415-ccdbMuta1 plasmid (Additional file 3: Figure S2A).
Two hundred fifty bp from each end of Muta1 were fused by overlap PCR to generate Muta1AR. The HIS3 fragment containing the yeast HIS3 coding sequence, HIS3 5’ and 3’ UTR, and HIS3 promoter was cloned from vector pGAL415-ccdb . The HIS3 fragment was then fused with 250 bp from each end of Muta1 through overlap PCR to generate Muta1HIS. Muta1NA1-5 elements were cloned directly from genomic DNA. All nonautonomous elements were inserted in the HpaI site of ADE2 for the exon excision assay or the XhoI site for 5’ UTR excision assay through homologous recombination in yeast as previously described . Donor site TSDs were introduced by adding corresponding TSD sequences in primers (Additional file 11: Table S3).
For Muta5 assay, plasmid pGAL415-ccdbMuta5 was constructed in the same way as pGAL415-ccdbMuta1, and Muta5NA was constructed by overlap PCR (Additional file 11: Table S3).
The yeast transposition assay using Saccharomyces cerevisiae strain DG2523 and the pWL89a vector was described previously [35, 36]. Transformation was performed using the Frozen-EZ Yeast Transformation kit (Zymo research). For excision assays, transformants were grown in 5 ml liquid media of CSM -leu-ura with 2% dextrose. After growth to saturation (36 h), cells were washed twice with 5 ml water, resuspended in 0.5 ml water and plated onto CSM -his-leu-ade with 2% galactose. Colonies were counted after incubation at 30°C for 15 days and viable counts were made by plating 100 μl of a 1 × 105 and 1 × 106 dilution on YPD plates.
For the reintegration assay, cells were grown to saturation in 5ml liquid CSM -leu-ura with 2% dextrose, cells were washed twice with 5 ml sterile water, resuspended in 0.5 ml water and plated onto CSM -leu-ura-ade with 2% galactose plate and CSM -his-leu + 5-FOA with 2% galactose plates. Colonies were counted after incubation at 30°C for 15 days, and viable counts were made by plating 100 μl of a 1 × 105 and 1 × 106 dilution on YPD plates. In another approach, individual Ade+ Muta1HIS excision revertant colonies isolated directly from plates of CSM -his-leu-ade with 2% galactose were streaked on CSM -his + 5-FOA plates to calculate the reintegration frequency.
Excision and reinsertion analysis
For footprint analysis, colony PCR was performed on ADE2 revertant colonies using primers (Table S3) flanking the insertion sites. PCR products were gel extracted (Zymoclean Gel DNA Recovery Kit) and sequenced. For reinsertion analysis, transposon display was conducted . Genomic DNA was extracted from revertant colonies using the Yeastar genomic DNA kit (Zymo research); DNA samples were digested by BfaI followed by adapter ligation. Pre-amplification and selective amplification were used to amplify the sequences between Muta1 TIR and the BfaI adapter sequence. Amplicons consisting of flanking sequences of the reinsertion sites and part of Muta1 TIR were resolved on a 4% agarose gel, and polymorphic fragments were recovered and sequenced. Flanking sequences were mapped to the yeast genome (S288C, http://yeastgenome.org/) and the reinsertion sites were determined with regard to the closest genomic features. The insertion site analysis figure was made using the program Pictogram http://genes.mit.edu/pictogram.html.
Mutagenesis of Muta1 transposase
Site-directed mutagenesis was used to generate mutant versions of Muta1 transposase. One pair of primers (Additional file 11: Table S3) was used for each mutation site, and pGAL415-ccdbMuta1 plasmid was used as template. PCR products were digested with Dpn1 to remove template, and the resulting plasmid was sequenced to confirm that mutations occurred as expected.
Mutator-like transposable elements
Polymerase chain reaction
Terminal inverted repeat
Target site duplication
We thank Jim Burnette, Lu Lu, Brad Cavinder, Peter Atkinson and Jason Stajich for their technical assistance and advice.
This research funded by a grant from the W.M. Keck Foundation to SRW.
Availability of data and materials
The materials and datasets in the current study are available from the corresponding author on request.
KL performed the sequence analysis, designed and carried out the yeast vector construction and transposition assays. SRW conceived the study, participated in its design and coordination and helped to draft the manuscript. Both authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Consent for publication
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Finnegan DJ. Eukaryotic transposable elements and genome evolution. Trends Genet. 1989;5:103–7.View ArticlePubMedGoogle Scholar
- Wicker T, Sabo F, Hua-Van A, Bennetzen JL, Capy P, Chalhoub B, Flavell A, Leroy P, Morgante M, Panaid O, Paux E, SanMiguel P, Schulman AH. A unified classification system for eukaryotic transposable elements. Nat Rev Genet. 2007;8:973–82.View ArticlePubMedGoogle Scholar
- Feschotte C, Pritham EJ. DNA transposons and the evolution of eukaryotic genomes. Annu Rev Genet. 2007;41:331–68.View ArticlePubMedPubMed CentralGoogle Scholar
- Bao W, Jurka MG, Kapitonov VV, Jurka J. New Superfamilies of Eukaryotic DNA Transposons and Their Internal Divisions. Mol Biol Evol. 2009;26(5):983–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Robertson DS. Characterization of a mutator system in maize. Mutat Res. 1978;51:21–8.View ArticleGoogle Scholar
- Lisch D. Regulation of the Mutator System of Transposons in Maize. Methods Mol Biol. 2013;1057:123–42.View ArticlePubMedGoogle Scholar
- Neuveglise C, Chalvet F, Wincker P, Gailardin C, Casaregola S. Mutator-like element in the yeast Yarrowia lipolytica displays multiple alternative splicings. Eukaryot Cell. 2005;4:615–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Pritham EJ, Feschotte C, Wessler SR. Un-expected diversity and differential success of DNA transposons in four species of entamoeba protozoans. Mol Biol Evoanl. 2005;22:1751–63.View ArticleGoogle Scholar
- Lopes FR, Silva JC, Benchimol M, Costa GG, Pereira GA, Carareto CM. The protist Trichomonas vaginalis harbors multiple lineages of transcriptionally active Mutator-like elements. BMC Genomics. 2009;10:330.View ArticlePubMedPubMed CentralGoogle Scholar
- Lisch D. Mutator transposons. Trends Plant Sci. 2002;7:498–504.View ArticlePubMedGoogle Scholar
- Jiang N, Bao Z, Zhang X, Eddy SR, Wessler SR. Pack-MULE transposable elements mediate gene evolution in plants. Nature. 2004;431:569–73.View ArticlePubMedGoogle Scholar
- Chalvet F, Grimaldi C, Kaper F, Langin T, Daboussi MJ. Hop, an active Mutator-like element in the genome of the fungus Fusarium oxysporum. Mol Biol Evol. 2003;20:1362–75.View ArticlePubMedGoogle Scholar
- Xu Z, Yan X, Maurais S, Fu H, O’Brien DG, Mottinger J, Dooner HK. Jittery, a Mutator distant relative with a paradoxical mobile behavior: Excision without reinsertion. Plant Cell. 2004;16:1105–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Li Y, Harris L, Dooner HK. TED, an autonomous and rare maize transposon of the Mutator superfamily with a high gametophytic excision frequency. Plant Cell. 2013;25:3251–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao D, Ferguson A, Jiang N. Transposition of a rice Mutator-like element in the yeast Saccharomyces cerevisiae. Plant Cell. 2015;27:132–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Babu MM, Iyer LM, Balaji S, Aravind L. The natural history of the WRKY-GCM1 zinc fingers and the relationship between transcription factors and transposons. Nucleic Acids Res. 2006;34(22):6505–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Nesmelova IV, Hackett PB. DDE transposases: Structural similarity and diversity. Adv Drug Deliv Rev. 2010;62(12):1187–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Eisen JA, Benito MI, Walbot V. Sequence similarity of putative transposases links the maize Mutator autonomous element and a group of bacterial insertion sequences. Nucleic Acids Res. 1994;22:2634–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Loessner I, Dietrich K, Dittrich D, Hacker J, Ziebuhr W. Transposase-dependent formation of circular IS256 derivatives in Staphylococcus epidermidis and Staphylococcus aureus. J Bacteriol. 2002;184:4709–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Yuan YW, Wessler SR. The catalytic domain of all eukaryotic cut-and-paste transposase superfamilies. Proc Natl Acad Sci U S A. 2011;108:7884–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Hickman AB, Ewis HE, Li X, Knapp JA, Laver T, Doss AL, Tolun G, Steven AC, Grishaev A, Bax A, Atkinson PW, Craig NL, Dyda F. Structural basis of hAT transposon end recognition by Hermes, an octameric DNA transposase from Musca domestica. Cell. 2014;158:353–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Calisher CH. Persistent Emergence of Dengue. Emerg Infect Dis. 2005;11:738–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Marchette NJ, Garcia R, Rudnick A. Isolation of Zika virus from Aedes aegypti mosquitoes in Malaysia. Am J Trop Med Hyg. 1969;18:411–5.PubMedGoogle Scholar
- Aviles G, Sabattini MS, Mitchell CJ. Peroral susceptibility of Aedes albifasciatus and Culex pipiens complex mosquitoes (Diptera: Culicidae) from Argentina to western equine encephalitis virus. Rev Saude Publica. 1990;24(4):265–9.View ArticlePubMedGoogle Scholar
- Lindsay SW, Wilkins HA, Zieler HA, Daly RJ, Petrarca V, Byass P. Ability of Anopheles gambiae mosquitoes to transmit malaria during the dry and wet seasons in an area of irrigated rice cultivation in The Gambia. J Trop Med Hyg. 1991;94(5):313–24.PubMedGoogle Scholar
- Fraser MJ. Insect Transgenesis: Current Applications and Future Prospects. Annu Rev Entomol. 2012;57:267–89.View ArticlePubMedGoogle Scholar
- O'Brochta DA, Sethuraman N, Wilson R, Hice RH, Pinkerton AC, Levesque CS, Bideshi DK, Jasinskiene N, Coates CJ, James AA, Lehane MJ, Atkinson PW. Gene vector and transposable element behavior in mosquitoes. J Exp Biol. 2003;206:3823–34.View ArticlePubMedGoogle Scholar
- Scali C, Nolan T, Sharakhov I, Sharakhova M, Crisanti A, Catteruccia F. Post-integration behavior of a Minos transposon in the malaria mosquito Anopheles stephensi. Mol Genet Genomics. 2007;278:575–84.View ArticlePubMedGoogle Scholar
- Nene V, et al. Genome sequence of Aedes aegypti, a major arbovirus vector. Science. 2007;316:1718–23.View ArticlePubMedGoogle Scholar
- Chen XG, et al. Genome sequence of the Asian Tiger mosquito, Aedes albopictus, reveals insights into its biology, genetics, and evolution. Proc Natl Acad Sci U S A. 2015;112:5907–15.View ArticleGoogle Scholar
- Arensburger P, Megy K, Waterhouse RM, Abrudan J, Amedeo P, Antelo B, Bartholomay L, Bidwell S, Caler E, Camara F. Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics. Science. 2010;330:86–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Holt RA, et al. The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002;298:129–49.View ArticlePubMedGoogle Scholar
- Han YJ, Burnette JM, Wessler SR. TARGeT: A web-based pipeline for retrieving and characterizing gene and transposable element families from genomic sequences. Nucleic Acids Res. 2009;37, e78.View ArticlePubMedPubMed CentralGoogle Scholar
- Bingham PE, Zachar Z. Retrotransposons and the FB transposon from Drosophila melanogaster. Mobile DNA Washington D.C: American Society for Microbiology Press. 1989;485–502.
- Yang G, Weil CF, Wessler SR. A rice Tc1/mariner-like element transposes in yeast. Plant Cell. 2006;18:2469–78.View ArticlePubMedPubMed CentralGoogle Scholar
- Weil CF, Kunze R. Transposition of maize Ac/Ds transposable elements in the yeast Saccharomyces cerevisiae. Nat Genet. 2000;26:187–90.View ArticlePubMedGoogle Scholar
- Hancock CN, Zhang F, Wessler SR. Transposition of the Tourist-MITE mPing in yeast: an assay that retains key features of catalysis by the class 2 PIF/Harbinger superfamily. Mob DNA. 2010;1:5.View ArticlePubMedPubMed CentralGoogle Scholar
- Mitra R, Fain-Thornton J, Craig NL. piggyBac can bypass DNA synthesis during cut and paste transposition. Embo J. 2008;27:1097–109.View ArticlePubMedPubMed CentralGoogle Scholar
- Biedler J, Qi Y, Holligan D, Della Torre A, Wessler SR, Tu Z. Transposable element (TE) display and rapid detection of TE insertion polymorphism in the Anopheles gambiae species complex. Insect Mol Biol. 2003;12(3):211–6.View ArticlePubMedGoogle Scholar
- Lynch M, Sung W, Morris K, Coffey N, Landry CR, Dopman EB, Dickinson WJ, Okamoto K, Kulkarni S, Hartl DL, Thomas WK. A genome-wide view of the spectrum of spontaneous mutations in yeast. Proc Natl Acad Sci U S A. 2008;8105(27):9272–7.View ArticleGoogle Scholar
- Sutton WD, Gerlach WL, Schwartz D, Peacock WJ. Molecular analysis of Ds controlling element mutations at the Adhl locus of maize. Science. 1983;223:1265–8.View ArticleGoogle Scholar
- Brillet B, Bigot Y, Auge-Gouillou C. Assembly of the Tc1 and mariner transposition initiation complexes depends on the origins of their transposase DNA binding domains. Genetica. 2007;130(2):105–20.View ArticlePubMedGoogle Scholar
- Zhou L, Mitra R, Atkinson PW, Hickman AB, Dyda F, Craig NL. Transposition of hAT elements links transposable elements and V(D)J recombination. Nature. 2004;432:995–1001.View ArticlePubMedGoogle Scholar
- Hennig S, Ziebuhr W. A transposase-independent mechanism gives rise to precise excision of IS256 from insertion sites in Staphylococcus epidermidis. J Bacteriol. 2008;190:1488–90.View ArticlePubMedGoogle Scholar
- Alberti S, Gitler AD, Lindquist S. A suite of gateway (R) cloning vectors for high-throughput genetic analysis in Saccharomyces cerevisiae. Yeast. 2007;24:913–9.View ArticlePubMedPubMed CentralGoogle Scholar