Previous Article |
Table of Contents
| Next Article
APPLIED PHYSICAL SCIENCES
Restriction mapping in nanofluidic devices



Departments of*Physics and
Molecular Biology, Princeton University, Princeton, NJ 08544
Contributed by Robert H. Austin, May 12, 2005
| Abstract |
|---|
|
|
|---|
We have performed restriction mapping of DNA molecules using restriction endonucleases in nanochannels with diameters of 100-200 nm. The location of the restriction reaction within the device is controlled by electrophoresis and diffusion of Mg2+ and EDTA. We have successfully used the restriction enzymes SmaI, SacI, and PacI, and have been able to measure the positions of restriction sites with a precision of
1.5 kbp in 1 min using single DNA molecules.
To date, the most powerful method for constructing restriction maps of long DNA molecules (100 kbp and above) is the one developed by Schwartz and coworkers (5, 6). Their technique consists of stretching the DNA on a surface to establish a one-to-one mapping between spatial and genomic position, initiating the restriction by exposing the DNA to restriction enzymes and Mg2+, and optically observing the lengths of the resulting fragments. An elegant study of the same basic approach showing separation of specific binding and induced cutting was published by Taylor et al. (7).
A fundamental principle in determining the error of a measurement, and a common strategy to reduce the experimental error, is to take multiple, statistically independent measurements of the same quantity. Note that the fixing of the molecule to the surface in Schwartz's approach prevents any fluctuations, and hence multiple molecules have to be observed to obtain statistically independent measurements of the same cut position. Moreover, genomic-length DNA molecules stretched on surfaces often exhibit inhomogeneous stretching, including breaks, and thus averaging over multiple molecules becomes imperative (6, 8). In contrast, a molecule stretched inside a nanochannel is not subject to any forces other than those causing lateral confinement and thus is able to thermally relax and fluctuate around an equilibrium conformation. The evaluation of a single molecule is thus sufficient if complete digestion is achieved. A detailed description of the statistics and dynamics of DNA molecules in nanochannels is published in ref. 9.
The main challenge in employing the concept of stretching and mapping DNA inside a closed fluidic system is to separate the steps of stretching and cutting. We have solved this problem by introducing the DNA electrophoretically into nanochannels and controlling the concentration of the enzyme cofactor Mg2+ in the device shown in Fig. 1. It consists of a microfluidic "loading" channel containing the DNA to be analyzed, the restriction enzyme, and EDTA, and a microfluidic "exit" channel containing Mg2+ and the restriction enzyme. The two microfluidic channels are linked by 10 nanofluidic channels of
100 nm in diameter. We expect that, in the absence of an applied voltage along the nanochannels, a Mg2+ gradient will be established as the result of Mg2+ and EDTA diffusion, and chelation of Mg2+ by EDTA. With an applied voltage, we expect the higher mobility of the Mg2+ ions to result in a constant Mg2+ concentration equal to that on the "exit" side. Note that the particular shape of the Mg2+ concentration is not crucial to the functioning of the device, only the facts that there is a concentration sufficient for restriction in the active region, i.e., the nanochannels, and a negligible concentration in the "loading" microchannel.
We present here restriction mapping of DNA molecules inside nanochannels controlled by this Mg2+-gradient, discuss the attainable resolution, and discuss how the cuts become visible.
| Methods |
|---|
|
|
|---|
1 µm deep and 100 µm wide were patterned by standard optical lithography and reactive ion etching using a CF4/H2 chemistry (10). We used 1.4-µm-thick AZ5214-E photoresist (Clariant, Somerville, NJ), exposed by using a mask aligner. The ratio of CF4 to H2 during the etching was 26.5-6, at a pressure of 35 mtorr (1 torr = 133 Pa) and a power density of 1 W/cm2. The microchannels were separated from each other by either 100 or 200 µm. Holes for accessing the microfluidic system were drilled by using dental sandblasting tool (Danville Engineering, San Ramon, CA). Nanochannels were prepared by focused ion milling after first coating with a 5-nm Au film (11). The gold was removed by using aqua regia, and the device was sealed with a quartz coverslip via quartz-quartz bonding (4).
Microscopy. We imaged DNA molecules and the Mg2+-sensitive dye using a 1.4-numerical aperture oil-immersion objective (Nikon) and an intensified charge-coupled device camera (Roper Scientific) under illumination with the 488 nm-line from an Ar/Kr ion laser. Images were collected by using 10-ms exposures spaced by
250 ms. The laser beam was shut between the frames of acquisition.
Magnesium Imaging. The "exit" side was filled with a solution containing 10 mM Tris, 1 mM DTT, and 100 µg/ml POP6 (Applied Biosystems). The "loading" side was filled with an identical solution except for the addition of 5 mM EDTA. Magnesium-acetate (5 mM) was placed into one reservoir on the "exit" side, and the Mg2+ was brought into the microchannel by applying 10 V across both microchannels, with no voltage across the nanochannels. During image acquisition, the microchannel voltage was 2 V, and the nanochannel voltage was 1 V. Images were acquired at a frame rate of 1 Hz.
Restriction Protocols. DNA was dyed with the bis-intercalating dye TOTO-1 (Molecular Probes) at a ratio of 1 dye molecule per 20 bp.
-DNA was obtained from New England Biolabs and used at a concentration of 0.5 µg/ml or 1 µg/ml. The basic buffer for restriction using SmaI, SacI, and PacI consisted of 20 mM Tris-acetate, 50 mM potassium-acetate, 1 mM DTT, 100 µg/ml BSA, and 100 µg/ml POP6. The restriction enzymes (New England Biolabs) were used with concentrations of
3 nM for SmaI,
3 nM for SacI, and
13 nM for PacI. We chose these concentrations because of the specific binding constant for SmaI of
109 M-1 (12). The solution that was filled into the "entry"-side also contained 5 mM EDTA. After filling of the device and inspection under the microscope, Mg2+ was brought into the "exit"-side microchannel by the method described in the Mg2+ imaging. During the restriction experiment, 2 V was applied over each microchannel. The potential of the two microfluidic systems was offset by 1 V for 100-µm-long nanochannels and 2 V for 200-µm-long nanochannels. Sealed chips could be used up to 5 days.
|
10 µg/ml with centrifugal filters (Microcon YM-100, Millipore). | Results and Discussion |
|---|
|
|
|---|
To perform real-time observations of restriction in nanochannels, the reaction has to be completed within seconds, to avoid photobleaching. This can be achieved by adjusting the enzyme concentration so that all recognition sites are occupied before they enter the channel. In that case, the rate-limiting step is the release of the DNA from the enzyme after cutting, or
1 s (14). Note that only a subset of restriction enzymes is able to bind to the correct sequence in the absence of divalent metal ions. Some of these are SmaI and EcoRI, for which the specific and nonspecific binding constants have been determined (12, 15). Kinetic tests of
-DNA restriction by SmaI using a moving Mg2+ front inside a microchannel has confirmed that the cutting occurs on a time scale shorter than 4 s with our combination of buffer and enzyme concentrations. We thus expect that in the case of SmaI, fully decorated DNA molecules enter the channel, and complete digestion of all accessible recognition sites occurs within a few seconds. If an enzyme does not exhibit specific binding in the absence of divalent metal ions, the enzyme needs to diffuse along the DNA molecule or bind to the DNA from solution before cutting. These additional steps will lead to somewhat longer times between entry of the DNA into the channel and cutting. However, the rate should rise with rising enzyme concentration.
|
-DNA (48.5 kbp) using the enzymes SmaI and SacI in nanochannels where the DNA was stretched to
40% and 30% of its full length, respectively (Fig. 3). Molecules were observed for
1 min, and both typical frames as well as time traces are shown. The latter is a stack of lines, each of which shows the intensity along the molecule for a single 10-ms frame. SmaI cuts
-DNA into fragments of 19.4, 12.2, 8.3, and 8.6 kbp, in that order. All four fragments were clearly observed. SacI cuts
-DNA into three fragments 22.6, 0.9, and 24.8 kbp long. Only the longer two fragments could be clearly observed. Fig. 3. demonstrates that cuts become visible because restriction fragments diffuse in the channels. We do not observe an obvious recoil of the ends at the cutting positions, probably because the thermalized polymer has a constant stretching independent of its length. Hence, the time of cutting can be inferred with a precision equal to the diffusion time of the fragment through a length corresponding to the resolution of the microscope. Note that the mild bleaching of DNA during the observation time in Fig. 3 may be either photobleaching or bleaching due to high Mg2+ concentrations. Thus, observation times longer than 3 min were not possible. However, one could easily extend the observation periods using lower Mg2+ concentrations and lower read-out times, possibly in connection with the use of an oxygen-scavenging system.
|
We have analyzed the size of the restriction fragments by fitting the intensity along the molecules to (4)
![]() | [1] |
Here, I is the intensity, the sum runs over all fragments, erf is the error function, w is the width of a segment, c its center,
is the sharpness of the molecules end, and b is a constant describing the background. Each experimental data frame was fit individually and thus, if frames are spaced by a time longer than the relaxation time for Brownian fluctuations, the measured lengths are statistically independent. The exact length of a fragment can hence in principle be determined from a single molecule by averaging over those independent measurements.
For small molecules and short exposure times (10 ms), we have found that
can become comparable to w, and hence w does not serve as a good measure of the amount of DNA in a single restriction fragment. Therefore, we determined the relative lengths of restriction fragments by integrating the fitted function and using the proportionality between total fluorescence intensity and length (17). The absolute number of base pairs was determined by first calculating the relative lengths of the restriction fragments and then multiplying it by the known length of the DNA. In principle, for well controlled nanochannels, the length of the DNA can be determined directly from the observed length of the whole DNA and the known channel dimensions.
Fig. 4 is a histogram of calculated cut positions of SmaI on
-DNA, constructed from all molecules with two or three observed cuts. There are three known restriction sites, and the positions derived from the histogram correspond to the positions in the sequence (Table 1). We believe that the incomplete digestion is due to the staining of the DNA with an intercalating dye, which is known to interfere with sequence-recognizing enzymes (18).
Furthermore, the recognition sequence for SmaI contains a known preferred binding site of the intercalating dye TOTO-1 (19). We find that Fig. 4 shows no signs of photoinduced or unspecific enzymatic cutting. Note also that cuts by SmaI on
-DNA appear within, at most, a few seconds (Fig. 3).
|
We noticed that some of the experimental values for the fragment lengths in Table 2 differ from the sequence data by more than the expected error of the measurement (the error of the mean was
0.3 kbp in this experiment). In particular, the 8.6-kbp fragment appears systematically underestimated, whereas the 12.2-kbp fragment was systematically overestimated. This misestimation probably stems from the combination of optical system response and the functional relationship used for fitting. In particular, the imaging system does not have a purely Gaussian response, but it also a slower-decaying component. The slowly decaying component causes a differently shaped molecule fluorescence profile depending on the length of the fragment. Shoulders that do not follow the error function model become particularly prominent for long molecules. These shoulders then lead to a greater apparent width of the longer restriction fragments and reduced apparent length for shorter molecules. Note also that for the same reason, the uncertainty for the position of the third cut is larger than that of the others.
|
|
-DNA in 100-nm-wide channels relaxes within 1 s or so. We thus expect the DNA to be partially or fully thermally relaxed at the time the cut occurs, and we expect no significant recoil of the ends due to tension along the DNA backbone. The gaps must thus arise mainly from diffusive motion of the fragments. The average relative diffusion coefficient of two fragments resulting from a cut using SacI was estimated from a fit to the graph of mean squared distance between the centers of mass versus time. We find a constant of 0.5 µm2/s for the 22- to 26-kbp-long fragments. This value has to be taken as an upper limit, because only molecules in which both segments did not touch after the initial cutting, and to which good curve fitting was possible for >20 s, were included. This restriction of data created average separations that overestimate the true ensemble. It appears reasonable to assume that the diffusion constant of DNA in nanochannels scales roughly linear with the length of the molecule. We thus anticipate that gap formation due to diffusion will be observable for 100- or 200-kbp-long DNA molecules within practicable times. We have not observed a passing-through of the long segments through each other, or of the short segments through the long one, but cannot exclude that the short segments exchanged their positions. We have explained the mechanism of stretching in the papers by Tegenfeldt et al. (4) and Reisner et al. (9), and in summary the stretching is reliant on semiflexibility, entropic effects, and self-exclusion. The self-exclusion places a large energy penalty on placing two polymer strands next to each other, and we expect that interpenetration by Brownian motion is unlikely. We note that the method for optical restriction mapping is also expected to be applicable to molecules beyond that length for the following reasons. The relaxation time of DNA confined to nanochannels is expected to scale with the square of the molecule length (9). Hence, very long molecules, which typically enter the nanochannels in an overstretched state, will not be thermally relaxed at the time of cutting. The resulting restriction fragments themselves will have a relaxation time that is only dependent on their own length and thus will relax much faster, opening visible gaps in the DNA and enabling observation of the cutting. We thus anticipate that the technique will be universally applicable, as long as the relaxation times of the fragments make time averaging feasible.
We return at this point to a comparison of our method to that of Schwartz et al. (5, 6). These authors have achieved a resolution of 1.1 kbp using 20-150 molecules (20). Alternatively, the error was given as 2.9%. We thus can see that the resolution attained by using a single molecule in our method is comparable. Note also that the error in our method is expected to scale with the square root of the fragment size. Thus, the technique should perform significantly better on long fragments. Note that a significant amount of noise in our device may stem from nonuniform channels, which can be reduced (4). Two key advantages of our technique are the easy integration into microfluidic systems able of handling microliters of fluid and the capability of single-molecule mapping. Isolated and purified DNA samples can be analyzed within 1 h, including chip wetting, sample loading, and mounting in a microscope.
In conclusion, we have demonstrated a method for the construction of ordered restriction maps using single DNA molecules stretched in nanofluidic channels, and integrated control of the enzymatic cofactor magnesium. We have used the method for the mapping of
-DNA and have shown a resolution of 1.5 kbp within 1 min. We believe the presented technology could be used to study the site-selective interaction between DNA and a range of metal-ion-induced enzymes.
| Acknowledgements |
|---|
| Footnotes |
|---|
Freely available online through the PNAS open access option.
We have observed that Mg2+ concentrations of 2 mM or more lead to a quenching of fluorescence from bis-intercalating DNA dyes. We believe this is due to removal of dye from the DNA. Thus, at sufficiently high Mg2+ concentrations, intercalating dyes such as TOTO-1 or Sybr Green do not interfere with sequence-specific enzymes, because the dyes are not effectively blocking the interaction anymore. ![]()
To whom correspondence should be addressed. E-mail: rha{at}suiling.princeton.edu.
© 2005 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
K. Jo, D. M. Dhingra, T. Odijk, J. J. de Pablo, M. D. Graham, R. Runnheim, D. Forrest, and D. C. Schwartz A single-molecule barcoding system using nanoslits for DNA analysis PNAS, February 20, 2007; 104(8): 2673 - 2678. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. T. Mannion, C. H. Reccius, J. D. Cross, and H. G. Craighead Conformational Analysis of Single DNA Molecules Undergoing Entropically Induced Motion in Nanochannels Biophys. J., June 15, 2006; 90(12): 4538 - 4545. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||