Long-read single-molecule maps of the functional methylome

Long-Read Structural and Epigenetic Profiling of a Kidney Tumor-Matched Sample with Nanopore Sequencing and Optical Genome Mapping

Carcinogenesis often involves significant alterations in the cancer genome architecture, marked by large structural and copy number variations (SVs and CNVs) that are difficult to capture with short-read sequencing. Traditionally, cytogenetic techniques are applied to detect such aberrations, but they are limited in resolution and do not cover features smaller than several hundred kilobases. Optical genome mapping and nanopore sequencing are attractive technologies that bridge this resolution gap and offer enhanced performance for cytogenetic applications. These methods profile native, individual DNA molecules, thus capturing epigenetic information. We applied both techniques to characterize a clear cell renal cell carcinoma (ccRCC) tumor's structural and copy number landscape, highlighting the relative strengths of each method in the context of variant size and average read length. Additionally, we assessed their utility for methylome and hydroxymethylome profiling, emphasizing differences in epigenetic analysis applicability.

OM2Seq: learning retrieval embeddings for optical genome mapping

Motivation

Genomics-based diagnostic methods that are quick, precise, and economical are essential for the advancement of precision medicine, with applications spanning the diagnosis of infectious diseases, cancer, and rare diseases. One technology that holds potential in this field is optical genome mapping (OGM), which is capable of detecting structural variations, epigenomic profiling, and microbial species identification. It is based on imaging of linearized DNA molecules that are stained with fluorescent labels, that are then aligned to a reference genome. However, the computational methods currently available for OGM fall short in terms of accuracy and computational speed.

Results

This work introduces OM2Seq, a new approach for the rapid and accurate mapping of DNA fragment images to a reference genome. Based on a Transformer-encoder architecture, OM2Seq is trained on acquired OGM data to efficiently encode DNA fragment images and reference genome segments to a common embedding space, which can be indexed and efficiently queried using a vector database. We show that OM2Seq significantly outperforms the baseline methods in both computational speed (by 2 orders of magnitude) and accuracy.

Availability and implementation

https://github.com/yevgenin/om2seq.

The expanding diagnostic toolbox for rare genetic diseases

Genomic technologies, such as targeted, exome and short-read genome sequencing approaches, have revolutionized the care of patients with rare genetic diseases. However, more than half of patients remain without a diagnosis. Emerging approaches from research-based settings such as long-read genome sequencing and optical genome mapping hold promise for improving the identification of disease-causal genetic variants. In addition, new omic technologies that measure the transcriptome, epigenome, proteome or metabolome are showing great potential for variant interpretation. As genetic testing options rapidly expand, the clinical community needs to be mindful of their individual strengths and limitations, as well as remaining challenges, to select the appropriate diagnostic test, correctly interpret results and drive innovation to address insufficiencies. If used effectively - through truly integrative multi-omics approaches and data sharing - the resulting large quantities of data from these established and emerging technologies will greatly improve the interpretative power of genetic and genomic diagnostics for rare diseases.

Benchmarking long-read genome sequence alignment tools for human genomics applications

Background

The utility of long-read genome sequencing platforms has been shown in many fields including whole genome assembly, metagenomics, and amplicon sequencing. Less clear is the applicability of long reads to reference-guided human genomics, which is the foundation of genomic medicine. Here, we benchmark available platform-agnostic alignment tools on datasets from nanopore and single-molecule real-time platforms to understand their suitability in producing a genome representation.

Results

For this study, we leveraged publicly-available data from sample NA12878 generated on Oxford Nanopore and sample NA24385 on Pacific Biosciences platforms. We employed state of the art sequence alignment tools including GraphMap2, long-read aligner (LRA), Minimap2, CoNvex Gap-cost alignMents for Long Reads (NGMLR), and Winnowmap2. Minimap2 and Winnowmap2 were computationally lightweight enough for use at scale, while GraphMap2 was not. NGMLR took a long time and required many resources, but produced alignments each time. LRA was fast, but only worked on Pacific Biosciences data. Each tool widely disagreed on which reads to leave unaligned, affecting the end genome coverage and the number of discoverable breakpoints. No alignment tool independently resolved all large structural variants (1,001-100,000 base pairs) present in the Database of Genome Variants (DGV) for sample NA12878 or the truthset for NA24385.

