Optical mapping of DNA: Single-molecule-based methods for mapping genomes

Sterically Enhanced Control of Enzyme-Assisted DNA Assembly

Traditional methods for the assembly of functionalised DNA structures, involving enzyme restriction and modification, present difficulties when working with small DNA fragments (<100 bp), in part due to a lack of control over enzymatic action during the DNA modification process. This limits the design flexibility and range of accessible DNA structures. Here, we show that these limitations can be overcome by introducing chemical modifications into the DNA that spatially restrict enzymatic activity. This approach, sterically controlled nuclease enhanced (SCoNE) DNA assembly, thereby circumvents the size limitations of conventional Gibson assembly (GA) and allows the preparation of well-defined, functionalised DNA structures with multiple probes for specific analytes, such as IL-6, procalcitonin (PCT), and a biotin reporter group. Notably, when using the same starting materials, conventional GA under typical conditions fails. We demonstrate successful analyte capture based on standard and modified sandwich ELISA and also show how the inclusion of biotin probes provides additional functionality for product isolation.

Remote Excitation of Tip-Enhanced Photoluminescence with a Parallel AgNW Coupler

Tip-enhanced photoluminescence (TEPL) microscopy allows for the correlation of scanning probe microscopic images and photoluminescent spectra at the nanoscale level in a similar way to tip-enhanced Raman scattering (TERS) microscopy. However, due to the higher cross-section of fluorescence compared to Raman scattering, the diffraction-limited background signal generated by far-field excitation is a limiting factor in the achievable spatial resolution of TEPL. Here, we demonstrate a way to overcome this drawback by using remote excitation TEPL (RE-TEPL). With this approach, the excitation and detection positions are spatially separated, minimizing the far-field contribution. Two probe designs are evaluated, both experimentally and via simulations. The first system consists of gold nanoparticles (AuNPs) through photoinduced deposition on a silver nanowire (AgNW), and the second system consists of two offset parallel AgNWs. This latter coupler system shows a higher coupling efficiency and is used to successfully demonstrate RE-TEPL spectral mapping on a MoSe2/WSe2 lateral heterostructure to reveal spatial heterogeneity at the heterojunction.

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.

Multi intercalators FRET enhanced detection of minute amounts of DNA

A novel approach is presented that increases sensitivity and specificity for detecting minimal traces of DNA in liquid and on solid samples. Förster Resonance Energy Transfer (FRET) from YOYO to Ethidium Bromide (EtBr) substantially increases the signal from DNA-bound EtBr highly enhancing sensitivity and specificity for DNA detection. The long fluorescence lifetime of the EtBr acceptor, when bound to DNA, allows for multi-pulse pumping with time gated (MPPTG) detection, which highly increases the detectable signal of DNA-bound EtBr. A straightforward spectra/image subtraction eliminates sample background and allows for a huge increase in the overall detection sensitivity. Using a combination of FRET and MPPTG detection an amount as small as 10 pg of DNA in a microliter sample can be detected without any additional sample purification/manipulation or use of amplification technologies. This amount of DNA is comparable to the DNA content of a one to two human cells. Such a detection method based on simple optics opens the potential for robust, highly sensitive DNA detection/imaging in the field, quick evaluation/sorting (i.e., triaging) of collected DNA samples, and can support various diagnostic assays.

Genomic Approaches for the Study of Flower Development in Floriculture Crops

The advances in genomics and bioinformatics have made possible the study in non-model plants of phenotypes associated to flower development. Floriculture crops are an interesting source of traits associated to flower development such as the transition between zygomorphic and actinomorphic flowers or the production of flowers with double and triple corollas. In this chapter, we summarize the material and methods for the use of floriculture crops to study flower development using genomic tools, from the sequencing and assembly of a reference genome to QTL and RNA-Seq analysis to search candidate genes associated to specific traits.

Complex genome assembly based on long-read sequencing

High-quality genome chromosome-scale sequences provide an important basis for genomics downstream analysis, especially the construction of haplotype-resolved and complete genomes, which plays a key role in genome annotation, mutation detection, evolutionary analysis, gene function research, comparative genomics and other aspects. However, genome-wide short-read sequencing is difficult to produce a complete genome in the face of a complex genome with high duplication and multiple heterozygosity. The emergence of long-read sequencing technology has greatly improved the integrity of complex genome assembly. We review a variety of computational methods for complex genome assembly and describe in detail the theories, innovations and shortcomings of collapsed, semi-collapsed and uncollapsed assemblers based on long reads. Among the three methods, uncollapsed assembly is the most correct and complete way to represent genomes. In addition, genome assembly is closely related to haplotype reconstruction, that is uncollapsed assembly realizes haplotype reconstruction, and haplotype reconstruction promotes uncollapsed assembly. We hope that gapless, telomere-to-telomere and accurate assembly of complex genomes can be truly routinely achieved using only a simple process or a single tool in the future.

Förster Resonance Energy Transfer-Enhanced Detection of Minute Amounts of DNA

This article presents a novel approach to increase the detection sensitivity of trace amounts of DNA in a sample by employing Förster resonance energy transfer (FRET) between intercalating dyes. Two intercalators that present efficient FRET were used to enhance sensitivity and improve specificity in detecting minute amounts of DNA. Comparison of steady-state acceptor emission spectra with and without the donor allows for simple and specific detection of DNA (acceptor bound to DNA) down to 100 pg/μL. When utilizing as an acceptor a dye with a significantly longer lifetime (e.g., ethidium bromide bound to DNA), multipulse pumping and time-gated detection enable imaging/visualization of picograms of DNA present in a microliter of an unprocessed sample or DNA collected on a swab or other substrate materials.

