A fast and scalable kymograph alignment algorithm for nanochannel-based optical DNA mappings

Random sampling of ligand arrangements on a one-dimensional lattice

We introduce a transfer-matrix-based sequential sampling scheme for generating random samples of ligand arrangements on one-dimensional templates. The number of ligand types is arbitrary, the binding constants can have positional dependence, and cooperativity parameters are included. From the random arrangements, any (linear or nonlinear) observable can be calculated using sample averaging. As an example case study, we investigate the competitive binding of three ligand types (the sequence-specific binder netropsin, YOYO-1, and ethidium bromide) to a DNA molecule. We also employ our random sampling method of ligands to determine the quality of synthetically generated DNA barcodes as a function of concentration of a ligand (e.g., netropsin) in optical DNA mapping (ODM) experiments. We provide publically available softwares, with a computational time that scales linearly with the lattice size, for generating random ligand arrangements and for generating synthetic barcodes.

Advances in optical mapping for genomic research

Recent advances in optical mapping have allowed the construction of improved genome assemblies with greater contiguity. Optical mapping also enables genome comparison and identification of large-scale structural variations. Association of these large-scale genomic features with biological functions is an important goal in plant and animal breeding and in medical research. Optical mapping has also been used in microbiology and still plays an important role in strain typing and epidemiological studies. Here, we review the development of optical mapping in recent decades to illustrate its importance in genomic research. We detail its applications and algorithms to show its specific advantages. Finally, we discuss the challenges required to facilitate the optimization of optical mapping and improve its future development and application.

Noise reduction in single time frame optical DNA maps

In optical DNA mapping technologies sequence-specific intensity variations (DNA barcodes) along stretched and stained DNA molecules are produced. These “fingerprints” of the underlying DNA sequence have a resolution of the order one kilobasepairs and the stretching of the DNA molecules are performed by surface adsorption or nano-channel setups. A post-processing challenge for nano-channel based methods, due to local and global random movement of the DNA molecule during imaging, is how to align different time frames in order to produce reproducible time-averaged DNA barcodes. The current solutions to this challenge are computationally rather slow. With high-throughput applications in mind, we here introduce a parameter-free method for filtering a single time frame noisy barcode (snap-shot optical map), measured in a fraction of a second. By using only a single time frame barcode we circumvent the need for post-processing alignment. We demonstrate that our method is successful at providing filtered barcodes which are less noisy and more similar to time averaged barcodes. The method is based on the application of a low-pass filter on a single noisy barcode using the width of the Point Spread Function of the system as a unique, and known, filtering parameter. We find that after applying our method, the Pearson correlation coefficient (a real number in the range from -1 to 1) between the single time-frame barcode and the time average of the aligned kymograph increases significantly, roughly by 0.2 on average. By comparing to a database of more than 3000 theoretical plasmid barcodes we show that the capabilities to identify plasmids is improved by filtering single time-frame barcodes compared to the unfiltered analogues. Since snap-shot experiments and computational time using our method both are less than a second, this study opens up for high throughput optical DNA mapping with improved reproducibility.

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.

Rapid identification of intact bacterial resistance plasmids via optical mapping of single DNA molecules

The rapid spread of antibiotic resistance – currently one of the greatest threats to human health according to WHO – is to a large extent enabled by plasmid-mediated horizontal transfer of resistance genes. Rapid identification and characterization of plasmids is thus important both for individual clinical outcomes and for epidemiological monitoring of antibiotic resistance. Toward this aim, we have developed an optical DNA mapping procedure where individual intact plasmids are elongated within nanofluidic channels and visualized through fluorescence microscopy, yielding barcodes that reflect the underlying sequence. The assay rapidly identifies plasmids through statistical comparisons with barcodes based on publicly available sequence repositories and also enables detection of structural variations. Since the assay yields holistic sequence information for individual intact plasmids, it is an ideal complement to next generation sequencing efforts which involve reassembly of sequence reads from fragmented DNA molecules. The assay should be applicable in microbiology labs around the world in applications ranging from fundamental plasmid biology to clinical epidemiology and diagnostics.

Visualizing the entire DNA from a chromosome in a single frame

The contiguity and phase of sequence information are intrinsic to obtain complete understanding of the genome and its relationship to phenotype. We report the fabrication and application of a novel nanochannel design that folds megabase lengths of genomic DNA into a systematic back-and-forth meandering path. Such meandering nanochannels enabled us to visualize the complete 5.7 Mbp (1 mm) stained DNA length of a Schizosaccharomyces pombe chromosome in a single frame of a CCD. We were able to hold the DNA in situ while implementing partial denaturation to obtain a barcode pattern that we could match to a reference map using the Poland-Scheraga model for DNA melting. The facility to compose such long linear lengths of genomic DNA in one field of view enabled us to directly visualize a repeat motif, count the repeat unit number, and chart its location in the genome by reference to unique barcode motifs found at measurable distances from the repeat. Meandering nanochannel dimensions can easily be tailored to human chromosome scales, which would enable the whole genome to be visualized in seconds.