Mojo Hand, a TALEN design tool for genome editing applications

Background Recent studies of transcription activator-like (TAL) effector domains fused to nucleases (TALENs) demonstrate enormous potential for genome editing. Effective design of TALENs requires a combination of selecting appropriate genetic features, finding pairs of binding sites based on a consensus sequence, and, in some cases, identifying endogenous restriction sites for downstream molecular genetic applications. Results We present the web-based program Mojo Hand for designing TAL and TALEN constructs for genome editing applications (http://www.talendesign.org). We describe the algorithm and its implementation. The features of Mojo Hand include (1) automatic download of genomic data from the National Center for Biotechnology Information, (2) analysis of any DNA sequence to reveal pairs of binding sites based on a user-defined template, (3) selection of restriction-enzyme recognition sites in the spacer between the TAL monomer binding sites including options for the selection of restriction enzyme suppliers, and (4) output files designed for subsequent TALEN construction using the Golden Gate assembly method. Conclusions Mojo Hand enables the rapid identification of TAL binding sites for use in TALEN design. The assembly of TALEN constructs, is also simplified by using the TAL-site prediction program in conjunction with a spreadsheet management aid of reagent concentrations and TALEN formulation. Mojo Hand enables scientists to more rapidly deploy TALENs for genome editing applications.

Background TAL domains exhibit programmable, sequence-specific binding to DNA, a feature that makes them a valuable addition to the tools of the molecular biologist. In particular, TAL domains may be used in combination with endonucleases to cause double-strand breaks which are exploited for genome editing, either by error-prone non-homologous end-joining repair of double-strand breaks or insertion of new sequence by homologous recombination. These exciting possibilities depend on the ability of a molecular biologist to design TAL binding sequences for specific genomic regions.
Sequence-specific DNA binding by TAL effectors is accomplished by individual sub-domains of 33-35 amino acids. These repeat variable di-residue (RVD) [1,2] domains contain a central pair of amino acids that determine the base to which it binds. A variety of these RVDs are found in nature, but artificial TAL effectors typically include: adenine=NI, cytosine=HD, guanine=NN, and thymine=NG. For example, these RVDs are used in the one-pot Golden Gate [3], FLASH [4], unit assembly [5], or iterative capped assembly [6] reactions to construct sequence-specific DNA binding proteins.
The TAL domain can bind nearly any DNA sequence. Early work on TAL effectors indicated a consensus sequence where a thymine must precede the binding site, followed by [ACG] and [CGT] [1][2][3][4]6]. These requirements are remarkably nonrestrictive, which makes the TAL proteins useful for targeting most genes and regulatory elements of sufficient length. Recent work, however, indicates that only the first of these consensus sequence rules appears to be a measurable constraint on TALEN design [7] (ALM, unpublished results), though even that may not be absolute when using appropriately designed N and C termini [6][7][8].
In the context of genome editing, software for designing binding sites for TALENs should be flexible and able to target both exons and introns. Also, detecting TALEN activity is often facilitated by restriction fragment length polymorphism (RFLP) assays (for example, [9]). Access to design software in a context suitable for use by molecular biologists and life scientists is also essential for TALEN use by the field.

Results and discussions
The Mojo Hand algorithm consists of multiple steps: user definition of sequence source, retrieval of non-FASTA sequence, identification of candidate TALEN sites, selection of restriction enzyme sites for TALEN efficiency analysis, and output to a text or .csv file. The output can be used directly or pasted into a supplied spreadsheet to ease in assembly of the designed TALENs. The process is depicted in Figure 1.

