纳米孔测序信号处理及其在DNA数据存储的应用

随着高通量测序技术的不断更新,可以在单个分子水平读取核苷酸序列的第三代测序技术迅速发展,纳米孔测序技术是其具有代表性的单分子测序技术,该技术通过检测DNA单链分子穿过纳米孔时引起的跨膜电流信号的变化,实现碱基识别.纳米孔测序仪在便携性、碱基读取速度、测序读段长度等方面较传统的第一代与第二代测序技术都有明显优势.随着纳米...     

单分子纳米孔测序技术及其应用研究进展

测序技术在通量和成本方面有了较大的改进,以单分子纳米孔测序技术为代表的第三代测序技术更是以其超长读长、实时检测和可以直接检测碱基甲基化修饰等优势在医学及生命科学等领域作出了较大贡献。文中就单分子纳米孔测序技术的原理进行了简要描述,并对其在临床、动物、植物、细菌及病毒等领域的应用和其未来的发展方向进行了讨论。     

Long-read whole-genome methylation patterning using enzymatic base conversion and nanopore sequencing

Long-read whole-genome sequencing analysis of DNA methylation would provide useful information on the chromosomal context of gene expression regulation. Here we describe the development of a method that improves the read length generated by using the bisulfite-sequencing-based approach. In this method, we combined recently developed enzymatic base conversion, where an unmethylated cytosine (C) should be converted to thymine (T), with nanopore sequencing. After methylation-sensitive base conversion, the sequencing library was constructed using long-range polymerase chain reaction. This type of analysis is possible using a minimum of 1 ng genomic DNA, and an N50 read length of 3.4–7.6 kb is achieved. To analyze the produced data, which contained a substantial number of base mismatches due to sequence conversion and an inaccurate base read of the nanopore sequencing, a new analytical pipeline was constructed. To demonstrate the performance of long-read methylation sequencing, breast cancer cell lines and clinical specimens were subjected to analysis, which revealed the chromosomal methylation context of key cancer-related genes, allele-specific methylated genes, and repetitive or deletion regions. This method should convert the intractable specimens for which the amount of available genomic DNA is limited to the tractable targets.     

Osmium-Based Pyrimidine Contrast Tags for Enhanced Nanopore-Based DNA Base Discrimination

Nanopores are a promising platform in next generation DNA sequencing. In this platform, an individual DNA strand is threaded into nanopore using an electric field, and enzyme-based ratcheting is used to move the strand through the detector. During this process the residual ion current through the pore is measured, which exhibits unique levels for different base combinations inside the pore. While this approach has shown great promise, accuracy is not optimal because the four bases are chemically comparable to one another, leading to small differences in current obstruction. Nucleobase-specific chemical tagging can be a viable approach to enhancing the contrast between different bases in the sequence. Herein we show that covalent modification of one or both of the pyrimidine bases by an osmium bipyridine complex leads to measureable differences in the blockade amplitudes of DNA molecules. We qualitatively determine the degree of osmylation of a DNA strand by passing it through a solid-state nanopore, and are thus able to gauge T and C base content. In addition, we show that osmium bipyridine reacts with dsDNA, leading to substantially different current blockade levels than exhibited for bare dsDNA. This work serves as a proof of principle for nanopore sequencing and mapping via base-specific DNA osmylation.     

Discrimination among individual Watson-Crick base pairs at the termini of single DNA hairpin molecules

Nanoscale α-hemolysin pores can be used to analyze individual DNA or RNA molecules. Serial examination of hundreds to thousands of molecules per minute is possible using ionic current impedance as the measured property. In a recent report, we showed that a nanopore device coupled with machine learning algorithms could automatically discriminate among the four combinations of Watson–Crick base pairs and their orientations at the ends of individual DNA hairpin molecules. Here we use kinetic analysis to demonstrate that ionic current signatures caused by these hairpin molecules depend on the number of hydrogen bonds within the terminal base pair, stacking between the terminal base pair and its nearest neighbor, and 5′ versus 3′ orientation of the terminal bases independent of their nearest neighbors. This report constitutes evidence that single Watson–Crick base pairs can be identified within individual unmodified DNA hairpin molecules based on their dynamic behavior in a nanoscale pore.     

Specific contacts between the FLP protein of the yeast 2-micron plasmid and its recombination site