Conclusions

These results suggest a combined approach is needed for LRS alignments for human genomics. Specifically, leveraging alignments from three tools will be more effective in generating a complete picture of genomic variability. It should be best practice to use an analysis pipeline that generates alignments with both Minimap2 and Winnowmap2 as they are lightweight and yield different views of the genome. Depending on the question at hand, the data available, and the time constraints, NGMLR and LRA are good options for a third tool. If computational resources and time are not a factor for a given case or experiment, NGMLR will provide another view, and another chance to resolve a case. LRA, while fast, did not work on the nanopore data for our cluster, but PacBio results were promising in that those computations completed faster than Minimap2. Due to its significant burden on computational resources and slow run time, Graphmap2 is not an ideal tool for exploration of a whole human genome generated on a long-read sequencing platform.

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.

Single-molecule approaches for DNA damage detection and repair: A focus on Repair Assisted Damage Detection (RADD)

The human genome is continually exposed to various stressors, which can result in DNA damage, mutations, and diseases. Among the different types of DNA damage, single-strand lesions are commonly induced by external stressors and metabolic processes. Accurate detection and quantification of DNA damage are crucial for understanding repair mechanisms, assessing environmental impacts, and evaluating response to therapy. However, traditional techniques have limitations in sensitivity and the ability to detect multiple types of damage. In recent years, single-molecule fluorescence approaches have emerged as powerful tools for precisely localizing and quantifying DNA damage. Repair Assisted Damage Detection (RADD) is a single-molecule technique that employs specific repair enzymes to excise damaged bases and incorporates fluorescently labeled nucleotides to visualize the damage. This technique provides valuable insights into repair efficiency and sequence-specific damage. In this review, we discuss the principles and applications of RADD assays, highlighting their potential for enhancing our understanding of DNA damage and repair processes.

Super-resolution imaging of linearized chromatin in tunable nanochannels

Nanofluidic linearization and optical mapping of naked DNA have been reported in the research literature, and implemented in commercial instruments. However, the resolution with which DNA features can be resolved is still inherently limited by both Brownian motion and diffraction-limited optics. Direct analysis of native chromatin is further hampered by difficulty in electrophoretic manipulation, which is routinely used for DNA analysis. This paper describes the development of a three-layer, tunable, nanochannel system that enables non-electrophoretic linearization and immobilization of native chromatin. Furthermore, through careful selection of self-blinking fluorescent dyes and the design of the nanochannel system, we achieve direct stochastic optical reconstruction microscopy (dSTORM) super-resolution imaging of the linearized chromatin. As an initial demonstration, rDNA chromatin extracted from Tetrahymena is analyzed by multi-color imaging of total DNA, newly synthesized DNA, and newly synthesized histone H3. Our analysis reveals a relatively even distribution of newly synthesized H3 across two halves of the rDNA chromatin with palindromic symmetry, supporting dispersive nucleosome segregation. As a proof-of-concept study, our work achieves super-resolution imaging of native chromatin fibers linearized and immobilized in tunable nanochannels. It opens up a new avenue for collecting long-range and high-resolution epigenetic information as well as genetic information.

Sequence-specific DNA labelling for fluorescence microscopy

The preservation of nucleus structure during microscopy imaging is a top priority for understanding chromatin organization, genome dynamics, and gene expression regulation. In this review, we summarize the sequence-specific DNA labelling methods that can be used for imaging in fixed and/or living cells without harsh treatment and DNA denaturation: (i) hairpin polyamides, (ii) triplex-forming oligonucleotides, (iii) dCas9 proteins, (iv) transcription activator-like effectors (TALEs) and (v) DNA methyltransferases (MTases). All these techniques are capable of identifying repetitive DNA loci and robust probes are available for telomeres and centromeres, but visualizing single-copy sequences is still challenging. In our futuristic vision, we see gradual replacement of the historically important fluorescence in situ hybridization (FISH) by less invasive and non-destructive methods compatible with live cell imaging. Combined with super-resolution fluorescence microscopy, these methods will open the possibility to look into unperturbed structure and dynamics of chromatin in living cells, tissues and whole organisms.

