Design of optimal labeling patterns for optical genome mapping via information theory.

MOTIVATION: Optical genome mapping (OGM) is a technique that extracts partial genomic information from optically imaged and linearized DNA fragments containing fluorescently labeled short sequence patterns. This information can be used for various genomic analyses and applications, such as the detection of structural variations and copy-number variations, epigenomic profiling, and microbial species identification. Currently, the choice of labeled patterns is based on the available biochemical methods and is not necessarily optimized for the application. RESULTS: In this work, we develop a model of OGM based on information theory, which enables the design of optimal labeling patterns for specific applications and target organism genomes. We validated the model through experimental OGM on human DNA and simulations on bacterial DNA. Our model predicts up to 10-fold improved accuracy by optimal choice of labeling patterns, which may guide future development of OGM biochemical labeling methods and significantly improve its accuracy and yield for applications such as epigenomic profiling and cultivation-free pathogen identification in clinical samples. AVAILABILITY AND IMPLEMENTATION: https://github.com/yevgenin/PatternCode.

DeepOM: single-molecule optical genome mapping via deep learning.

MOTIVATION: Efficient tapping into genomic information from a single microscopic image of an intact DNA molecule is an outstanding challenge and its solution will open new frontiers in molecular diagnostics. Here, a new computational method for optical genome mapping utilizing deep learning is presented, termed DeepOM. Utilization of a convolutional neural network, trained on simulated images of labeled DNA molecules, improves the success rate in the alignment of DNA images to genomic references. RESULTS: The method is evaluated on acquired images of human DNA molecules stretched in nano-channels. The accuracy of the method is benchmarked against state-of-the-art commercial software Bionano Solve. The results show a significant advantage in alignment success rate for molecules shorter than 50 kb. DeepOM improves the yield, sensitivity, and throughput of optical genome mapping experiments in applications of human genomics and microbiology. AVAILABILITY AND IMPLEMENTATION: The source code for the presented method is publicly available at https://github.com/yevgenin/DeepOM.

Chemoenzymatic labeling of DNA methylation patterns for single-molecule epigenetic mapping.

DNA methylation, specifically, methylation of cytosine (C) nucleotides at the 5-carbon position (5-mC), is the most studied and significant epigenetic modification. Here we developed a chemoenzymatic procedure to fluorescently label non-methylated cytosines in CpG context, allowing epigenetic profiling of single DNA molecules spanning hundreds of thousands of base pairs. We used a CpG methyltransferase with a synthetic S-adenosyl-l-methionine cofactor analog to transfer an azide to cytosines instead of the natural methyl group. A fluorophore was then clicked onto the DNA, reporting on the amount and position of non-methylated CpGs. We found that labeling efficiency was increased up to 2-fold by the addition of a nucleosidase, presumably by degrading the inactive by-product of the cofactor after labeling, preventing its inhibitory effect. We used the method to determine the decline in global DNA methylation in a chronic lymphocytic leukemia patient and then performed whole-genome methylation mapping of the model plant Arabidopsis thaliana. Our genome maps show high concordance with published bisulfite sequencing methylation maps. Although mapping resolution is limited by optical detection to 500-1000 bp, the labeled DNA molecules produced by this approach are hundreds of thousands of base pairs long, allowing access to long repetitive and structurally variable genomic regions.

Chemo-Enzymatic Fluorescence Labeling Of Genomic DNA For Simultaneous Detection Of Global 5-Methylcytosine And 5-Hydroxymethylcytosine.

5-Methylcytosine and 5-hydroxymethylcytosine are epigenetic modifications involved in gene regulation and cancer. We present a new, simple, and high-throughput platform for multi-color epigenetic analysis. The novelty of our approach is the ability to multiplex methylation and de-methylation signals in the same assay. We utilize an engineered methyltransferase enzyme that recognizes and labels all unmodified CpG sites with a fluorescent cofactor. In combination with the already established labeling of the de-methylation mark 5-hydroxymethylcytosine via enzymatic glycosylation, we obtained a robust platform for simultaneous epigenetic analysis of these marks. We assessed the global epigenetic levels in multiple samples of colorectal cancer and observed a 3.5-fold reduction in 5hmC levels but no change in DNA methylation levels between sick and healthy individuals. We also measured epigenetic modifications in chronic lymphocytic leukemia and observed a decrease in both modification levels (5-hydroxymethylcytosine: whole blood 30 %; peripheral blood mononuclear cells (PBMCs) 40 %. 5-methylcytosine: whole blood 53 %; PBMCs 48 %). Our findings propose using a simple blood test as a viable method for analysis, simplifying sample handling in diagnostics. Importantly, our results highlight the assay's potential for epigenetic evaluation of clinical samples, benefiting research and patient management.