Contact points between the FLP protein of the yeast 2-micron plasmid and its recombination site have been defined. Important features of the region previously defined as the minimal recombination site in vitro include a pair of 13-base pair inverted repeats separated by an 8-base pair spacer. The two FLP protein-binding sites within this region are 12 base pairs in length. In each case they include the internal 11 base pairs of one of the 13-base pair repeats, as well as the adjacent base pair within the spacer. The internal 6 base pairs within the spacer are not involved in binding or recognition by FLP protein. When the size of the spacer is increased or decreased by one base pair, the distance between the cleavage points is also increased or decreased correspondingly by one base pair. Points of cleavage are unaffected by changes in the spacer sequence. Specific contact points involving purine residues, identified by methylation protection and recombination interference experiments, are located in both the major and minor grooves of the DNA. Additional contact points between FLP protein and phosphate groups in the phosphate-deoxyribose backbone are clustered near the cleavage sites.     

Automated forward and reverse ratcheting of DNA in a nanopore at 5-Å precision

An emerging DNA sequencing technique uses protein or solid-state pores to analyze individual strands as they are driven in single-file order past a nanoscale sensor. However, uncontrolled electrophoresis of DNA through these nanopores is too fast for accurate base reads. Here, we describe forward and reverse ratcheting of DNA templates through the α-hemolysin nanopore controlled by phi29 DNA polymerase without the need for active voltage control. DNA strands were ratcheted through the pore at median rates of 2.5-40 nucleotides per second and were examined at one nucleotide spatial precision in real time. Up to 500 molecules were processed at ～130 molecules per hour through one pore. The probability of a registry error (an insertion or deletion) at individual positions during one pass along the template strand ranged from 10% to 24.5% without optimization. This strategy facilitates multiple reads of individual strands and is transferable to other nanopore devices for implementation of DNA sequence analysis.     

A parallel graph decomposition algorithm for DNA sequencing with nanopores

MOTIVATION: With the potential availability of nanopore devices that can sense the bases of translocating single-stranded DNA (ssDNA), it is likely that 'reads' of length approximately 10(5) will be available in large numbers and at high speed. We address the problem of complete DNA sequencing using such reads.We assume that approximately 10(2) copies of a DNA sequence are split into single strands that break into randomly sized pieces as they translocate the nanopore in arbitrary orientations. The nanopore senses and reports each individual base that passes through, but all information about orientation and complementarity of the ssDNA subsequences is lost. Random errors (both biological and transduction) in the reads create further complications. RESULTS: We have developed an algorithm that addresses these issues. It can be considered an extreme variation of the well-known Eulerian path approach. It searches over a space of de Bruijn graphs until it finds one in which (a) the impact of errors is eliminated and (b) both possible orientations of the two ssDNA sequences can be identified separately and unambiguously.Our algorithm is able to correctly reconstruct real DNA sequences of the order of 10(6) bases (e.g. the bacterium Mycoplasma pneumoniae) from simulated erroneous reads on a modest workstation in about 1 h. We describe, and give measured timings of, a parallel implementation of this algorithm on the Cray Multithreaded Architecture (MTA-2) supercomputer, whose architecture is ideally suited to this 'unstructured' problem. Our parallel implementation is crucial to the problem of rapidly sequencing long DNA sequences and also to the situation where multiple nanopores are used to obtain a high-bandwidth stream of reads.     

The impact of the number of layers of the graphene nanopores and the electrical field on ssDNA translocation

Graphene-based nanopore devices hold great promise for the next generation DNA sequencing because graphene is atomically thin which is extremely important for single base recognition. To understand the fundamental details of DNA translocation through a graphene nanopore, in this work, molecular dynamics simulations of ssDNA translocation through the nanopore were performed to trace the nucleobase trajectories and to investigate the impact of the number of layers of the graphene membrane and the electrical field on ssDNA translocation. We found that the velocity of ssDNA translocation was speeded up with the higher bias voltage, and the two-layered and five-layered graphene membrane with 1.0-nm diameter circular nanopore could discern different DNA strand by the translocation time.     

DNA base pair hybridization and water-mediated metastable structures studied by molecular dynamics simulations

The base pair hybridization of a DNA segment was studied using molecular dynamics simulation. The results show the obvious correlation between the probability of successful hybridization and the accessible surface area to water of two successive base pairs, including the unpaired base pair adjacent to paired base pair and this adjacent paired base pair. Importantly, two metastable structures in an A-T base pair were discovered by the analysis of the free energy landscape. Both structures involved addition of a water molecule to the linkage between the two nucleobases in one base pair. The existence of the metastable structures provide potential barriers to the Watson-Crick base pair, and numerical simulations show that those potential barriers can be surmounted by thermal fluctuations at higher temperatures. These studies contribute an important step toward the understanding of the mechanism in DNA hybridization, particularly the effect of temperature on DNA hybridization and polymerase chain reaction. These observations are expected to be helpful for facilitating experimental bio/nanotechnology designs involving fast hybridization.     