Synthesis of S-Adenosyl-L-Methionine Analogs with Extended Transferable Groups for Methyltransferase-Directed Labeling of DNA and RNA

S-Adenosyl-L-methionine (AdoMet) is a ubiquitous methyl donor for a variety of biological methylation reactions catalyzed by methyltransferases (MTases). AdoMet analogs with extended propargylic chains replacing the sulfonium-bound methyl group can serve as surrogate cofactors for many DNA and RNA MTases, enabling covalent derivatization and subsequent labeling of their cognate target sites in DNA or RNA. Although AdoMet analogs with saturated aliphatic chains are less popular than propargylic ones, they can be useful for dedicated studies that require certain chemical derivatization. Here we describe synthetic procedures for the preparation of two AdoMet analogs, one with a transferable 6-azidohex-2-ynyl group (carrying an activating C≡C triple bond and a terminal azide functionality), and the other one with a transferable ethyl-2,2,2-d₃ group (an isotope-labeled aliphatic moiety). Our synthetic approach is based on direct chemoselective alkylation of S-adenosyl-L-homocysteine at sulfur with a corresponding nosylate or triflate, respectively, under acidic conditions. We also describe synthetic routes to 6-azidohex-2-yn-1-ol and conversion of the alcohols to corresponding nosylate and triflate alkylators. Using these protocols, the synthetic AdoMet analogs can be prepared within 1 to 2 weeks. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Synthesis of 6-azidohex-2-yn-1-ol Basic Protocol 2: Synthesis of 4-nitrobenzenesulfonate Basic Protocol 3: Synthesis of trifluoromethanesulfonates Basic Protocol 4: S-Alkylation of AdoHcy with sulfonates Basic Protocol 5: Purification and characterization of AdoMet analogs.

Dam Assisted Fluorescent Tagging of Chromatin Accessibility (DAFCA) for Optical Genome Mapping in Nanochannel Arrays

Proteins and enzymes in the cell nucleus require physical access to their DNA target sites in order to perform genomic tasks such as gene activation and transcription. Hence, chromatin accessibility is a central regulator of gene expression, and its genomic profile holds essential information on the cell type and state. We utilized the E. coli Dam methyltransferase in combination with a fluorescent cofactor analogue to generate fluorescent tags in accessible DNA regions within the cell nucleus. The accessible portions of the genome are then detected by single-molecule optical genome mapping in nanochannel arrays. This method allowed us to characterize long-range structural variations and their associated chromatin structure. We show the ability to create whole-genome, allele-specific chromatin accessibility maps composed of long DNA molecules extended in silicon nanochannels.

Methylation of the 4q35 D4Z4 repeat defines disease status in facioscapulohumeral muscular dystrophy

Genetic diagnosis of facioscapulohumeral muscular dystrophy (FSHD) remains a challenge in clinical practice as it cannot be detected by standard sequencing methods despite being the third most common muscular dystrophy. The conventional diagnostic strategy addresses the known genetic parameters of FSHD: the required presence of a permissive haplotype, a size reduction of the D4Z4 repeat of chromosome 4q35 (defining FSHD1) or a pathogenic variant in an epigenetic suppressor gene (consistent with FSHD2). Incomplete penetrance and epistatic effects of the underlying genetic parameters as well as epigenetic parameters (D4Z4 methylation) pose challenges to diagnostic accuracy and hinder prediction of clinical severity. In order to circumvent the known limitations of conventional diagnostics and to complement genetic parameters with epigenetic ones, we developed and validated a multistage diagnostic workflow that consists of a haplotype analysis and a high-throughput methylation profile analysis (FSHD-MPA). FSHD-MPA determines the average global methylation level of the D4Z4 repeat array as well as the regional methylation of the most distal repeat unit by combining bisulphite conversion with next-generation sequencing and a bioinformatics pipeline and uses these as diagnostic parameters. We applied the diagnostic workflow to a cohort of 148 patients and compared the epigenetic parameters based on FSHD-MPA to genetic parameters of conventional genetic testing. In addition, we studied the correlation of repeat length and methylation level within the most distal repeat unit with age-corrected clinical severity and age at disease onset in FSHD patients. The results of our study show that FSHD-MPA is a powerful tool to accurately determine the epigenetic parameters of FSHD, allowing discrimination between FSHD patients and healthy individuals, while simultaneously distinguishing FSHD1 and FSHD2. The strong correlation between methylation level and clinical severity indicates that the methylation level determined by FSHD-MPA accounts for differences in disease severity among individuals with similar genetic parameters. Thus, our findings further confirm that epigenetic parameters rather than genetic parameters represent FSHD disease status and may serve as a valuable biomarker for disease status.

