Optical Restriction Mapping Research Articles

OPTIMA: sensitive and accurate whole-genome alignment of error-prone genomic maps by combinatorial indexing and technology-agnostic statistical analysis.

BackgroundResolution of complex repeat structures and rearrangements in the assembly and analysis of large eukaryotic genomes is often aided by a combination of high-throughput sequencing and genome-mapping technologies (for example, optical restriction mapping). In particular, mapping technologies can generate sparse maps of large DNA fragments (150 kilo base pairs (kbp) to 2 Mbp) and thus provide a unique source of information for disambiguating complex rearrangements in cancer genomes. Despite their utility, combining high-throughput sequencing and mapping technologies has been challenging because of the lack of efficient and sensitive map-alignment algorithms for robustly aligning error-prone maps to sequences.ResultsWe introduce a novel seed-and-extend glocal (short for global-local) alignment method, OPTIMA (and a sliding-window extension for overlap alignment, OPTIMA-Overlap), which is the first to create indexes for continuous-valued mapping data while accounting for mapping errors. We also present a novel statistical model, agnostic with respect to technology-dependent error rates, for conservatively evaluating the significance of alignments without relying on expensive permutation-based tests.ConclusionsWe show that OPTIMA and OPTIMA-Overlap outperform other state-of-the-art approaches (1.6−2 times more sensitive) and are more efficient (170−200 %) and precise in their alignments (nearly 99 % precision). These advantages are independent of the quality of the data, suggesting that our indexing approach and statistical evaluation are robust, provide improved sensitivity and guarantee high precision.Electronic supplementary materialThe online version of this article (doi:10.1186/s13742-016-0110-0) contains supplementary material, which is available to authorized users.

The Release 6 reference sequence of the Drosophila melanogaster genome

Drosophila melanogaster plays an important role in molecular, genetic, and genomic studies of heredity, development, metabolism, behavior, and human disease. The initial reference genome sequence reported more than a decade ago had a profound impact on progress in Drosophila research, and improving the accuracy and completeness of this sequence continues to be important to further progress. We previously described improvement of the 117-Mb sequence in the euchromatic portion of the genome and 21 Mb in the heterochromatic portion, using a whole-genome shotgun assembly, BAC physical mapping, and clone-based finishing. Here, we report an improved reference sequence of the single-copy and middle-repetitive regions of the genome, produced using cytogenetic mapping to mitotic and polytene chromosomes, clone-based finishing and BAC fingerprint verification, ordering of scaffolds by alignment to cDNA sequences, incorporation of other map and sequence data, and validation by whole-genome optical restriction mapping. These data substantially improve the accuracy and completeness of the reference sequence and the order and orientation of sequence scaffolds into chromosome arm assemblies. Representation of the Y chromosome and other heterochromatic regions is particularly improved. The new 143.9-Mb reference sequence, designated Release 6, effectively exhausts clone-based technologies for mapping and sequencing. Highly repeat-rich regions, including large satellite blocks and functional elements such as the ribosomal RNA genes and the centromeres, are largely inaccessible to current sequencing and assembly methods and remain poorly represented. Further significant improvements will require sequencing technologies that do not depend on molecular cloning and that produce very long reads.

Optical Whole-Genome Restriction Mapping as a Tool for Rapidly Distinguishing and Identifying Bacterial Contaminants in Clinical Samples.

Optical restriction genome mapping is a technology in which a genome is linearized on a surface and digested with specific restriction enzymes, giving an arrangement of the genome with gaps whose order and size are unique for a given organism. Current applications of this technology include assisting with the correct scaffolding and ordering of genomes in conjunction with whole-genome sequencing, observation of genetic drift and evolution using comparative genomics and epidemiological monitoring of the spread of infections. Here, we investigated the suitability of genome mapping for use in clinical labs as a potential diagnostic tool. Using whole genome mapping, we investigated the basic performance of the technology for identifying two bacteria of interest for food-safety (Lactobacilli spp. and Enterohemorrhagic Escherichia coli). We further evaluated the performance for identifying multiple organisms from both simple and complex mixtures. We were able to successfully generate optical restriction maps of four Lactobacillus species as well as a strain of Enterohemorrhagic Escherichia coli from within a mixed solution, each distinguished using a common compatible restriction enzyme. Finally, we demonstrated that optical restriction maps were successfully obtained and the correct organism identified within a clinical matrix. With additional development, whole genome mapping may be a useful clinical tool for rapid invitro diagnostics.