Sequence input
Entering the sequence to be targeted by the TAL effector may be accomplished through automatic download of the gene through the NCBI Gene or Nucleotide databases or FASTA-formatted text file. Users specify the unique identifier for the gene of interest. Then the exons and introns of which are retrieved using E-Utilities [10], the API for Pub-Med and other NCBI databases.
To process subsequences (exon or intron) in genes downloaded from NCBI Gene, Mojo Hand requires at least one mRNA, CDS, misc RNA, or exon feature. No subsequence information is used when users enter an arbitrary sequence through FASTA files or Nucleotide requests. The Gene record must also include genomic location information so that the correct genomic sequence may be requested from the NCBI Nucleotide database. In a random sample of genes from diverse organisms, we found that 80% of genes were precisely located. Of these, we noted that 63% had at least one mRNA record, 89% had at least one CDS record, and 4% had at least one miscellaneous RNA record (Table 1). Only 5% of the genes lacked all subsequence features (below). In some cases, features for several genes (other reading frames, especially on the opposite strand) are annotated but Mojo Hand can determine which feature to use based on the gene symbol and its aliases. However, when multiple mRNA features are included for the same gene, there is no obvious way to determine which is most appropriate, so Mojo Hand selects the first feature found unless the user specifies a preference by index (or by transcript identifier). Similarly, because CDS records lack transcript or other unique tags, the user may only specify an index preference. Exon features are the smallest defined feature used as a group. Therefore, no further subsequence selection can be made.
Many genes are not fully annotated and may not have annotated mRNA features. Mojo Hand prefers mRNA features, but if none are found, CDS and miscellaneous RNA records may be used. If none of these subsequence features exist, the entire sequence is used. The mRNA record is given the highest priority because it gives the user flexibility in selecting any exon. However, users may also set the priority to CDS records (−−cds-index=1) that do not contain promoter and other regulatory elements that may complicate results.  Figure 1 Flowchart of Mojo Hand. Inputs are the Entrez Gene identifier, a nucleotide fragment, or a text file containing the sequence in FASTA format. The output is a report containing possible binding sites, suitable restriction sites, and RVDs. This output can be further processed with the included spreadsheet, which produces a customized recipe for each TALEN. Frequency of feature count for 2052 randomly selected genes. Gene annotations often do not include mRNA records, but most (95%) have at least one CDS record. Miscellaneous RNA features are relatively rare. Genes may also contain multiple subsequence features that include variations in exon locations, identical reports from multiple authors, or genes in other reading frames. This table only includes genes for the~80% of genes for which precise genomic location.
Users may select exons or introns for analysis or use the entire sequence. Depending on the input, some flanking sequence may be included for further analysis. The flanking sequence is desirable because some exons are short compared to the length of the TAL binding site. Flanking sequence may also be needed for selecting a restriction enzyme that cuts only once in an amplicon of specific length, especially for binding sites that occur near the ends of an exon or for very short exons as illustrated in Figure 2.

NCBI access
Sequence data was downloaded using ESummary and EFetch -URL-based methods of requesting information from the various NCBI database. The gene of interest (e.g. gene ID 567858) is requested by the unique identifier in the Gene database first. Mojo Hand then constructs a request for a detailed record based on the genomic location as defined by the RefSeq accession number and indices of the beginning and end of the sequence. The source of these 3 parameters can be automatically determined from the XML output of an ESummary request from a Gene database or manual entry at the command line. For example, Gene 567858 can be requested by Gene Identifier or its location NC 007118.5, 21501610-21527471. This identifier and range define the gene of interest and may be used if the unique identifier in Gene does not produce the desired result. Mojo Hand may also be used to find binding sites in sequences that are not designated as a gene by NCBI. Mojo Hand then requests the detailed record from Nucleotide in XML format using EFetch, captures the beginning and end points of the exons from the mRNA section, and stores the genomic sequence for later analysis. The forward or reverse strand is requested based on the beginning and end indices mentioned above. Because any number of mRNA features may be available, Mojo Hand parses the XML to find those features that are designated with the same symbol as the gene of interest. This procedure distinguishes between the gene of interest and other genes encoded in another reading frame. Some records (e.g. gene ID 32619) have many mRNA entries for the gene of interest. Since there is no automatic way to determine which entry is the most appropriate, the first entry is used. In cases of manual entry of gene location, no symbol is available and the first mRNA is used, and a warning issued. Troublesome cases may be handled by downloading the gene, extracting the exons manually and using the file input mode.
The length of the requested sequence is modified to include a user-defined number of bases upstream and downstream of the gene. If the gene is very near the beginning of the contig, some of the sequence may be undefined and is filled with placeholders. Likewise, if the requested sequence falls near the end of the contig such that the trailing flanking sequence extends beyond its end (e.g. gene ID 802118), placeholders are added.
Records containing long sequences (e.g. gene ID 19091, 1.2 million base pair (bp)) are processed somewhat differently. The XML output is not processed beyond the end of the GSeq feature table field.