Computational Methods for Single-cell DNA Methylome Analysis

Dissecting intercellular epigenetic differences is key to understanding tissue heterogeneity. Recent advances in single-cell DNA methylome profiling have presented opportunities to resolve this heterogeneity at the maximum resolution. While these advances enable us to explore frontiers of chromatin biology and better understand cell lineage relationships, they pose new challenges in data processing and interpretation. This review surveys the current state of computational tools developed for single-cell DNA methylome data analysis. We discuss critical components of single-cell DNA methylome data analysis, including data preprocessing, quality control, imputation, dimensionality reduction, cell clustering, supervised cell annotation, cell lineage reconstruction, gene activity scoring, and integration with transcriptome data. We also highlight unique aspects of single-cell DNA methylome data analysis and discuss how techniques common to other single-cell omics data analyses can be adapted to analyze DNA methylomes. Finally, we discuss existing challenges and opportunities for future development.

Testing for Facioscapulohumeral Muscular Dystrophy with Optical Genome Mapping

The introduction of optical genome mapping has improved time constraints and a lack of specificity from previous methodologies when performing genome-wide analyses of samples. Optical genome mapping allows for the detection of structural variations, aberrations, and functionality traits from a single stained molecule of DNA. Though the preparation time is increased compared to previously utilized visualization techniques, optical genome mapping significantly reduces the time needed for analysis. Specifically, individual disease pipelines have been developed to rapidly analyze prepared samples. One of these diseases, Facioscapulohumeral Muscular Dystrophy (FSHD), is detected through quantification of the D4Z4 repeat array on chromosome 4q35. Optical genome mapping, with the ability to enumerate the repeats of the D4Z4 array, has demonstrated the capability to precisely diagnose FSHD. In this protocol, the preparation of samples and subsequent loading and analysis in an optical genome mapping system is discussed for the detection and analysis of FSHD. These methods should prove highly useful in FSHD analyses and beyond with the development of further disease analysis pipelines within the instrument. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Genomic DNA isolation, labeling, and staining Basic Protocol 2: Mapping and analysis with the Bionano Saphyr® system.

MGMT and Whole-Genome DNA Methylation Impacts on Diagnosis, Prognosis and Therapy of Glioblastoma Multiforme

Epigenetic changes in DNA methylation contribute to the development of many diseases, including cancer. In glioblastoma multiforme, the most prevalent primary brain cancer and an incurable tumor with a median survival time of 15 months, a single epigenetic modification, the methylation of the O⁶-Methylguanine-DNA Methyltransferase (MGMT) gene, is a valid biomarker for predicting response to therapy with alkylating agents and also, independently, prognosis. More recently, the progress from single gene to whole-genome analysis of DNA methylation has allowed a better subclassification of glioblastomas. Here, we review the clinically relevant information that can be obtained by studying MGMT gene and whole-genome DNA methylation changes in glioblastomas, also highlighting benefits, including those of liquid biopsy, and pitfalls of the different detection methods. Finally, we discuss how changes in DNA methylation, especially in glioblastomas bearing mutations in the Isocitrate Dehydrogenase (IDH) 1 and 2 genes, can be exploited as targets for tailoring therapy.

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.

Single-molecule analysis of genome-wide DNA methylation by fiber FISH coupled with atomic force microscopy

