Genomic Encryption of Digital Data Stored in Synthetic DNA

Today, we can read human genomes and store digital data robustly in synthetic DNA. Herein, we report a strategy to intertwine these two technologies to enable the secure storage of valuable information in synthetic DNA, protected with personalized keys. We show that genetic short tandem repeats (STRs) contain sufficient entropy to generate strong encryption keys, and that only one technology, DNA sequencing, is required to simultaneously read the key and the data. Using this approach, we experimentally generated 80 bit strong keys from human DNA, and used such a key to encrypt 17 kB of digital information stored in synthetic DNA. Finally, the decrypted information was recovered perfectly from a single massively parallel sequencing run.     

Large-scale pyrosequencing of synthetic DNA: a comparison with results from Sanger dideoxy sequencing

Pyrosequencing is a relatively recent method for sequencing short stretches of DNA. Because both Pyrosequencing and Sanger dideoxy sequencing were recently used to characterize and validate DNA molecular barcodes in a large yeast gene-deletion project, a meta-analysis of those data allow an excellent and timely opportunity for evaluating Pyrosequencing against the current gold standard, Sanger dideoxy sequencing. Starting with yeast genomic DNA, parallel PCR amplification methods were used to prepared 4747 short barcode-containing constructs from 6000 Saccharomyces cerevisiae gene-deletion strains. Pyrosequencing was optimized for average read lengths of 25-30 bases, which included in each case a 20-mer barcode sequence. Results were compared with sequence data obtained by the standard Sanger dideoxy chain termination method. In most cases, sequences obtained by Pyrosequencing and Sanger dideoxy sequencing were of comparable accuracy, and the overall rate of failure was similar. The DNA in the barcodes is derived from synthetic oligonucleotide sequences that were inserted into yeast-deletion-strain genomic DNA by homologous recombination and represents the most significant amount of DNA from a synthetic source that has been sequenced to date. Although more automation and quality control measures are needed, Pyrosequencing was shown to be a fast and convenient method for determining short stretches of DNA sequence.     

Copy number variants calling for single cell sequencing data by multi-constrained optimization

Variations in DNA copy number carry important information on genome evolution and regulation of DNA replication in cancer cells. The rapid development of single-cell sequencing technology allows one to explore gene expression heterogeneity among single-cells, thus providing important cancer cell evolution information. Single-cell DNA/RNA sequencing data usually have low genome coverage, which requires an extra step of amplification to accumulate enough samples. However, such amplification will introduce large bias and makes bioinformatics analysis challenging. Accurately modeling the distribution of sequencing data and effectively suppressing the bias influence is the key to success variations analysis.Recent advances demonstrate the technical noises by amplification are more likely to follow negative binomial distribution, a special case of Poisson distribution. Thus, we tackle the problem CNV detection by formulating it into a quadratic optimization problem involving two constraints, in which the underling signals are corrupted by Poisson distributed noises. By imposing the constraints of sparsity and smoothness, the reconstructed read depth signals from single-cell sequencing data are anticipated to fit the CNVs patterns more accurately. An efficient numerical solution based on the classical alternating direction minimization method (ADMM) is tailored to solve the proposed model. We demonstrate the advantages of the proposed method using both synthetic and empirical single-cell sequencing data. Our experimental results demonstrate that the proposed method achieves excellent performance and high promise of success with single-cell sequencing data.     

Identification of methylated regions with peak search based on Poisson model from massively parallel methylated DNA immunoprecipitation-sequencing data

DNA methylation is one of the most important epigenetic modification types, which plays a critical role in gene expression. High efficient surveying of whole genome DNA methylation has been aims of many researchers for long. Recently, the rapidly developed massively parallel DNA‐sequencing technologies open the floodgates to vast volumes of sequence data, enabling a paradigm shift in profiling the whole genome methylation. Here, we describe a strategy, combining methylated DNA immunoprecipitation sequencing with peak search to identify methylated regions on a whole‐genome scale. Massively parallel methylated DNA immunoprecipitation sequencing combined with methylation DNA immunoprecipitation was adopted to obtain methylated DNA sequence data from human leukemia cell line K562, and the methylated regions were identified by peak search based on Poison model. From our result, 140 958 non‐overlapping methylated regions have been identified in the whole genome. Also, the credibility of result has been proved by its strong correlation with bisulfite‐sequencing data (Pearson R²=0.92). It suggests that this method provides a reliable and high‐throughput strategy for whole genome methylation identification.     

Contaminant-induced current decline in capillary array electrophoresis