Combining dense and sparse labeling in optical DNA mapping

Optical DNA mapping (ODM) is based on fluorescent labeling, stretching and imaging of single DNA molecules to obtain sequence-specific fluorescence profiles, DNA barcodes. These barcodes can be mapped to theoretical counterparts obtained from DNA reference sequences, which in turn allow for DNA identification in complex samples and for detecting structural changes in individual DNA molecules. There are several types of DNA labeling schemes for ODM and for each labeling type one or several types of match scoring methods are used. By combining the information from multiple labeling schemes one can potentially improve mapping confidence; however, combining match scores from different labeling assays has not been implemented yet. In this study, we introduce two theoretical methods for dealing with analysis of DNA molecules with multiple label types. In our first method, we convert the alignment scores, given as output from the different assays, into p-values using carefully crafted null models. We then combine the p-values for different label types using standard methods to obtain a combined match score and an associated combined p-value. In the second method, we use a block bootstrap approach to check for the uniqueness of a match to a database for all barcodes matching with a combined p-value below a predefined threshold. For obtaining experimental dual-labeled DNA barcodes, we introduce a novel assay where we cut plasmid DNA molecules from bacteria with restriction enzymes and the cut sites serve as sequence-specific markers, which together with barcodes obtained using the established competitive binding labeling method, form a dual-labeled barcode. All experimental data in this study originates from this assay, but we point out that our theoretical framework can be used to combine data from all kinds of available optical DNA mapping assays. We test our multiple labeling frameworks on barcodes from two different plasmids and synthetically generated barcodes (combined competitive-binding- and nick-labeling). It is demonstrated that by simultaneously using the information from all label types, we can substantially increase the significance when we match experimental barcodes to a database consisting of theoretical barcodes for all sequenced plasmids.

Assemblies of the genomes of parasitic wasps using meta-assembly and scaffolding with genetic linkage

Safe, effective biological-control introductions against invasive pests depend on narrowly host-specific natural enemies with the ability to adapt to a changing environment. As part of a project on the genetic architectures of these traits, we assembled and annotated the genomes of two aphid parasitoids, Aphelinus atriplicis and Aphelinus certus. We report here several assemblies of A. atriplicis made with Illumina and PacBio data, which we combined into a meta-assembly. We scaffolded the meta-assembly with markers from a genetic map of hybrids between A. atriplicis and A. certus. We used this genetic-linkage scaffolded (GLS) assembly of A. atriplicis to scaffold a de novo assembly of A. certus. The de novo assemblies of A. atriplicis differed in contiguity, and the meta-assembly of these assemblies was more contiguous than the best de novo assembly. Scaffolding with genetic-linkage data allowed chromosomal-level assembly of the A. atriplicis genome and scaffolding a de novo assembly of A. certus with this GLS assembly, greatly increased the contiguity of the A. certus assembly to the point where it was also at the chromosomal-level. However, completeness of the A. atriplicis assembly, as measured by percent complete, single-copy BUSCO hymenopteran genes, varied little among de novo assemblies and was not increased by meta-assembly or genetic scaffolding. Furthermore, the greater contiguity of the meta-assembly and GLS assembly had little or no effect on the numbers of genes identified, the proportions with homologs or functional annotations. Increased contiguity of the A. certus assembly provided modest improvement in assembly completeness, as measured by percent complete, single-copy BUSCO hymenopteran genes. The total genic sequence increased, and while the number of genes declined, gene length increased, which together suggest greater accuracy of gene models. More contiguous assemblies provide uses other than gene annotation, for example, identifying the genes associated with quantitative trait loci and understanding of chromosomal rearrangements associated with speciation.

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.

Electrical DNA Sequence Mapping Using Oligodeoxynucleotide Labels and Nanopores

Identifying DNA species is crucial for diagnostics. For DNA identification, single-molecule DNA sequence mapping is an alternative to DNA sequencing toward fast point-of-care testing, which traditionally relies on targeting and labeling DNA sequences with fluorescent labels and readout using optical imaging methods. A nanopore is a promising sensor as a complement to optical mapping with advantages of electric measurement suitable for portable devices and potential for high resolution. Here, we demonstrate a high-resolution nanopore-based DNA sequence mapping by labeling specific short sequence motifs with oligodeoxynucleotides (ODNs) using DNA methyltransferase (MTase) and detecting them using nanopores. We successfully detected ODNs down to the size of 11 nucleotides without introducing extra reporters and resolved neighboring sites with a distance of 141 bp (∼48 nm) on a single DNA molecule. To accurately locate the sequence motif positions on DNA, a nanopore data analysis method is proposed by considering DNA velocity change through nanopores and using ensemble statistics to translate the time-dependent signals to the location information. Our platform enables high-resolution detection of small labels on DNA and high-accuracy localization of them for DNA species identification in an all-electrical format. The method presents an alternative to optical techniques relying on fluorescent labels and is promising for miniature-scale integration for diagnostic applications.

CRISPR-based diagnostics for detection of pathogens

The improved sensitivity and superior specificity associated with the use of molecular assays has improved the fate of disease diagnosis by bestowing the clinicians with outcomes that are both rapid and precise. In recent years, CRISPR has made considerable progress in in vitro diagnostic platform which has paved its way for developing rapid and sensitive CRISPR-based diagnostic tools. Improved perception and better understanding of diverse CRISPR-Cas systems has broadened the reach of CRISPR applications for not just early detection of pathogens but also for early onset of diseases such as cancer. The inherent allele specificity of CRISPR is the predominant reason for its application in designing a diagnostic-tool that is field-deployable, portable, sensitive, specific and rapid. In this chapter, we highlight various CRISPR-based diagnostic platforms, its applications, challenges and future prospects of the CRISPR-Cas system.