A straightforward single-molecule approach was developed for identifying whole-genome DNA methylation through fiber-FISH coupled with AFM. This method has advantages of low DNA input, reproduction, long reads and low cost.

Facioscapulohumeral muscular dystrophy type 2: an update on the clinical, genetic, and molecular findings

Facioscapulohumeral muscular dystrophy (FSHD) is a common genetic disease of the skeletal muscle with a characteristic pattern of weakness. Facioscapulohumeral muscular dystrophy type 2 (FSHD2) accounts for approximately 5% of all cases of FSHD and describes patients without a D4Z4 repeat contraction on chromosome 4. Phenotypically FSHD2 shows virtually no difference from FSHD1 and both forms of FSHD arise via a common downstream mechanism of epigenetic derepression of the transcription factor DUX4 in skeletal muscle cells. This results in expression of DUX4 and target genes leading to skeletal muscle toxicity. Over the past decade, major progress has been made in our understanding of the genetic and epigenetic architecture that underlies FSHD2 pathogenesis, as well as the clinical manifestations and disease progression. These include the finding that FSHD2 is a digenic disease and that mutations in the genes SMCHD1, DNMT3B, and more recently LRIF1, can cause FSHD2. FSHD2 is complex and it is important that clinicians keep abreast of recent developments; this review aims to serve as an update of the clinical, genetic, and molecular research into this condition.

Towards improved genetic diagnosis of human differences of sex development

Despite being collectively among the most frequent congenital developmental conditions worldwide, differences of sex development (DSD) lack recognition and research funding. As a result, what constitutes optimal management remains uncertain. Identification of the individual conditions under the DSD umbrella is challenging and molecular genetic diagnosis is frequently not achieved, which has psychosocial and health-related repercussions for patients and their families. New genomic approaches have the potential to resolve this impasse through better detection of protein-coding variants and ascertainment of under-recognized aetiology, such as mosaic, structural, non-coding or epigenetic variants. Ultimately, it is hoped that better outcomes data, improved understanding of the molecular causes and greater public awareness will bring an end to the stigma often associated with DSD.

Precise Epigenetic Analysis Using Targeted Bisulfite Genomic Sequencing Distinguishes FSHD1, FSHD2, and Healthy Subjects

The true prevalence of facioscapulohumeral muscular dystrophy (FSHD) is unknown due to difficulties with accurate clinical evaluation and the complexities of current genetic diagnostics. Interestingly, all forms of FSHD are linked to epigenetic changes in the chromosome 4q35 D4Z4 macrosatellite, suggesting that epigenetic analysis could provide an avenue for sequence-based FSHD diagnostics. However, studies assessing DNA methylation at the FSHD locus have produced conflicting results; thus, the utility of this technique as an FSHD diagnostic remains controversial. Here, we critically compared two protocols for epigenetic analysis of the FSHD region using bisulfite genomic sequencing: Jones et al., that contends to be individually diagnostic for FSHD1 and FSHD2, and Gaillard et al., that can identify some changes in DNA methylation levels between groups of clinically affected FSHD and healthy subjects, but is not individually diagnostic for any form of FSHD. We performed both sets of assays on the same genetically confirmed samples and showed that this discrepancy was due strictly to differences in amplicon specificity. We propose that the epigenetic status of the FSHD-associated D4Z4 arrays, when accurately assessed, is a diagnostic for genetic FSHD and can readily distinguish between healthy, FSHD1 and FSHD2. Thus, epigenetic diagnosis of FSHD, which can be performed on saliva DNA, will greatly increase accessibility to FSHD diagnostics for populations around the world.

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.

Multicolor Whole-Genome Mapping in Nanochannels for Genetic Analysis