Long-read single-molecule maps of the functional methylome.

We report on the development of a methylation analysis workflow for optical detection of fluorescent methylation profiles along chromosomal DNA molecules. In combination with Bionano Genomics genome mapping technology, these profiles provide a hybrid genetic/epigenetic genome-wide map composed of DNA molecules spanning hundreds of kilobase pairs. The method provides kilobase pair-scale genomic methylation patterns comparable to whole-genome bisulfite sequencing (WGBS) along genes and regulatory elements. These long single-molecule reads allow for methylation variation calling and analysis of large structural aberrations such as pathogenic macrosatellite arrays not accessible to single-cell second-generation sequencing. The method is applied here to study facioscapulohumeral muscular dystrophy (FSHD), simultaneously recording the haplotype, copy number, and methylation status of the disease-associated, highly repetitive locus on Chromosome 4q.

Long reads capture simultaneous enhancer-promoter methylation status for cell-type deconvolution.

MOTIVATION: While promoter methylation is associated with reinforcing fundamental tissue identities, the methylation status of distant enhancers was shown by genome-wide association studies to be a powerful determinant of cell-state and cancer. With recent availability of long reads that report on the methylation status of enhancer-promoter pairs on the same molecule, we hypothesized that probing these pairs on the single-molecule level may serve the basis for detection of rare cancerous transformations in a given cell population. We explore various analysis approaches for deconvolving cell-type mixtures based on their genome-wide enhancer-promoter methylation profiles. RESULTS: To evaluate our hypothesis we examine long-read optical methylome data for the GM12878 cell line and myoblast cell lines from two donors. We identified over 100 000 enhancer-promoter pairs that co-exist on at least 30 individual DNA molecules. We developed a detailed methodology for mixture deconvolution and applied it to estimate the proportional cell compositions in synthetic mixtures. Analysis of promoter methylation, as well as enhancer-promoter pairwise methylation, resulted in very accurate estimates. In addition, we show that pairwise methylation analysis can be generalized from deconvolving different cell types to subtle scenarios where one wishes to resolve different cell populations of the same cell-type. AVAILABILITY AND IMPLEMENTATION: The code used in this work to analyze single-molecule Bionano Genomics optical maps is available via the GitHub repository https://github.com/ebensteinLab/Single_molecule_methylation_in_EP.

Enzyme-free optical DNA mapping of the human genome using competitive binding.

Optical DNA mapping (ODM) allows visualization of long-range sequence information along single DNA molecules. The data can for example be used for detecting long range structural variations, for aiding DNA sequence assembly of complex genomes and for mapping epigenetic marks and DNA damage across the genome. ODM traditionally utilizes sequence specific marks based on nicking enzymes, combined with a DNA stain, YOYO-1, for detection of the DNA contour. Here we use a competitive binding approach, based on YOYO-1 and netropsin, which highlights the contour of the DNA molecules, while simultaneously creating a continuous sequence specific pattern, based on the AT/GC variation along the detected molecule. We demonstrate and validate competitive-binding-based ODM using bacterial artificial chromosomes (BACs) derived from the human genome and then turn to DNA extracted from white blood cells. We generalize our findings with in-silico simulations that show that we can map a vast majority of the human genome. Finally, we demonstrate the possibility of combining competitive binding with enzymatic labeling by mapping DNA damage sites induced by the cytotoxic drug etoposide to the human genome. Overall, we demonstrate that competitive-binding-based ODM has the potential to be used both as a standalone assay for studies of the human genome, as well as in combination with enzymatic approaches, some of which are already commercialized.

5-Hydroxymethylcytosine as a clinical biomarker: Fluorescence-based assay for high-throughput epigenetic quantification in human tissues.