DNA Manipulation and Single-Molecule Imaging

DNA replication, repair, and recombination in the cell play a significant role in the regulation of the inheritance, maintenance, and transfer of genetic information. To elucidate the biomolecular mechanism in the cell, some molecular models of DNA replication, repair, and recombination have been proposed. These biological studies have been conducted using bulk assays, such as gel electrophoresis. Because in bulk assays, several millions of biomolecules are subjected to analysis, the results of the biological analysis only reveal the average behavior of a large number of biomolecules. Therefore, revealing the elementary biological processes of a protein acting on DNA (e.g., the binding of protein to DNA, DNA synthesis, the pause of DNA synthesis, and the release of protein from DNA) is difficult. Single-molecule imaging allows the analysis of the dynamic behaviors of individual biomolecules that are hidden during bulk experiments. Thus, the methods for single-molecule imaging have provided new insights into almost all of the aspects of the elementary processes of DNA replication, repair, and recombination. However, in an aqueous solution, DNA molecules are in a randomly coiled state. Thus, the manipulation of the physical form of the single DNA molecules is important. In this review, we provide an overview of the unique studies on DNA manipulation and single-molecule imaging to analyze the dynamic interaction between DNA and protein.

Aligning optical maps to de Bruijn graphs

MOTIVATION

Optical maps are high-resolution restriction maps (Rmaps) that give a unique numeric representation to a genome. Used in concert with sequence reads, they provide a useful tool for genome assembly and for discovering structural variations and rearrangements. Although they have been a regular feature of modern genome assembly projects, optical maps have been mainly used in post-processing step and not in the genome assembly process itself. Several methods have been proposed for pairwise alignment of single molecule optical maps-called Rmaps, or for aligning optical maps to assembled reads. However, the problem of aligning an Rmap to a graph representing the sequence data of the same genome has not been studied before. Such an alignment provides a mapping between two sets of data: optical maps and sequence data which will facilitate the usage of optical maps in the sequence assembly step itself.

RESULTS

We define the problem of aligning an Rmap to a de Bruijn graph and present the first algorithm for solving this problem which is based on a seed-and-extend approach. We demonstrate that our method is capable of aligning 73% of Rmaps generated from the Escherichia coli genome to the de Bruijn graph constructed from short reads generated from the same genome. We validate the alignments and show that our method achieves an accuracy of 99.6%. We also show that our method scales to larger genomes. In particular, we show that 76% of Rmaps can be aligned to the de Bruijn graph in the case of human data.

AVAILABILITY AND IMPLEMENTATION

The software for aligning optical maps to de Bruijn graph, omGraph is written in C++ and is publicly available under GNU General Public License at https://github.com/kingufl/omGraph.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cultivation-Free Typing of Bacteria Using Optical DNA Mapping

A variety of pathogenic bacteria can infect humans, and rapid species identification is crucial for the correct treatment. However, the identification process can often be time-consuming and depend on the cultivation of the bacterial pathogen(s). Here, we present a stand-alone, enzyme-free, optical DNA mapping assay capable of species identification by matching the intensity profiles of large DNA molecules to a database of fully assembled bacterial genomes (>10 000). The assay includes a new data analysis strategy as well as a general DNA extraction protocol for both Gram-negative and Gram-positive bacteria. We demonstrate that the assay is capable of identifying bacteria directly from uncultured clinical urine samples, as well as in mixtures, with the potential to be discriminative even at the subspecies level. We foresee that the assay has applications both within research laboratories and in clinical settings, where the time-consuming step of cultivation can be minimized or even completely avoided.

Identifying microbial species by single-molecule DNA optical mapping and resampling statistics

Single-molecule DNA mapping has the potential to serve as a powerful complement to high-throughput sequencing in metagenomic analysis. Offering longer read lengths and forgoing the need for complex library preparation and amplification, mapping stands to provide an unbiased view into the composition of complex viromes and/or microbiomes. To fully enable mapping-based metagenomics, sensitivity and specificity of DNA map analysis and identification need to be improved. Using detailed simulations and experimental data, we first demonstrate how fluorescence imaging of surface stretched, sequence specifically labeled DNA fragments can yield highly sensitive identification of targets. Second, a new analysis technique is introduced to increase specificity of the analysis, allowing even closely related species to be resolved. Third, we show how an increase in resolution improves sensitivity. Finally, we demonstrate that these methods are capable of identifying species with long genomes such as bacteria with high sensitivity.

DNA barcodes for rapid, whole genome, single-molecule analyses

We report an approach for visualizing DNA sequence and using these ‘DNA barcodes’ to search complex mixtures of genomic material for DNA molecules of interest. We demonstrate three applications of this methodology; identifying specific molecules of interest from a dataset containing gigabasepairs of genome; identification of a bacterium from such a dataset and, finally, by locating infecting virus molecules in a background of human genomic material. As a result of the dense fluorescent labelling of the DNA, individual barcodes of the order 40 kb pairs in length can be reliably identified. This means DNA can be prepared for imaging using standard handling and purification techniques. The recorded dataset provides stable physical and electronic records of the total genomic content of a sample that can be readily searched for a molecule or region of interest.

Fluorescent SAM analogues for methyltransferase based DNA labeling

