Email updates

Keep up to date with the latest news and content from BMC Research Notes and BioMed Central.

Open Access Technical Note

A combinatorial approach to the restriction of a mouse genome

Leonid V Bystrykh

Author Affiliations

Laboratory of Ageing Biology and Stem Cells, European Research Institute for the Biology of Ageing, University Medical Center Groningen, University of Groningen, Antonius Deusinglaan 1, 9700 AD, Groningen, The Netherlands

BMC Research Notes 2013, 6:284  doi:10.1186/1756-0500-6-284


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1756-0500/6/284


Received:10 September 2012
Accepted:16 July 2013
Published:22 July 2013

© 2013 Bystrykh; licensee 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.

Abstract

Background

A fragmentation of genomic DNA by restriction digestion is a popular step in many applications. Usually attention is paid to the expected average size of the DNA fragments. Another important parameter, randomness of restriction, is regularly implied but rarely verified. This parameter is crucial to the expectation, that either all fragments made by restriction will be suitable for the method of choice, or only a fraction of those will be effectively used by the method. If only a fraction of the fragments are used, we often should know whether the used fragments are representative of the whole genome. With a modern knowledge of mouse, human and many other genomes, frequencies and distributions of restriction sites and sizes of corresponding DNA fragments can be analyzed in silico. In this manuscript, the mouse genome was systematically scanned for frequencies of complementary 4-base long palindromes.

Findings and conclusions

The study revealed substantial heterogeneity in distribution of those sites genome-wide. Only few palindromes showed close to random pattern of distribution. Overall, the distribution of frequencies for most palindromes is much wider than expected by random occurrence. In practical terms, accessibility of genome upon restriction can be improved by a selective combination of restrictases using a few combinatorial rules. It is recommended to mix at least 3 restrictases, their recognition sequences (palindrome) should be the least similar to each other. Principles of the optimization and optimal combinations of restrictases are provided.

Keywords:
DNA fragmentation; Genome restriction; Randomness of genome

Background

Fragmentation of genomic DNA is very usual step in many protocols aimed at genomic analysis. For instance, studies of DNA methylation employ a genomic restriction protocol with a pair of enzymes, one of which is sensitive to the methyl cytosine in the context of a CG di-nucleotide [1-3]. Telomere length can be measured by use of a combination of for instance two restriction enzymes [4] with the expectation that the longest fragment in the digest will be the telomeric region. A reduced representation approach employs the idea that a randomly selected fraction of DNA fragments adequately represent the whole genome, yet the size of the genomic sample can be substantially reduced [5-8]. A strategy of restriction site associated markers (RAD) for high throughput genotyping is another example of this concept [9]. A few methods aim for maximal accessibility of the genome fragments for detection of some particular genomic elements, for instance retro elements [10-13] or retroviral vectors integrating into the genome [14-16].

In all those kinds of applications it is important to know the cost of deviation from the expected pattern of genomic fragmentation. Issues of genomic bias from randomness can be largely ignored in case of, for instance, restriction-based analysis of DNA methylation: they are exclusively aimed at CpG analysis rather than a randomly sampled genome, the restriction schemes used work well with CpG rich loci but miss large parts of CpG poor regions. In the case of reduced representation, biased genome sampling might affect the calculated genome diversity and phylogeny as a consequence [17]. In the case of genome accessibility, the requirements of the genomic restriction are the highest due to a risk of missing the genomic element in question and drawing incorrect conclusions as a consequence. In the search for endogenous retro elements, we deal with hundreds of those elements across the entire genome, therefore missing a fraction of those might be acceptable. In the case of an analysis of retroviral integration sites, a relatively small number of integrations per cell are generally expected. Therefore the risk and the cost of missing those few might be very high. Yet a common practice of using particular restriction schemes is highly empirical. For a long time, only one restriction enzyme was used for the recovery of retroviral integration sites [18]. It wasn’t until 2007 that Harkey et al., [14] confronted this practice and initiated a search for the optimal combination of restriction enzymes, since a single enzyme restriction was claimed to be ineffective. Ideally all fragments should be as close as possible to the preferred average. An average fragment size should fit the requirements of the assay. Note that studies performed on the analysis of the endogenous retro elements use combinations of 6-base enzyme restrictases [12,13], which are generally very rare cutters falling below 1 cut per 1 kb. Those authors claim a rather broad window of detection, namely up to 4 kb long fragments, therefore the restriction they used would serve the purpose. In a retroviral integration site analysis, the distance of the integrated vector to the restriction site is recommended to be in a range of 20–200 bp [15]. In practice this means that restriction with an enzyme which has an average frequency in the range 5–50 cuts per 1 kb will serve the purpose. Frequent 4-base cutters formally suit those requirements on average. However, it is equally important to ensure that restriction is positionally close to random, as defined for instance in [19]. Until now, this point was rarely expressed or taken in consideration. Since DNA is highly non-random and variable in local frequencies of particular nucleotides is highly unlikely to expect random distribution of restriction sites for any particular restrictase across genome. In the field of retroviral integration site analysis, authors have attempted to circumvent the genome accessibility problem by combining 4-base recognizing restriction enzymes either empirically [16,20], or using a combination of the most frequent cutters [15]. Those combinations are claimed to reach >95% of the whole mouse genome. Although increasing the average restriction frequency might be helpful, this alone is not sufficient to ensure improvement along the entire genomic sequence. It is important to verify that genomic loci void of some particular restriction site will be cut by another enzyme added to the combination and eventually reduce the variation of restricted fragment sizes to the minimum. Alternatively, if the added enzyme would cut at the same loci as another enzyme from the selected enzyme mixture, such a combination will probably be rather counter-effective. It is important therefore to check whether the spectrum of restriction fragments used fits the requirements of the analysis. This inevitably leads to the question how close the distribution of any enzyme or a combination of enzymes is to the random distribution (since an option of making equal fragments is not realistic). If it is too wide, how can it be reduced in practical terms? To illustrate this, all possible complementary palindromes were scanned along the mouse genome, those sites were analyzed for frequencies and randomness of distribution both for each site separately and in combinations. The analysis revealed that distribution of all 4-base long restriction sites is wider than expected by random. A few simple rules can be followed to achieve combinations of restriction sites which fit to a random model. None of the combinations could create a narrower distribution of restriction sites than expected by random.