High-throughput capillary array electrophoresis (CAE) instruments for DNA sequencing suffer to varying degrees from read length degradation associated with electrophoretic current decline and inhibition or delay in the arrival of fragments at the detector. This effect is known to be associated with residual amounts of large, slow-moving fragments of template or genomic DNA carried through from sample preparation and sequencing reactions. Here, we investigate the creation and expansion of an ionic depletion region induced by overloading the capillary with low-mobility DNA fragments, and the effect of growth of this region on electrophoresis run failure. Slow-moving fragments are analytically and experimentally shown to reduce the ionic concentration of the downstream electrolyte. With injection of large fragments beyond a threshold quantity, the anode-side boundary of the nascent depletion region begins to propagate toward the anode at a rate faster than the contaminant DNA migration. Under such conditions, the depletion region expands, the capillary current declines dramatically, and the electrophoresis run yields a short read length or fails completely.     

Genomsequenzierer. Die dritte Generation

Although the living cell can read its genetic information from a single molecule of DNA, researchers have used vast numbers of identical copies in traditional sequencing methods. Dramatic improvements in sequencing technology only became available after the first sequencing of the human genome was finished and led to a rapid drop in costs by five orders of magnitude. Instruments of the third generation, now reaching the market, will for the first time read single molecules and reduce costs even further. The implications for medicine and society will be profound.     

A theoretical study of the possible use of electroosmotic flow to extend the read length of DNA sequencing by end-labeled free solution electrophoresis

End-labeled free solution electrophoresis (ELFSE) provides a means of separating DNA with free-solution CE, eliminating the need for gels and polymer solutions which increase the run time and can be difficult to load into a capillary. In free-solution electrophoresis, DNA is normally free-draining and all fragments reach the detector at the same time, whereas ELFSE uses an uncharged label molecule attached to each DNA fragment in order to render the electrophoretic mobility size-dependent. With ELFSE, however, the larger molecules are not separated enough (limiting the read length in the case of ssDNA sequencing) while the smaller ones are overseparated; the larger ones are too fast while the shorter ones are too slow, which is the opposite of traditional gel-based methods. In this article, we show how an EOF could be used to overcome these problems and extend the DNA sequencing read length of ELFSE. This counterflow would allow the larger, previously unresolved molecules more time to separate and thereby increase the read length. Through our theoretical investigation, we predict that an EOF mobility of approximately the same magnitude as that of unlabeled DNA would provide the best results for the regime where all molecules move in the same direction. Even better resolution would be possible for smaller values of EOF which allow different directions of migration; however, the migration times then would become too large. The flow would need to be well controlled since the gain in read length decreases as the magnitude of the counterflow increases; an EOF mobility double that of unlabeled DNA would no longer increase the read length, although ELFSE would still benefit from a reduction in migration time.     

Dendritic Quaternary-Encoded Oligourethanes for Data Encryption

The use of sequence-defined digital polymers for data storage and encryption has received increasing attention due to their precision structures similar to natural biomacromolecules (e.g., DNA) but increased stability. However, the rapid development of sequencing techniques raises the concern of information leakage. Herein, dendritic quaternary-encoded oligourethanes bearing a photoresponsive trigger, self-immolative backbones, and a mass spectrometry tag of PEG dendron have been developed for data encryption. Although the sequence information in linear analogs can be readily deciphered by mass spectrometry, sequencing of dendritic oligourethanes cannot be achieved by either primary MS or tandem MS/MS owing to the unique spatial conformation. Intriguingly, the fragmentation pathways of a quaternary dendrimer under MS/MS conditions can be converted to 2772-bit 2D matrices with ≈1.98×10⁸⁷ permutations, serving as high-strength encryption keys for highly reliable data encryption.     

Dissecting the differences between the alpha and beta anomers of the oxidative DNA lesion FaPydG

The oxidative DNA lesion, FaPydG rapidly anomerizes to form a mixture of the alpha and beta anomer. To investigate the mutagenic potential of both forms, we prepared stabilized bioisosteric analogues of both configurational isomers and incorporated them into oligonucleotides. These were subsequently used for thermodynamic melting-point studies and for primer-extension experiments. While the beta compound, in agreement with earlier data, prefers cytidine as the pairing partner, the alpha compound is not able form a stable base pair with any natural base. In primer-extension studies with the high-fidelity polymerase Bst Pol I, the polymerase was able to read through the lesion. The beta compound showed no strong mutagenic potential. The alpha compound, in contrast, strongly destabilized DNA duplexes and also blocked all of the tested DNA polymerases, including two low-fidelity polymerases of the Y-family.     

Information Coding in a Reconfigurable DNA Origami Domino Array