In this work, the preparation of new S-adenosyl-l-methionine (SAM) analogues for sequence specific DNA labeling is evaluated. Fluorescent cofactors were synthesized and their applicability in methyltransferase based optical mapping is demonstrated.

A chromosome-level genome assembly of Pyropia haitanensis (Bangiales, Rhodophyta)

Pyropia haitanensis (Bangiales, Rhodophyta), a major economically important marine crop, is also considered as an ideal research model of Rhodophyta to address several major biological questions such as sexual reproduction and adaptation to intertidal abiotic stresses. However, comparative genomic analysis to decipher the underlying molecular mechanisms is hindered by the lack of high-quality genome information. Therefore, we integrated sequencing data from Illumina short-read sequencing, PacBio single-molecule sequencing and BioNano optical genome mapping. The assembled genome was approximately 53.3 Mb with an average GC% of 67.9%. The contig N50 and scaffold N50 were 510.3 kb and 5.8 Mb, respectively. Additionally, 10 superscaffolds representing 80.9% of the total assembly (42.7 Mb) were anchored and orientated to the 5 linkage groups based on markers and genetic distance; this outcome is consistent with the karyotype of five chromosomes (n = 5) based on cytological observation in P. haitanensis. Approximately 9.6% and 14.6% of the genomic region were interspersed repeat and tandem repeat elements, respectively. Based on full-length transcriptome data generated by PacBio, 10,903 protein-coding genes were identified. The construction of a genome-wide phylogenetic tree demonstrated that the divergence time of P. haitanensis and Porphyra umbilicalis was ~204.4 Ma. Interspecies comparison revealed that 493 gene families were expanded and that 449 were contracted in the P. haitanensis genome compared with those in the Po. umbilicalis genome. The genome identified is of great value for further research on the genome evolution of red algae and genetic adaptation to intertidal stresses.

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.

Modelling BioNano optical data and simulation study of genome map assembly

Motivation

The launch of the BioNano next-generation mapping system has greatly enhanced the performance of physical map construction, thus rapidly expanding the application of optical mapping in genome research. Data biases have profound implications for downstream applications. However, very little is known about the properties and biases of BioNano data, and the very factors that contribute to whole-genome optical map assembly.

Results

We generated BioNano molecule data from eight organisms with diverse base compositions. We first characterized the properties/biases of BioNano molecule data, i.e. molecule length distribution, false labelling signal, variation of optical resolution and coverage distribution bias, and their inducing factors such as chimeric molecules, fragile sites and DNA molecule stretching. Second, we developed the BioNano Molecule SIMulator (BMSIM), a novel computer simulation program for optical data. BMSIM, is of great use for future genome mapping projects. Third, we evaluated the experimental variables that impact whole-genome optical map assembly. Specifically, the effects of coverage depth, molecule length, false-positive and false-negative labelling signals, chimeric molecules and nicking enzyme and nick site density were investigated. Our simulation study provides the empirical findings on how to control experimental variables and gauge analytical parameters to maximize benefit and minimize cost on whole-genome optical map assembly.

Availability and implementation

BMSIM is freely available on: https://github.com/pingchen09990102/BMSIM.

Supplementary information

Supplementary data are available at Bioinformatics online.

A general strategy for direct, enzyme-catalyzed conjugation of functional compounds to DNA

The methyltransferase enzymes can be applied to deliver a range of modifications to pre-determined sites on large DNA molecules with exceptional specificity and efficiency. To date, however, a limited number of modifications have been delivered in this way because of the complex chemical synthesis that is needed to produce a cofactor analogue carrying a specific function, such as a fluorophore. Here, we describe a method for the direct transfer of a series of functional compounds (seven fluorescent dyes, biotin and polyethylene glycol) to the DNA duplex. Our approach uses a functional cofactor analogue, whose final preparative step is performed alongiside the DNA modification reaction in a single pot, with no purification needed. We show that fluorophore conjugation efficiency in these mixtures is significantly improved compared to two-step labeling approaches. Our experiments highlight the remarkable malleability and selectivity of the methyltransferases tested. Additional analysis using high resolution localization of the fluorophore distribution indicates that target sites for the methyltransferase are predominantly labeled on a single strand of their palindromic site and that a small and randomly-distributed probability of off-site labeling exists.

Overview of Next-Generation Sequencing Technologies

High throughput DNA sequencing methodology (next generation sequencing; NGS) has rapidly evolved over the past 15 years and new methods are continually being commercialized. As the technology develops, so do increases in the number of corresponding applications for basic and applied science. The purpose of this review is to provide a compendium of NGS methodologies and associated applications. Each brief discussion is followed by web links to the manufacturer and/or web-based visualizations. Keyword searches, such as with Google, may also provide helpful internet links and information. © 2018 by John Wiley & Sons, Inc.

A Near Chromosome Assembly of the Dromedary Camel Genome

The dromedary camel is an economically and socially important species of livestock in many parts of the world, being used for transport and the production of milk and meat. Much like cattle and horses, the camel may be found in industrial farming conditions as well as used in sporting. Camel racing is a multi-million dollar industry, with some specimens being valued at upward of 9.5 million USD. Despite its apparent value to humans, the dromedary camel is a neglected species in genomics. While cattle and other domesticated species have had much attention in terms of genome assembly, the camel has only been assembled to scaffold level, which does not give a clear indication of the order or chromosomal location of sequenced fragments. In this study, the Reference Assistant Chromosome Assembly (RACA) algorithm was implemented to use read-pair information of camel scaffolds, aligned with the cattle and human genomes in order to organize and orient these scaffolds in a near-chromosome level assembly. This method generated 72 large size fragments (N50 54.36 Mb). These predicted chromosome fragments (PCFs) were then compared with comparative maps of camel and cytogenetic map of alpaca chromosomes, allowing us to further upgrade the assembly. This dromedary camel assembly will be an invaluable tool to verify future camel assemblies generated with chromatin conformation or/and long read technologies. This study provides the first near-chromosome assembly of the dromedary camel, thus adding this economically important species to a growing pool of knowledge regarding the genome structure of domesticated livestock.

