- Research article
- Open Access
A high-throughput screen identifying sequence and promiscuity characteristics of the loxP spacer region in Cre-mediated recombination
© Missirlis et al; licensee BioMed Central Ltd. 2006
- Received: 21 December 2005
- Accepted: 04 April 2006
- Published: 04 April 2006
Cre-loxP recombination refers to the process of site-specific recombination mediated by two loxP sequences and the Cre recombinase protein. Transgenic experiments exploit integrative recombination, where a donor plasmid carrying a loxP site and DNA of interest integrate into a recipient loxP site in a target genome. Unfortunately, integrative recombination is highly inefficient because the insert is flanked by two loxP sites, which themselves become targets for Cre and lead to subsequent excision of the insert. A small number of mutations have been discovered in parts of the loxP sequence, specifically the spacer and inverted repeat segments, that increase the efficiency of integrative recombination. In this study we introduce a high-throughput in vitro assay to rapidly detect novel loxP spacer mutants and describe the sequence characteristics of successful recombinants.
We created synthetic loxP oligonucleotides that contained a combination of inverted repeat mutations (the lox 66 and lox 71 mutations) and mutant spacer sequences, degenerate at 6 of the 8 positions. After in vitro Cre recombination, 3,124 recombinant clones were identified by sequencing. Included in this set were 31 unique, novel, self-recombining sequences. Using network visualization tools, we recognized 12 spacer sets with restricted promiscuity. We observed that increased guanine content at all spacer positions save for position 8 resulted in increased recombination. Interestingly, recombination between identical spacers was not preferred over non-identical spacers. We also identified a set of 16 pairs of loxP spacers that reacted at least twice with another spacer, but not themselves. Further, neither the wild-type P1 phage loxP sequence nor any of the known loxP spacer mutants appeared to be kinetically favoured by Cre recombinase.
This study approached loxP spacer mutant screening in an unbiased manner, assuming nothing about candidate loxP sites save for the conserved 4 and 5 spacer positions. Candidate sites were free to recombine with any other sequence in the pool of all possible sites. The subset of loxP sites identified here are candidates for in vivo serial recombination as they have already demonstrated limited promiscuity with other loxP spacer and stability in the presence of Cre.
- Space Sequence
- Inverted Repeat
- loxP Site
- Recombination Reaction
- loxP Sequence
Published mutant loxP sites
Left inverted repeat sequence
Right inverted repeat sequence
lox 511 
lox 5171 
lox 2272 
lox 71 [12, 13]
lox 66 [12, 13]
Two classes of variant loxP sites are available to engineer stable Cre-loxP integrative recombination. Both exploit sequence mutations in the Cre recognition sequence, either within the 8 bp spacer region or the 13-bp inverted repeats. Spacer mutants such as lox 511, lox 5171 , lox 2272 , m2, m3, m7, and m11  (Table 1) recombine readily with themselves but have a markedly reduced rate of recombination with the wild-type site. This class of mutants has been exploited for DNA insertion by Recombinase Mediated Cassette Exchange (RMCE) [9–11]. Inverted repeat mutants represent the second class available and contain altered bases in the left inverted repeat (LE mutant) or the right inverted repeat (RE mutant). The LE mutant, lox 71, has 5 bp on the 5' end of the left inverted repeat that are changed from the wild type sequence to TACCG [12, 13] (Table 1). Similarly, the RE mutant, lox 66, has the five 3'-most bases changed to CGGTA (Table 1). Inverted repeat mutants are used for integrating plasmid inserts into chromosomal DNA with the LE mutant designated as the "target" chromosomal loxP site into which the "donor" RE mutant recombines. Post-recombination, loxP sites are located in cis, flanking the inserted segment. The mechanism of recombination is such that post-recombination one loxP site is a double mutant (containing both the LE and RE inverted repeat mutations) and the other is wild type [14–16]. The double mutant is sufficiently different from the wild-type site that it is unrecognized by Cre recombinase and the inserted segment is not excised. Recently, spacer and inverted repeat mutants have been combined to increase the specificity and stability of integrative recombination [8, 17].
In past experiments, novel spacer mutants were discovered by mutating a suite of bases [6, 7] or generating a set of potential spacer mutants and testing recombination between these spacers with the wild-type loxP site . Lee and Saito used an in vitro assay that evaluated the recombination efficiency of 24 spacers with 1 bp substitutions and 30 spacers with 2 bp substitutions from the sequence of the wild-type loxP . Their data suggested that homology was required at positions 2–5 and positions 6–7 for efficient strand exchange and resolution of the Holiday junction whereas positions 1 and 8 had relaxed homology requirements. They concluded that homology was essential to achieve recombination rates between mutant loxP spacers comparable to that of the wild-type sequence. Their success with the lox 2272 mutant suggested that positions 2 and 7 were particularly important in blocking promiscuous recombination. Langer et al. focused on the in vivo functionality of novel loxP sites in bacterial and eukaryotic systems . In their work, a suite of random loxP sites were generated by utilizing a template loxP oligonucleotide with degenerate spacers at positions 1,2,3,6,7,8. They discovered four spacers with vigorous homologous recombination and limited recombination with the wild-type sequence (m2, m3, m7, m11). Only the m2 mutant was tested and found to be functional for RMCE in bacteria and Cre-expressing eukaryotic cells (HEK-293). Other studies have evaluated the in vivo heterospecificity of the wild-type spacer versus other published loxP sequences. Siegel et al. tested heterospecific recombination between five published loxP spacer mutants and found that wild-type/FAS, wild-type/2272 and 2272/5171 pairings were incompatible while limited compatibility was found between wild-type/511 sites . However, another study found equal recombination efficiency into a BAC between the wild-type and lox 511 spacers . Recently, a phage P1 transduction assay recapitulated the sensitivity and exclusivity of the lox 2272 and lox 5171 spacer mutants . However, to the best of our knowledge, no studies have investigated if the wild-type spacer is the most efficient sequence for Cre-mediated recombination.
Given that spacer and inverted repeat mutants have been used together successfully, it may be possible to introduce numerous DNA segments into a given target molecule, chromosome or genome if a sufficient number of non-promiscuous LE/RE-spacer mutants can be identified. In the present study we have developed and executed a high-throughput in vitro strategy to rapidly identify novel LE/RE-spacer mutants and determine their specificity.
The recombination reaction
Oligonucleotide sequences (LE, RE) and primers used in the high-throughput LE/RE-spacer mutant loxP screen
5' < CTGAGCTCGGGTGCCTGTTACCGTTCGTATANNNTANNNTATACGAAGTTATGACTGCACCCATCAGAGC > 3'
5' < TGACGAATTCCTGAGCTCGGGTGCCTGT > 3'
5' < GCTAAAGCTTGCTCTGATGGGTGCAGTC > 3'
5' < TGTCTCTGGCACGCTAGGATAACTTCGTATANNNTANNNTATACGAACGGTATATGGCCTTCTGTAGCGG > 3'
5' < GTACGCGGCCGCTGTCTCTGGCACGCTAGG > 3'
5' < GTACGCGGCCGCCCGCTACAGAAGGCCATA > 3'
Sequencing and analyzing the successful recombinants
The desired 1.8 kbp recombination products were amplified by PCR using the RE_NotI forward and reverse primers (Table 2) then digested with NotI, agarose gel purified and re-circularized with T4 DNA ligase to generate a library of paired loxP recombination products in pUC19 vector (Figure 1). These plasmids were transformed into DH10B cells, grown overnight, and plated on solid media. Each individual colony (clone) represented a distinct, successful, recombination reaction between two loxP spacer sequences. A total of 5,670 clones were picked, grown overnight, and plasmid DNA was isolated and sequenced with M13 Reverse sequencing primer using previously described methods (Figure 1) . Of these clones, 4,992 yielded successful sequence (See Additional files 1 &2).
There were 3,124 reverse strand sequence reads from successful recombination reactions that were used for further analysis. However, these sequences could not be analyzed as is. First, each spacer was reverse complemented to facilitate comparisons with published loxP spacers as most spacers are published in the positive strand orientation [6, 7]. Since each sequence read represented the final product of recombination, the spacer sequence of the original input LE and RE oligonucleotides had to be inferred based on the published location of the scissile bonds  and mechanism of recombination (Figure 3) [14–16]. For each recombination reaction, the input LE oligonucleotide spacer was defined as loxP spacer position 1 from the LE/RE double inverted repeat mutant and positions 2–8 from the spacer with wild-type inverted repeats (Figure 3). Accordingly, the input RE oligonucleotide spacer was defined as position 1 from the wild-type inverted repeat loxP spacer sequence and positions 2–8 of the LE/RE double inverted repeat mutant (Figure 3).
Total reactions and mutations per position in the 8 bp loxP spacer region organized by mismatch Type
A Perl script (see Additional file 3) was run on these sequences to transform them to their inferred original format and all intermediates based on the reaction mechanisms (Figure 3) (see Additional file 4). There were a total of 47 Type I reactions, 303 Type II reactions, 284 Type III reactions and 2,490 Type IV reactions with an average number of mismatches/spacer of 0, 1.5, 2.05 and 3.56 respectively (Table 3). The minimum number of mismatches varied by type: zero in Type I or perfect matches, 1–2 in Type II, 1–4 in Type III, and 2–6 in the Type IV. Type II reactions had an approximately equal number of single and double mismatches (150 and 153 respectively). Type III reactions had 94 single mismatches, 102 double mismatches, 60 triple and 19 quadruple mismatches. Lastly, Type IV reactions had 379, 878, 785, 368, and 80 of two, three, four, five and six bp mismatches, respectively. The observation of 80 Type IV reactions with six mismatches reflects the ability of the loxP spacer to tolerate numerous mutations and continue to recombine.
Sequence characteristics of the Cre-loxP spacer region
Base frequency (%) by position of each input loxP spacer sequence that successfully recombined
Novel recombining LE/RE-spacer mutants and their promiscuity profile
Non-self recombining spacers
Selecting candidate spacers
The present study approaches loxP spacer mutant screening in an unbiased manner, assuming nothing about candidate loxP sites save for the conserved 4 and 5 spacer positions. Candidate sites were free to recombine with any other sequence in the pool of all possible sites. The present strategy also tested the general promiscuity of loxP spacer sequences. Cre-loxP spacer crosstalk was directly assayed since any two spacer sequences could interact in the Cre reaction. This revealed the promiscuity profile of each spacer sequence and identified the spacer sequences that were favoured by the Cre recombinase protein. By applying profiling algorithms typically used to analyze transcription factor binding sites , we have generated the sequence profile of the consensus, that is, the spacer most favoured by Cre recombinase (Figure 5). The most notable feature of this "sequence logo" is the abundance of guanine residues. The reason for the strong bias towards guanine content remains unclear.
A number of assumptions about Cre-loxP recombination were directly challenged by the present study, the first being that Cre-loxP recombination between identical spacer sequences is favoured over recombination between non-identical spacer sequences [6, 7]. Cre-loxP site promiscuity is not a new observation; for example initial studies identified another loxP recombination site in the E. coli chromosome, loxB . However, since most applications of Cre-loxP focus on pre-defined spacers there has been limited attention paid to promiscuity. Our data showed that for each spacer, self-recombination was not significantly greater than recombination with other spacers. Second, it has been assumed that the wild-type loxP spacer sequence is a favoured substrate for Cre recombinase. Although this assumption has never been explicitly stated, the majority of screens to identify novel spacers focus on point mutations to the wild-type spacer sequence. In the present study the wild-type spacer was detected in only 1 recombination reaction. This was a non-self recombination and not once did this sequence undergo self recombination (see Additional file 4). Our data show that GGGTAGGN is the sequence most favoured by Cre recombinase (Figure 5). Third, Lee and Saito argued that domains exist in the spacer region which cannot tolerate mutations due the requirements of strand exchange and resolution of the Holiday junction in the loxP intermediate . Lee and Saito described three spacer domains, two in which homology was required for exclusive self-recombination. We did not observe this in our data as all possible positions in the spacer sequence could tolerate mismatches and proceed to recombine in our in vitro assay. The observation of 80 reactions in our data that had all six possible degenerate spacers mismatched suggests there are few absolute spacer homology requirements required for Cre-mediated recombination (at least in an in vitro setting). Finally, it has been assumed that directed, non-promiscuous recombination is best mediated by identical spacer mutants. In fact, none of our self-recombining spacer mutants that reacted more than once were non-promiscuous. Further, we identified 16 pairs of non-self recombining spacers; some of which recombined at levels similar to the top self-recombiners (Table 5). For example, the non-self recombining spacers TGGTAGGC and TGGTAGGT reacted together four times which is identical to the number of GGGTAGGG self-recombinations (Table 5, Figure 7). All spacers, save for GGGTAGGG, reacted with equal probability with a suite of spacer sequences however in most cases there was one spacer that had a higher number of reactions. These results question the necessity, or preference, of the Cre recombinase to react with homologous identical sequences. In future Cre-loxP applications it may prove reasonable to select any two sites that show specificity for each other, regardless of whether or not they share identical spacer sequence.
Several published spacer mutants have been described in the literature as self-recombining and are non-promiscuous (at least with the wild-type loxP spacer sequence). Some sites are extremely non-reactive with other loxP sites, namely lox 2272 , m3, and m7 . We screened our data set for published sequences (Table 1) to see if these trends were recapitulated. The wild-type spacer was detected in one reaction but no other published spacers were detected at all. We did detect a number of the spacers in the Lee and Saito data  but only one of their spacers, AGGTATGC or lox 23, was detected in our self-recombining data. These results do not challenge the utility of published spacer mutants but indicate that published spacer sequences are potentially not kinetically favoured by Cre recombinase relative to other spacers. Additionally, there has been speculation in the literature that spacers at positions 2 and 7 are essential for exclusive recombination, at least between loxP mutants and the wild-type spacer [7, 8]. In our set of self-recombining sites we detected recombination at both positions 2 and 7 indicating that mismatches at these positions resulted in spacer promiscuity (see Additional file 4). In fact, all position, but not all nucleotide, mismatch combinations were detected in the complete dataset. However, our results cannot rule out the possibility that certain combinations of specific nucleotide mismatches can block promiscuous recombination of self-recombining spacers.
When using inverted repeat mutants the recombination reaction yields two different tandem loxP sites; one site with double (LE and RE) inverted repeat mutants and the other with wild-type sequence. These products are assumed to be resistant to subsequent Cre-mediated recombination. We screened all of our sequences for evidence of failure of the inverted repeat mutants to block further rounds of recombination. We found only three examples of additional recombination within our set of 3,124 successful reactions. The fact that even the most prolific spacer sequences such as GGGTAGGC did not react with the double mutants suggested that products were stable post-recombination.
There is a need for functional genomics tools that can expand the number of transgenic DNA segments that can be specifically incorporated into a target sequence. For example, synthetic biology applications require a mechanism to stitch together disparate DNA molecules from various sources in an ordered, controlled manner. Serializing RMCE or insertional recombination via inverted repeats has been limited by the small number of stable, non-promiscuous loxP sites identified to date. Furthermore, previous attempts to discover novel loxP spacers have focused on mutants slightly different than, and non-reactive with, the wild-type sequence. In this work we deploy a high-throughput pipeline that rapidly identifies novel loxP spacers, regardless of sequence composition, with stability post-recombination and limited promiscuity. The set of 12 spacer sequences detected with these properties have an increased probability of successfully tolerating multiple loxP recombinations in vivo and, therefore, these sites have potential utility.
In this study we introduce a novel, high-throughput assay to rapidly screen for successful recombination between loxP spacers using standard molecular biology techniques and high throughput sequencing. Using the results from 3,124 successful reactions, we provide the first profile and comprehensive analysis of sequence spacers favouring Cre-loxP recombination. Increasing the guanine content at all spacer positions save for position 8 resulted in increased recombination reactions. Interestingly, self-recombination is not preferred over non-self recombination and the collection of known loxP spacer mutants do not appear to be favoured by Cre recombinase. We have identified novel self and non-self recombining loxP spacer sites that expand the repertoire of sites available for further study and potential in vivo recombination applications. Further, given that these sites have been shown to be functional in combination with inverted repeat mutations, they should be useful for mediating direct serial insertions into a given target construct, chromosome or genome.
Generating two pools of LE/RE-spacer mutants
Two oligonucleotides were designed that contained loxP sites with six degenerate spacer nucleotides (positions 1,2,3,6,7,8) and two central fixed spacer nucleotides (4,5) (Table 2). One of the two oligonucleotide pools (LE) contained the lox 71 left arm mutant sequence and the other (RE) contained the lox 66 right arm mutant sequence . These oligonucleotides were synthesized (Invitrogen), purified by polyacrylamide gel electorophoresis, and double-stranded product was generated by PCR.
Generating a library of circularized LE-spacer mutants
The amplified, double stranded LE-spacer fragment was double digested with EcoRI and HindIII, agarose gel purified, and ligated to EcoRI/HindIII linearized and agarose gel purified pUC19 vector to create a library of LE mutant sequences (Figure 1). Ligation products were transformed into DH10B cells (Invitrogen) and amplified on 2xYT solid media. Supercoiled plasmid DNA was isolated by alkaline lysis and then subsequently agarose gel purified.
The recombination reaction
A 20 μL in vitro recombination reaction was set up with approximately 300 ng of supercoiled LE/pUC19 plasmid, 30 ng of double stranded RE PCR product, 200 ng of CRE Recombinase (BD Bioscience) and 1x (final concentration) Cre Recombinase Buffer (BD Biosciences). The reaction was incubated overnight at room temperature and the desired 1.8 kbp recombination products were agarose gel purified and amplified by PCR using the RE_NotI_Forward and RE_NotI_Reverse primers (Table 2).
Sequencing recombination reaction products
The amplified recombination products were digested with NotI, agarose gel purified and re-circularized with T4 DNA ligase to generate a library of paired loxP recombination products in pUC19 vector. A total of 5,670 clones were sequenced with M13 Reverse sequencing primer using previously described methods . Vector sequences were removed and reads were quality trimmed to the Q20 standard .
Identifying successful recombinants
Sequence data was analyzed using Perl v5.6.1 and MySQL v10.8 Distribution 3.23.21-beta (for pc-linux-gnu). A Perl script termed pair_count4.pl (See Additional file 3) identified recombinants within sequence traces that contained two loxP sites and applied the recombination model (Figure 3) to infer the sequence of the input oligonucleotides. Successful recombinants contained serial loxP sites with one containing the lox 66 and lox 71 inverted repeat mutants and another with the wild-type inverted repeats. The inverted repeats flanked an 8 bp spacer region. All analysis details and scripts are provided online (see Additional file 5).
Properties of recombinant spacers
Spacer sequences from the sequencing reads and the inferred input sequences were stored in a MySQL database (see Additional file 5). From these tables, all unique input spacers were identified. For each unique spacer, the number of self-recombining reactions, non-self recombining reactions, and total interactions were stored in MySQL tables. All successful recombinant spacer sequences were used to build a sequence logo representing the preferred nucleotide for Cre recombinase at each position in the spacer. A local Java implementation of the "sequence logo" algorithm  was used to generate a graphic representing the Cre recombinase binding site profile. Information content of each base is defined as I = sum × frequency (nt) × log(2) (4 × frequency (nt)). Therefore the height (H) of each base in the logo is: H = frequency (nt) × I. The effect of varying bases at each particular position in the spacer was evaluated by a Perl script that calculated the percent of the maximum recombination achieved by each polymorphism (observed/total) (see Additional file 5). This was calculated for all single nucleotide combinations flanking the axis of symmetry (positions 4,5) that were present in our data. Cytoscape, a software package for visualizing biomolecular interaction networks , was used to map out the promiscuity profile of each self-and non-self recombining spacer.
We thank Dr. Gregg Morin for helpful discussions. RAH is a Michael Smith Foundation for Health Research Scholar. This work was funded by the British Columbia Cancer Agency.
- Hoess RH, Ziese M, Sternberg N: P1 site-specific recombination: nucleotide sequence of the recombining sites. Proc Natl Acad Sci U S A. 1982, 79: 3398-3402.PubMedPubMed CentralView ArticleGoogle Scholar
- Sauer B: Functional expression of the cre-lox site-specific recombination system in the yeast Saccharomyces cerevisiae. Mol Cell Biol. 1987, 7: 2087-2096.PubMedPubMed CentralView ArticleGoogle Scholar
- Sauer B, Henderson N: Site-specific DNA recombination in mammalian cells by the Cre recombinase of bacteriophage P1. Proc Natl Acad Sci U S A. 1988, 85: 5166-5170.PubMedPubMed CentralView ArticleGoogle Scholar
- Sauer B, Henderson N: Cre-stimulated recombination at loxP-containing DNA sequences placed into the mammalian genome. Nucleic Acids Res. 1989, 17: 147-161.PubMedPubMed CentralView ArticleGoogle Scholar
- Sauer B, Henderson N: Targeted insertion of exogenous DNA into the eukaryotic genome by the Cre recombinase. New Biol. 1990, 2: 441-449.PubMedGoogle Scholar
- Hoess RH, Wierzbicki A, Abremski K: The role of the loxP spacer region in P1 site-specific recombination. Nucleic Acids Res. 1986, 14: 2287-2300.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee G, Saito I: Role of nucleotide sequences of loxP spacer region in Cre-mediated recombination. Gene. 1998, 216: 55-65.PubMedView ArticleGoogle Scholar
- Langer SJ, Ghafoori AP, Byrd M, Leinwand L: A genetic screen identifies novel non-compatible loxP sites. Nucleic Acids Res. 2002, 30: 3067-3077.PubMedPubMed CentralView ArticleGoogle Scholar
- Seibler J, Bode J: Double-reciprocal crossover mediated by FLP-recombinase: a concept and an assay. Biochemistry. 1997, 36: 1740-1747.PubMedView ArticleGoogle Scholar
- Schlake T, Bode J: Use of mutated FLP recognition target (FRT) sites for the exchange of expression cassettes at defined chromosomal loci. Biochemistry. 1994, 33: 12746-12751.PubMedView ArticleGoogle Scholar
- Baer A, Bode J: Coping with kinetic and thermodynamic barriers: RMCE, an efficient strategy for the targeted integration of transgenes. Curr Opin Biotechnol. 2001, 12: 473-480.PubMedView ArticleGoogle Scholar
- Albert H, Dale EC, Lee E, Ow DW: Site-specific integration of DNA into wild-type and mutant lox sites placed in the plant genome. Plant J. 1995, 7: 649-659.PubMedView ArticleGoogle Scholar
- Araki K, Araki M, Yamamura K: Targeted integration of DNA using mutant lox sites in embryonic stem cells. Nucleic Acids Res. 1997, 25: 868-872.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Duyne GD: A structural view of cre-loxp site-specific recombination. Annu Rev Biophys Biomol Struct. 2001, 30: 87-104.PubMedView ArticleGoogle Scholar
- Lee L, Sadowski PD: Strand selection by the tyrosine recombinases. Prog Nucleic Acid Res Mol Biol. 2005, 80: 1-42.PubMedView ArticleGoogle Scholar
- Lee L, Sadowski PD: Sequence of the loxP site determines the order of strand exchange by the Cre recombinase. J Mol Biol. 2003, 326: 397-412.PubMedView ArticleGoogle Scholar
- Araki K, Araki M, Yamamura K: Site-directed integration of the cre gene mediated by Cre recombinase using a combination of mutant lox sites. Nucleic Acids Res. 2002, 30: e103-PubMedPubMed CentralView ArticleGoogle Scholar
- Siegel RW, Jain R, Bradbury A: Using an in vivo phagemid system to identify non-compatible loxP sequences. FEBS Lett. 2001, 505: 467-473.PubMedView ArticleGoogle Scholar
- Wang Z, Engler P, Longacre A, Storb U: An efficient method for high-fidelity BAC/PAC retrofitting with a selectable marker for mammalian cell transfection. Genome Res. 2001, 11: 137-142.PubMedPubMed CentralView ArticleGoogle Scholar
- Chatterjee PK, Shakes LA, Srivastava DK, Garland DM, Harewood KR, Moore KJ, Coren JS: Mutually exclusive recombination of wild-type and mutant loxP sites in vivo facilitates transposon-mediated deletions from both ends of genomic DNA in PACs. Nucleic Acids Res. 2004, 32: 5668-5676.PubMedPubMed CentralView ArticleGoogle Scholar
- Smailus DE, Marziali A, Dextras P, Marra MA, Holt RA: Simple, robust methods for high-throughput nanoliter-scale DNA sequencing. Genome Res. 2005Google Scholar
- Hoess RH, Abremski K: Mechanism of strand cleavage and exchange in the Cre-lox site-specific recombination system. J Mol Biol. 1985, 181: 351-362.PubMedView ArticleGoogle Scholar
- Klug SJ, Famulok M: All you wanted to know about SELEX. Mol Biol Rep. 1994, 20: 97-107.PubMedView ArticleGoogle Scholar
- Schneider TD, Stephens RM: Sequence logos: a new way to display consensus sequences. Nucleic Acids Res. 1990, 18: 6097-6100.PubMedPubMed CentralView ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504.PubMedPubMed CentralView ArticleGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.PubMedView ArticleGoogle Scholar
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.