DNA nanostructures with programmable nanoscale patterns has been achieved in the past decades, and molecular information coding (MIC) on those designed nanostructures has gained increasing attention for information security. However, achieving steganography and cryptography synchronously on DNA nanostructures remains a challenge. Herein, we demonstrated MIC in a reconfigurable DNA origami domino array (DODA), which can reconfigure intrinsic patterns but keep the DODA outline the same for steganography. When a set of keys (DNA strands) are added, the cryptographic data can be translated into visible patterns within DODA. More complex cryptography with the ASCII code within a programmable 6×6 lattice is demonstrated to demosntrate the versatility of MIC in the DODA. Furthermore, an anti‐counterfeiting approach based on conformational transformation‐mediated toehold strand displacement reaction is designed to protect MIC from decoding and falsification.     

Robust Chemical Preservation of Digital Information on DNA in Silica with Error‐Correcting Codes

Information, such as text printed on paper or images projected onto microfilm, can survive for over 500 years. However, the storage of digital information for time frames exceeding 50 years is challenging. Here we show that digital information can be stored on DNA and recovered without errors for considerably longer time frames. To allow for the perfect recovery of the information, we encapsulate the DNA in an inorganic matrix, and employ error‐correcting codes to correct storage‐related errors. Specifically, we translated 83 kB of information to 4991 DNA segments, each 158 nucleotides long, which were encapsulated in silica. Accelerated aging experiments were performed to measure DNA decay kinetics, which show that data can be archived on DNA for millennia under a wide range of conditions. The original information could be recovered error free, even after treating the DNA in silica at 70 °C for one week. This is thermally equivalent to storing information on DNA in central Europe for 2000 years.     

Improvement of base-calling in multilane automated DNA sequencing by use of electrophoretic calibration standards, data linearization, and trace alignment

We present a new method for the linearization and alignment of data traces generated by multilane automated DNA sequencing instruments. Application of this method to data generated with the Visible Genetics Open Gene DNA sequencing system (using MicroCel 700 gel cassettes, with a 25 cm separation distance) allows read lengths of > 1,000 nucleotides to be routinely obtained with high confidence and > 97% accuracy. This represents an increase of 10-15% in average read length, relative to data from this system that have not been processed in the fashion described herein. Most importantly, the linearization and alignment method allows usable sequence to be obtained from a fraction of 10-15% of data sets which, because of original trace misalignment problems, would otherwise have to be discarded. Our method involves adding electrophoretic calibration standards to the DNA sequencing fragments. The calibration standards are labeled with a dye that differs spectrally from the dye attached to the sequencing fragments. The calibration standards are identical in all the lanes. Analysis of the mobilities of the calibration standards allows correction for both systematic and random variation of electrophoretic properties between gel lanes. We have successfully used this method with two-dye and three-dye DNA sequencing instruments.     

毛细管阵列电泳与规模化DNA测序

 根据 10 0 80份基因组DNA测序的结果 ,讨论了毛细管阵列电泳测序方法的技术特点 ,并对影响测序结果的一些因素进行了分析。在此基础上与平板凝胶电泳方法进行比较 ,显示了毛细管阵列电泳的优点。同时也对大规模测序技术环节之间的协调进行了探讨 。     

Length-based encoding of binary data in DNA

We developed a system to encode digital information in DNA polymers based on the partial restriction digest (PRD). Our encoding method relies on the length of the fragments obtained by the PRD rather than the actual content of the nucleotide sequence, thus eliminating the need for expensive sequencing machinery. In this letter, we report on the encoding of 12 bits of data in a DNA fragment of 110 nucleotides and the process of recovering the data.     

High-speed DNA sequencing by tube-based capillary electrophoresis

We assessed the feasibility of high-speed DNA sequencing by tube-based capillary electrophoresis (TCE) with electrokinetic sample injections. We developed a water-circulated TCE system to control the capillary temperature precisely. Using this system and a ready-made sieving matrix at 50 degrees C, single-stranded DNA size marker fragments were separated at various pairs of the electric field strength, E, of 128-480 V/cm and the capillary effective length, L, of 100-360 mm. Assuming the read length (RL) is the fragment size at which the peak width equals the peak interval per base in obtained electropherograms, we estimated the values of RL (E, L), the RL at the pair (E, L). The points in ELz-space, (E, L, RL(E, L)), form a curved surface expressed by z = RL(E, L). Analyzing the contour lines of this curved surface, we determined the pairs of E and L providing target RLs of 300-500 bases within a minimum time. At a pair optimized for a 500-base RL (330 V/cm, 200 mm), one-color sequencing fragments were successfully separated up to 529 bases within 9.6 min. These results demonstrate that high-speed DNA sequencing comparable with that obtained by microfabricated chip-based capillary electrophoresis (MCE) can be achieved with TCE, which is more suitable in automation than MCE.     

Massively parallel display of genomic DNA fragments by rolling-circle amplification and strand displacement amplification on chip