Shining a Spotlight on DNA: Single-Molecule Methods to Visualise DNA

The ability to watch single molecules of DNA has revolutionised how we study biological transactions concerning nucleic acids. Many strategies have been developed to manipulate DNA molecules to investigate mechanical properties, dynamics and protein–DNA interactions. Imaging methods using small molecules and protein-based probes to visualise DNA have propelled our understanding of complex biochemical reactions involving DNA. This review focuses on summarising some of the methodological developments made to visualise individual DNA molecules and discusses how these probes have been used in single-molecule biophysical assays.

Advances in detection and quantification of methylcytosine and its derivatives

Methylation of the fifth carbon atom in cytosine is an epigenetic modification of deoxyribonucleic acid that plays important roles in numerous cellular processes and disease pathogenesis. Three additional states of cytosine, that is, 5-hydroxymethylcytosine, 5-formylcytosine and 5-carboxylcytosine, have been identified and associated with the diagnosis and/or prognosis of diseases. However, accurate measurement of those intermediates is a challenge since their global levels are relatively low. A number of innovative methods have been developed to detect and quantify these compounds in biological samples, such as blood, tissue and urine, etc. This review focuses on recent advancement in detection and quantification of four cytosine modifications, based on which, the development, diagnosis, and prognosis of diseases could be monitored through non-invasive procedures.

Entropic trap purification of long DNA

Long-read genomic applications, such as genome mapping in nanochannels, require long DNA that is free of small-DNA impurities. We have developed a chip-based system based on entropic trapping that can simultaneously concentrate and purify a long DNA sample under the alternating application of an applied pressure (for sample injection) and an electric field (for filtration and concentration). In contrast, short DNA tends to pass through the filter owing to its comparatively weak entropic penalty for entering the nanoslit. The single-stage prototype developed here, which operates in a continuous pulsatile manner, achieves selectivities of up to 3.5 for λ-phage DNA (48.5 kilobase pairs) compared to a 2 kilobase pair standard based on experimental data for the fraction filtered using pure samples of each species. The device is fabricated in fused silica using standard clean-room methods, making it compatible for integration with long-read genomics technologies.

By land, air, and sea: hemipteran diversity through the genomic lens

Thanks to a recent spate of sequencing projects, the Hemiptera are the first hemimetabolous insect order to achieve a critical mass of species with sequenced genomes, establishing the basis for comparative genomics of the bugs. However, as the most speciose hemimetabolous order, there is still a vast swathe of the hemipteran phylogeny that awaits genomic representation across subterranean, terrestrial, and aquatic habitats, and with lineage-specific and developmentally plastic cases of both wing polyphenisms and flightlessness. In this review, we highlight opportunities for taxonomic sampling beyond obvious pest species candidates, motivated by intriguing biological features of certain groups as well as the rich research tradition of ecological, physiological, developmental, and particularly cytogenetic investigation that spans the diversity of the Hemiptera.

From Short Reads to Chromosome-Scale Genome Assemblies

A high-quality, annotated genome assembly is the foundation for many downstream studies. However, obtaining such an assembly is a complex, reiterative process that requires the assimilation of high-quality data and combines different approaches and data types. While some software packages incorporating multiple steps of genome assembly are commercially available, they may not be flexible enough to be routinely applied to all organisms, particularly to nonmodel species such as pathogenic oomycetes and fungi. If researchers understand and apply the most appropriate, currently available tools for each step, it is possible to customize parameters and optimize results for their organism of study. Based on our experience of de novo assembly and annotation of several oomycete species, this chapter provides a modular workflow from processing of raw reads, to initial assembly generation, through optimization, chromosome-scale scaffolding and annotation, outlining input and output data as well as examples and alternative software used for each step. The accompanying Notes provide background information for each step as well as alternative options. The final result of this workflow could be an annotated, high-quality, validated, chromosome-scale assembly or a draft assembly of sufficient quality to meet specific needs of a project.

Enrichment of megabase-sized DNA molecules for single-molecule optical mapping and next-generation sequencing

Next-generation sequencing (NGS) has caused a revolution, yet left a gap: long-range genetic information from native, non-amplified DNA fragments is unavailable. It might be obtained by optical mapping of megabase-sized DNA molecules. Frequently only a specific genomic region is of interest, so here we introduce a method for selection and enrichment of megabase-sized DNA molecules intended for single-molecule optical mapping: DNA from a human cell line is digested by the NotI rare-cutting enzyme and size-selected by pulsed-field gel electrophoresis. For demonstration, more than 600 sub-megabase- to megabase-sized DNA molecules were recovered from the gel and analysed by denaturation-renaturation optical mapping. Size-selected molecules from the same gel were sequenced by NGS. The optically mapped molecules and the NGS reads showed enrichment from regions defined by NotI restriction sites. We demonstrate that the unannotated genome can be characterized in a locus-specific manner via molecules partially overlapping with the annotated genome. The method is a promising tool for investigation of structural variants in enriched human genomic regions for both research and diagnostic purposes. Our enrichment method could potentially work with other genomes or target specified regions by applying other genomic editing tools, such as the CRISPR/Cas9 system.