Epigenetic transformations may provide early indicators for cancer and other disease. Specifically, the amount of genomic 5-hydroxymethylcytosine (5-hmC) was shown to be globally reduced in a wide range of cancers. The integration of this global biomarker into diagnostic workflows is hampered by the limitations of current 5-hmC quantification methods. Here we present and validate a fluorescence-based platform for high-throughput and cost-effective quantification of global genomic 5-hmC levels. We utilized the assay to characterize cancerous tissues based on their 5-hmC content, and observed a pronounced reduction in 5-hmC level in various cancer types. We present data for glioblastoma, colorectal cancer, multiple myeloma, chronic lymphocytic leukemia and pancreatic cancer, compared to corresponding controls. Potentially, the technique could also be used to follow response to treatment for personalized treatment selection. We present initial proof-of-concept data for treatment of familial adenomatous polyposis.

Rapid Quantification of Oxidation and UV Induced DNA Damage by Repair Assisted Damage Detection-(Rapid RADD).

Knowing the amount and type of DNA damage is of great significance for a broad range of clinical and research applications. However, existing methods are either lacking in their ability to distinguish between types of DNA damage or limited in their sensitivity and reproducibility. The method described herein enables rapid and robust quantification of type-specific single-strand DNA damage. The method is based on repair-assisted damage detection (RADD) by which fluorescent nucleotides are incorporated into DNA damage sites using type-specific repair enzymes. Up to 90 DNA samples are then deposited on a multiwell glass slide, and analyzed by a conventional slide scanner for quantification of DNA damage levels. Accurate and sensitive measurements of oxidative or UV-induced DNA damage levels and repair kinetics are presented for both in vitro and in vivo models.

Selective nanopore sequencing of human BRCA1 by Cas9-assisted targeting of chromosome segments (CATCH).

Next generation sequencing (NGS) is challenged by structural and copy number variations larger than the typical read length of several hundred bases. Third-generation sequencing platforms such as single-molecule real-time (SMRT) and nanopore sequencing provide longer reads and are able to characterize variations that are undetected in NGS data. Nevertheless, these technologies suffer from inherent low throughput which prohibits deep sequencing at reasonable cost without target enrichment. Here, we optimized Cas9-Assisted Targeting of CHromosome segments (CATCH) for nanopore sequencing of the breast cancer gene BRCA1. A 200 kb target containing the 80 kb BRCA1 gene body and its flanking regions was isolated intact from primary human peripheral blood cells, allowing long-range amplification and long-read nanopore sequencing. The target was enriched 237-fold and sequenced at up to 70× coverage on a single flow-cell. Overall performance and single-nucleotide polymorphism (SNP) calling were directly compared to Illumina sequencing of the same enriched sample, highlighting the benefits of CATCH for targeted sequencing. The CATCH enrichment scheme only requires knowledge of the target flanking sequence for Cas9 cleavage while providing contiguous data across both coding and non-coding sequence and holds promise for characterization of complex disease-related or highly variable genomic regions.

Irys Extract.

Summary: Irys Extract is a software tool for generating genomic information from data collected by the BioNano Genomics Irys platform. The tool allows the user easy access to the raw data in the form of cropped images and genetically aligned intensity profiles. The latter are also made compatible with the BED format for using with popular genomic browsers such as the UCSC Genome Browser. Availability and implementation: Irys Extract has been developed in Matlab R2015a, it was tested to work with IrysView 2.4.0.15879 and AutoDetect 2.1.4.9159, and it currently runs under Microsoft Windows operating systems (7-10). Irys Extract can be downloaded alongside its manual and a demo dataset at http://www.nanobiophotonix.com and https://sites.google.com/site/raniarielly/. Contact: uv@post.tau.ac.il. Supplementary information: Supplementary data are available at Bioinformatics online.

Microfluidic DNA combing for parallel single-molecule analysis.