Massively parallel genomic DNA fragments display on chip plays a key role in the new generation DNA sequencing. Here, we developed a new technology to display the parallel genomic DNA fragment massively based on two-step reaction with Ф29 DNA polymerase. The genomic DNA fragments were firstly amplified by rolling-circle amplification (RCA) reaction in liquid phase, and then amplified further on the chip by the strand displacement of Ф29 DNA polymerase. In our experiments, through DNA colonies produced by two-step amplification reaction T7 genomic DNA fragments are displayed massively and parallely on the chip, which has been verified through hybridizing the probe labeled with fluorescence or extension reaction with fluorescent-dNTP. The significant difference of fluourescence signals between background and displayed DNA fragments could be obtained. Our results show that the method has good reproducibility in experiments, which may be hopeful to serve the high-throughput sequencing.     

Low viscosity DNA sieving matrices for capillary electrophoresis

The rapid development of DNA capillary electrophoresis (CE) technology has increased the demand of new low viscosity sieving matrices with high separation capacity. The high throughput, resolution and automatic operation of CE systems have stimulated the application of the technique to different kinds of DNA analysis, including DNA sequencing, separation of restriction fragments, PCR products and synthetic oligonucleotides. In addition specific methods for PCR-based mutation assays for the study of known and unknown point mutations have been developed for use in CE. The key component for a large scale application of CE to DNA analysis is the availability of appropriate sieving matrices. This article gives an overview of the linear polymers used as DNA separation matrices with particular emphasis on the polymers that combine high sieving capacity, low viscosity and chemical resistance.     

Separation efficiency of free‐solution conjugated electrophoresis with drag‐tags incorporating a synthetic amino acid

DNA sequencing or separation by conventional capillary electrophoresis with a polymer matrix has some inherent drawbacks, such as the expense of polymer matrix and limitations in sequencing read length. As DNA fragments have a linear charge‐to‐friction ratio in free solution, DNA fragments cannot be separated by size. However, size‐based separation of DNA is possible in free‐solution conjugate electrophoresis (FSCE) if a “drag‐tag” is attached to DNA fragments because the tag breaks the linear charge‐to‐friction scaling. Although several previous studies have demonstrated the feasibility of DNA separation by free‐solution conjugated electrophoresis, generation of a monodisperse drag‐tag and identification of a strong, site‐specific conjugation method between a DNA fragment and a drag‐tag are challenges that still remain. In this study, we demonstrate an efficient FSCE method by conjugating a biologically synthesized elastin‐like polypeptide (ELP) and green fluorescent protein (GFP) to DNA fragments. In addition, to produce strong and site‐specific conjugation, a methionine residue in drag‐tags is replaced with homopropargylglycine (Hpg), which can be conjugated specifically to a DNA fragment with an azide site.     

Functional Nanopores Enabled with DNA

Membrane-spanning nanopores are used in label-free single-molecule sensing and next-generation portable nucleic acid sequencing, and as powerful research tools in biology, biophysics, and synthetic biology. Naturally occurring protein and peptide pores, as well as synthetic inorganic nanopores, are used in these applications, with their limitations. The structural and functional repertoire of nanopores can be considerably expanded by functionalising existing pores with DNA strands and by creating an entirely new class of nanopores with DNA nanotechnology. This review outlines progress in this area of functional DNA nanopores and outlines developments to open up new applications.     

Bisulfite-free and single-nucleotide resolution sequencing of DNA epigenetic modification of 5-hydroxymethylcytosine using engineered deaminase

The discovery of 5-hydroxymethylcytosine (5hmC) in mammalian genomes is a landmark in epigenomics study. Similar to 5-methylcytosine (5mC), 5hmC is viewed as a critical epigenetic modification. Deciphering the functions of 5hmC necessitates the location analysis of 5hmC in genomes. Here, we proposed an engineered deaminase-mediated sequencing (EDM-seq) method for the quantitative detection of 5hmC in DNA at single-nucleotide resolution. This method capitalizes on the engineered human apolipoprotein B mRNA-editing catalytic polypeptide-like 3A (A3A) protein to produce differential deamination activity toward cytosine, 5mC, and 5hmC. In EDM-seq, the engineered A3A (eA3A) protein can deaminate C and 5mC but not 5hmC. The original C and 5mC in DNA are deaminated by eA3A to form U and T, both of which are read as T during sequencing, while 5hmC is resistant to deamination by eA3A and is still read as C during sequencing. Therefore, the remaining C in the sequence manifests the original 5hmC. By EDM-seq, we achieved the quantitative detection of 5hmC in genomic DNA of lung cancer tissue. The EDM-seq method is bisulfite-free and does not require DNA glycosylation or chemical treatment, which offers a valuable tool for the straightforward and quantitative detection of 5hmC in DNA at single-nucleotide resolution.

In EDM-seq, the original C and 5mC in DNA are deaminated by eA3A to form U and T, both of which are read as T during sequencing. While the 5hmC is resistant to deamination by eA3A and is still read as C during sequencing.