Analysis of structural variations (SVs) is important to understand mutations underlying genetic disorders and pathogenic conditions. However, characterizing SVs using short-read, high-throughput sequencing technology is difficult. Although long-read sequencing technologies are being increasingly employed in characterizing SVs, their low throughput and high costs discourage widespread adoption. Sequence motif-based optical mapping in nanochannels is useful in whole-genome mapping and SV detection, but it is not possible to precisely locate the breakpoints or estimate the copy numbers. We present here a universal multicolor mapping strategy in nanochannels combining conventional sequence-motif labeling system with Cas9-mediated target-specific labeling of any 20-base sequences (20mers) to create custom labels and detect new features. The sequence motifs are labeled with green fluorophores and the 20mers are labeled with red fluorophores. Using this strategy, it is possible to not only detect the SVs but also utilize custom labels to interrogate the features not accessible to motif-labeling, locate breakpoints, and precisely estimate copy numbers of genomic repeats. We validated our approach by quantifying the D4Z4 copy numbers, a known biomarker for facioscapulohumeral muscular dystrophy (FSHD) and estimating the telomere length, a clinical biomarker for assessing disease risk factors in aging-related diseases and malignant cancers. We also demonstrate the application of our methodology in discovering transposable long non-interspersed Elements 1 (LINE-1) insertions across the whole genome.

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.

nanotatoR: a tool for enhanced annotation of genomic structural variants

BACKGROUND

Whole genome sequencing is effective at identification of small variants, but because it is based on short reads, assessment of structural variants (SVs) is limited. The advent of Optical Genome Mapping (OGM), which utilizes long fluorescently labeled DNA molecules for de novo genome assembly and SV calling, has allowed for increased sensitivity and specificity in SV detection. However, compared to small variant annotation tools, OGM-based SV annotation software has seen little development, and currently available SV annotation tools do not provide sufficient information for determination of variant pathogenicity.

RESULTS

We developed an R-based package, nanotatoR, which provides comprehensive annotation as a tool for SV classification. nanotatoR uses both external (DGV; DECIPHER; Bionano Genomics BNDB) and internal (user-defined) databases to estimate SV frequency. Human genome reference GRCh37/38-based BED files are used to annotate SVs with overlapping, upstream, and downstream genes. Overlap percentages and distances for nearest genes are calculated and can be used for filtration. A primary gene list is extracted from public databases based on the patient's phenotype and used to filter genes overlapping SVs, providing the analyst with an easy way to prioritize variants. If available, expression of overlapping or nearby genes of interest is extracted (e.g. from an RNA-Seq dataset, allowing the user to assess the effects of SVs on the transcriptome). Most quality-control filtration parameters are customizable by the user. The output is given in an Excel file format, subdivided into multiple sheets based on SV type and inheritance pattern (INDELs, inversions, translocations, de novo, etc.). nanotatoR passed all quality and run time criteria of Bioconductor, where it was accepted in the April 2019 release. We evaluated nanotatoR's annotation capabilities using publicly available reference datasets: the singleton sample NA12878, mapped with two types of enzyme labeling, and the NA24143 trio. nanotatoR was also able to accurately filter the known pathogenic variants in a cohort of patients with Duchenne Muscular Dystrophy for which we had previously demonstrated the diagnostic ability of OGM.

CONCLUSIONS

The extensive annotation enables users to rapidly identify potential pathogenic SVs, a critical step toward use of OGM in the clinical setting.

Label as you fold: methyltransferase-assisted functionalization of DNA nanostructures

Non-DNA labels are key components for the construction of functional DNA nanostructures. Here, we present a method to graft covalent labels onto DNA origami nanostructures in an enzymatic one-pot reaction. The DNA methyltransferase M.TaqI labels the DNA nanostructures with azide groups, which serve as universal attachment points via click chemistry. Direct labeling with fluorescent dyes is also demonstrated. The procedure yields structures with high fluorescence intensities and narrow intensity distributions. In combination with UV crosslinking it enables the creation of temperature-stable, intense fluorescent beacons.

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.

Quantifying DNA damage induced by ionizing radiation and hyperthermia using single DNA molecule imaging