Ionic Current-Based Mapping of Short Sequence Motifs in Single DNA Molecules Using Solid-State Nanopores

Nanopore sensors show great potential for rapid, single-molecule determination of DNA sequence information. Here, we develop an ionic current-based method for determining the positions of short sequence motifs in double-stranded DNA molecules with solid-state nanopores. Using the DNA-methyltransferase M.TaqI and a biotinylated S-adenosyl-l-methionine cofactor analogue we create covalently attached biotin labels at 5'-TCGA-3' sequence motifs. Monovalent streptavidin is then added to bind to the biotinylated sites giving rise to additional current blockade signals when the DNA passes through a conical quartz nanopore. We determine the relationship between translocation time and position along the DNA contour and find a minimum resolvable distance between two labeled sites of ∼200 bp. We then characterize a variety of DNA molecules by determining the positions of bound streptavidin and show that two short genomes can be simultaneously detected in a mixture. Our method provides a simple, generic single-molecule detection platform enabling DNA characterization in an electrical format suited for portable devices for potential diagnostic applications.

Noise reduction in single time frame optical DNA maps

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

Applications of optical DNA mapping in microbiology

Optical mapping (OM) has been used in microbiology for the past 20 years, initially as a technique to facilitate DNA sequence-based studies; however, with decreases in DNA sequencing costs and increases in sequence output from automated sequencing platforms, OM has grown into an important auxiliary tool for genome assembly and comparison. Currently, there are a number of new and exciting applications for OM in the field of microbiology, including investigation of disease outbreaks, identification of specific genes of clinical and/or epidemiological relevance, and the possibility of single-cell analysis when combined with cell-sorting approaches. In addition, designing lab-on-a-chip systems based on OM is now feasible and will allow the integrated and automated microbiological analysis of biological fluids. Here, we review the basic technology of OM, detail the current state of the art of the field, and look ahead to possible future developments in OM technology for microbiological applications.

Whole-Genome Restriction Mapping by "Subhaploid"-Based RAD Sequencing: An Efficient and Flexible Approach for Physical Mapping and Genome Scaffolding

Assembly of complex genomes using short reads remains a major challenge, which usually yields highly fragmented assemblies. Generation of ultradense linkage maps is promising for anchoring such assemblies, but traditional linkage mapping methods are hindered by the infrequency and unevenness of meiotic recombination that limit attainable map resolution. Here we develop a sequencing-based "in vitro" linkage mapping approach (called RadMap), where chromosome breakage and segregation are realized by generating hundreds of "subhaploid" fosmid/bacterial-artificial-chromosome clone pools, and by restriction site-associated DNA sequencing of these clone pools to produce an ultradense whole-genome restriction map to facilitate genome scaffolding. A bootstrap-based minimum spanning tree algorithm is developed for grouping and ordering of genome-wide markers and is implemented in a user-friendly, integrated software package (AMMO). We perform extensive analyses to validate the power and accuracy of our approach in the model plant Arabidopsis thaliana and human. We also demonstrate the utility of RadMap for enhancing the contiguity of a variety of whole-genome shotgun assemblies generated using either short Illumina reads (300 bp) or long PacBio reads (6-14 kb), with up to 15-fold improvement of N50 (∼816 kb-3.7 Mb) and high scaffolding accuracy (98.1-98.5%). RadMap outperforms BioNano and Hi-C when input assembly is highly fragmented (contig N50 = 54 kb). RadMap can capture wide-range contiguity information and provide an efficient and flexible tool for high-resolution physical mapping and scaffolding of highly fragmented assemblies.

Methyltransferase-Directed Labeling of Biomolecules and its Applications

Methyltransferases (MTases) form a large family of enzymes that methylate a diverse set of targets, ranging from the three major biopolymers to small molecules. Most of these MTases use the cofactor S-adenosyl-l-Methionine (AdoMet) as a methyl source. In recent years, there have been significant efforts toward the development of AdoMet analogues with the aim of transferring moieties other than simple methyl groups. Two major classes of AdoMet analogues currently exist: doubly-activated molecules and aziridine based molecules, each of which employs a different approach to achieve transalkylation rather than transmethylation. In this review, we discuss the various strategies for labelling and functionalizing biomolecules using AdoMet-dependent MTases and AdoMet analogues. We cover the synthetic routes to AdoMet analogues, their stability in biological environments and their application in transalkylation reactions. Finally, some perspectives are presented for the potential use of AdoMet analogues in biology research, (epi)genetics and nanotechnology.

Optical DNA mapping in nanofluidic devices: principles and applications

Optical DNA mapping has over the last decade emerged as a very powerful tool for obtaining long range sequence information from single DNA molecules. In optical DNA mapping, intact large single DNA molecules are labeled, stretched out, and imaged using a fluorescence microscope. This means that sequence information ranging over hundreds of kilobasepairs (kbp) can be obtained in one single image. Nanochannels offer homogeneous and efficient stretching of DNA that is crucial to maximize the information that can be obtained from optical DNA maps. In this review, we highlight progress in the field of optical DNA mapping in nanochannels. We discuss the different protocols for sequence specific labeling and divide them into two main categories, enzymatic labeling and affinity-based labeling. Examples are highlighted where optical DNA mapping is used to gain information on length scales that would be inaccessible with traditional techniques. Enzymatic labeling has been commercialized and is mainly used in human genetics and assembly of complex genomes, while the affinity-based methods have primarily been applied in bacteriology, for example for rapid analysis of plasmids encoding antibiotic resistance. Next, we highlight how the design of nanofluidic channels can been altered in order to obtain the desired information and discuss how recent advances in the field make it possible to retrieve information beyond DNA sequence. In the outlook section, we discuss future directions of optical DNA mapping, such as fully integrated devices and portable microscopes.