DNA combing is a widely used method for stretching and immobilising DNA molecules on a surface. Fluorescent labelling of genomic information enables high-resolution optical analysis of DNA at the single-molecule level. Despite its simplicity, the application of DNA combing in diagnostic workflows is still limited, mainly due to difficulties in analysing multiple small-volume DNA samples in parallel. Here, we report a simple and versatile microfluidic DNA combing technology (µDC), which allows manipulating, stretching and imaging of multiple, microliter scale DNA samples by employing a manifold of parallel microfluidic channels. Using DNA molecules with repetitive units as molecular rulers, we demonstrate that the µDC technology allows uniform stretching of DNA molecules. The stretching ratio remains consistent along individual molecules as well as between different molecules in the various channels, allowing simultaneous quantitative analysis of different samples loaded into parallel channels. Furthermore, we demonstrate the application of µDC to characterise UVB-induced DNA damage levels in human embryonic kidney cells and the spatial correlation between DNA damage sites. Our results point out the potential application of µDC for quantitative and comparative single-molecule studies of genomic features. The extremely simple design of µDC makes it suitable for integration into other microfluidic platforms to facilitate high-throughput DNA analysis in biological research and medical point-of-care applications.

Sizing femtogram amounts of dsDNA by single-molecule counting.

Modern molecular-biology applications raise renewed interest in sizing minute-amounts of DNA. In this work we utilize single-molecule imaging with in situ size calibration to accurately analyze the size and mass distribution of DNA samples. We exploit the correlation between DNA length and its fluorescence intensity after staining in order to assess the length of individual DNA fragments by fluorescence microscopy. Synthetic reference DNA standards are added to the investigated sample before staining and serve as internal size calibrators, supporting a robust assay for accurate DNA sizing. Our results demonstrate the ability to reconstruct the exact length distribution in a complex DNA sample by sizing a subset containing only femtogram amounts of DNA, thus, outperforming microfluidic gel electrophoresis which is the currently accepted gold standard. This assay may find useful applications for genetic analysis where the exact size distribution of DNA molecules is critical and the availability of genetic material is limited.

Bacteriophage strain typing by rapid single molecule analysis.

Rapid characterization of unknown biological samples is under the focus of many current studies. Here we report a method for screening of biological samples by optical mapping of their DNA. We use a novel, one-step chemo-enzymatic reaction to covalently bind fluorophores to DNA at the four-base recognition sites of a DNA methyltransferase. Due to the diffraction limit of light, the dense distribution of labels results in a continuous fluorescent signal along the DNA. The amplitude modulations (AM) of the fluorescence intensity along the stretched DNA molecules exhibit a unique molecular fingerprint that can be used for identification. We show that this labelling scheme is highly informative, allowing accurate genotyping. We demonstrate the method by labelling the genomes of λ and T7 bacteriophages, resulting in a consistent, unique AM profile for each genome. These profiles are also successfully used for identification of the phages from a background phage library. Our method may provide a facile route for screening and typing of various organisms and has potential applications in metagenomics studies of various ecosystems.

Spatiotemporal manipulation of retinoic acid activity in zebrafish hindbrain development via photo-isomerization.

All-trans retinoic acid (RA) is a key player in many developmental pathways. Most methods used to study its effects in development involve continuous all-trans RA activation by incubation in a solution of all-trans RA or by implanting all-trans RA-soaked beads at desired locations in the embryo. Here we show that the UV-driven photo-isomerization of 13-cis RA to the trans-isomer (and vice versa) can be used to non-invasively and quantitatively control the concentration of all-trans RA in a developing embryo in time and space. This facilitates the global or local perturbation of developmental pathways with a pulse of all-trans RA of known concentration or its inactivation by UV illumination. In zebrafish embryos in which endogenous synthesis of all-trans RA is impaired, incubation for as little as 5 minutes in 1 nM all-trans RA (a pulse) or 5 nM 13-cis RA followed by 1-minute UV illumination is sufficient to rescue the development of the hindbrain if performed no later than bud stage. However, if subsequent to this all-trans RA pulse the embryo is illuminated (no later than bud stage) for 1 minute with UV light (to isomerize, i.e. deactivate, all-trans RA), the rescue of hindbrain development is impaired. This suggests that all-trans RA is sequestered in embryos that have been transiently exposed to it. Using 13-cis RA isomerization with UV light, we further show that local illumination at bud stage of the head region (but not the tail) is sufficient to rescue hindbrain formation in embryos whose all-trans RA synthetic pathway has been impaired.

One-Pot Chemoenzymatic Cascade for Labeling of the Epigenetic Marker 5-Hydroxymethylcytosine.