Efficient and scalable scaffolding using optical restriction maps.

In the next generation sequencing techniques millions of short reads are produced from a genomic sequence at a single run. The chances of low read coverage to some regions of the sequence are very high. The reads are short and very large in number. Due to erroneous base calling, there could be errors in the reads. As a consequence, sequence assemblers often fail to sequence an entire DNA molecule and instead output a set of overlapping segments that together represent a consensus region of the DNA. This set of overlapping segments are collectively called contigs in the literature. The final step of the sequencing process, called scaffolding, is to assemble the contigs into a correct order. Scaffolding techniques typically exploit additional information such as mate-pairs, pair-ends, or optical restriction maps. In this paper we introduce a series of novel algorithms for scaffolding that exploit optical restriction maps (ORMs). Simulation results show that our algorithms are indeed reliable, scalable, and efficient compared to the best known algorithms in the literature.

Potential impact on kidney infection: a whole-genome analysis of Leptospira santarosai serovar Shermani

Leptospira santarosai serovar Shermani is the most frequently encountered serovar, and it causes leptospirosis and tubulointerstitial nephritis in Taiwan. This study aims to complete the genome sequence of L. santarosai serovar Shermani and analyze the transcriptional responses of L. santarosai serovar Shermani to renal tubular cells. To assemble this highly repetitive genome, we combined reads that were generated from four next-generation sequencing platforms by using hybrid assembly approaches to finish two-chromosome contiguous sequences without gaps by validating the data with optical restriction maps and Sanger sequencing. Whole-genome comparison studies revealed a 28-kb region containing genes that encode transposases and hypothetical proteins in L. santarosai serovar Shermani, but this region is absent in other pathogenic Leptospira spp. We found that lipoprotein gene expression in both L. santarosai serovar Shermani and L. interrogans serovar Copenhageni were upregulated upon interaction with renal tubular cells, and LSS19962, a L. santarosai serovar Shermani-specific gene within a 28-kb region that encodes hypothetical proteins, was upregulated in L. santarosai serovar Shermani-infected renal tubular cells. Lipoprotein expression during leptospiral infection might facilitate the interactions of leptospires within kidneys. The availability of the whole-genome sequence of L. santarosai serovar Shermani would make it the first completed sequence of this species, and its comparison with that of other Leptospira spp. may provide invaluable information for further studies in leptospiral pathogenesis.

Exploring both sequence detection and restriction endonuclease cleavage kinetics by recognition site via single-molecule microfluidic trapping.

We demonstrate the feasibility of a single-molecule microfluidic approach to both sequence detection and obtaining kinetic information for restriction endonucleases on dsDNA. In this method, a microfluidic stagnation point flow is designed to trap, hold, and linearize double-stranded (ds) genomic DNA to which a restriction endonuclease has been pre-bound sequence-specifically. By introducing the cofactor magnesium, we determine the binding location of the enzyme by the cleavage process of dsDNA as in optical restriction mapping, however here the DNA need not be immobilized on a surface. We note that no special labeling of the enzyme is required, which makes it simpler than our previous scheme using stagnation point flows for sequence detection. Our accuracy in determining the location of the recognition site is comparable to or better than other single molecule techniques due to the fidelity with which we can control the linearization of the DNA molecules. In addition, since the cleavage process can be followed in real time, information about the cleavage kinetics, and subtle differences in binding and cleavage frequencies among the recognition sites, may also be obtained. Data for the five recognition sites for the type II restriction endonuclease EcoRI on λ-DNA are presented as a model system. While the roles of the varying fluid velocity and tension along the chain backbone on the measured kinetics remain to be determined, we believe this new method holds promise for a broad range of studies of DNA-protein interactions, including the kinetics of other DNA cleavage processes, the dissociation of a restriction enzyme from the cleaved substrate, and other macromolecular cleavage processes.