Microfluidic DNA Stretching Device for Single-Molecule Diagnostics

The method described here enables the automatic stretching and patterning of single DNA molecules onto a solid surface. It does not require chemical modification of the DNA or surface modification of the substrate. To detect a signal variation caused by sequence-specific dye binding or partial melting, it is crucial that the DNA molecules are arrayed in a parallel direction inside the narrow microscopic field. The method uses zigzag-shaped microgrooves in a densely-arranged molecular patterning apparatus in a microfluidic channel. By syringing through the microchannel, over 1500 DNA molecules can be arrayed simultaneously in the microgrooves. It will therefore serve as a template preparation for DNA molecular diagnosis by high-resolution imaging.

Genetic Drivers of Multidrug Resistance in Candida glabrata

Both the incidence of invasive fungal infections and rates of multidrug resistance associated with fungal pathogen Candida glabrata have increased in recent years. In this perspective, we will discuss the mechanisms underlying the capacity of C. glabrata to rapidly develop resistance to multiple drug classes, including triazoles and echinocandins. We will focus on the extensive genetic diversity among clinical isolates of C. glabrata, which likely enables this yeast to survive multiple stressors, such as immune pressure and antifungal exposure. In particular, over half of C. glabrata clinical strains collected from U.S. and non-U.S. sites have mutations in the DNA mismatch repair gene MSH2, leading to a mutator phenotype and increased frequencies of drug-resistant mutants in vitro. Furthermore, recent studies and data presented here document extensive chromosomal rearrangements among C. glabrata strains, resulting in a large number of distinct karyotypes within a single species. By analyzing clonal, serial isolates derived from individual patients treated with antifungal drugs, we were able to document chromosomal changes occurring in C. glabrata in vivo during the course of antifungal treatment. Interestingly, we also show that both MSH2 genotypes and chromosomal patterns cluster consistently into specific strain types, indicating that C. glabrata has a complex population structure where genomic variants arise, perhaps during the process of adaptation to environmental changes, and persist over time.

Direct identification of antibiotic resistance genes on single plasmid molecules using CRISPR/Cas9 in combination with optical DNA mapping

Bacterial plasmids are extensively involved in the rapid global spread of antibiotic resistance. We here present an assay, based on optical DNA mapping of single plasmids in nanofluidic channels, which provides detailed information about the plasmids present in a bacterial isolate. In a single experiment, we obtain the number of different plasmids in the sample, the size of each plasmid, an optical barcode that can be used to identify and trace the plasmid of interest and information about which plasmid that carries a specific resistance gene. Gene identification is done using CRISPR/Cas9 loaded with a guide-RNA (gRNA) complementary to the gene of interest that linearizes the circular plasmids at a specific location that is identified using the optical DNA maps. We demonstrate the principle on clinically relevant extended spectrum beta-lactamase (ESBL) producing isolates. We discuss how the gRNA sequence can be varied to obtain the desired information. The gRNA can either be very specific to identify a homogeneous group of genes or general to detect several groups of genes at the same time. Finally, we demonstrate an example where we use a combination of two gRNA sequences to identify carbapenemase-encoding genes in two previously not characterized clinical bacterial samples.

Genome Mapping in Plant Comparative Genomics

Genome mapping produces fingerprints of DNA sequences to construct a physical map of the whole genome. It provides contiguous, long-range information that complements and, in some cases, replaces sequencing data. Recent advances in genome-mapping technology will better allow researchers to detect large (>1kbp) structural variations between plant genomes. Some molecular and informatics complications need to be overcome for this novel technology to achieve its full utility. This technology will be useful for understanding phenotype responses due to DNA rearrangements and will yield insights into genome evolution, particularly in polyploids. In this review, we outline recent advances in genome-mapping technology, including the processes required for data collection and analysis, and applications in plant comparative genomics.

Generalized recovery algorithm for 3D super-resolution microscopy using rotating point spread functions

Super-resolution microscopy with phase masks is a promising technique for 3D imaging and tracking. Due to the complexity of the resultant point spread functions, generalized recovery algorithms are still missing. We introduce a 3D super-resolution recovery algorithm that works for a variety of phase masks generating 3D point spread functions. A fast deconvolution process generates initial guesses, which are further refined by least squares fitting. Overfitting is suppressed using a machine learning determined threshold. Preliminary results on experimental data show that our algorithm can be used to super-localize 3D adsorption events within a porous polymer film and is useful for evaluating potential phase masks. Finally, we demonstrate that parallel computation on graphics processing units can reduce the processing time required for 3D recovery. Simulations reveal that, through desktop parallelization, the ultimate limit of real-time processing is possible. Our program is the first open source recovery program for generalized 3D recovery using rotating point spread functions.

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

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

Rapid Tracing of Resistance Plasmids in a Nosocomial Outbreak Using Optical DNA Mapping

Resistance to life-saving antibiotics increases rapidly worldwide, and multiresistant bacteria have become a global threat to human health. Presently, the most serious threat is the increasing spread of Enterobacteriaceae carrying genes coding for extended spectrum β-lactamases (ESBL) and carbapenemases on highly mobile plasmids. We here demonstrate how optical DNA maps of single plasmids can be used as fingerprints to trace plasmids, for example, during resistance outbreaks. We use the assay to demonstrate a potential transmission route of an ESBL-carrying plasmid between bacterial strains/species and between patients, during a polyclonal outbreak at a neonatal ward at Sahlgrenska University Hospital (Gothenburg, Sweden). Our results demonstrate that optical DNA mapping is an easy and rapid method for detecting the spread of plasmids mediating resistance. With the increasing prevalence of multiresistant bacteria, diagnostic tools that can aid in solving ongoing routes of transmission, in particular in hospital settings, will be of paramount importance.