Ionizing radiation (IR) is a common mode of cancer therapy, where DNA damage is the major reason of cell death. Here, we use an assay based on fluorescence imaging of single damaged DNA molecules isolated from radiated lymphocytes, to quantify IR induced DNA damage. The assay uses a cocktail of DNA-repair enzymes that recognizes and excises DNA lesions and then a polymerase and a ligase incorporate fluorescent nucleotides at the damage sites, resulting in a fluorescent “spot” at each site. The individual fluorescent spots can then be counted along single stretched DNA molecules and the global level of DNA damage can be quantified. Our results demonstrate that inclusion of the human apurinic/apyrimidinic endonuclease 1 (APE1) in the enzyme cocktail increases the sensitivity of the assay for detection of IR induced damage significantly. This optimized assay also allowed detection of a cooperative increase in DNA damage when IR was combined with mild hyperthermia, which is sometimes used as an adjuvant in IR therapy. Finally, we discuss how the method may be used to identify patients that are sensitive to IR and other types of DNA damaging agents.

Long-read sequencing to understand genome biology and cell function

Determining the sequence of DNA and RNA molecules has a huge impact on the understanding of cell biology and function. Recent advancements in next-generation short-read sequencing (NGS) technologies, drops in cost and a resolution down to the single-cell level shaped our current view on genome structure and function. Third-generation sequencing (TGS) methods further complete the knowledge about these processes based on long reads and the ability to analyze DNA or RNA at single molecule level. Long-read sequencing provides additional possibilities to study genome architecture and the composition of highly complex regions and to determine epigenetic modifications of nucleotide bases at a genome-wide level. We discuss the principles and advancements of long-read sequencing and its applications in genome biology.

Single-molecule analysis of nucleic acid biomarkers - A review

Nucleic acids are important biomarkers for disease detection, monitoring, and treatment. Advances in technologies for nucleic acid analysis have enabled discovery and clinical implementation of nucleic acid biomarkers. However, challenges remain with technologies for nucleic acid analysis, thereby limiting the use of nucleic acid biomarkers in certain contexts. Here, we review single-molecule technologies for nucleic acid analysis that can be used to overcome these challenges. We first discuss the various types of nucleic acid biomarkers important for clinical applications and conventional technologies for nucleic acid analysis. We then discuss technologies for single-molecule in vitro and in situ analysis of nucleic acid biomarkers. Finally, we discuss other ultra-sensitive techniques for nucleic acid biomarker detection.

Structural variation in the sequencing era

Identifying structural variation (SV) is essential for genome interpretation but has been historically difficult due to limitations inherent to available genome technologies. Detection methods that use ensemble algorithms and emerging sequencing technologies have enabled the discovery of thousands of SVs, uncovering information about their ubiquity, relationship to disease and possible effects on biological mechanisms. Given the variability in SV type and size, along with unique detection biases of emerging genomic platforms, multiplatform discovery is necessary to resolve the full spectrum of variation. Here, we review modern approaches for investigating SVs and proffer that, moving forwards, studies integrating biological information with detection will be necessary to comprehensively understand the impact of SV in the human genome.

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.

Simultaneous detection of multiple DNA damage types by multi-colour fluorescent labelling

Herein we present an assay allowing concurrent detection of oxidative DNA damage and photoproducts. We apply DNA repair enzymes specific for each lesion type to incorporate spectrally distinct fluorescent nucleotides, enabling simultaneous quantification of the lesions on individual DNA molecules. We follow the repair of both damage types in skin cells exposed to artificial sunlight.

Global modulation in DNA epigenetics during pro-inflammatory macrophage activation

DNA methylation patterns create distinct gene-expression profiles. These patterns are maintained after cell division, thus enabling the differentiation and maintenance of multiple cell types from the same genome sequence. The advantage of this mechanism for transcriptional control is that chemical-encoding allows to rapidly establish new epigenetic patterns 'on-demand' through enzymatic methylation and demethylation of DNA. Here we show that this feature is associated with the fast response of macrophages during their pro-inflammatory activation. By using a combination of mass spectroscopy and single-molecule imaging to quantify global epigenetic changes in the genomes of primary macrophages, we followed three distinct DNA marks (methylated, hydroxymethylated and unmethylated), involved in establishing new DNA methylation patterns during pro-inflammatory activation. The observed epigenetic modulation together with gene-expression data generated for the involved enzymatic machinery may suggest that de-methylation upon LPS-activation starts with oxidation of methylated CpGs, followed by excision-repair of these oxidized bases and their replacement with unmodified cytosine.