Measurements of DNA barcode label separations in nanochannels from time-series data

Single-molecule optical genome mapping in nanochannels: multidisciplinarity at the nanoscale

The human genome contains multiple layers of information that extend beyond the genetic sequence. In fact, identical genetics do not necessarily yield identical phenotypes as evident for the case of two different cell types in the human body. The great variation in structure and function displayed by cells with identical genetic background is attributed to additional genomic information content. This includes large-scale genetic aberrations, as well as diverse epigenetic patterns that are crucial for regulating specific cell functions. These genetic and epigenetic patterns operate in concert in order to maintain specific cellular functions in health and disease. Single-molecule optical genome mapping is a high-throughput genome analysis method that is based on imaging long chromosomal fragments stretched in nanochannel arrays. The access to long DNA molecules coupled with fluorescent tagging of various genomic information presents a unique opportunity to study genetic and epigenetic patterns in the genome at a single-molecule level over large genomic distances. Optical mapping entwines synergistically chemical, physical, and computational advancements, to uncover invaluable biological insights, inaccessible by sequencing technologies. Here we describe the method's basic principles of operation, and review the various available mechanisms to fluorescently tag genomic information. We present some of the recent biological and clinical impact enabled by optical mapping and present recent approaches for increasing the method's resolution and accuracy. Finally, we discuss how multiple layers of genomic information may be mapped simultaneously on the same DNA molecule, thus paving the way for characterizing multiple genomic observables on individual DNA molecules.

On-Chip Stretching, Sorting, and Electro-Optical Nanopore Sensing of Ultralong Human Genomic DNA

Solid-state nanopore sensing of ultralong genomic DNA molecules has remained challenging, as the DNA must be controllably delivered by its leading end for efficient entry into the nanopore. Herein, we introduce a nanopore sensor device designed for electro-optical detection and sorting of ultralong (300+ kilobase pair) genomic DNA. The fluidic device, fabricated in-silicon and anodically bonded to glass, uses pressure-induced flow and an embedded pillar array for controllable DNA stretching and delivery. Extremely low concentrations (50 fM) and sample volumes (∼1 μL) of DNA can be processed. The low height profile of the device permits high numerical aperture, high magnification imaging of DNA molecules, which remain in focus over extended distances. We demonstrate selective DNA sorting based on sequence-specific nick translation labeling and imaging at high camera frame rates. Nanopores are fabricated directly in the assembled device by laser etching. We show that uncoiling and stretching of the ultralong DNA molecules permits efficient nanopore capture and threading, which is simultaneously and synchronously imaged and electrically measured. Furthermore, our technique provides key insights into the translocation behavior of ultralong DNA and promotes the development of all-in-one micro/nanofluidic platforms for nanopore sensing of biomolecules.

Extension distribution for DNA confined in a nanochannel near the Odijk regime

DNA confinement in a nanochannel typically is understood via mapping to the confinement of an equivalent neutral polymer by hard walls. This model has proven to be effective for confinement in relatively large channels where hairpin formation is frequent. An analysis of existing experimental data for Escherichia coli DNA extension in channels smaller than the persistence length, combined with an additional dataset for λ-DNA confined in a 34 nm wide channel, reveals a breakdown in this approach as the channel size approaches the Odijk regime of strong confinement. In particular, the predicted extension distribution obtained from the asymptotic solution to the weakly correlated telegraph model for a confined wormlike chain deviates significantly from the experimental distribution obtained for DNA confinement in the 34 nm channel, and the discrepancy cannot be resolved by treating the alignment fluctuations or the effective channel size as fitting parameters. We posit that the DNA-wall electrostatic interactions, which are sensible throughout a significant fraction of the channel cross section in the Odijk regime, are the source of the disagreement between theory and experiment. Dimensional analysis of the wormlike chain propagator in channel confinement reveals the importance of a dimensionless parameter, reflecting the magnitude of the DNA-wall electrostatic interactions relative to thermal energy, which has not been considered explicitly in the prevailing theories for DNA confinement in a nanochannel.

Measuring the wall depletion length of nanoconfined DNA

Efforts to study the polymer physics of DNA confined in nanochannels have been stymied by a lack of consensus regarding its wall depletion length. We have measured this quantity in 38 nm wide, square silicon dioxide nanochannels for five different ionic strengths between 15 mM and 75 mM. Experiments used the Bionano Genomics Irys platform for massively parallel data acquisition, attenuating the effect of the sequence-dependent persistence length and finite-length effects by using nick-labeled E. coli genomic DNA with contour length separations of at least 30 µm (88 325 base pairs) between nick pairs. Over 5 × 106 measurements of the fractional extension were obtained from 39 291 labeled DNA molecules. Analyzing the stretching via Odijk's theory for a strongly confined wormlike chain yielded a linear relationship between the depletion length and the Debye length. This simple linear fit to the experimental data exhibits the same qualitative trend as previously defined analytical models for the depletion length but now quantitatively captures the experimental data.

Distribution of label spacings for genome mapping in nanochannels