Parallel Arrays of Geometric Nanowells for Assembling Curtains of DNA with Controlled Lateral Dispersion

The analysis of individual molecules is evolving into an important tool for biological research, and presents conceptually new ways of approaching experimental design strategies. However, more robust methods are required if these technologies are to be made broadly available to the biological research community. To help achieve this goal we have combined nanofabrication techniques with single-molecule optical microscopy for assembling and visualizing curtains comprised of thousands of individual DNA molecules organized at engineered diffusion barriers on a lipid bilayer-coated surface. Here we present an important extension of this technology that implements geometric barrier patterns comprised of thousands of nanoscale wells that can be loaded with single molecules of DNA. We show that these geometric nanowells can be used to precisely control the lateral distribution of the individual DNA molecules within curtains assembled along the edges of the engineered barrier patterns. The individual molecules making up the DNA curtain can be separated from one another by a user-defined distance dictated by the dimensions of the nanowells. We demonstrate the broader utility of these patterned DNA curtains in a novel, real time restriction assay that we refer to as dynamic optical restriction mapping, which can be used to rapidly identify entire sets of cleavage sites within a large DNA molecule.

The genome of the diatom Thalassiosira pseudonana: ecology, evolution, and metabolism.

Diatoms are unicellular algae with plastids acquired by secondary endosymbiosis. They are responsible for approximately 20% of global carbon fixation. We report the 34 million-base pair draft nuclear genome of the marine diatom Thalassiosira pseudonana and its 129 thousand-base pair plastid and 44 thousand-base pair mitochondrial genomes. Sequence and optical restriction mapping revealed 24 diploid nuclear chromosomes. We identified novel genes for silicic acid transport and formation of silica-based cell walls, high-affinity iron uptake, biosynthetic enzymes for several types of polyunsaturated fatty acids, use of a range of nitrogenous compounds, and a complete urea cycle, all attributes that allow diatoms to prosper in aquatic environments.

A new method for straightening DNA molecules for optical restriction mapping.

We have developed an improved method of straightening DNA molecules for use in optical restriction mapping. The DNA was straightened on 3-aminopropyltriethoxysilane-coated glass slides using surface tension generated by a moving meniscus. In our method the meniscus motion was controlled mechanically, which provides advantages of speed and uniformity of the straightened molecules. Variation in the affinity of the silanized surfaces for DNA was compensated by precoating the slide with single-stranded non-target blocking DNA. A small amount of MgCl2 added to the DNA suspension increased the DNA-surface affinity and was necessary for efficient restriction enzyme digestion of the straightened surface-bound DNA. By adjusting the amounts of blocking DNA and MgCl2, we prepared slides that contained many straight parallel DNA molecules. Straightened lambda phage DNA (48 kb) bound to a slide surface was digested by EcoRI restriction endonuclease, and the resulting restriction fragments were imaged by fluorescence microscopy using a CCD camera. The observed fragment lengths showed excellent agreement with their predicted lengths.

Identification of seven new human MHC class I region genes around the HLA-F locus.

Using cDNA hybridization selection techniques, we identified seven new genes in a 280 kilobase YAC covering the HLA-F locus. The new genes were mapped back to the YAC by a combination of optical restriction mapping and pulse field gel electrophoresis. Northern analysis of individual clones demonstrated the presence of either different mRNA sizes or different expression patterns. Two of the cDNA clones were expressed only in lymphoid cell lines: one in Jurkat cells (T cell) and another in JY cells (B cell). All the genes lacked sequence similarity to any known classical and non-classical major histocompatibility complex (MHC) class I genes, indicating that the MHC class I region has more functions than anticipated. Of the seven new genes, one is highly similar (97%) to mouse 60S ribosomal protein, and another is homologous to diubiquitin proteins. Of the two G-coupled receptor-like cDNAs, one was fully sequenced and found to be an olfactory receptor-like gene. The study strengthens evidence that the MHC complex not only plays a key role in the immune system, but also contributes to non-immunological functions.