The epigenetic DNA modification 5-hydroxymethylcytosine (5-hmC) is important for the regulation of gene expression during development and in tumorigenesis. 5-hmC can be selectively glycosylated by T4 ß-glucosyltransferase (ß-GT); introduction of an azide on the attached sugar provides a chemical handle for isolation or fluorescent tagging of 5-hmC residues by click chemistry. This approach has not been broadly adopted because of the challenging synthesis and limited commercial availability of the glycosylation substrate, 6-deoxy-6-azido-α-D-glucopyranoside. We report the enzyme-assisted synthesis of this precursor by the uridylyltransferase from Pasteurella multocida (PmGlmU). We were able to directly label 5-hmC in genomic DNA by an enzymatic cascade involving successive action of PmGlmU and ß-GT. This is a facile and cost-effective one-pot chemoenzymatic methodology for 5-hmC analysis.

Lighting up individual DNA damage sites by in vitro repair synthesis.

DNA damage and repair are linked to fundamental biological processes such as metabolism, disease, and aging. Single-strand lesions are the most abundant form of DNA damage; however, methods for characterizing these damage lesions are lacking. To avoid double-strand breaks and genomic instability, DNA damage is constantly repaired by efficient enzymatic machinery. We take advantage of this natural process and harness the repair capacity of a bacterial enzymatic cocktail to repair damaged DNA in vitro and incorporate fluorescent nucleotides into damage sites as part of the repair process. We use single-molecule imaging to detect individual damage sites in genomic DNA samples. When the labeled DNA is extended on a microscope slide, damage sites are visualized as fluorescent spots along the DNA contour, and the extent of damage is easily quantified. We demonstrate the ability to quantitatively follow the damage dose response to different damaging agents as well as repair dynamics in response to UV irradiation in several cell types. Finally, we show the modularity of this single-molecule approach by labeling DNA damage in conjunction with 5-hydroxymethylcytosine in genomic DNA extracted from mouse brain tissue.

Enzyme-responsive amphiphilic PEG-dendron hybrids and their assembly into smart micellar nanocarriers.

Enzyme-responsive micelles have great potential as drug delivery platforms due to the high selectivity of the activating enzymes. Here we report a highly modular design for the efficient and simple synthesis of amphiphilic block copolymers based on a linear hydrophilic polyethyleneglycol (PEG) and an enzyme-responsive hydrophobic dendron. These amphiphilic hybrids self-assemble in water into micellar nanocontainers that can disassemble and release encapsulated molecular cargo upon enzymatic activation. The utilization of monodisperse dendrons as the stimuli-responsive block enabled a detailed kinetic study of the molecular mechanism of the enzymatically triggered disassembly. The modularity of these PEG-dendron hybrids allows control over the disassembly rate of the formed micelles by simply tuning the PEG length. Such smart amphiphilic hybrids could potentially be applied for the fabrication of nanocarriers with adjustable release rates for delivery applications.

Spectroscopic quantification of 5-hydroxymethylcytosine in genomic DNA.

5-Hydroxymethylcytosine (5hmC), a modified form of the DNA base cytosine, is an important epigenetic mark linked to regulation of gene expression in development, and tumorigenesis. We have developed a spectroscopic method for a global quantification of 5hmC in genomic DNA. The assay is performed within a multiwell plate, which allows simultaneous recording of up to 350 samples. Our quantification procedure of 5hmC is direct, simple, and rapid. It relies on a two-step protocol that consists of enzymatic glucosylation of 5hmC with an azide-modified glucose, followed by a "click reaction" with an alkyne-fluorescent tag. The fluorescence intensity recorded from the DNA sample is proportional to its 5hmC content and can be quantified by a simple plate reader measurement. This labeling technique is specific and highly sensitive, allowing detection of 5hmC down to 0.002% of the total nucleotides. Our results reveal significant variations in the 5hmC content obtained from different mouse tissues, in agreement with previously reported data.

Programmable comprehensive controller for multi-color 3D confocal spinning-disk image scanning microscope.

This paper introduces a compact, portable, and highly accurate triggering control system for a 3D confocal spinning-disk image scanning microscope (CSD-ISM). Building upon on our previously published research, we expanded the hardware of the controller and synchronized it with a sub-micron translator which scans the object in the z-direction. As well as expanding the hardware, the software also was extended from previously published work similarly as it is stated for hardware while allowing full control over the 3D movement. We showed a clear and smooth 3D image made up of a collection of 2D images at different heights.