Findings

Fitness to the Poisson distribution

At present detailed description of mouse and other genomes can be found at several WEB sites, such as NCBI (http://www.ncbi.nlm.nih.gov/genome webcite), Ensembl (http://www.ensembl.org/index.html webcite), or UCSC (http://genome.ucsc.edu webcite). Since the mouse genome is almost completely assembled, frequencies of all possible restriction sites can be mapped with a high precision. This work is mainly focused on an analysis of 4-base long palindromic sequences. Naturally 24 = 16 combinations of complementary palindromes are possible since the last two bases always depend on the combination of the 1st two bases. All chromosomes were scanned for each palindrome, and frequencies of all palindromes per 1 kb non-overlapping windows along the entire chromosome were recorded and further analyzed. Each restriction site showed a rather variable average, varying from 8.25 sites/kb for ATAT on chromosome X to 0.06 sites/kb for CGCG on chromosome Y. For each palindrome, estimated frequencies of cuts were compared with the random hypothesis (Poisson distribution). From the results shown on a Figure 1A is clear that AT-rich palindromes are the most frequent cutters and they deviate the most from the random distribution. Balanced palindromes (those which contain all 4 bases) tend to fit the best to the random test across all frequencies of restrictions, CG-rich palindromes show reasonable randomness but low frequencies of occurrence. The way experimentally found frequencies deviate from the Poisson model is very uniform: Poisson systematically underestimates both low and high frequencies of the distribution (Figure 2A), in other words: the real distribution is systematically wider than expected by random. This deviation can be better visualized in a kind of a “lasso” plot: by plotting observed and predicted frequencies of cuts against each other (Figure 2B). Ideally if observed and predicted frequencies are equal, such a plot should generate a straight line. Instead, it shows a “lasso”-like curve. Importantly, in the case of low cutters, the shoulders of the distribution are smaller, which can cause a better fit to the Poisson when low efficiency cutters are compared to high efficiency enzymes (Figure 1A). To test this, the bin size was varied from 1 kb to 2, 5, 10, 25 kb respectively and the fitness of the lower frequency palindromes to the random model was repeatedly measured (Figure 1B). By increasing bin size from 1 kb to 8 kb (and consequently increasing the average frequency per bin), P values from an F-test became slightly worse for balanced, non-CG containing palindromes (GATC, CTAG, TGCA). Among unbalanced palindromes, a profound change in F-test P value was recorded. In fact the difference between AT and GC rich sites appeared to be due to the differences in frequencies only. When the bin size was increased, CG rich palindromes showed very comparable frequencies per bin and randomness values as found for AT-rich palindromes. To conclude: CG-rich and AT-rich palindromes, although different in frequencies, are not different in their randomness, rather the F-test is sensitive to the frequencies instead.

thumbnailFigure 1. Frequencies of occurrence and randomness of 4-base palindromes across mouse genome. 2D plot of randomness, measured as an F-test P value versus average frequencies of tetra-palindromes on all mouse chromosomes (A) and trends of randomness for non-AT rich palindromes upon changing bin size and average frequencies per bin therefore (B). Note that GATC, CTAG, TGCA remain quite stable, whereas CG-containing palindromes are moving towards the area of AT-rich palindromes.

thumbnailFigure 2. Example of deviation from randomness for ACGT and CCGG. Broad shoulders detected upon changing bin size and average frequency per bin (A), non-linearity in expected versus observed frequencies plot (B). Note that linearity of such a plot became worse upon increased bin size and corresponding frequency of restriction.

Similarity metric by di-nucleotide count

In addition to classification of palindromes as AT-rich, GC-rich and balanced, it is useful to introduce more a precise measure of similarity or difference between palindromes. Since all of them are composed of complementary pairs of bases, a comparison of di-nucleotide composition in each palindrome is a naturally suitable metric. Each palindrome can contain up to 3 different di-nucleotides, for instance the AGCT palindrome contains AG, GC, and CT di-nucleotides. As is shown in Figure 3A, many palindromes have variable similarity to others by having 2 to 3 di-nucleotides in common. For each of 16 possible palindromes, there is at least one highly similar palindrome with 2 di-nucleotides in common.

thumbnailFigure 3. Relationship between similarity of palindromes and correlation in their restriction frequencies. A, Similarity network between all complementary palindromes. All edges connect palindromes similar by 2 to 3 di-nucleotide combinations. B, Pearson correlation frequencies for all pairs of palindromes on chr1. C, Top- and bottom-12 correlations from all pairs on chr1. Note, that the highest positive correlations are found between similar palindromes, whereas negatively correlating palindromes are of an opposite type.

The same metric was applied for a list of currently known 6-base restriction sites (Table 1). Important difference between 4-base cutters and 6-base cutters is that whereas tetra-palindromes fall naturally in 3 different categories, AT-rich, CG- rich and balanced, hexa-palindromes cannot be fully balanced. Instead, all variety of 6-base sites (totally 55 known restriction sites) falls into 4 groups: AT-rich (e.g. TTTAAA), CG-rich (like CGGCCG), but also AT-enriched (CTTAAG) and GC enriched (CTCGAG).

Table 1. Pearson-based combinations of known 6-mer restriction enzymes

Correlations between palindrome frequencies along the chromosome

Cross correlation of all tetramer palindromes across mouse chromosome 1 revealed that an average of those correlations slightly shifted into negative region (Figure 3B), the positive shoulder is longer than the negative one. Note that a correlation is analyzed through tens of thousands data points and even small correlation values are highly significant. For instance, a correlation of 0.05 along one of the smallest chromosome 19 collects approximately 60000 data points and a corresponding p-value < 0.0001. A selection of the lowest 10 and the highest 10 correlations (Figure 3C), which constitutes 19% of all correlated pairs, appeared to be enriched with unbalanced palindromes (note that the number of balanced and unbalanced palindromes is equal). Of those, palindromes similar in their di-nucleotide composition and type of enrichment (GC or AT) show significant positive correlations. The complementary palindromes from opposite enrichment groups tend to show negative correlation. Analysis of correlations was repeated for all chromosomes. Similar results were obtained for the same palindrome along the whole mouse genome (Figure 4, Additional file 1: Table S1). Note that only X and Y chromosomes deviated most drastically from the common trend whereas all autosomes showed highly similar results. In practical terms, this finding ensures that analysis of one autosome is already well representative of the whole mouse genome. It saves computational time as well as simplifying visualization of the data. The trends found by di-nucleotide context and palindrome enrichment type could be further followed and analyzed in more detail during the selection of a series of best anti-correlating palindromes. Obviously, the purpose of such selection is to test the possibility of finding combinations of various restriction sites where a desired average frequency is achieved together with minimal bias and maximally possible randomness of distribution.

thumbnailFigure 4. Correlation between all 4-base palindromes on all mouse chromosomes. The most deviating from the common trend were ChrX (red line) and ChrY (blue line).

Additional file 1: Table S1. Pearson correlation data for all tetra palindromes all mouse chromosomes. Table S2: Selection of all palindromic combinations using 3 different approaches. Table S3: Pearson correlation data for all tetra palindromes mouse chromosome 1. Table S4: Commercial names for all 4- and 6- base restriction endonucleases. Additional files. The Supplementary files S1, S2, S3, S4 are deposited in an open access repository https://easy.dans.knaw.nl/ui/datasets/id/easy-dataset:53177 webcite. http://www.persistent-identifier.nl?identifier=urn:nbn:nl:ui:13-d9zq-rm webcite.

Format: ZIP Size: 196KB Download fileOpen Data

A negative correlation strategy

We can try to begin with each of the 14 palindromes for which restriction enzymes are available and follow the most negatively correlating site as the best candidate for a mixture. At this stage ATAT and TATA palindromes were excluded from analysis since they do not represent any known restriction enzyme. An algorithm will stepwise select for the most negatively correlating palindrome and consecutive count of combined frequencies of already selected palindromes. An example of such a strategy is shown in Figure 5 for the selected sequences AATT, GGCC, AGCT, and CATG. A network graph shows that the total number of negatively correlating palindromes gradually decreases upon addition of the next palindrome to the mixture. Whereas AATT is negatively correlating with 8 different palindromes, addition of the third palindrome, AGCT shows only two candidates, ACGT and CATG. Upon addition of the last palindrome, CATG, no more candidates to the mixture were revealed (in a threshold < −0.05). Note that the correlation value between the selected group and every subsequently added palindrome is decreasing (and the cumulative correlation value is decreasing as well, Figure 5B). At the same time, randomness (fitness to the Poisson distribution) of the combined palindrome frequencies is asymptotically increasing, as shown in the “lasso” plot (Figure 5C), namely a plot of observed frequencies versus expected gradually became linear. A full screen of all possible starting palindromes revealed that, in fact, good randomness is already achieved in a mixture of 3 palindromes. Further additions are useful mainly to satisfy average restriction frequency.

thumbnailFigure 5. Example of a combination by negative correlation strategy. A, Choices of negatively correlating palindromes (with a threshold < −0.05) along consequently made choices, AATT, GGCC, AGCT, CATG. B, Parameters for the selected group gradually reaching their maximum/minimum values. C, a “lasso” plot shows improved linearity upon gradual addition of each subsequent palindrome to the mixture.

When starting with different palindromes, out of 14 trials (with all possible palindromes as the 1st choice) only 10 unique combinations were found (for details see Additional file 1: Table S2). In this series there is a strong effect of the most frequent and most biased palindromes, like TTAA, AATT, which always require CG rich palindromes for compensation, the most frequent compensating palindrome was found to be GGCC. A network of 1st 2 choices made by the negative correlation strategy illustrates this trend (Figure 6A) and summarizes a “popularity contest” among different palindromes. Figure 6B shows the total frequency of palindromes in a complete series of 5 out of 14. Similarly, AATT, AGCT, CATG and GGCC appeared to be the most frequent choices for a Pearson-based strategy.

thumbnailFigure 6. Typical choices and total numbers of unique sets in different strategies of combining restriction enzymes. A, a network representation of the first two choices made by the most negative correlation strategy when using all possible palindromes. Note that AATT and GGCC are the most frequent choices. B, a histogram count of all used palindromes using different strategies. Note that Pearson and Poisson strategies used 14 trials in total, the combinatorial approach used 60 trials (total unique combinations).

During this analysis, it was found that there is an approximate solution for the laborious recalculation of correlation coefficients upon addition of a new palindrome to the mixture (it was actually done for all discussed examples). In fact the next best palindrome can be found using the primary paired correlation table by estimating the weighted (by restriction frequency) average between two previously added components (group or single palindrome, more details in Additional file 1: Table S1, page example), the palindrome with the most negative correlation to the group can be selected as the next best candidate to the mixture.

Use of an F-test for finding palindrome combinations

Considering such a small number of found combinations listed above we could try some alternative strategies. For instance, a search by Pearson correlations works well with biased distributions, but it is probably less sensitive when dealing with more balanced and randomized palindromes. With this idea in mind, the F-test was tried instead. Screening through all starting palindromes showed good agreement with AT rich and balanced palindromes. One new combination with a good randomness was revealed (Additional file 1: Table S2). However, this algorithm appeared to be completely nonfunctional in the case of low frequency GC rich palindromes, such as CGCG, TCGA, ACGT and GCGC. Apparently such enrichment occurs due to the problem of the F-test estimation (mentioned earlier in this text): with increased frequency of the restriction, the F-test value became worse, which hinders the search for the optimal randomization. This approach is inferior to the selection by negative correlation and can be recommended for balanced and frequent palindromes only.

Comparison with combinatorial model

It was already mentioned above (Figure 3C) that the resulting combinations based on a Pearson correlation revealed interesting combinatorial rules: similar palindromes will be naturally excluded and the most dissimilar palindromes will be selected (as they are negatively correlating). We can follow this trend in more detail by counting cumulative frequencies of all possible di-nucleotides within previously found combinations of palindromes (Additional file 1: Table S2). Note that since each palindrome can contain 3 different di-nucleotides at maximum the non-redundant set of restriction sites (no repeated use of di-nucleotides is allowed) should not be higher than 5. Although a maximum of 15 di-nucleotide combinations out of 16 possible could be reached, the results show that maximal coverage of possible di-nucleotides in a set of 15 palindromes is 14. Out of 8568 combinations screened, a list of 60 combinations containing 14 different di-nucleotides was generated. A variation of total restriction frequencies varied from 1 to 19 per 1 kb bin size. Selective testing for randomness with the F-test demonstrated that these are quite credible combinations quite closely resembling the sets generated by the negative correlation strategy (Figure 7). Overall, the cloud of all combinations found is broad in terms of an average frequency and fitness to the random distribution. Out of 60 total combinations, all except CGCG and GCGC were equally represented in these series (Figure 6B). A selection by Pearson correlation is localized with few more examples found by combinatorial approach. Both of them are substantially better than combinations found using F-test fitness to Poisson distribution. Only one alternative combination with high total frequency of restriction was found in a set of F-test generated series, the same combination was the best regarding randomization of di-nucleotides.

thumbnailFigure 7. Combinations of palindromes generated by different approaches. Pink dots, selection by negative Pearson correlation; green, selection by F-test; blue, all combinations by 5 from 14 palindromes with 14 di-nucleotides used, red, combinations from literature:1, Harkey et al. [14], 2, Bystrykh et al. [21], 3, Biasco et al. [20], 4, Gabriel et al. [15], which is identical to simplified Wang et al. [16]. Details in a Additional file 1: Table S1.

Previously suggested combinations

As already mentioned before, a use of a single cut for studies of retro elements is problematic, since neither frequency nor accessibility of fragments are at their optimum. In the case of retroviral integration site analysis, one restriction enzyme, cutting at TTAA was used for quite some time. It is only since 2007 [14], when a mixture of enzymes was suggested in order to maximize genome accessibility. Figure 7 and Additional file 1: Table S2 demonstrate average frequencies and randomness of a few published combinations. A mix by Wang et al. [16] is unnecessarily redundant. After removing redundancy it became identical to a mix suggested by Gabriel et al. [15], who were aiming on the maximal restriction frequency ignoring a di-nucleotide bias in selecting two AT-rich enzymes, one balanced (also AT containing) and no GC rich enzymes. As a consequence, such a mixture is highly non-random and the spread of the cut fragment size will be unnecessarily broad: it will systematically undercut GC-rich regions, many of which are situated in promoter regions, which are potential targets for integration of retroviral vectors. Similarly, an empirical set by Biasco et al. [20] is almost the same in average frequency as given by Harkey et al., (2007). However, the difference in randomness of those mixtures is dramatic. Extensive restriction of AT rich regions by mixing the most frequent restrictases faces a risk of generating too small fragments (primarily in AT rich regions) to be used for site identification. On the contrary, the addition of GC-rich restrictases, albeit with a minimal increment in average frequency, will substantially improve randomness of the cut and therefore will improve the overall chances of detecting integrated vectors. Note that a mixture suggested in our recent review [21] was purely based on the combinatorial principle before the current analysis was done.

Use of 6-base restrictases

Whereas 4-base restrictases show frequencies of restrictions in a range from 0.05 to 8 cuts per 1 kb, 6- base cutters are approximately 10 times less frequent (see data for chr1 in Additional file 1: Table S3 for details). Therefore their use must be justified after critical evaluation of the method of choice. General trends among 6-base cutters are similar to those found and described above for tetra-palindromes. AT-rich restriction sites are the most frequent, balanced 6-cutters do not exist. Instead AT- or GC- enriched palindromes show intermediate frequencies, and GC-rich sites are the least frequent along the genome. If two restrictases would be combined, the combinatorial rule still stands: AT-rich palindromes a best matching to GC-rich (as can be revealed by negative correlations). All other palindromes, however, behave less predictable (regarding combinatorial rule) and more elaborate concepts should be developed. For practical purpose, if any selected combinations will be implemented, names of known 4- and 6-base restriction endonucleases and corresponding palindromes are provided in a Additional file 1: Table S4.

Discussion

Fragmentation of genomic DNA by restriction remains a routine element in modern genomic studies. Unfortunately, the selection of particular restriction enzymes for the application is often empirical and not subjected to any optimization or validation. An issue of prediction and distribution of nucleotide dimers, restriction sites and small words was previously discussed in the literature. Originally prediction of restriction sites was based on a single-base probabilistic rule [22]. In the case of microbial genomes, near to random distributions were suggested for restriction sites [19] or small words [23]. However, a bias in di-nucleotide frequencies [24] and genomic k-mers [25] were also documented in large genomes. Therefore the assumption of randomness should be made carefully for analysis of short DNA sequences, including restriction sites. The concern that restriction with a particular enzyme or mixture of enzymes might substantially deviate from random did not spread far outside the field of theoretical genome biology. For instance, in studies of genomic retro elements, restriction of DNA is performed with empirical combinations of 6-base restrictases [12,13], and no issue of randomness is addressed. Note that 6-base cutters show considerably lower frequencies compared to 4-cutters (see some examples in [21]) and therefore inferior in frequencies to the examples of combinations given in this paper. In fact use of 6-base cutters for fragmentation of large genomes should be theoretically justified, otherwise not recommended. In the field of population genetics, the problem of biased sampling and its possible drawbacks is already introduced [17]. In a recent publication it was anticipated that the performance of genome sequencing itself depends on representative and random fragmentation [26]. An issue of non-optimal fragment size is one of many problems recognized in chromosome conformation capture protocols [27,28]. Yet, in many of those publications authors have used 1 or 2 empirically selected 6- or 4-base cutters to create genomic fragments. In similar studies [6,7] only a single 4-base restriction enzyme was used for genome digestion. Van Tassel et al., [6] selected their 4-base cutter by computer analysis to minimize the appearance of repeated elements in their virtual fragments. This selection makes sense to improve efficiency of mapping. However, as was already mentioned here, a single cut strategy is always inferior to the combination of restrictases regarding its randomness. Some genome restriction based protocols have already reached clinical studies, for instance gene therapy treated patients. In those studies, the problem of accessibility has been already identified and discussed [14,15]. However authors still focus their attention on an average restriction frequency without much consideration of the randomness. As a result, authors almost empirically suggested mixtures of enzymes with low to medium randomness (examples are shown in Figure 7). Note that those few suggested combinations are further replicated in more recent clinical studies. To conclude, the problem of a proper fragmentation of genomic DNA and its verification for randomness remains an unexplored option for further improvement. As it is demonstrated here, none of the restriction sites are fully random. However, some restrictases are more random than others. Combinations of restrictases presented in this paper provide the best random distribution across the mouse genome.

The F-test was used here for the assessment of randomness, which although superior to the Chi-squared test [29], should be used with care. It is sensitive to the average frequency of the distribution. Nevertheless the F-test served the purpose quite well, namely it helped to compare different restrictases, and also helped to establish how combinations of restriction enzymes could improve the randomness of genomic restriction.

Although the combinatorial principles revealed might look puzzling and unexpected, the explanation is quite simple. The reason similar palindromes correlate is because they correlate to the regions enriched in subsequences which they contain. If they contain similar subsequences, for instance similar di-nucleotides, they will tend to be positively correlating. For the same reason, complementary di-nucleotides (AT vs GC, AA vs GG etc.) will tend to show negative correlation. This combinatorial rule is essentially based on Markov chains models, it is well reproducible in a different genomic context (see consistency of correlations across autosomes), it is also resistant to SNP variations and other minor mutations in a population of considered species. In the case of balanced palindromes, this rule will be more problematic, because 4 bases cannot be enriched all at the same time: any enrichment also means loss of complementary bases or di-nucleotides. Consequently it will hamper enrichment of corresponding subsequences. Although the combinatorial principle was found to be useful as a guideline, some deviations from it were also evident. For instance, CG-containing palindromes tend to correlate with each other more than would be expected from similarity measurements. Mutually exclusive pairwise palindromes with relatively low frequency might co-occur in a bigger set of palindromes due to the dominating effect of more frequent and biased palindrome.

There can be many factors indirectly causing non-randomness of restriction sites along the genome. It might reflect an evolutionary expansion of DNA by self-duplication, selective pressure in coding regions and in CpG islands for CG- rich palindromes. Some evolutionary pressure can be expected for ATAT and TATA sequences due to their similarities to TATA boxes. It is, however, difficult to expect direct selective pressure on the distribution of restriction sites themselves (unless a clear functional demonstration of those sequences will be provided). Instead, an effect of local enrichment for sub-words, for instance di-nucleotides, presented in GC or AT rich palindromes is most likely. This selectivity of genomic landscape fundamentally breaks the principle of equal chance along the genome and this might be a primary cause of deviation from the Poisson model. Analytically, the Poisson model (see also Methods) can be adjusted to such situations, if λ is the average frequency of an event and k is the expected exact frequency, we can consider λ as being a function of k (equation can be found in the Methods). The notion that the genomic landscape is primarily responsible for non-randomness of restriction sites is an important example how genomic non-randomness in general can be interpreted beside other existing explanations. As an example, Falconer et al. [30] found a wider distribution of segregating DNA template strands than predicted by a random distribution (very similar to the case of non-randomness found in this paper). Such non-randomness was interpreted as selective process, disregarding the alternative option of having non-random (selective) genomic background for the fairly random process of template segregation. Those two options, not one, should be considered for further experimental verifications.

During preparation of our recent review [21] on a subject of clonal analysis of hematopoietic cells, we briefly suggested the combinatorial principle of mixing restrictases for optimal genome fragmentation, which was hypothetical at that time. This paper provides the necessary background to this concept in sufficient detail and demonstrates how this principle helps in analysis and finding the best combinations of restrictases. Although mainly tetra-nucleotide palindromes were studied here, the strategy described can be applied to restrictases of any length.

Conclusions

Enzymatic restriction is a very frequent step in fragmentation of genomic DNA. A question of the randomness of restriction is often expected or implied, but rarely verified. A genome-wide restriction analysis of 4-base cutters revealed significant deviation from randomness for practically all tested restriction enzymes. The randomness can be improved by combining separate restrictases in the mixture. A combinatorial approach is probably the simplest but most effective principle to achieve such improvement.

Methods

For the mouse genome, the Bioconductor (version 2.9) metadata package Bsgenome.Mmusculus.UCSC.mm10 for R version 2.15 was used which consists of the mmu10 assembly from UCSC (based on C57BL/6 J mouse strain). This package was used in conjunction with R (version 2.14.1) to scan the genome for the presence of all 16 possible complement palindromes in non-overlapping bins of 1 kb in size. All bins with N < 500 were included in the resulting table. For each complement palindrome, fitness to the Poisson distribution per chromosome was calculated using a custom Python script and Gnumeric spreadsheet (in Linux). For a random model, a Poisson distribution was used as following:

<a onClick="popup('http://www.biomedcentral.com/1756-0500/6/284/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1756-0500/6/284/mathml/M1">View MathML</a>

Where λ is observed average, k is exact number of events expected at given average, f(k;λ) is a probability of such event.

Randomness of restriction was tested by comparing for observed and Poisson predicted data sets using the F-test as commonly described [27]. A p-value for the two-tailed hypothesis test comparing the variances of two populations, in our case palindrome frequencies and Poisson distribution, was calculated using the built in F-test function in the Gnumeric spreadsheet program. When used in the text, this is referred as an F-test P value or F-test P on the plot. Pearson correlation was calculated along the entire chromosome of choice using a custom script in python, using a double-pass strategy and the equation;

<a onClick="popup('http://www.biomedcentral.com/1756-0500/6/284/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1756-0500/6/284/mathml/M2">View MathML</a>

The probability of correlation was assessed using the t-test transformation.

<a onClick="popup('http://www.biomedcentral.com/1756-0500/6/284/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1756-0500/6/284/mathml/M3">View MathML</a>

and an online probability calculator http://www.statstodo.com/TTest_Pgm.php webcite.

Through the text all palindromes are mentioned by their DNA sequence. Names of the corresponding restriction endonucleases are provided in the Additional file 1: Table S4.

Competing interests

Author declares no competing interests.

Acknowledgments

Author thank Erik Zwart for help in downloading of mouse genome and its virtual restriction, Rupert Overall for reading the manuscript and useful remarks. The project was partially supported by HEM-ID grant (http://www.hemid.eu) and NIRM grant (http://www.nirm.nl).

References

  1. Yegnasubramanian S, Lin X, Haffner MC, DeMarzo AM, Nelson WG: Combination of methylated-DNA precipitation and methylation-sensitive restriction enzymes (COMPARE-MS) for the rapid, sensitive and quantitative detection of DNA methylation.

    Nucleic Acids Res 2006, 34(3):e19. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Khulan B, Thompson RF, Ye K, Fazzari MJ, Suzuki M, Stasiek E, Figueroa ME, Glass JL, Chen Q, Montagna C, Hatchwell E, Selzer RR, Richmond TA, Green RD, Melnick A, Greally JM: Comparative isoschizomer profiling of cytosine methylation: the HELP assay.

    Genome Res 2006, 16(8):1046-55. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Kitamura E, Igarashi J, Morohashi A, Hida N, Oinuma T, Nemoto N, Song F, Ghosh S, Held WA, Yoshida-Noro C, Nagase H: Analysis of tissue-specific differentially methylated regions (TDMs) in humans.

    Genomics 2007, 89(3):326-37. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Aviv A, Hunt SC, Lin J, Cao X, Kimura M, Blackburn E: Impartial comparative analysis of measurement of leukocyte telomere length/DNA content by Southern blots and qPCR.

    Nucleic Acids Res 2011, 39(20):e134. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Baird NA, Etter PD, Atwood TS, Currey MC, Shiver AL, Lewis ZA, Selker EU, Cresko WA, Johnson EA: Rapid SNP discovery and genetic mapping using sequenced RAD markers.

    PLoS One 2008, 3(10):e3376. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Van Tassell CP, Smith TP, Matukumalli LK, Taylor JF, Schnabel RD, Lawley CT, Haudenschild CD, Moore SS, Warren WC, Sonstegard TS: SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries.

    Nat Methods 2008, 5(3):247-52. PubMed Abstract | Publisher Full Text OpenURL

  7. Luca F, Hudson RR, Witonsky DB, Di Rienzo A: A reduced representation approach to population genetic analyses and applications to human evolution.

    Genome Res 2011, 21(7):1087-98. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Monson-Miller J, Sanchez-Mendez DC, Fass J, Henry IM, Tai TH, Comai L: Reference genome-independent assessment of mutation density using restriction enzyme-phased sequencing.

    BMC Genomics 2012, 13:72. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  9. Miller MR, Dunham JP, Amores A, Cresko WA, Johnson EA: Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers.

    Genome Res 2007, 17(2):240-8. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Ishihara H, Tanaka I: Detection and cloning of unique integration sites of retrotransposon, intracisternal A-particle element in the genome of acute myeloid leukemia cells in mice.

    FEBS Lett 1997, 418(1–2):205-9. PubMed Abstract | Publisher Full Text OpenURL

  11. Ribet D, Dewannieux M, Heidmann T: An active murine transposon family pair: retrotransposition of "master" MusD copies and ETn trans-mobilization.

    Genome Res 2004, 14(11):2261-7. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Takabatake T, Ishihara H, Ohmachi Y, Tanaka I, Nakamura MM, Fujikawa K, Hirouchi T, Kakinuma S, Shimada Y, Oghiso Y, Tanaka K: Microarray-based global mapping of integration sites for the retrotransposon, intracisternal A-particle, in the mouse genome.

    Nucleic Acids Res 2008, 36(10):e59. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. An W, Davis ES, Thompson TL, O'Donnell KA, Lee CY, Boeke JD: Plug and play modular strategies for synthetic retrotransposons.

    Methods 2009, 49(3):227-35. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Harkey MA, Kaul R, Jacobs MA, Kurre P, Bovee D, Levy R, Blau CA: Multiarm high-throughput integration site detection: limitations of LAM-PCR technology and optimization for clonal analysis.

    Stem Cells Dev 2007, 16(3):381-92. PubMed Abstract | Publisher Full Text OpenURL

  15. Gabriel R, Eckenberg R, Paruzynski A, Bartholomae CC, Nowrouzi A, Arens A, Howe SJ, Recchia A, Cattoglio C, Wang W, Faber K, Schwarzwaelder K, Kirsten R, Deichmann A, Ball CR, Balaggan KS, Yáñez-Muñoz RJ, Ali RR, Gaspar HB, Biasco L, Aiuti A, Cesana D, Montini E, Naldini L, Cohen-Haguenauer O, Mavilio F, Thrasher AJ, Glimm H, Von Kalle C, Saurin W, Schmidt M: Comprehensive genomic access to vector integration in clinical gene therapy.

    Nat Med 2009, 15(12):1431-6. PubMed Abstract | Publisher Full Text OpenURL

  16. Wang GP, Berry CC, Malani N, Leboulch P, Fischer A, Hacein-Bey-Abina S, Cavazzana-Calvo M, Bushman FD: Dynamics of gene-modified progenitor cells analyzed by tracking retroviral integration sites in a human SCID-X1 gene therapy trial.

    Blood 2010, 115(22):4356-66. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Pérez-Enciso M, Ferretti L: Massive parallel sequencing in animal genetics: wherefroms and wheretos.

    Anim Genet 2010, 41(6):561-9. PubMed Abstract | Publisher Full Text OpenURL

  18. Schmidt M, Glimm H, Lemke N, Muessig A, Speckmann C, Haas S, Zickler P, Hoffmann G, Von Kalle C: A model for the detection of clonality in marked hematopoietic stem cells.

    Ann N Y Acad Sci 2001, 938:146-55. PubMed Abstract | Publisher Full Text OpenURL

  19. Churchill GA, Daniels DL, Waterman MS: The distribution of restriction enzyme sites in Escherichia coli.

    Nucleic Acids Res 1990, 18(3):589-97. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Biasco L, Ambrosi A, Pellin D, Bartholomae C, Brigida I, Roncarolo MG, Di Serio C, Von Kalle C, Schmidt M, Aiuti A: Integration profile of retroviral vector in gene therapy treated patients is cell-specific according to gene expression and chromatin conformation of target cell.

    EMBO Mol Med 2011, 3(2):89-101. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Bystrykh LV, Verovskaya E, Zwart E, Broekhuis M, De Haan G: Counting stem cells: methodological constraints.

    Nat Methods 2012, 9(6):567-74. PubMed Abstract | Publisher Full Text OpenURL

  22. Wijsman EM: Optimizing selection of restriction enzymes in the search for DNA variants.

    Nucleic Acids Res 1984, 12(23):9209-26. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Fofanov Y, Luo Y, Katili C, Wang J, Belosludtsev Y, Powdrill T, Belapurkar C, Fofanov V, Li TB, Chumakov S, Pettitt BM: How independent are the appearances of n-mers in different genomes?

    Bioinformatics 2004, 20(15):2421-8. PubMed Abstract | Publisher Full Text OpenURL

  24. Karlin S, Mrázek J: Compositional differences within and between eukaryotic genomes.

    Proc Natl Acad Sci USA 1997, 94(19):10227-32. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Chor B, Horn D, Goldman N, Levy Y, Massingham T: Genomic DNA k-mer spectra: models and modalities.

    Genome Biol 2009, 10(10):R108. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  26. Young AL, Abaan HO, Zerbino D, Mullikin JC, Birney E, Margulies EH: A new strategy for genome assembly using short sequence reads and reduced representation libraries.

    Genome Res 2010, 20(2):249-56. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Dekker J: The three ‘C’ s of chromosome conformation capture: controls, controls, controls.

    Nat Methods 2006, 3(1):17-21. PubMed Abstract | Publisher Full Text OpenURL

  28. Cournac A, Marie-Nelly H, Marbouty M, Koszul R, Mozziconacci J: Normalizatio An issue of not optimal fragment An issue of not optimal fragment n of a chromosomal contact map.

    BMC Genomics 2012, 13:436. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. Altman DG: Practical statistics for medical research. Boca Raton: Chapman & Hall/CRC; 1999. OpenURL

  30. Falconer E, Chavez EA, Henderson A, Poon SS, McKinney S, Brown L, Huntsman DG, Lansdorp PM: Identification of sister chromatids by DNA template strand sequences.

    Nature 2010, 463(7277):93-7. PubMed Abstract | Publisher Full Text OpenURL