Identification of binding sites
We designed Mojo Hand to identify TAL binding sites based on a user-defined template. We initially used the template sequence of Ts[ACG] [CGT].*Te. The notation s and e indicate the start and end of the binding sequence, bases enclosed in brackets represent a choice, and .* indicate zero or more bases for which there is no preference. So this template will identity TAL binding sites that start with a T, are followed by A, C, or G (not a T), and are followed by a C, G or T (not an A). Another template, based on the work of Sun and coworkers [7], is Ts.*e. In practice, we have found no substantial functional constraints besides the initial 5' T bp, so we now use this as the default template sequence parameter for Mojo Hand.
In addition to the constraints described above, the user specifies a range of allowed lengths for the binding sites and the spacer. The algorithm generates the entire set of possible binding sites based on these lengths. Each candidate binding site is then filtered based on TAL-site restrictions. The candidate sites are made by iterating the possible lengths (Figure 2a When the length of an exon is long, execution time may be extensive and the number of TALEN sites becomes unusably large. Therefore, for long exons, some binding sites are skipped by incrementing the beginning of the first TAL site by values larger than unity. For exons longer than 1000 bp, the increment is 10; for exons longer than 5000, the increment is 20. Using this skipping method and the classical consensus sequence, we observed that many binding sites were found in our test set (35 genes, listed below). Every possible binding site may be obtained by submitting fragments of length 999 base pairs or less as a FASTA file or NCBI Nucleotide request.

Restriction enzyme analysis
The biological effect of TALEN activity may be observed in several ways, including restriction fragment length polymorphism (RFLP), sequencing, and phenotype. If RFLP will be used as an evaluation approach, a restriction site should be present within the space between the TAL binding sites. In conjunction with double-strand break repair by error-prone non-homologous end-joining (NHEJ), the nearby restriction site is often disrupted. We subjected candidate binding site pairs and spacers to further analysis to find those candidates with unique restriction-enzyme binding sites within their spacer region. If the user requests the single TAL binding site option, no enzyme analysis is performed.
We used REBASE [11], a database of restriction enzymes hosted by New England Biolabs (NEB), as the default restriction enzyme database collection within Mojo Hand. We use only the subset of enzymes that are commercially available. Future releases of the REBASE database may be downloaded and used in the place of the version distributed with this manuscript. Custom inhouse databases may also be used if the format matches that of REBASE. We also used results published by NEB [12] to determine which enzymes were compatible with several PCR buffers (standard, Thermopol, Phusion, and Crimson). These scores are rescaled on the range [0,9] and displayed with each enzyme. Formal and prototype names (first described enzyme of a particular family) are displayed for ease of use. Mojo Hand permits selection of enzyme site screening by vendor(s). Narrowing the selection decreases computational load thereby decreasing processing time for TALEN site selection using Mojo Hand.

Mojo Hand output
The output of Mojo Hand is either a flat text file or HTML that includes the sequences to which the TAL proteins will eventually bind, RVDs, and enzymes that cut in each spacer (Figure 3a) or a downloadable comma-separated values (.csv) file that can be imported into common programs like Microsoft Excel (Figure 3b). The binding sites are denoted TAL1 and TAL2, as illustrated in Figure 2b. The first site and spacer are considered the forward direction, regardless of the strand in the NCBI Nucleotide database. The reverse complement of the second binding site is also provided to show the nucleotide sequence to which the TALEN protein will bind.
The gene symbol is indicated at the beginning of each TALEN site. The notation symbol E# indicates which exon each binding site was found within. The prefix E indicates that the region of interest was an exon; I for intron, and A when no subsequence region was available and the entire sequence was used.
The position and length of each binding site are provided. These values are relative to the beginning of the subsequence fragment (exon or intron) and include the user-defined short flanking sequence. The notation #-#(#) indicates the first and last index of the binding site, with the length in parentheses.
The restriction sites for RFLP are listed for each TALEN site. The enzyme name is listed with its prototype. The compatibility of the enzymes with full-strength PCR buffer [12] is listed, rescaled so that it falls on the range [0, 9], with 0 representing no activity. The position of the first base of the enzyme's recognition sequence is given for nucleases that cut only once in the amplicon (binding site with default of 150 flanking bp on either end). The minimum distance to the second cut site is configurable (default = 80 bp). If an enzyme cuts in two positions that are at least this minimum distance apart, the second cut is given relative to the first restriction site. The indexing is based on the beginning of the subsequence fragment, including the long flanking length. If using the Mojo Hand web service, enzyme restriction site matches are highlighted in black in the spacer while reverse complement restriction site matches are highlighted in red.

TALEN recipe generator
Commonly used TALEN assembly protocols involve a large number of plasmids [3]. Mojo Hand includes a spreadsheet that aids in TALEN formulation that may be used in conjunction with the Golden Gate method and the TALEN DNA kit from Addgene (TALEN Kit #1000000016; Cambridge, MA, USA) [3,9,13]. Mojo Hand outputs tabdelimited or CSV downloadable RVD sequences, which can be transferred into the spreadsheet using the clipboard. The spreadsheet then produces recipes that facilitate the molecular assembly of TALENs.

Verification of Mojo Hand and measurement of TALEN activity
To prospectively test the Mojo Hand software, we designed TALENs targeting zebrafish flt3 sequence using Mojo Hand. A 450 bp genomic fragment including Exon 1 was used as input. An Exon 1 targeting TALEN pair with 15-RVD binding sites including a 5' upstream T nucleotide at both arms was chosen, which has a 16 bp spacer and a unique NsiI site in the spacer (Figures 3  and 4). The TALEN constructs were subsequently synthesized with the Golden Gate Method [3] using the GoldyTALEN scaffold [9,13]. In vitro transcribed mRNAs (mMessage mMachine T3 kit, Life Technology, Carlsbad, CA, US) encoding the TALEN pair were injected into the cytoplasm of 1-cell stage wild-type zebrafish. At 48 hours-post-fertilization (hpf ), genomic DNA was extracted from single and groups of 10 phenotypically normal embryos. The in vivo activity of the TALEN was screened by RFLP assay as previously described [9]. In embryos injected with the TALEN mRNA, uncut PCR product was detected representing the loss of NsiI site through TALEN-induced small deletion, with the percentage of converted chromosomes measured to be 84 ± 5% (average of three different experiments). These results demonstrated a successful TALEN design with Mojo Hand. We have conducted a separate prospective analysis of 12 more TALENs made using Mojo Hand, with each pair garnering over 20% activity and the set averaging over 50% converted chromosomes in this zebrafish embryo assay (data not shown). The record includes the binding sequences in proper orientation, the RVDs (tab delimited in normal output for easy transfer to a spreadsheet) and the restriction enzymes that could be used. Note that this example has only a single cutting site so the absolute position of the beginning of the restriction site is given but the column for relative position of the second cut is empty. b) Screen-shot of downloaded spreadsheet in .CSV format of the same example.

Source code and Web service
Mojo Hand is available as a web service at www.talendesign. org. The site allows access to the program without the trouble of installation and with the ease of a familiar interface. Point-of-use help is available for each field. The source code and spreadsheet are also available for non-commercial use with applicable license.

Conclusions
We developed Mojo Hand based on an initial training set of 35 genes, further tested the program with users from 3 separate laboratories, and finally conducted a prospective study of over a dozen TALEN pairs. During the initial vetting process, we showed that the correct sequences were downloaded from NCBI databases using NCBI nr/ntBLAST [14]. We manually confirmed that the Mojo Hand predicted binding sites were within the expected exon of the correct gene in all cases. We also compared Mojo Hand output to manually retrieved GenBank records and verified several binding sites to ensure that they occurred in the correct exon and that the prefix requirements (a 5' thymine, in most cases) were satisfied. Our test set was constructed so that different numbers (0-3) of subsequence features were present, allowing us to assess how Mojo Hand prioritizes mRNA, CDS, and misc RNA features. We included genes with multiple aliases and features labeled with an alias to test if all appropriate subsequence records were found. Mojo Hand complements previously described software such as the ISU TALEN Targeter [3,15] because Mojo Hand addresses the difficulty of downloading the sequence and extracting exons (or introns) based on annotation from GenBank. The ISU TALEN Targeter currently only accepts sequences entered in a text box or FASTA file upload. Mojo Hand can also screen for possible TALEN sites using more extensive databases of restriction enzymes (REBASE) rather than the NEB database that was recently added to the ISU tool program. Mojo Hand also provides a spreadsheet that bridges the gap between the RVD output and the bench. The spreadsheet produces recipes for individual TALENs that take into account local reagent concentrations.
We also compared our work to idTALE, a web service provided by King Abdullah University of Science and Technology [16]. idTALE allows users to provide sequence directly or by Ensembl gene identifier. Genes are, however, restricted to just a few species (A. thaliana, P. patens, C. elegans, D. melanogaster, S. cerevisiae), and no restriction analysis is done. Mojo Hand appears more flexible because any gene or sequence can be entered, any consensus sequence may be used, and restriction analysis is available.
Beyond the single RVD to nucleotide recognition cipher of TALENs, other interactions that affect TALEN activity appear to be minor. However, factors that potentially affect TALEN efficiency continue to be investigated. New information regarding TALEN design may require rapid change of current TALEN design software. Therefore, Mojo Hand has been designed to permit user-defined adjustments. Beyond TALEN design based on binding alone, Mojo Hand provides an integrated a way to download the exons or introns of any gene and to filter the results based on restriction enzyme recognition sites.
We have designed Mojo Hand to be as general as possible, but there are several limitations. The annotated features in GenBank can vary based on what is known about a particular gene, so Mojo Hand may not be able to download certain predicted genes. We found that a significant proportion of randomly selected genes in NCBI have no subsequence records (5% , Table 1), and many (20%) have no genomic location. Also, we limit our automatic download features to work only on genes in the GenBank and NCBI Nucleotide databases, which excludes much of the non-coding, non-repetitive regions of the genome. Users may manually enter sequence data to overcome this limitation.