Whole-genome optical mapping reveals a mis-assembly between two rRNA operons of Corynebacterium pseudotuberculosis strain 1002

Background

Studies have detected mis-assemblies in genomes of the species Corynebacterium pseudotuberculosis. These new discover have been possible due to the evolution of the Next-Generation Sequencing platforms, which have provided sequencing with accuracy and reduced costs. In addition, the improving of techniques for construction of high accuracy genomic maps, for example, Whole-genome mapping (WGM) (OpGen Inc), have allow high-resolution assembly that can detect large rearrangements.

Results

In this work, we present the resequencing of Corynebacterium pseudotuberculosis strain 1002 (Cp1002). Cp1002 was the first strain of this species sequenced in Brazil, and its genome has been used as model for several studies in silico of caseous lymphadenitis disease. The sequencing was performed using the platform Ion PGM and fragment library (200 bp kit). A restriction map was constructed, using the technique of WGM with the enzyme KpnI. After the new assembly process, using WGM as scaffolder, we detected a large inversion with size bigger than one-half of genome. A specific analysis using BLAST and NR database shows that the inversion occurs between two homology RNA ribosomal regions.

Conclusion

In conclusion, the results showed by WGM could be used to detect mismatches in assemblies, providing genomic maps with high resolution and allow assemblies with more accuracy and completeness. The new assembly of C. pseudotuberculosis was deposited in GenBank under the accession no. CP012837.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2673-7) contains supplementary material, which is available to authorized users.

How to determine local stretching and tension in a flow-stretched DNA molecule

We determine the nonuniform stretching of and tension in a mega base pairs-long fragment of deoxyribonucleic acid (DNA) that is flow stretched in a nanofluidic chip. We use no markers, do not know the contour length of the DNA, and do not have the full DNA molecule inside our field of view. Instead, we analyze the transverse thermal motion of the DNA. Tension at the center of the DNA adds up to 16 pN, giving almost fully stretched DNA. This method was devised for optical mapping of DNA, specifically, DNA denaturation patterns. It may be useful also for other studies, e.g., DNA-protein interactions, specifically, their tension dependence. Generally, wherever long strands of DNA-e.g., native DNA extracted from human cells or bacteria-must be stretched with ease for inspection, this method applies.

Improvement of the banana "Musa acuminata" reference sequence using NGS data and semi-automated bioinformatics methods

Background

Recent advances in genomics indicate functional significance of a majority of genome sequences and their long range interactions. As a detailed examination of genome organization and function requires very high quality genome sequence, the objective of this study was to improve reference genome assembly of banana (Musa acuminata).

Results

We have developed a modular bioinformatics pipeline to improve genome sequence assemblies, which can handle various types of data. The pipeline comprises several semi-automated tools. However, unlike classical automated tools that are based on global parameters, the semi-automated tools proposed an expert mode for a user who can decide on suggested improvements through local compromises. The pipeline was used to improve the draft genome sequence of Musa acuminata. Genotyping by sequencing (GBS) of a segregating population and paired-end sequencing were used to detect and correct scaffold misassemblies. Long insert size paired-end reads identified scaffold junctions and fusions missed by automated assembly methods. GBS markers were used to anchor scaffolds to pseudo-molecules with a new bioinformatics approach that avoids the tedious step of marker ordering during genetic map construction. Furthermore, a genome map was constructed and used to assemble scaffolds into super scaffolds. Finally, a consensus gene annotation was projected on the new assembly from two pre-existing annotations. This approach reduced the total Musa scaffold number from 7513 to 1532 (i.e. by 80 %), with an N50 that increased from 1.3 Mb (65 scaffolds) to 3.0 Mb (26 scaffolds). 89.5 % of the assembly was anchored to the 11 Musa chromosomes compared to the previous 70 %. Unknown sites (N) were reduced from 17.3 to 10.0 %.

Conclusion

The release of the Musa acuminata reference genome version 2 provides a platform for detailed analysis of banana genome variation, function and evolution. Bioinformatics tools developed in this work can be used to improve genome sequence assemblies in other species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2579-4) contains supplementary material, which is available to authorized users.

Effects of genome structure variation, homeologous genes and repetitive DNA on polyploid crop research in the age of genomics

Compared to diploid species, allopolyploid crop species possess more complex genomes, higher productivity, and greater adaptability to changing environments. Next generation sequencing techniques have produced high-density genetic maps, whole genome sequences, transcriptomes and epigenomes for important polyploid crops. However, several problems interfere with the full application of next generation sequencing techniques to these crops. Firstly, different types of genomic variation affect sequence assembly and QTL mapping. Secondly, duplicated or homoeologous genes can diverge in function and then lead to emergence of many minor QTL, which increases difficulties in fine mapping, cloning and marker assisted selection. Thirdly, repetitive DNA sequences arising in polyploid crop genomes also impact sequence assembly, and are increasingly being shown to produce small RNAs to regulate gene expression and hence phenotypic traits. We propose that these three key features should be considered together when analyzing polyploid crop genomes. It is apparent that dissection of genomic structural variation, elucidation of the function and mechanism of interaction of homoeologous genes, and investigation of the de novo roles of repeat sequences in agronomic traits are necessary for genomics-based crop breeding in polyploids.