In genome mapping experiments, long DNA molecules are stretched by confining them to very narrow channels, so that the locations of sequence-specific fluorescent labels along the channel axis provide large-scale genomic information. It is difficult, however, to make the channels narrow enough so that the DNA molecule is fully stretched. In practice, its conformations may form hairpins that change the spacings between internal segments of the DNA molecule, and thus the label locations along the channel axis. Here, we describe a theory for the distribution of label spacings that explains the heavy tails observed in distributions of label spacings in genome mapping experiments.

Evaluation of Blob Theory for the Diffusion of DNA in Nanochannels

We have measured the diffusivity of λ-DNA molecules in approximately square nanochannels with effective sizes ranging from 117 nm to 260 nm at moderate ionic strength. The experimental results do not agree with the non-draining scaling predicted by blob theory. Rather, the data are consistent with the predictions of previous simulations of the Kirkwood diffusivity of a discrete wormlike chain model, without the need for any fitting parameters.

Odijk excluded volume interactions during the unfolding of DNA confined in a nanochannel

We report experimental data on the unfolding of human and E. coli genomic DNA molecules shortly after injection into a 45 nm nanochannel. The unfolding dynamics are deterministic, consistent with previous experiments and modeling in larger channels, and do not depend on the biological origin of the DNA. The measured entropic unfolding force per friction per unit contour length agrees with that predicted by combining the Odijk excluded volume with numerical calculations of the Kirkwood diffusivity of confined DNA. The time scale emerging from our analysis has implications for genome mapping in nanochannels, especially as the technology moves towards longer DNA, by setting a lower bound for the delay time before making a measurement.

Electrophoretic stretching and imaging of single native chromatin fibers in nanoslits

Stretching single chromosomal DNA fibers in nanofluidic devices has become a valuable tool for studying the genome and more recently the epigenome. Although nanofluidic technology has been extensively used in single molecular DNA analysis, compared to bare DNA, much less work has been done to elongate chromatin, and only a few studies utilize more biologically relevant samples such as native eukaryotic chromatin. Here, we provide a method for stretching and imaging individual chromatin fibers within a micro- and nanofluidic device. This device was used to electrophoretically stretch and image single native chromatin fibers extracted from human cancer cells (HeLa cells) by attaching the chromatin to microspheres held at the entrance of a nanoslit. To further demonstrate the potential of this device in epigenetics, histone modification H3k79me2 was optically detected by fluorescence microscopy.

Wall depletion length of a channel-confined polymer

Numerous experiments have taken advantage of DNA as a model system to test theories for a channel-confined polymer. A tacit assumption in analyzing these data is the existence of a well-defined depletion length characterizing DNA-wall interactions such that the experimental system (a polyelectrolyte in a channel with charged walls) can be mapped to the theoretical model (a neutral polymer with hard walls). We test this assumption using pruned-enriched Rosenbluth method (PERM) simulations of a DNA-like semiflexible polymer confined in a tube. The polymer-wall interactions are modeled by augmenting a hard wall interaction with an exponentially decaying, repulsive soft potential. The free energy, mean span, and variance in the mean span obtained in the presence of a soft wall potential are compared to equivalent simulations in the absence of the soft wall potential to determine the depletion length. We find that the mean span and variance about the mean span have the same depletion length for all soft potentials we tested. In contrast, the depletion length for the confinement free energy approaches that for the mean span only when depletion length no longer depends on channel size. The results have implications for the interpretation of DNA confinement experiments under low ionic strengths.

Optical DNA mapping in nanofluidic devices: principles and applications

Optical DNA mapping has over the last decade emerged as a very powerful tool for obtaining long range sequence information from single DNA molecules. In optical DNA mapping, intact large single DNA molecules are labeled, stretched out, and imaged using a fluorescence microscope. This means that sequence information ranging over hundreds of kilobasepairs (kbp) can be obtained in one single image. Nanochannels offer homogeneous and efficient stretching of DNA that is crucial to maximize the information that can be obtained from optical DNA maps. In this review, we highlight progress in the field of optical DNA mapping in nanochannels. We discuss the different protocols for sequence specific labeling and divide them into two main categories, enzymatic labeling and affinity-based labeling. Examples are highlighted where optical DNA mapping is used to gain information on length scales that would be inaccessible with traditional techniques. Enzymatic labeling has been commercialized and is mainly used in human genetics and assembly of complex genomes, while the affinity-based methods have primarily been applied in bacteriology, for example for rapid analysis of plasmids encoding antibiotic resistance. Next, we highlight how the design of nanofluidic channels can been altered in order to obtain the desired information and discuss how recent advances in the field make it possible to retrieve information beyond DNA sequence. In the outlook section, we discuss future directions of optical DNA mapping, such as fully integrated devices and portable microscopes.

Super-Resolution Genome Mapping in Silicon Nanochannels