Structure of the DNA interstrand cross-link of 4,5',8-trimethylpsoralen.

4,5',8-Trimethylpsoralen (TMP) cross-links a 5' TpA or a 5' ApT site by photoreacting with one thymine moiety in each DNA strand. We are interested in whether psoralen interstrand cross-links all share one structure or whether there are significant differences. In this paper, we employed a rapid method for probing the structure of the cross-link by making a series of TMP cross-linked duplexes containing specific base-pair mismatches. The relative stability provided by a base pair can be correlated with neighboring base pairs by comparing the extents of gel retardation when base-pair mismatches happen in each position. From our studies, we infer that with respect to the furan-side strand, the 5'T.A base pair of the two T.A base pairs in the TpA site is not hydrogen bonded. Immediately on each side of the cross-linked TpA site is a highly stabilized base pair. Next, a region of decreased stability occurs in each arm of a cross-linked duplex and these base pairs of least stability are located farther away from the cross-linked thymines as the lengths of the arms of the cross-linked helix increase. Finally, even in 7 M urea at 49 degrees C the cross-linked helix is hydrogen bonded at both ends of a duplex of 22 base pairs. We propose that the structures of interstrand cross-links in DNA vary appreciably with the DNA sequence, the length of the DNA duplex, and the structures of the DNA cross-linking agents.     

The flexibility of A-form DNA.

We have determined the rise per base pair and persistence length of A-form DNA in trifluoroethanol solutions for fragments 350-900 base pairs in length that best describe rotational diffusion coefficients determined by transient electric birefringence. The 2.6 A spacing between base pairs found in crystal and fiber A-form structures is preserved in solution. The persistence length is about 1500 A, or about three times longer than for B-form DNA. There is no apparent electrostatic contribution to the persistence length in the salt concentration range 0.2-2.0 mM Na cacodylate. This suggests an even closer association between DNA and its neutralizing counterions than predicted by condensation theory, perhaps due to a sheath of trifluoroethanol excluded water surrounding the A-form helix.     

The Flexibility of A-Form DNA

Abstract

We have determined the rise per base pair and persistence length of A-form DNA in trifluoroethanol solutions for fragments 350–900 base pairs in length that best describe rotational diffusion coefficients determined by transient electric birefringence. The 2.6 A spacing between base pairs found in crystal and fiber A-form structures is preserved in solution. The persistence length is about 1500 A, or about three times longer than for B-form DNA. There is no apparent electrostatic contribution to the persistence length in the salt concentration range 0.2–2.0 mM Na cacodylate. This suggests an even closer association between DNA and its neutralizing counterions than predicted by condensation theory, perhaps due to a sheath of trifluoroethanol excluded water surrounding the A-form helix.     

纳米孔测序技术应用于食品微生物检测的可行性分析

近些年来DNA测序技术发展迅速,已经从第一代生化测序发展到第三代单分子测序。作为第三代测序技术中的一种不同于当前流行的其他测序技术,纳米孔测序技术是基于电信号的一种物理方法测序。许多研究者通常将高通量测序技术应用于食品微生物的研究,但是将纳米孔测序技术应用于食品中微生物的检测却鲜有报道。Oxford Nanopore Technologies(牛津纳米孔科技公司)研发的DNA测序仪MinION,是世界首例用于商业测序的纳米孔测序仪,经过不断完善,近年来MinION在DNA测序中被广泛应用。MinION 测序一次需要的DNA量约1μg,其标准识别速度为一秒钟识别250个碱基,平均读长可至13kb~20kb,测序准确率可以达到98%。纳米孔测序的高识别速度和高准确率,完全满足快速检测的要求,将其应用于食品中微生物检测是完全可行的。     

207 Nanopore immobilization of DNA polymerase enhances single-molecule sequencing