Optical genome mapping in nanochannels is a powerful genetic analysis method, complementary to deoxyribonucleic acid (DNA) sequencing. The method is based on detecting a pattern of fluorescent labels attached along individual DNA molecules. When such molecules are extended in nanochannels, the labels create a fluorescent genetic barcode that is used for mapping the DNA molecule to its genomic locus and identifying large-scale variation from the genome reference. Mapping resolution is currently limited by two main factors: the optical diffraction limit and the thermal fluctuations of DNA molecules suspended in the nanochannels. Here, we utilize single-molecule tracking and super-resolution localization in order to improve the mapping accuracy and resolving power of this genome mapping technique and achieve a 15-fold increase in resolving power compared to currently practiced methods. We took advantage of a naturally occurring genetic repeat array and labeled each repeat with custom-designed Trolox conjugated fluorophores for enhanced photostability. This model system allowed us to acquire extremely long image sequences of the equally spaced fluorescent markers along DNA molecules, enabling detailed characterization of nanoconfined DNA dynamics and quantitative comparison to the Odijk theory for confined polymer chains. We present a simple method to overcome the thermal fluctuations in the nanochannels and exploit single-step photobleaching to resolve subdiffraction spaced fluorescent markers along fluctuating DNA molecules with ∼100 bp resolution. In addition, we show how time-averaging over just ∼50 frames of 40 ms enhances mapping accuracy, improves mapping P-value scores by 3 orders of magnitude compared to nonaveraged alignment, and provides a significant advantage for analyzing structural variations between DNA molecules with similar sequence composition.

Dynamic simulations show repeated narrowing maximizes DNA linearization in elastomeric nanochannels

This paper uses computer simulations to reveal unprecedented details about linearization of deoxyribonucleic acid (DNA) inside dynamic nanochannels that can be repeatedly widened and narrowed. We first analyze the effect of rate of channel narrowing on DNA linearization dynamics. Quick (∼0.1 s) narrowing of nanoscale channels results in rapid overstretching of the semi-flexible chain followed by a slower (∼0.1-10 s) relaxation to an equilibrium extension. Two phenomena that induce linearization during channel narrowing, namely, elongational-flow and confinement, occur simultaneously, regardless of narrowing speed. Interestingly, although elongational flow is a minimum at the mid-point of the channel and increases towards the two ends, neither the linearization dynamics nor the degree of DNA extension varies significantly with the center-of-mass of the polymer projected on the channel axis. We also noticed that there was a significant difference in time to reach the equilibrium length, as well as the degree of DNA linearization at short times, depending on the initial conformation of the biopolymer. Based on these observations, we tested a novel linearization protocol where the channels are narrowed and widened repeatedly, allowing DNA to explore multiple conformations. Repeated narrowing and widening, something uniquely enabled by the elastomeric nanochannels, significantly decrease the time to reach the equilibrium-level of stretch when performed within periods comparable to the chain relaxation time and more effectively untangle chains into more linearized biopolymers.

Modeling the relaxation of internal DNA segments during genome mapping in nanochannels

We have developed a multi-scale model describing the dynamics of internal segments of DNA in nanochannels used for genome mapping. In addition to the channel geometry, the model takes as its inputs the DNA properties in free solution (persistence length, effective width, molecular weight, and segmental hydrodynamic radius) and buffer properties (temperature and viscosity). Using pruned-enriched Rosenbluth simulations of a discrete wormlike chain model with circa 10 base pair resolution and a numerical solution for the hydrodynamic interactions in confinement, we convert these experimentally available inputs into the necessary parameters for a one-dimensional, Rouse-like model of the confined chain. The resulting coarse-grained model resolves the DNA at a length scale of approximately 6 kilobase pairs in the absence of any global hairpin folds, and is readily studied using a normal-mode analysis or Brownian dynamics simulations. The Rouse-like model successfully reproduces both the trends and order of magnitude of the relaxation time of the distance between labeled segments of DNA obtained in experiments. The model also provides insights that are not readily accessible from experiments, such as the role of the molecular weight of the DNA and location of the labeled segments that impact the statistical models used to construct genome maps from data acquired in nanochannels. The multi-scale approach used here, while focused towards a technologically relevant scenario, is readily adapted to other channel sizes and polymers.

Measurements of DNA barcode label separations in nanochannels from time-series data

We analyzed time-series data for fluctuations of intramolecular segments of barcoded E. coli genomic DNA molecules confined in nanochannels with sizes near the persistence length of DNA. These dynamic data allowed us to measure the probability distribution governing the distance between labels on the DNA backbone, which is a key input into the alignment methods used for genome mapping in nanochannels. Importantly, this dynamic method does not require alignment of the barcode to the reference genome, thereby removing a source of potential systematic error in a previous study of this type. The results thus obtained support previous evidence for a left-skewed probability density for the distance between labels, albeit at a lower magnitude of skewness. We further show that the majority of large fluctuations between labels are short-lived events, which sheds further light upon the success of the linearized DNA genome mapping technique. This time-resolved data analysis will improve existing genome map alignment algorithms, and the overall idea of using dynamic data could potentially improve the accuracy of genome mapping, especially for complex heterogeneous samples such as cancer cells.