A zero-mode waveguide (ZMW) is a nanoscale optical waveguide driven at a frequency below its cut-off. In this mode, the electric field, instead of traveling down the axis of the conducting cavity, decays exponentially. By fabricating waveguides with sub-wavelength diameters and illuminating them with laser light, the electric field in the waveguide is confined enough to enable single-molecule optical detection at micromolar concentration [1]. Immobilizing single DNA polymerases in ZMWs and using special phosphate-fluorescently labeled dNTPs form the basis for single-molecule real-time DNA sequencing, one of the most promising next-generation sequencing platforms [2]. In this method, the polymerase replicates the sample DNA, and as it incorporates new bases into the product strand, the labeled dNTPs emit a burst of light before the phosphate is cleaved off. The sequence of colors corresponds to the DNA sequence (see Figure 1 below from Eid et al., 2009). Because the ZMW aperture’s diameter is sub-diffraction-limit, it is impossible to optically distinguish one polymerase in a ZMW from two. Having only one polymerase in each waveguide is critical to sequencing accuracy. In its present state, experimenters use diffusion to fill ZMWs with polymerases, resulting in a Poisson distribution for filling ZMWs, and consequently a theoretical limit of 36.8% of ZMWs having only one polymerase [2]. We achieve full polymerase occupancy of ZMWs by fabricating the structures on an ultrathin silicon nitride membrane and drilling a nanopore at the base of each waveguide with an ion beam. A short DNA fragment with biotin on either end is conjugated to a streptavidin and then drawn into the nanopore with a voltage bias. There is then a free biotin at the base of the ZMW. A polymerase–streptavidin complex can diffuse into the ZMW and bind to the exposed biotin. Because the nanopore is too small to fit more than one molecule, only one ZMW will bind to a biotin in the nanopore. Upon flushing the ZMW chamber, the biotin-bound polymerase will remain trapped in the pore, and only a single polymerase will remain at the base of each waveguide.      

Investigation of DNA dynamics and drug-DNA interaction by steady state fluorescence anisotropy.

We have used steady-state fluorescence polarization anisotropy (FPA) of ethidium probe molecules bound to DNA to investigate DNA-DNA interactions and the effect of high densities of intercalating drugs on the internal motions of DNA responsible for depolarization of the ethidium fluorescence. To calibrate the method, we examined the effect of DNA length on (FPA) using DNA varying in size from 10-150 base pair. The association of approximately 30 base pair DNA at high concentrations was then detected by its effect on (FPA). With sample concentrations approaching those commonly used in various physical experiments (NMR, Raman) significant DNA-DNA interactions are observed. With high molecular weight DNA (greater than 500 base pair), the limiting value of the (FPA) (0.23) is due to internal motions of the DNA (and bound chromophores). The (FPA) of ethidium probe molecules (1 drug/200 base pair) is unaffected by the addition of high levels (1 drug/2 base pair) proflavine. This indicates that either the elastic properties of DNA are unaffected by high densities of intercalated drug or that the depolarization of the ethidium fluorescence is due to highly localized motions of the base pairs that are unperturbed by binding of drugs at neighboring sites.     

S-DNA is oversupercoiled DNA with 1.94-2.19 A rise per base pair along duplex axis

Supercoiled pGEMEX DNA with length of 3993 nucleotides was immobilized on four substrates (freshly cleaved mica, standard amino mica, modified amino mica with increased and decreased surface charge density compared with standard amino mica) and it was visualized by atomic force microscopy (AFM) in air. Plectonomically supercoiled DNA molecules as well as single molecules with extremely high level of compaction (i.e. molecules with significantly higher superhelix density values on comparison with previously experimentally measured and theoretically investigated ones) were visualized on modified amino mica which was characterized by increased surface charge density. Distance between base pairs along duplex axis was determined by measurements of contour length of single oversupercoiled DNA molecules. Determined rise per base pair was varied from 1.94 to 2.19 A. These compressed supercoiled DNA molecules like a spring with decreased rise/base pair on comparison with well-known DNA forms were called new DNA form--S-DNA. A model of S-DNA was built. Formation of the S-DNA molecules was suggested to be an intermediate stage on the compaction of the single supercoiled DNA molecules up to the spheroids and minitoroids. Oversupercoiling and further compression of the supercoiled DNA molecules was shown to cause by high surface charge density of amino mica on which DNA molecules were immobilized.     

The flexibility of alternating dA-dT sequences

The flexibility of alternating poly (dA-dT) has been investigated by the technique of transient electric dichroism. Rotational relaxation times, which are very sensitive to changes in the end-to-end length of flexible polymers, are determined from the field free dichroism decay curves of four, well defined fragments of poly (dA-dT) ranging in size from 136 to 270 base pairs. Persistence lengths, calculated from the results of Hagerman and Zimm (Biopolymers (1981) 29, 1481-1502), are in the range 200-250 A. This makes alternating dA-dT sequences about twice as flexible as naturally occurring, "random" sequence DNA. Considering a bend around a nucleosome, for example, this difference in persistence length translates to an energy difference between poly (dA-dT) and random sequence DNA of 0.17 kT/base pair or 1 kcal per 10 base pair stretch. This energy difference is sufficiently large to suggest that dA-dT sequences could serve as markers in DNA packaging, for example, at sites where DNA must tightly bend to accommodate structures.